不能R中

Question

我尝试使用不复制弄清楚双因素方差分析方差分析（）和似乎无法得到我尝试了这种方式的正确程序（在不同的岗位）。

manova in R error message: length of 'dimnames' [1] not equal to array extent

文本问题与踏板旋转和初始速度是加速度的预测因素有关。

这里是数据：

acc <- data.frame(Degrees = c("d5","d8","d10"), 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)) 

ACC

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 

不知道下一步怎么办。

Answer 1

传统的答案是：

macc <- melt(acc, id.var="Degrees") 
lm(value ~ Degrees + variable, macc) 
anova(lm(value ~ Degrees + variable, macc)) 

而所有剩下的正在建设中的结果的正确描述。 （请注意，我使用“+”而不是“*”）。

anova(lm(value ~ Degrees * variable, macc)) 

你可能已经编写任何一个或者度或MPH变量既作为数字和得到的不饱和型号：施工时饱和模型（一个无残留物）使用交互模式，当你得到一个近乎完美的答案。但它仍然会增加描述结果的复杂性。

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)) 
macc <- melt(acc, id.var="Degrees") 
anova(lm(value ~ Degrees * variable, macc)) 

使用sub从字符变量中删除“MPH”。我认为有必要使用as.numeric(as.character())作为我认为是因素变量，但sub操作显然剥离了因子属性，仅使用as.numeric就足够了。

macc$speed <- as.numeric(sub("MPH", "", macc$variable)) 
anova(lm(value ~ Degrees + speed, macc)) 
# output omitted 
anova(lm(value ~ Degrees * speed, macc)) 
#------------------- 
    Analysis of Variance Table 

Response: value 
       Df Sum Sq Mean Sq F value Pr(>F)  
Degrees  1 0.016827 0.016827 16.904 0.003384 ** 
speed   1 0.048735 0.048735 48.959 0.000113 *** 
Degrees:speed 1 0.006367 0.006367 6.396 0.035309 * 
Residuals  8 0.007963 0.000995      
--- 
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1