Тестирование случайного эффекта с использованием RLRsim приводит к случайным эффектам сообщения об ошибке, не зависящим - PullRequest
Новая
       9

Тестирование случайного эффекта с использованием RLRsim приводит к случайным эффектам сообщения об ошибке, не зависящим

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

Проблема как написано в заголовке.

Вопрос в том, знает ли кто-нибудь решение проблемы?Поскольку набор данных содержит данные парного исследования, случайные эффекты не являются независимыми.Тем не менее, я хочу проверить их на значимость.

набор данных и, если предпочтительный r-скрипт replex_RLRsim можно найти здесь 

#Replex for RLRsim:

##reduced model-specific data-set:
datafclr <-read.csv("datafclr.csv", header = TRUE, sep = ",", dec = ".", fill = TRUE)

##Testing female random intercept:

#The full-model maf is already prepared --> tim2!

tim2 <- lme(fixed=EERTmn ~ male + female + 
              (male:time7c)          + (female:time7c)          +                           
              (male:IERT_Cp)         + (female:IERT_Cp)         + 
              (male:IERT_Cp_Partner) + (female:IERT_Cp_Partner)-1,

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

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

mA.f <-tim2

#control:

summary(mA.f)

#The model m0 containing only the other random effect not to be tested needs
#to be created first:

m0.f <- lme(fixed=EERTmn ~ male + female + 
              (male:time7c)          + (female:time7c)          +                           
              (male:IERT_Cp)         + (female:IERT_Cp)         + 
              (male:IERT_Cp_Partner) + (female:IERT_Cp_Partner)-1,

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

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

#control:

summary(m0.f)

#The model m containing only the effect to be tested also:

m.f <- lme(fixed=EERTmn ~ male + female + 
             (male:time7c)          + (female:time7c)          +                           
             (male:IERT_Cp)         + (female:IERT_Cp)         + 
             (male:IERT_Cp_Partner) + (female:IERT_Cp_Partner)-1,

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

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

#control:
summary(m.f)

#Now test the female random intercept for significance:

#basic pattern of the cmd: exactRLRT(m, mA, m0)

#inserted:

exactRLRT(m.f, mA.f, m0.f)

#Output:

#Error in exactRLRT(m.f, mA.f, m0.f) : 
#  Random effects not independent - covariance(s) set to 0 under H0.

#exactRLRT can only test a single variance.
#In addition: Warning message:
#  In extract.lmeDesign(m) : Removing incomplete cases from supplied data.

Большое спасибо заранее!

Бест, Патрик

...