Nothing
dgh <-
function(x,A,B,g,h,log=FALSE) {
sapply(x,function(x) {
x=(x-A)/B
if(g==0 && h==0) v=stats::dnorm(x)/B
else if(h==0) {
if((x<=-1/g && g>0) || (x>=-1/g && g<0)) y=NaN
else y=log(1+g*x)/g
v=stats::dnorm(y)/(B*exp(g*y))
}
else if(h > 0) {
y=gh_inv(x,g,h)
v=stats::dnorm(y)/(B*deriv_gh(y,g,h))
}
else {
base::stop("negative kurtosis parameter")
}
if(log==FALSE) v
else log(v)
})
}
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.