library(Amelia)
amelia(a3,idvars="tr",m=5,noms="d30",ords=c("ai","gcm","gce","gcv"))->a3.impu
#plot(a3.impu,overimpute=T)--OK
a3.imput<-a3.impu$imputations[1:5]
library(MatchIt)
cem("tr",datalist=a3.imput,drop="bac",data=a3)->imp.match
imp.match
##z.o<-zelig(gcsmottak~ais+as.numeric(bac)+as.numeric(age),data=list(en),model="ls")..........cannot
match using zelig
out <- att(imp.match, gcso ~ as.numeric(bac)+as.numeric(age), data=a3)
##results in the following:
> out
Linear regression model on CEM matched data:
Error in print.cem.att(list(mult = list(list(att.model = c(14.5304012233110, :
subscript out of bounds
Any ideas?
Also-Is there an easy way to extract the matched cases from imp.match?
Matchit has the match.data option, but as far as I understand I cannot
use imputed data as input to matchit.
Regards,
M
-
cem Mailing List, served by Harvard-MIT Data Center
Send messages: cem(a)lists.gking.harvard.edu
[un]subscribe Options: http://lists.gking.harvard.edu/?info=cem
More information on cem: http://gking.harvard.edu/cem