# Load the data
data("EEGwide", package = "MANOVA.RM")
vars <- colnames(EEGwide)[1:6]
# Part the data into six groups of sex and diagnosis
X_list <- list(t(EEGwide[EEGwide$sex=="M" & EEGwide$diagnosis=="AD",vars]),
t(EEGwide[EEGwide$sex=="M" & EEGwide$diagnosis=="MCI",vars]))
nv <- unlist(lapply(X_list, ncol))
set.seed(31415)
test_combined(X_list, nv)
Run the code above in your browser using DataLab