R

Question

我想表明，由于物种种类，不同物种的种子显示出不同的长度。 对于每个物种，我有几棵树，每棵树都有几粒种子。

使用R，我做了一个方差分析：

summary(aov(Length ~ Species))

然而，审稿人注意到了一个独立的问题，因为种子可能​​来自同一棵树。 （这确实是一个真正的问题！）

要回答这个问题，我认为我应该做一个嵌套的ANOVA。是吗？

但是，有很多方法可以编写代码：

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 

这确实意味着物种对种子长度有显着影响，是吗？

非常感谢你的帮助!!

穆里尔

Answer 1

有关R中嵌套ANOVA的一些示例以及对混合模型的一些了解，请参阅here。

我将在R中安装软件包lme4，执行?lmer，并查看提供的页面上的“混合和多级模型”部分。也许这是一种更好的数据处理方法。