In statlab/permuter: permuter: Permutation Tests

Vignette Info

This example comes from the textbook section 8.3 (p 397-398). It involves testing for stochastic ordering of time curves.

library(permuter)
data(rats)

g <- rats[, 1]
B <- 100
data <- rats[, -c(1, 2)]
p <- dim(data)[2]
P <- array(0, dim = c((B + 1), p))

for (j in 1:p) {
    P[, j] <- stoch.ord(data[, j], g, alt = -1, B = B)
    cat(j, "variable/", p, "\n")
}


colnames(P) <- colnames(data)
P.FWE <- FWE.minP_old(P)
res <- cbind(P[1, ], P.FWE)
colnames(res) <- c("Raw-p", "Adj-p")
res

statlab/permuter documentation built on May 30, 2019, 9:45 a.m.