problèmes de syntaxe MANOVA.RM dans R

Question

J'ai une trame de données composée de trois variables de réponse continues, deux variables et un sujet ID

library(MANOVA.RM) 
y1<-sample(1:150, 100, replace=T) 
y2<-sample(1:150, 100, replace=T) 
y3<-sample(1:150, 100, replace=T) 
x1<-sample(x=c("S1", "S2"), size=100, replace=T, prob=rep(1/2,2)) 
x2<-sample(x=c("T1", "T2"), size=100, replace=T, prob=rep(1/2,2)) 
id <- seq (from=1, to=100, by=1) 

data <- data.frame(y1, y2, y3, x1, x2, id) 

Je suis en train d'exécuter un MANOVA mais ai des problèmes avec la syntaxe autour de l'attribution de variables de réponse.

Le fichier d'aide pour MANOVA.RM utilise la syntaxe suivante comme exemple

EEG_MANOVA <- MANOVA(resp ~ sex * diagnosis, data = EEG, subject = "id", resampling = "paramBS", iter = 1000, alpha = 0.05)

Cependant, le côté gauche de la formule (ie « resp ») ne se réfère à une seule variable de réponse, tandis que je veux faire référence à une matrice de variables (dans ce cas y1, y2, and y3)

Comment puis-je modifier la syntaxe pour inclure trois variables de réponse?

Answer 1

Il devrait fonctionner si vous liez les colonnes de réponse comme (cela fonctionne dans manova paquet):

test <- manova(cbind(y1,y2,y3) ~ x1, data = data, subject = "id")