Monte Carlo 的预测分布

Predictive distribution with Monte Carlo

我想获得每个预测点的分布,即因变量 y 的每个值。在因变量 y 中,我有 199 个观察值,对于每个观察值,我想要一个 distribution.So 我的输出应该是一个 199*100 的矩阵。我怎样才能得到它?我每次都尝试了 T=100.But 的 Monte Carlo 模拟,它是一个 199 的向量而不是矩阵 199*100.R-Code is here

你可以这样做

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

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