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

Question

我想找出我在这道题中哪里出错了？

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) 出错：未找到对象 'speed'

提前致谢

Answer 1

我在调用您的代码的以下部分时收到您发布的错误消息：

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

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

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

然后我必须更改下一行，因为在使用连接时出现错误。我把它改成这个并且它有效：

plot(post, col="red")

最后一行对我来说没有错误。

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

How do i sample the posterior distribution for the regression coefficients?

r

mcmc