TWOpovPV:

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

Usage

TWOpovPV(x, y, alpha = 0.05, CN = FALSE)

Arguments

x
y
alpha
CN

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, y, alpha = 0.05, CN = FALSE) 
{
    alph <- c(1:99)/100
    for (i in 1:99) {
        irem <- i
        chkit <- TWOpov(x, y, alpha = alph[i], CN = CN)$ci
        if (sign(chkit[1] * chkit[2]) == 1) 
            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]
            alph <- c(1:99)/100
            for (i in 1:99) {
                irem <- i
                chkit <- TWOpov(x, y, alpha = alph[i], CN = CN)$ci
                if (sign(chkit[1] * chkit[2]) == 1) 
                  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 <- TWOpov(x, y, alpha = alph[i], CN = CN)$ci
            if (sign(chkit[1] * chkit[2]) == 1) 
                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 <- TWOpov(x, y, alpha = alph[i], CN = CN)$ci
            if (sign(chkit[1] * chkit[2]) == 1) 
                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 <- TWOpov(x, y, alpha = alph[i], CN = CN)$ci
            if (sign(chkit[1] * chkit[2]) == 1) 
                break
        }
    }
    res = TWOpov(x, y, alpha = alpha, CN = CN)
    list(p.value = p.value, est.rho1 = res$est.rho1, est.rho2 = res$est.rho2, 
        ci = res$ci)
  }

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