蒙特卡罗的预测分布

Question

我想获得每个预测点的分布，即因变量y的每个值。在因变量y中，我有199个观测值，对于每个观测值，我想要一个分布。所以我的输出应该是199 * 100的矩阵。我怎么能得到它？我尝试了T = 100的蒙特卡罗模拟。但每次，它的矢量为199而不是矩阵199 * 100. R-Code is here

Answer 1

您可以执行类似

的操作

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

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