标签： r
我在R中使用MCMCpack为我的数据提供“贝叶斯”线性回归模型（在使用PostProbMod函数比较一组可能模型的贝叶斯因子之后）。假设我选择的模型是Y = X3 + X7 + X9。
现在，我希望能够计算新观测值的估计Y值的概率分布。也就是说，我想将X3，X7和X9的新值输入到模型中，并根据每组自变量查看模型预测的概率分布图。
我是否可以利用MCMCpack提供的有关X3，X7和X9的回归系数分位数的信息来制作这些图？
...谢谢