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)
)
result <- DERLS(data, rho = 0.01, lambda = 0.95, R = 3, nb = 50)
head(result$Yhat) # inspect imputed Y
result$betahat # inspect estimated coefficients
Run the code above in your browser using DataLab