Мне нужно включить еще один случайный эффект в мою модель lme
. Однако эти 2 случайных эффекта не связаны, поэтому не могут быть вложенными.
Я не уверен, как написать код для этого. Я пробовал следующее:
modelA <- lme(FT~ Init.Age + Sex + Status, random= (1|Patient) + (1|Site), data = data, na.action = na.omit)
modelB <- lme(FT~ Init.Age + Sex + Status, random= list (~1|Patient, ~1|Site), data = data, na.action = na.omit)
работает только модель B. Но это дает мне точно такие же результаты, как random= ~1|Patient
, и я думаю, что это не может быть правильным. Результаты следующие:
Спасибо за вашу помощь! лил