Question

我的问题是关于以下 YouTube 视频 https://www.youtube.com/watch?v=0F0QoMCSKJ4 （贝叶斯统计介绍，第 1 部分：StataCorp LLC 的基本概念）在上面的视频中 4.22 分钟后验 = Beta(34,36) P(0.4 < θ < 0.499) = 0.434 还可以从视频中找到此计算的附加屏幕截图。当我尝试复制上述内容时 a) 在 R

pbeta(0.499,34,36)-pbeta(0.4,34,36) 
= 0.5140403

b) 在 Excel 2013 中

=BETA.DIST(0.499,34,36,1) - BETA.DIST(0.4,34,36,1) 
= 0.588551 - 0.074511
= 0.51404

我想知道为什么我得到的结果与视频中显示的结果不同。提前致谢，

Stata area calculation

Answer 1

所以你注意到这里有些不对劲是非常正确的。虽然他们为 alpha / (alpha + beta) = 34/(34+36) = 0.4875 的后验分布提供的平均值是正确的，但他们错误地在 (0.4,0.499) 区间内取概率时没有更新先验。也就是说，他们计算的 0.434 值基于 beta(30,30) 分布。在 stata 中，累积 Beta 分布的特征为 ibeta(alpha,beta,theta)。这就是他们计算 P(0.4 < theta < 0.499) = 0.434 的方法：

disp ibeta(30,30,0.499) - ibeta(30,30,0.4)
.43429354

有关 Beta 发行版背后的直觉的精彩对话，请参阅问题 here。

使用 Stata 进行贝叶斯统计

1 个答案: