p<-5
nj<-c(16,8,12,10,13)
y11<-round(rmvnorm(nj[1],mean = xbar71, sigma = s1))
y21<-round(rmvnorm(nj[2],mean = xbar72, sigma = s2))
y31<-round(rmvnorm(nj[3],mean = xbar73, sigma = s3))
y41<-round(rmvnorm(nj[4],mean = xbar74, sigma = s4))
y51<-round(rmvnorm(nj[5],mean = xbar75, sigma = s5))
dfMan1<-data.frame(Y=rbind(y11,y21,y31,y41,y51), IV=factor(rep(1:5, nj)))
manres1<-manova(cbind(Y.1, Y.2)~ IV, data=dfMan1)
manres1
Call:
manova(cbind(Y.1, Y.2) ~ IV, data = dfMan1)
Terms:
IV Residuals
Y.1 4.511039 15.997436
Y.2 0 0
Deg. of Freedom 4 54
Residual standard errors: 0.5442874 0
Предполагаемые эффекты могут быть несбалансированными
summary(manres1)
Ошибка в qr.default (D% % ss [[nt]]% % D, tol = tol): NA / NaN / Inf при вызове сторонней функции (arg 1)
На самом деле пропущенных значений нет.