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





Вы могли бы сделать что-то вроде
mc_pred <- matrix(nrow = 199, ncol = 100)
for (i in 1:100) {
mc_pred[, i] <- predict(model, x_test)
}