1 |
r1 |
|
n1 |
|
r2 |
|
n2 |
|
x |
|
y |
|
alpha |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (r1 = sum(x), n1 = length(x), r2 = sum(y), n2 = length(y),
x = NA, y = NA, alpha = 0.05)
{
cr <- qchisq(1 - alpha, 1)
p1 <- r1/n1
p2 <- r2/n2
a <- p1 + p2
b <- p1 - p2
u <- 0.25 * (1/n1 + 1/n2)
v <- 0.25 * (1/n1 - 1/n2)
V <- u * ((2 - a) * a - b^2) + 2 * v * (1 - a) * b
A <- sqrt(cr * (V + cr * u^2 * (2 - a) * a + cr * v^2 * (1 -
a)^2))
B <- (b + cr * v * (1 - a))/(1 + cr * u)
ci <- NA
ci[1] <- B - A/(1 + cr * u)
ci[2] <- B + A/(1 + cr * u)
list(ci = ci)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.