acbinomciv2:

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

Usage

acbinomciv2(x = sum(y), nn = length(y), y = NULL, n = NA, alpha = 0.05, nullval = 0.5)

Arguments

x
nn
y
n
alpha
nullval

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 (x = sum(y), nn = length(y), y = NULL, n = NA, alpha = 0.05, 
    nullval = 0.5) 
{
    res = acbinomci(x = x, nn = nn, y = y, alpha = alpha)
    ci = res$ci
    alph <- c(1:99)/100
    for (i in 1:99) {
        irem <- i
        chkit <- acbinomci(x = x, nn = nn, y = y, alpha = alph[i])$ci
        if (chkit[1] > nullval || chkit[2] < nullval) 
            break
    }
    p.value <- irem/100
    if (p.value <= 0.1) {
        iup <- (irem + 1)/100
        alph <- seq(0.001, iup, 0.001)
        for (i in 1:length(alph)) {
            p.value <- alph[i]
            chkit <- acbinomci(x = x, nn = nn, y = y, alpha = alph[i])$ci
            if (chkit[1] > nullval || chkit[2] < nullval) 
                break
        }
    }
    if (p.value <= 0.001) {
        alph <- seq(1e-04, 0.001, 1e-04)
        for (i in 1:length(alph)) {
            p.value <- alph[i]
            chkit <- acbinomci(x = x, nn = nn, y = y, alpha = alph[i])$ci
            if (chkit[1] > nullval || chkit[2] < nullval) 
                break
        }
    }
    list(n = nn, phat = res$phat, ci = res$ci, p.value = p.value)
  }

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