Я пытаюсь запустить модель смешанного эффекта, которая состоит из 2 переменных с фиксированным эффектом, где первый фиксированный эффект имеет два уровня, а другой представляет непрерывные данные. Непрерывные данные - это pH, и у меня есть 12 уникальных значений, которые связаны с различными участками моего блока лечения. Поскольку я использую непрерывные данные, использование среды cld (emmeans ()) дает мне только среднее значение pH из 12 уникальных значений, поэтому я заинтересован в использовании emtrends () для разделения pH по отдельным значениям, а неединственное усредненное значение pH. Я пытаюсь запустить функцию emtrends (), чтобы определить, значительно ли отличаются уклоны друг от друга
Я уже пробовал код, который следует позже с этим комментарием. Когда я запускаю свою модель, мы обнаруживаем, что pH значительно влияет на переменную отклика. Однако, когда я использую функцию emtrends (), чтобы понять, как pH влияет на переменную ответа, я получаю следующую ошибку: «Ошибка в emm_basis.merMod (модель, attr (данные,« термины »), RG @ model).info $ xlev,: аргумент "options" отсутствует, без значения по умолчанию "
Пример кадра данных:
site species nitrogen_treat pH response_var
A B nit_added 3.4 36.8
A B nit_added 3.4 35.7
B A no_nitrogen 5.6 32.1
A C no_nitrogen 5.6 30.1
B D no_nitrogen 4.3 30.2
C C nit_added 3.1 37.2
C A nit_added 7.2 10.2
Итак, мой код для функции lmer выглядит следующим образом:
lm <- lmer(response_var ~ nitrogen_treat * pH +
(1 | site) + (1 | species))
, который я затем использовал, чтобы попробовать свою функцию em_trends (), которая не работает:
emtrends(lm, -1, var = "pH")
Я также экспериментировал с обменом "-1 "с ~ азотной обработкой * pH и всеми возможными комбинациями, но я постоянно получаю сообщение об ошибке, о котором я упоминал выше. Буду признателен за понимание того, почему я получаю сообщения об ошибках, которые я получаю!