This example comes from the textbook section 10.5.3 (p 459-462).
library(permuter) data(monachus) p <- dim(monachus)[2] - 1 n <- table(monachus[, 1]) n contr <- rep(c(1/sum(n[2]), -1/sum(n[1])), n[2:1]) B <- 100 T <- array(0, dim = c((B + 1), p)) g <- monachus[, 1] data <- monachus[, -1] T[1, ] <- t(data) %*% contr for (bb in 2:(B + 1)) { data.star <- data[sample(1:sum(n)), ] T[bb, ] <- t(data.star) %*% contr } P <- t2p_old(abs(T)) colnames(P) <- colnames(data) cat("Rough p-values \n") t(t(P[1, ]))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.