# Example data
set.seed(123)
n <- 1000 # Number of rows
p <- 5 # Number of independent variables
data <- data.frame(Y = rnorm(n), X1 = rnorm(n), X2 = rnorm(n))
data[sample(n, 100), 1] <- NA # Introduce missing response values
# Call fiMI function
result <- fiMI(data, R = 10, n = n, M = 20)
# View results
print(result$Yhat) # Predicted response values
Run the code above in your browser using DataLab