Повторные меры ANOVA в R?

Question

Я смотрю на средний домашний диапазон на двух участках (один подвергся восстановлению среды обитания, а другой - экспериментальный контроль) в течение трех этапов процесса восстановления (до, во время и через два года после). Я хочу увидеть, отличаются ли различия в среднем размере домашнего диапазона между сайтами и периодами. Исходя из наличия двух категориальных переменных (место и период), я предполагаю, что это будет сделано с помощью повторных измерений ANOVA? Мне нужно было посмотреть, какой код будет использоваться, поскольку я никогда раньше не делал ANOVA в R.

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) для обсуждения анова против регрессии.

код будет выглядеть примерно так:

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

Первый бит этого кода просто хранит эту линейную модель (lm) как объект R, который мы назвали lm1.

тогда вы можете сделать:

summary(lm1)

Здесь будут рассмотрены эффекты Участка (восстановление среды обитания по сравнению с контролем), Период (до, во время и после), а также взаимодействие между ними.

Есть много постов о интерпретации этих итоговых результатов. Я разместил некоторые ниже. Первый может быть полезен, если вы не знаете, как интерпретировать условия взаимодействия:

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