imp <- mice(nhanes2, m = 2, print = FALSE, seed = 14221)
# descriptive statistics
getfit(with(imp, table(hyp, age)))
# model fitting and testing
fit1 <- with(imp, lm(bmi ~ age + hyp + chl))
fit2 <- with(imp, glm(hyp ~ age + chl, family = binomial))
fit3 <- with(imp, anova(lm(bmi ~ age + chl)))
Run the code above in your browser using DataLab