Question

以下脚本

#!/usr/bin/Rscript --vanilla
x <- c(4.5,6.4,7.2,6.7,8.8,7.8,9.6,7.0,5.9,6.8,5.7,5.2)
fertilizer<- factor(c('A','A','A','A','B','B','B','B','C','C','C','C'))
crop <- factor(c('I','II','III','IV','I','II','III','IV','I','II','III','IV'))
av <- aov(x~fertilizer*crop)
summary(av)

产量

                Df  Sum Sq Mean Sq
fertilizer       2 13.6800  6.8400
crop             3  2.8200  0.9400
fertilizer:crop  6  6.5800  1.0967

对于其他数据，aov通常会给出F统计量和相关的p值。导致R省略多汁部分的数据有什么不对/特殊？

Answer 1

你应该在公式中使用+代替*吗？

> summary(aov(x~fertilizer + crop))
            Df  Sum Sq Mean Sq F value  Pr(>F)  
fertilizer   2 13.6800  6.8400  6.2371 0.03426 *
crop         3  2.8200  0.9400  0.8571 0.51218  
Residuals    6  6.5800  1.0967                  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

我如何让aov向我展示F统计量和p值？

1 个答案: