bi2CR:

In musto101/wilcox_R: What the package does (short line)

Usage

bi2CR(k1, n1, k2, n2, alpha = 0.05, fin = 100, xlab = "p1", ylab = "p2")

Arguments

k1
n1
k2
n2
alpha
fin
xlab
ylab

Examples

##---- 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 (k1, n1, k2, n2, alpha = 0.05, fin = 100, xlab = "p1", 
    ylab = "p2") 
{
    nom = 1 - alpha
    x = rep(0:fin/fin, (fin + 1))
    y = rep(0:fin/fin, rep((fin + 1), (fin + 1)))
    indik = rep(0, (fin + 1)^2)
    for (i1 in 0:fin) {
        p1 = i1/fin
        pr1 = dbinom(0:n1, n1, p1)
        for (i2 in 0:fin) {
            p2 = i2/fin
            pr2 = dbinom(0:n2, n2, p2)
            jpr = rep(pr1, n2 + 1) * rep(pr2, rep(n1 + 1, n2 + 
                1))
            obs = dbinom(k1, n1, p1) * dbinom(k2, n2, p2)
            sumpr = sum(jpr[jpr > obs])
            if (sumpr < nom) {
                k = (i1 + 1) + (fin + 1) * i2
                indik[k] = indik[k] + 1
            }
        }
    }
    plot(x[indik == 1], y[indik == 1], xlim = c(0, 1), ylim = c(0, 
        1), col = "green", pch = 15, xlab = xlab, ylab = ylab)
    points(k1/n1, k2/n2, pch = 19)
  }

musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.