bsqrmbt:

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

Usage

bsqrmbt(x, y, alpha = 0.05, bend = 1.28, nboot = 599, SEED = TRUE)

Arguments

x
y
alpha
bend
nboot
SEED

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, bend = 1.28, nboot = 599, SEED = TRUE) 
{
    if (SEED) 
        set.seed(2)
    T <- bsqrm(x, y, alpha, bend)$teststat
    TT <- 0
    bsqrmbt <- numeric(2)
    xone <- onestep(x, bend = bend)
    yone <- onestep(y, bend = bend)
    for (j in 1:nboot) {
        xx <- (sample(x, length(x), replace = TRUE) - xone)
        yy <- (sample(y, length(y), replace = TRUE) - yone)
        TT[j] <- bsqrm(xx, yy, alpha, bend)$teststat
    }
    TT <- sort(TT)
    bott <- round(alpha * nboot) + 1
    bsqrmbt <- TT[nboot - bott]
    pv = mean(T <= TT)
    list(critval = bsqrmbt, teststat = TRUE, p.value = pv)
  }

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