Как рассчитать все парные сравнения основных факторов в пакете afex с помощью функции mixed ()? - PullRequest
Новая
       109

Как рассчитать все парные сравнения основных факторов в пакете afex с помощью функции mixed ()?

0 голосов
/ 29 марта 2019

Я пытаюсь смоделировать свою модель смешанных эффектов, используя пакет afex и функцию mixed () (в R), но получаю только общие p-значения для каждого фактора. Как я могу рассчитать парные сравнения основных эффектов? В частности, я тестировал каждое животное 3 раза и мне интересно, в чем различия между тестами.

Вот мой код: 

m1<-mixed(stand~trt*test+(1|calf), 
          method="KR",
          check_contrasts = TRUE,
          data 
print(m1) #summary(m1) only reports correlations of the fixed effects

Выход: Таблица Anova для смешанной модели (тесты типа 3, метод KR) Модель: стенд ~ trt * test + (1 | теленок) Данные: данные Эффект df F p.value 1 трт 1, 22 0,78 .39 2 тест 2, 44 7,93 ** .001 3 трт: тест 2, 44 1,19 .31

К сожалению, я не знал, как подготовить пример данных. Итак, я скопировал часть своих данных. Надеюсь, это поможет.

подставка для теста для теленка trt groom

1 30918, мать 5 275
2 30919 OF контроля 58 313
3 30920 Мама 31 397
4 30921 ОФМ 0 239
5 30922 контроля 16 674
6 30923 контроля 1 434

...