powered by
Hotelling(Y, X)
manova
Y1<-rnorm(100,15,3) Y2<-rnorm(100,17,2) X<-factor(c(rep(1,50),rep(2,50))) Hotelling(cbind(Y1,Y2),X) anova(lm(cbind(Y1,Y2)~X))
Run the code above in your browser using DataLab