bdms1:

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

Usage

bdms1(x, con)

Arguments

x
con

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, con) 
{
    pool <- x[[1]]
    JK <- length(x)
    for (j in 2:JK) pool <- c(pool, x[[j]])
    N <- length(pool)
    rval <- rank(pool)
    rvec <- list()
    up <- length(x[[1]])
    rvec[[1]] <- rval[1:up]
    rbar <- mean(rvec[[1]])
    nvec <- length(rvec[[1]])
    for (j in 2:JK) {
        down <- up + 1
        up <- down + length(x[[j]]) - 1
        rvec[[j]] <- rval[down:up]
        nvec[j] <- length(rvec[[j]])
        rbar[j] <- mean(rvec[[j]])
    }
    phat <- (rbar - 0.5)/N
    phat <- as.matrix(phat)
    svec <- NA
    for (j in 1:JK) svec[j] <- sum((rvec[[j]] - rbar[j])^2)/(nvec[j] - 
        1)
    svec <- svec/N^2
    VN <- N * diag(svec/nvec)
    top <- con[1, 1] * sum(diag(VN))
    Ftest <- N * (t(phat) %*% con %*% phat)/top
    nu1 <- con[1, 1]^2 * sum(diag(VN))^2/sum(diag(con %*% VN %*% 
        con %*% VN))
    lam <- diag(1/(nvec - 1))
    nu2 <- sum(diag(VN))^2/sum(diag(VN %*% VN %*% lam))
    sig <- 1 - pf(Ftest, nu1, nu2)
    list(F = Ftest, nu1 = nu1, nu2 = nu2, q.hat = phat, p.value = sig)
  }

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