rqfitpv:

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

Usage

rqfitpv(x, y, alpha = 0.05, nullval = rep(0, ncol(cbind(x, y))), qval = 0.5, xout = FALSE, outfun = out, ...)

Arguments

x
y
alpha
nullval
qval
xout
outfun
...

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, nullval = rep(0, ncol(cbind(x, 
    y))), qval = 0.5, xout = FALSE, outfun = out, ...) 
{
    stop("This function has been removed. Use qregci")
    x <- as.matrix(x)
    p <- ncol(x)
    np <- p + 1
    output <- matrix(NA, ncol = 4, nrow = np)
    dimnames(output) <- list(NULL, c("Param.", "ci.low", "ci.up", 
        "p.value"))
    for (j in 1:np) {
        output[j, 1] <- j - 1
        ci <- rqfit(x, y, qval = qval, xout = xout, outfun = outfun, 
            ...)$ci[j, ]
        output[j, 2] <- ci[1]
        output[j, 3] <- ci[2]
        alph <- c(1:99)/100
        for (i in 1:99) {
            irem <- i
            chkit <- rqfit(x, y, alpha = alph[i], qval = qval, 
                xout = xout, outfun = outfun, ...)$ci[j, ]
            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 <- rqfit(x, y, alpha = alph[i], qval = qval, 
                  xout = xout, outfun = outfun, ...)$ci[j, ]
                if (chkit[1] > nullval || chkit[2] < nullval) 
                  break
            }
        }
        output[j, 4] <- p.value
    }
    output
  }

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