如何对回归系数的后验分布进行采样？

【问题标题】：How do i sample the posterior distribution for the regression coefficients?如何对回归系数的后验分布进行采样？
【发布时间】：2020-11-25 20:33:58
【问题描述】：

我想弄清楚我在这个问题上哪里出错了？

library(MCMCpack)
library(mcmcplots)
X=c(66, 51, 60, 48, 48, 34, 48, 46, 51, 58, 51, 62, 43, 42, 44, 50, 57, 52, 40, 42)
Y=c(75.8, 56.8, 68.0, 55.4, 56.4, 44.2, 57.4, 53.8, 63.8, 67.4, 58.8, 66.6, 55.4, 49.6, 59.2, 58.0, 64.6, 58.6, 54.0, 57.6)
data<-list(X,Y)
post= MCMCregress(dist~speed, data=data)
summary(post)
plot(post, col=c("red"))
denplot(post)

我不断收到错误消息-

eval 中的错误（predvars、data、env）：找不到对象“速度”

提前谢谢你

【问题讨论】：

标签： r mcmc

【解决方案1】：

当我调用您的以下代码部分时，我收到了您发布的错误：

post= MCMCregress(dist~speed, data=data)

这是因为在 data 中，您已经指定了名称 X 和 Y。当我将其更改为以下内容时，我可以成功运行该行：

post= MCMCregress(Y~X, data=data)

然后我必须更改下一行，因为在使用连接时出现错误。我把它改成这样，它就可以工作了：

plot(post, col="red")

最后一行对我有用，没有错误。

【讨论】：

谢谢，我明白了。最好的问候