如何在LME中指定1级随机效应

Question

使用下面的代码在两级之间估计随机截距。

但是，令人感兴趣的是，在截距范围内，截距在整个时间范围内是否存在可变性。

谁能给我一个提示，说明如何在公式中指定随机效果？

非常感谢。

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)

最好

帕特里克