R中的嵌套Anova

n53p2ov0  于 2023-03-15  发布在  其他
关注(0)|答案(1)|浏览(164)

我想展示的是,不同物种的种子显示出不同的长度,这是由于物种的因素,对于每一个物种,我有几棵树,对于每一棵树,我有几颗测量过的种子。
我用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

这确实意味着物种对种子长度有显著的影响,对吗?
非常感谢你的帮助!!
穆里尔

ubbxdtey

ubbxdtey1#

有关R中嵌套ANOVA的一些示例以及对混合模型的一些了解,请参见here
我会安装lme4包,在R中运行?lmer,然后查看提供的页面上的“混合和多级模型”部分,也许这是处理数据的更好方法。

相关问题