[R] MICE Package and LMER

KDT dkadengye at gmail.com
Mon Jun 14 16:58:08 CEST 2010


Hi R users,

I am estimating a multilevel model using lmer.  My dataset has missing
values and I am using MICE package to make Multiple Imputations. 
Everything works good until i reach the POOLING stage using the pool()
function. I am able to get a summary of the pooled fixed effects but not the
random effects. No errors or warnings are given to me.
I checked the "help" file in R and the developers of MICE noted that "The
function pools also estimates obtained with lme() and lmer(), BUT only the
fixed part of the model."

Does anyone have any ideas on how I can get a summary of pooled random
effects? 

Below is my code
imp<-mice(mydata,m=3, imputationMethod =c("","","","logreg"),maxit=2, pri=F) 
model <- with(data=imp, lmer(miss~ sex + age + (1|id) + (1|sch),
family=binomial(link="logit")))
result<-pool(model)
summary(result)

Thanks
Trevor
-- 
View this message in context: http://r.789695.n4.nabble.com/MICE-Package-and-LMER-tp2254504p2254504.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list