olsLmcp:

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

Usage

olsLmcp(x, y, xout = TRUE, outfun = outpro, ISO = FALSE, STAND = TRUE, alpha = 0.05, pr = TRUE, BLO = FALSE, HC3 = FALSE, ...)

Arguments

x
y
xout
outfun
ISO
STAND
alpha
pr
BLO
HC3
...

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, xout = TRUE, outfun = outpro, ISO = FALSE, STAND = TRUE, 
    alpha = 0.05, pr = TRUE, BLO = FALSE, HC3 = FALSE, ...) 
{
    if (!is.list(x)) 
        stop("Argument x should have list mode")
    if (!is.list(y)) 
        stop("Argument y should have list mode")
    J = length(x)
    x = lapply(x, as.matrix)
    pchk = lapply(x, ncol)
    temp = matl(pchk)
    if (var(as.vector(temp)) != 0) 
        stop("Something is wrong. \nNumber of covariates differs among the groups being compared")
    nv = NULL
    p = ncol(x[[1]])
    p1 = p + 1
    for (j in 1:J) {
        xy = elimna(cbind(x[[j]], y[[j]]))
        x[[j]] = xy[, 1:p]
        y[[j]] = xy[, p1]
        x[[j]] = as.matrix(x[[j]])
        nv = c(nv, nrow(x[[j]]))
    }
    nv.keep = nv
    critrad = NULL
    if (xout) {
        temp = lapply(x, outfun, plotit = FALSE, STAND = STAND, 
            ...)
        for (j in 1:J) {
            x[[j]] = x[[j]][temp[[j]]$keep, ]
            y[[j]] = y[[j]][temp[[j]]$keep]
            nv.keep[j] = length(y[[j]])
        }
    }
    tot = (J^2 - J)/2
    dvec <- alpha/c(1:tot)
    outl = list()
    nr = tot * p1
    outp = matrix(NA, ncol = 5, nrow = tot)
    x = lapply(x, as.matrix)
    xx = list()
    yy = list()
    iall = 0
    ivp = c(1, tot) - tot
    i = 0
    for (j in 1:J) {
        for (k in 1:J) {
            if (j < k) {
                i = i + 1
                xx[[1]] = x[[j]]
                xx[[2]] = x[[k]]
                yy[[1]] = y[[j]]
                yy[[2]] = y[[k]]
                if (!ISO) 
                  all = ols1way(xx, yy, HC3 = HC3)
                if (ISO) 
                  all = ols1wayISO(xx, yy, HC3 = HC3)
                temp = all$adjusted.p.value
                if (is.null(temp)) 
                  temp = all$p.value
                outp[i, 1] = j
                outp[i, 2] = k
                outp[i, 3] = temp
            }
        }
        temp2 <- order(0 - outp[, 3])
        icc = c(1:tot)
        icc[temp2] = dvec
        outp[, 4] = icc
    }
    flag = (outp[, 3] <= outp[, 4])
    outp[, 5] = rep(0, tot)
    outp[flag, 5] = 1
    dimnames(outp) = list(NULL, c("Group", "Group", "p.value", 
        "p.crit", "sig"))
    list(n = nv, n.keep = nv.keep, output = outp)
  }

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