# Load the data
data("EEGwide", package = "MANOVA.RM")
# Select only the males with the diagnosis AD
X <- as.matrix(EEGwide[EEGwide$sex == "W" & EEGwide$diagnosis == "AD",
c("brainrate_temporal", "brainrate_frontal",
"brainrate_central","complexity_temporal",
"complexity_frontal", "complexity_central")])
test_covariance_structure(X = X, structure = "diagonal", method = "MC")
Run the code above in your browser using DataLab