Displaying 1 result from an estimated 1 matches for "grd3mod".
2012 Mar 30
3
pooling in MICE
...vice 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....