set.seed(123)
n <- 60
data <- data.frame(
Y = c(rnorm(n - 10), rep(NA, 10)), # 50 observed + 10 missing
X1 = rnorm(n),
X2 = rnorm(n)
)
res <- DfiMI(data, R = 3, M = 5)
head(res$Yhat) # inspect imputed Y
res$betahat # inspect coefficients
Run the code above in your browser using DataLab