Nothing
"epi.bohning" <- function(obs, exp, alpha = 0.05){
J <- length(obs)
smr <- obs / exp
smr.bar <- sum(smr, na.rm = TRUE) / J
# Bohning's test:
top <- (1 / (J - 1)) * sum(((obs - (smr.bar * exp))^2 / exp), na.rm = TRUE) - smr.bar
bottom <- sqrt((2 * smr.bar) / (J - 1))
bohning <- top / bottom
p <- 1 - pnorm(q = bohning, mean = 0, sd = 1)
rval <- as.data.frame(cbind(test.statistic = bohning, p.value = p))
return(rval)
}
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.