我想展示的是,不同物种的种子显示出不同的长度,这是由于物种的因素,对于每一个物种,我有几棵树,对于每一棵树,我有几颗测量过的种子。
我用R做了方差分析:
summary(aov(Length ~ Species))
然而,审阅者注意到了一个独立性的问题,因为种子可能从同一棵树提供。(而这确实是一个真实的的问题!)
为了回答这个问题,我想我应该做一个嵌套方差分析,对吗?
但是,有很多方法可以编写代码:
summary(aov(Length ~ Species*Tree))
summary(aov(Length ~ Tree*Species))
summary(aov(Length ~ Species/Tree))
summary(aov(Length ~ Species+Error(Tree)))
我相信这是列出的最后一种可能性,它将允许我表明种子的长度因物种而异,并考虑到种子可能来自同一棵树。
你能证实吗?
当我运行该命令时,我获得以下内容:
Error: Tree
Df Sum Sq Mean Sq F value Pr(>F)
Species 12 320.6 26.715 14.98 4.96e-15 ***
Residuals 71 126.6 1.784
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Error: Within
Df Sum Sq Mean Sq F value Pr(>F)
Residuals 1541 11.92 0.007733
这确实意味着物种对种子长度有显著的影响,对吗?
非常感谢你的帮助!!
穆里尔
1条答案
按热度按时间ubbxdtey1#
有关R中嵌套ANOVA的一些示例以及对混合模型的一些了解,请参见here。
我会安装
lme4
包,在R中运行?lmer
,然后查看提供的页面上的“混合和多级模型”部分,也许这是处理数据的更好方法。