代码之家  ›  专栏  ›  技术社区  ›  unutbu

如何让aov显示F统计量和p值?

  •  1
  • unutbu  · 技术社区  · 16 年前

    下面的脚本

    #!/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

    1 回复  |  直到 16 年前
        1
  •  4
  •   Tristan    16 年前

    你应该在公式中用+而不是*吗?

    > 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