在 R 中重复测量的 ANCOVA

机器算法验证 r 安乔娃
2022-04-12 14:37:07

我是 R 编程的新手,并且一直在尝试执行重复测量 ANCOVA,其中温度为因变量,站点为自变量,日期为协变量,年份为重复测量。

我的数据集包含 20 天、2 个不同年份的 4 个地点的温度。地点、日期和年份似乎对温度有显着影响,但地点和年份或地点和日期之间没有显着的相互作用。我最初执行了两个独立的 ANCOVA,但我的合著者希望它们在一个模型中。据我所知,没有像 ANOVA 那样进行重复测量 ANCOVA 的简单方法。

我尝试过使用: lm(formula = Temp~Site*Date + Year),但这不只是将 year 作为不影响模型的交互吗?

您可以提供的任何帮助将不胜感激!

3个回答

这里有一个关于这个主题的教程列表:

http://www.r-statistics.com/2010/04/repeated-measures-anova-with-r-tutorials/

祝你好运。

样本数据?

set.seed(5)
d <- expand.grid(Site=LETTERS[1:4], Date=1:20, Year=factor(1:2))
d$Temp <- round(rnorm(nrow(d), mean=60, sd=15))

也许你可以试试

lm(formula = Temp~Site*(Date + Year))

这样,您将与站点进行两次交互,并且日期和年份之间将没有交互