MANOVA中的R错误消息： 'dimnames'[1]不等于阵列程度的长度

Question

试图对这个数据执行MANOVA：

创建data.frame：

acc <- data.frame(Degrees = c("5","8","10"), MPH10=c(0.35, 0.37, 0.32), 
MPH25=c(0.19, 0.28, 0.30), MPH40=c(0.14, 0.19, 0.29), MPH55=c(0.10, 0.19, 0.23)) 

校验数据.frame：

acc 
    Degrees MPH10 MPH25 MPH40 MPH55 
1  5 0.35 0.19 0.14 0.10 
2  8 0.37 0.28 0.19 0.19 
3  10 0.32 0.30 0.29 0.23 

我输入：

acc_manova <- manova(cbind(MPH10,MPH25,MPH40,MPH55) ~ Degrees, data = acc) 

然后运行它：

acc_manova 

我得到一个错误信息：

Call: 
    manova(cbind(MPH10, MPH25, MPH40, MPH55) ~ as.factor(Degrees), 
    data = acc) 

Terms: 
Error in dimnames(tmp) <- list(c(rn, "Deg. of Freedom"), nmeffect) : 
    length of 'dimnames' [1] not equal to array extent 

所以我想它与度列的名称做：D05，D08，D10，所以我放弃了d和0名持有人。有相同的错误信息

然后我添加as.factor（度），再次运行acc_manova，并提出了相同的错误。

对此的任何想法？

Answer 1

您的学位栏不是数字，而是一个因子（分类数据）。将因子更改为数字可解决您的问题：

acc$Degrees = as.numeric(acc$Degrees) 
acc_manova <- manova(cbind(MPH10,MPH25,MPH40,MPH55) ~ Degrees, data = acc) 
> acc_manova 
Call: 
    manova(cbind(MPH10, MPH25, MPH40, MPH55) ~ Degrees, data = acc) 

Terms: 
        Degrees Residuals 
resp 1   1.2500e-03 1.6667e-05 
resp 2   0.000200000 0.006666667 
resp 3   0.005000000 0.006666667 
resp 4   0.000800000 0.008066667 
Deg. of Freedom   1   1 

Residual standard error: 0.004082483 0.08164966 0.08164966 0.08981462 
Estimated effects may be unbalanced