несколько случайных эффектов в NLME - PullRequest
0 голосов
/ 03 октября 2019

Мне нужно включить еще один случайный эффект в мою модель 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, и я думаю, что это не может быть правильным. Результаты следующие: enter image description here

Спасибо за вашу помощь! лил

...