Hi folks, I am using r to run multiple imputation on propensity score matching, followed by my outcome regression model which is a HLM model. I read a thread: https://stats.stackexchange.com/questions/257789/manipulating-data-for-propensity-score-matching-following-multiple-imputation-wi that talks about the use of package mira to get the average regression coefficients from the imputed datasets but the package doesn't work in my occasion. Anyone knows how I could get the pooled regression estimate if the outcome regression is not a simple regression but a complicated one such as HLM? Thanks.
Hi Ted, the package mira doesn't work because mira cannot handle the pooled estimation of regression coefficients from HLM model. The package mira works fine only with simple regression models. In order to use mira, the coef() and vcov() must be available but in my case, they are not.
- Badges
- Science topic
- Similar topics
- Mathematics
- Statistics
- Statistical Software
- R Statistical Package