Nothing
`dgpd` <-
function(x, m = 0, lambda = 1, xi = 0)
{
k <- xi
if (!SHAPE.XI) k <- -xi
if ((length(m)!=length(lambda))|(length(m)!=length(xi))|(length(xi)!=length(lambda)))
stop("m, lambda and xi should have the same lengths")
n <- length(m)
if ((n > 1) & (length(x) != n))
stop("When vectors of lengths > 1, m, lambda and xi should have the same length as x")
val <- rep(Inf, length(x))
if (n==1)
{
if (k==0)
{
val <- exp(-(x - m)/lambda)/lambda
val[x < m] <- 0
}
else
{
val <- 1. + k * (x - m)/lambda
if (k>0)
{
val[x<=m] <- 0
val[x>m] <- val[x>m]^(-1/k-1)/lambda
}
else
{
val[x<=m | x >= m-lambda/k] <- 0
val[x>m & x < m-lambda/k] <- val[x>m & x < m-lambda/k]^(-1/k-1)/lambda
}
}
}
else
{
K0 <- (k==0)
KPOS <- (k>0)
KNEG <- (k<0)
val[K0] <- exp(-(x[K0] - m[K0])/lambda[K0])/lambda[K0]
val[(x < m) & K0] <- 0
val[!K0] <- 1. + k[!K0] * (x[!K0] - m[!K0])/lambda[!K0]
val[KPOS & (x<=m)] <- 0
val[KPOS & (x>m)] <- val[(!K0) & (x>m)]^(-1/k[(!K0) & (x>m)]-1)/lambda[(!K0) & (x>m)]
val[KNEG &(x<=m | x >= m-lambda/k)] <- 0
val[KNEG & x>m & x < m-lambda/k] <- val[KNEG & x>m & x < m-lambda/k]^(-1/k[KNEG & x>m & x < m-lambda/k]-1)/lambda[KNEG & x>m & x < m-lambda/k]
}
val
}
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.