Nothing
pval.logRR <- function(f0 = 0, f1 = 0, e0, n0, e1, n1,
alternative = "two.sided", RR.null = 1){
e0.m <- e0 + f0
ne0.m <- n0 - e0.m
e1.m <- e1 + f1
ne1.m <- n1 - e1.m
if((abs(e0.m) < 1e-6 & abs(e1.m) < 1e-6) |
(abs(ne0.m) < 1e-6 & abs(ne1.m) < 1e-6)){
pval <- 1
}else{
if(abs(e0.m) < 1e-6 | abs(ne0.m) < 1e-6 |
abs(e1.m) < 1e-6 | abs(ne1.m) < 1e-6){
e0.m <- e0.m + 0.5
ne0.m <- ne0.m + 0.5
e1.m <- e1.m + 0.5
ne1.m <- ne1.m + 0.5
}
logRR <- log(e1.m/(e1.m + ne1.m)) - log(e0.m/(e0.m + ne0.m))
se.logRR <- sqrt(1/e0.m + 1/e1.m -
1/(e0.m + ne0.m) - 1/(e1.m + ne1.m))
if(alternative == "two.sided"){
pval <- 2*pnorm(-abs(logRR - log(RR.null))/se.logRR)
}
if(alternative == "one.sided"){
pval <- pnorm(-abs(logRR - log(RR.null))/se.logRR)
}
}
return(pval)
}
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.