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
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
