R中的重復測量方差分析？

Question

我在恢復過程的三個階段（之前，期間和之后兩年）觀察兩個站點（一個經歷了棲息地恢復，另一個是實驗控制）的平均家庭范圍大小。 我想看看不同站點和時段的平均家庭范圍大小的差異是否不同。 基於具有兩個分類變量（站點和時段），我假設這將使用重復測量ANOVA來完成？ 我需要查看將使用哪些代碼，因為我之前從未在R中進行過ANOVA。

rm (list = ls())
hrdata=read.csv(xxx)
hrdata

Answer 1

我認為你可以用線性模型做到這一點，但請參閱（ https://stats.stackexchange.com/questions/20002/regression-vs-anova-discrepancy-aov-vs-lm-in-r ）進行討論anova vs回歸。

代碼看起來像這樣：

lm1 <- lm(HRS ~ Site * Period, data=hrdata)

這段代碼的第一部分就是將這個線性模型（ lm ）存儲為R對象，我們將其命名為lm1 。

然后你可以這樣做：

summary(lm1)

這將考察Site（棲息地恢復與控制），Period（之前，期間和之后）以及兩者之間的相互作用的影響。

有很多關於解釋這些摘要結果的帖子。 我在下面貼了一些。 如果您不確定如何解釋交互術語，則第一個可能很有用：

https://stats.stackexchange.com/questions/56784/how-to-interpret-the-interaction-term-in-lm-formula-in-r

https://stats.stackexchange.com/questions/59250/how-to-interpret-the-output-of-the-summary-method-for-an-lm-object-in-r

https://stats.stackexchange.com/questions/115304/interpreting-output-from-anova-when-using-lm-as-input