Как указать случайные эффекты уровня 1 в lme - PullRequest
Новая
       59

Как указать случайные эффекты уровня 1 в lme

0 голосов
/ 22 декабря 2018

Используя приведенный ниже код, случайные перехваты оцениваются на уровне между диадами 2.

Однако, что интересно знать, так это то, есть ли изменчивость в перехвате во времени на уровне внутри диад.

Может кто-нибудь дать мне подсказку о том, как указать это в формуле для случайных эффектов?

Заранее большое спасибо. 

tim_mod2bin <- lme(fixed=EERTmn ~ male + female + 
                    (male:time7c)          + (female:time7c)          +                           
                    (male:IERT_Cp)         + (female:IERT_Cp)         + 
                    (male:IERT_Cp_Partner) + (female:IERT_Cp_Partner) +
                                                                                        #final model!
                    (male:PHQ9sB)         + (female:PHQ9sB_Partner)  +                                                  

                    (male:IERT_Cp_Partner:PHQ9sB)  + (female:IERT_Cp_Partner:PHQ9sB_Partner) -1,  

               control=list(maxIter=100000), data=full37finalr,                                                               

               random=~male + female -1| dyade, correlation=corAR1(), na.action=na.omit)

Best,

Patrick

...