在R中的重复测量方差分析中指定Error()项


10

我在定义R中的双向重复测量方差分析的误差项时遇到问题。我的数据包括沿从树中提取的芯线的三个径向位​​置(内部,中间和外部)的木材密度估计。一共有20种树,每种树有6个个体,每个树有两个核心。

为了测试径向位置对木材密度的影响,我使用以下两种方差分析模型,其中的误差项解释了个体之间的差异:

radpos.aov <- aov(WD ~ Species*Radialposition + Error(Individual), data=Radpos)

但是,我不确定我对错误术语的说明是否足够。我还应该考虑核心内的可变性吗?对我来说,这种可变性与径向位置是相同的,而径向位置是我感兴趣的主要因素。

尽管我花了一些时间阅读有关在“重复测量ANOVA”中指定错误项的信息,但是在实际指定错误项方面仍然存在问题。我将对此有所帮助。

Answers:


8

这将是

 radpos.aov <- aov(WD ~ Species*Radialposition + Error(Individual/(Radialposition)), data=Radpos)
 summary(radpos.aov, type=3)

这说明存在的主题错误Radialposition。如果您还有其他主题内因素,则将它们与Radialposition误差分母一起放入(在交互中)+ Error(Individual/(Radialpostion*wiFactorA))。这就是我的理解。GLM如果您没有丢失的数据,这与SPSS的重复测量相符。


1
+1。请参阅此处的详细备忘单:dwoll.de/rexrepos/posts/anovaMixed.html。同样在这里:personality-project.org/r/r.guide/r.anova.html
变形虫
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.