r - 蒙特卡洛预测分布

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)
}