HI Nicole, Homework? No, it is a part of analysis for my paper. Any idea? ya
At 2012-03-30 23:32:21,"Nicole Marie Ford" <nmf...@uwm.edu> wrote: >ya, > >is this homework? > > >----- Original Message ----- >From: "ya" <xinxi...@163.com> >To: r-help@r-project.org >Sent: Friday, March 30, 2012 9:04:59 AM >Subject: [R] pooling in MICE > >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. [[alternative HTML version deleted]] ______________________________________________ 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.
