Displaying 1 result from an estimated 1 matches for "erq2tnc".
2007 Mar 02
1
Mitools and lmer
...t; #Read data
> data.dir<-system.file("dta",package="mitools")
> files.imp<-imputationList(lapply(list.files(data.dir,
> pattern="imp.\\.dta", full=TRUE), read.dta))
>
> #estimate model over each imputed dataset
> model0<-with(files.imp,lmer( erq2tnc ~1+trt2+nash+wash+male+coh2+coh3+(1 |
> sitebeth)))
> #extract betas and standard errors
> betas<-MIextract(model0,fun=coef)
> vars<-MIextract(model0,fun=vcov)
> #Combine the results
> summary(MIcombine(betas,vars))
>
> Error in cbar + results[[i]] : non-numeric argume...