inst/examples/R/data/missing/norm/example_amelia.R

In sigbertklinke/shinyExample: Create and run simple (interactive) programs as Shiny apps

# run example_mar.R before
library("Amelia")
library("mitools")
# run NORM
aobj <- amelia(xmar, noms=c("chas","rad"))
# compute linear regressions
models <- lapply(aobj$imputations, function(x) {
  lm(medv~lstat, data=x) 
})
# look at one model
summary(models[[1]])
# extract 
beta <- MIextract(models, fun="coef")
vcov <- MIextract(models, fun="vcov")
summary(MIcombine(beta, vcov))
summary(MIcombine(models))