regci:

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

Usage

regci(x, y, regfun = tsreg, nboot = 599, alpha = 0.05, SEED = TRUE, pr = TRUE, null.val = NULL, xout = FALSE, outfun = outpro, plotit = FALSE, xlab = "Predictor 1", ylab = "Predictor 2", ...)

Arguments

x
y
regfun
nboot
alpha
SEED
pr
null.val
xout
outfun
plotit
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 (x, y, regfun = tsreg, nboot = 599, alpha = 0.05, SEED = TRUE, 
    pr = TRUE, null.val = NULL, xout = FALSE, outfun = outpro, 
    plotit = FALSE, xlab = "Predictor 1", ylab = "Predictor 2", 
    ...) 
{
    x <- as.matrix(x)
    p1 <- ncol(x) + 1
    p <- ncol(x)
    xy <- cbind(x, y)
    xy <- elimna(xy)
    x <- xy[, 1:p]
    y <- xy[, p1]
    nrem = length(y)
    estit = regfun(x, y, xout = xout, ...)$coef
    if (xout) {
        if (pr) 
            print("Default for outfun is now outpro, not out")
        m <- cbind(x, y)
        flag <- outfun(x, plotit = FALSE, ...)$keep
        m <- m[flag, ]
        x <- m[, 1:p]
        y <- m[, p1]
    }
    if (is.null(null.val)) 
        null.val = rep(0, p1)
    flagF = identical(regfun, tsreg)
    if (flagF) {
        if (pr) {
            if (sum(duplicated(y) > 0)) 
                print("Duplicate values detected; tshdreg might have more power than tsreg")
        }
    }
    nv = length(y)
    x <- as.matrix(x)
    if (SEED) 
        set.seed(2)
    if (pr) 
        print("Taking bootstrap samples. Please wait.")
    data <- matrix(sample(length(y), size = length(y) * nboot, 
        replace = TRUE), nrow = nboot)
    bvec <- apply(data, 1, regboot, x, y, regfun, xout = FALSE, 
        ...)
    regci <- matrix(0, p1, 5)
    vlabs = "Intercept"
    for (j in 2:p1) vlabs[j] = paste("Slope", j - 1)
    dimnames(regci) <- list(vlabs, c("ci.low", "ci.up", "Estimate", 
        "S.E.", "p-value"))
    ilow <- round((alpha/2) * nboot)
    ihi <- nboot - ilow
    ilow <- ilow + 1
    se <- NA
    pvec <- NA
    for (i in 1:p1) {
        bsort <- sort(bvec[i, ])
        pvec[i] <- (sum(bvec[i, ] < null.val[i]) + 0.5 * sum(bvec[i, 
            ] == null.val[i]))/nboot
        if (pvec[i] > 0.5) 
            pvec[i] <- 1 - pvec[i]
        regci[i, 1] <- bsort[ilow]
        regci[i, 2] <- bsort[ihi]
        se[i] <- sqrt(var(bvec[i, ]))
    }
    if (p1 == 3) {
        if (plotit) {
            plot(bvec[2, ], bvec[3, ], xlab = xlab, ylab = ylab)
        }
    }
    regci[, 3] = estit
    pvec <- 2 * pvec
    regci[, 4] = se
    regci[, 5] = pvec
    list(regci = regci, n = nrem, n.keep = nv)
  }

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