Nothing
fitpot<-function (data, threshold = NA, nextremes = NA, evi = NA)
{
Call <- match.call()
data <- as.numeric(data)
data <- data[!is.na(data)]
data <- sort(data)
n <- length(data)
if (!is.na(nextremes) && !is.na(threshold))
stop("Enter EITHER a threshold or the number of upper extremes")
if (is.na(nextremes) && is.na(threshold)) {
threshold <- min(data)
}
else {
if (!is.na(nextremes))
threshold <- rev(data)[as.numeric(nextremes)]
}
xdat <- data[data >= threshold] - threshold
tes <- ifelse(Tm(xdat, nsim = 0, m = 5, omit = 8)["evi"] >
0.1, T, F)
prob <- length(xdat)/n
aux <- egpd(xdat, evi = evi, heavy = tes)
aux2 <- egpd(xdat, evi=NA, heavy=tes)
if (is.na(evi))
evi <- aux$evi
var.evi <- aux2$evi.var
psi <- aux$psi
var.psi <- aux2$psi.var
out <- list()
out$coeff <- c(evi = evi, psi = psi, threshold = threshold, prob = prob)
class(out) <- "fitpot"
attr(out, "Call") <- Call
attr(out, "evi.sd") <- sqrt(var.evi)
attr(out, "psi.sd") <- sqrt(var.psi)
names(attr(out, "evi.sd")) <- "evi"
names(attr(out, "psi.sd")) <- "psi"
attr(out, "df.residual") <- length(xdat)-1
return(out)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.