Прогнозирующее распределение с Монте-Карло

Question

Я хочу получить распределение по каждой точке прогнозирования, то есть для каждого значения зависимой переменной y. В зависимой переменной y у меня есть 199 наблюдений, и для каждого наблюдения я хочу распределение. Так что мой вывод должен быть матрицей 199 * 100. Как я могу получить это? Я попробовал симуляцию Монте-Карло для T = 100. Но каждый раз, это вектор 199 вместо матрицы 199 * 100. R-код здесь

Answer 1

Вы могли бы сделать что-то вроде

mc_pred <- matrix(nrow = 199, ncol = 100)

for (i in 1:100) {
    mc_pred[, i] <- predict(model, x_test)
}