我在恢复过程的三个阶段（之前，期间和之后两年）观察两个站点（一个经历了栖息地恢复，另一个是实验控制）的平均家庭范围大小。我想看看不同站点和时段的平均家庭范围大小的差异是否不同。基于具有两个分类变量（站点和时段），我假设这将使用重复测量ANOVA来完成？我需要查看将使用哪些代码，因为我之前从未在R中进行过ANOVA。
rm (list = ls())
hrdata=read.csv(xxx)
hrdata
我认为你可以用线性模型做到这一点，但请参阅（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（之前，期间和之后）以及两者之间的相互作用的影响。
有很多关于解释这些摘要结果的帖子。我在下面贴了一些。如果您不确定如何解释交互术语，则第一个可能很有用：