variance.test.twosample.independent.fx <- function(
fx
,data
,...
) {
summary.out <- summary.impl(fx, data = data, stat.n = T, stat.var = T, format.generate.cellcodes = T)
if (nrow(summary.out) == 2) {
variance.test.twosample.independent.simple(
sample.variance.g1 = summary.out$var[1]
,sample.size.g1 = summary.out$n[1]
,sample.variance.g2 = summary.out$var[2]
,sample.size.g2 = summary.out$n[2]
,...
)
} else if (nrow(summary.out) > 2) {
cmbn <- combn(summary.out$cell.code, 2)
ret <- list()
for (i in 1:ncol(cmbn)) {
g1 <- cmbn[1,i]
g2 <- cmbn[2,i]
idx_1 <- which(summary.out$cell.code == g1)
idx_2 <- which(summary.out$cell.code == g2)
ret[[paste(g1,"vs.",g2)]] <- variance.test.twosample.independent.simple(
sample.variance.g1 = summary.out$var[idx_1]
,sample.size.g1 = summary.out$n[idx_1]
,sample.variance.g2 = summary.out$var[idx_2]
,sample.size.g2 = summary.out$n[idx_2]
,...
)
ret[[paste(g1,"vs.",g2)]]$data.name <- paste(g1,"vs.",g2)
}
ret
} else {
stop("Need to provide at least 2 groups.")
}
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.