#r #anova
#r #anova
Вопрос:
Я смотрю на средний размер домашнего ареала на двух участках (один из которых подвергся восстановлению среды обитания, а другой является экспериментальным контролем) в течение трех этапов процесса восстановления (до, во время и через два года после).). Я хочу посмотреть, различаются ли различия в среднем размере домашнего диапазона на разных сайтах и в разные периоды. Основываясь на наличии двух категориальных переменных (site и period), я предполагаю, что это будет сделано с использованием повторяющихся мер ANOVA? Мне нужно было посмотреть, какой код будет использоваться, поскольку я никогда раньше не делал ANOVA в R.
rm (list = ls())
hrdata=read.csv(xxx)
hrdata
Ответ №1:
Я думаю, вы могли бы сделать это с помощью линейной модели, но смотрите (https://stats.stackexchange.com/questions/20002/regression-vs-anova-discrepancy-aov-vs-lm-in-r ) для обсуждения anova против регрессии.
код будет выглядеть примерно так:
lm1 <- lm(HRS ~ Site * Period, data=hrdata)
Первый бит этого кода просто сохраняет эту линейную модель ( lm ) как объект R, который мы назвали lm1 .
тогда вы можете сделать:
summary(lm1)
Это будет рассматривать влияние сайта (восстановление среды обитания против контроля), периода (до, во время и после) и взаимодействия между ними.
Есть много сообщений об интерпретации этих итоговых результатов. Я опубликовал некоторые из них ниже. Этот первый может быть полезен, если вы не уверены, как интерпретировать условия взаимодействия: