Nothing
pgh <-
function(q,A,B,g,h,log.p=FALSE) {
sapply(q,function(q) {
if(g==0 && h==0) l=stats::pnorm(q,A,B)
else if(h==0) {
q=(q-A)/B
if((q<=-1/g && g>0) || (q>=-1/g && g<0)) l=NaN
else {
l=stats::pnorm(log(1+g*q)/g)
}
}
else if(h>0) {
q=(q-A)/B
z=gh_inv(q,g,h)
l=stats::pnorm(z)
}
else {
base::stop("negative kurtosis parameter")
}
if(log.p==FALSE) l
else log(l)
})
}
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.