ancovaG:

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

Usage

ancovaG(x1, y1, x2, y2, fr1 = 1, fr2 = 1, tr = 0.2, alpha = 0.05, plotit = FALSE, pts = NULL, sm = FALSE, pr = TRUE, xout = FALSE, outfun = out, test = medpb2, ...)

Arguments

x1
y1
x2
y2
fr1
fr2
tr
alpha
plotit
pts
sm
pr
xout
outfun
test
...

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 (x1, y1, x2, y2, fr1 = 1, fr2 = 1, tr = 0.2, alpha = 0.05, 
    plotit = FALSE, pts = NULL, sm = FALSE, pr = TRUE, xout = FALSE, 
    outfun = out, test = medpb2, ...) 
{
    xy = elimna(cbind(x1, y1))
    x1 = xy[, 1]
    y1 = xy[, 2]
    xy = elimna(cbind(x2, y2))
    x2 = xy[, 1]
    y2 = xy[, 2]
    output = list()
    if (is.null(pts[1])) {
        mat <- matrix(NA, 5, 3)
        dimnames(mat) <- list(NULL, c("X", "n1", "n2"))
        npt <- 5
        isub <- c(1:5)
        xorder <- order(x1)
        y1 <- y1[xorder]
        x1 <- x1[xorder]
        xorder <- order(x2)
        y2 <- y2[xorder]
        x2 <- x2[xorder]
        n1 <- 1
        n2 <- 1
        vecn <- 1
        for (i in 1:length(x1)) n1[i] <- length(y1[near(x1, x1[i], 
            fr1)])
        for (i in 1:length(x1)) n2[i] <- length(y2[near(x2, x1[i], 
            fr2)])
        for (i in 1:length(x1)) vecn[i] <- min(n1[i], n2[i])
        sub <- c(1:length(x1))
        isub[1] <- min(sub[vecn >= 12])
        isub[5] <- max(sub[vecn >= 12])
        isub[3] <- floor((isub[1] + isub[5])/2)
        isub[2] <- floor((isub[1] + isub[3])/2)
        isub[4] <- floor((isub[3] + isub[5])/2)
        for (i in 1:5) {
            mat[i, 1] = x1[isub[i]]
            g1 <- y1[near(x1, x1[isub[i]], fr1)]
            g2 <- y2[near(x2, x1[isub[i]], fr2)]
            g1 <- g1[!is.na(g1)]
            g2 <- g2[!is.na(g2)]
            mat[i, 2] = length(g1)
            mat[i, 3] = length(g2)
            output[[i]] <- test(g1, g2, ...)
        }
    }
    if (!is.null(pts[1])) {
        mat <- matrix(NA, length(pts), 3)
        dimnames(mat) <- list(NULL, c("X", "n1", "n2"))
        n1 <- 1
        n2 <- 1
        vecn <- 1
        for (i in 1:length(pts)) {
            n1[i] <- length(y1[near(x1, pts[i], fr1)])
            n2[i] <- length(y2[near(x2, pts[i], fr2)])
        }
        for (i in 1:length(pts)) {
            mat[i, 1] = pts[i]
            g1 <- y1[near(x1, pts[i], fr1)]
            g2 <- y2[near(x2, pts[i], fr2)]
            g1 <- g1[!is.na(g1)]
            g2 <- g2[!is.na(g2)]
            mat[i, 2] = length(g1)
            mat[i, 3] = length(g2)
            output[[i]] <- test(g1, g2, ...)
        }
    }
    if (plotit) 
        runmean2g(x1, y1, x2, y2, fr = fr1, est = mean, tr = tr, 
            sm = sm, xout = xout, outfun = outfun, ...)
    list(mat, output)
  }

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