# Some utility functions
alpha <- function(conf.level){
c(1 - conf.level, 1 + conf.level) / 2
}
errorNA <- function(var){
deparse(substitute(var))
not.ok <- anyNA(var)
if(not.ok) stop("missing values in '",deparse(substitute(var)),"'") else (invisible(NULL))
}
# g
`edgeworth.skew` <- function(x,w,t){
v <- sum(x*w^2)
b1 <- (sum(x*w^3) / (v^(3/2)))^2
b2 <- sum((x*w^4)*(1+3*x))/(v^2)
sqrt(b1)/6 + ((b2-3)/8 +15*b1/72) -
(sqrt(b1)/6)*t^2 - ((b2-3)/24 + 4*b1/36)*t^3 -
(b1/72)*t^5
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.