使用R语言MASS包的manova函数进行双因素多元方差分析,一共有18个因变量,在summary结果的时候,报错:“Error in summary.manova(fit) : residuals have rank 17 < 18”,去掉一个因变量就可以成功输出summary结果,请问这个报错的原因是什么呢?
是MASS包的,代码非常简单,data为原数据,X1~X18为18个因变量,position和class为双因素:
X <- select(data,x1:x18,position,class)
class <- factor(X$class)
position<- factor(X$position)
fit1<-manova(X~class*position)