如何讓aov向我展示F統計量和p值？

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