Nothing
## ----step0-----------------------------------------------------
options(width = 65)
suppressMessages(library(mi))
data(nlsyV, package = "mi")
## ----step1-----------------------------------------------------
mdf <- missing_data.frame(nlsyV)
## ----step1.5---------------------------------------------------
show(mdf) # momrace is guessed to be ordered
## ---- step2----------------------------------------------------
mdf <- change(mdf, y = c("income", "momrace"), what = "type",
to = c("non", "un"))
show(mdf)
## ---- step3----------------------------------------------------
summary(mdf)
image(mdf)
hist(mdf)
## ---- step4----------------------------------------------------
rm(nlsyV) # good to remove large unnecessary objects to save RAM
options(mc.cores = 2)
imputations <- mi(mdf, n.iter = 30, n.chains = 4, max.minutes = 20)
show(imputations)
## ---- step5A---------------------------------------------------
round(mipply(imputations, mean, to.matrix = TRUE), 3)
Rhats(imputations)
## ---- step5B---------------------------------------------------
imputations <- mi(imputations, n.iter = 5)
## ---- step6----------------------------------------------------
plot(imputations)
plot(imputations, y = c("ppvtr.36", "momrace"))
hist(imputations)
image(imputations)
summary(imputations)
## ---- step7----------------------------------------------------
analysis <- pool(ppvtr.36 ~ first + b.marr + income + momage + momed + momrace,
data = imputations, m = 5)
display(analysis)
## ---- step8----------------------------------------------------
dfs <- complete(imputations, m = 2)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.