我有一个贝叶斯代码用于估算RStan中的模型（有点复杂的模型）。在估计模型后，我从后验分布中获取了500个参数集样本，以模拟未来（接下来的1个月）模型中的数据。最后，我采用预测值的方法（每个时间点500个预测值），然后将它们与实际观察结果（用图表）进行比较。
我的问题是 - 如何计算预测值的这些均值的区间？
示例：从后验分布中绘制参数样本后，我使用这些参数模拟了模型中的变量X1，X2，.... X30：
X1= (33,25,10,19,25)
X2= (11,10,15,13.5,17)
.......
X30= (40,33.3,50,29,45.1)
现在我找到了意思（X1）; mean（X2）; .... mean（X30）并根据时间绘制它们。我想找到这些方法的间隔。
答案 0 :(得分：3)
每个X1 ... X30的示例后抽取：
m <- matrix(data = rep(rnorm(500, 0, 1), 30), nrow = 500, ncol = 30, byrow = FALSE)
colnames(m) <- paste0("X", 1:30)
获得这些抽奖的可靠间隔：
ci <- apply(m, 2, quantile, probs = c(0.1, 0.9))
colnames(ci) <- paste0("X", 1:30)