未来预测值均值的区间估计

Question

我在RStan中有一个用来估计模型(有点复杂的模型)的贝叶斯代码。在估计了模型之后，我从后验分布中提取了500个参数集样本，以模拟未来(下一个月)的模型数据。最后，我取预测值的平均值(每个时间点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)

现在我找到了mean(X1)；mean(X2);....mean(X30)，并将它们与时间作图。我想找出这些方法的间隔。

Answer 1

每个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)