[R] pooling in MICE

ya xinxi813 at 163.com
Fri Mar 30 16:04:59 CEST 2012


Hi everyone,

Does anyone here has experience using MICE to impute missing value? I am 
having problem to pool the imputed dataset for a MANOVA test, could you 
give me some advice please?
Here is my code:
 > library(mice)
 > 
grd3dat=subset(paper2,Control==1&Grade_1==3,select=c(Boy,BVCategoryT1,reclassify,meanseap_1,meandep_1,meananx_1,meanper_1,meanrelma_1,meandiff_1,meanascon_1,meanclasc_1,meanschoc_1))
 > grd3imp=mice(grd3dat,m=20,seed=123,print=F)
 > grd3mod=with(grd3imp,lm(cbind(meanseap_1,meandep_1,meananx_1)~Boy*BVCategoryT1))
 > summary(pool(grd3mod))
Error in table[, 8] <- x$nmis[names(x$qbar)] :
   replacement has length zero

 > pool(grd3mod)
Call: pool(object = grd3mod)
Pooled coefficients:
  [1]  2.51538922  0.08041626 -0.02266421  0.33964482 -0.38691307 -0.28778570 -0.10405965  0.12206124  0.86739342 -0.17195284  0.04784570 -0.38342776
[13] -0.07785758  0.27914375  0.09205530 -0.09185746  2.04934495 -0.34833583 -0.02972867 -0.65564544 -0.33033151  0.29788500  0.11309193  0.32081963
Fraction of information about the coefficients missing due to nonresponse:
  [1] 0.2113750 0.2137368 0.3849378 0.1885110 0.2056638 0.4089043 0.1960244 0.1116587 0.2015711 0.1960397 0.3745362 0.2108851 0.2545268 0.4241958 0.2146440
[16] 0.2520252 0.2380274 0.2131245 0.3006583 0.2464491 0.2915607 0.2173454 0.1970864 0.1422042
Thank you very much.

YA



More information about the R-help mailing list