之前的贝叶斯击球率

机器算法验证 贝叶斯 事先的
2022-01-29 10:01:04

我想问一个问题,这个问题的灵感来自对有关 beta 分布直觉的查询的出色回答。我想更好地理解击球率的先验分布的推导。看起来大卫正在从均值和范围中退出参数。

假设均值为且标准差为 ,您能否通过求解以下两个方程 来回退0.270.18αβ

αα+β=0.27αβ(α+β)2(α+β+1)=0.182

2个回答

请注意:

αβ(α+β)2=(αα+β)(1αα+β)

这意味着方差因此可以用均值表示为

σ2=μ(1μ)α+β+1

如果您想要的标准差(方差),只需计算:.27.18.0324

α+β=μ(1μ)σ21=.27(1.27).03241=5.083333

现在您知道总数了,很简单:αβ

α=μ(α+β)=.275.083333=1.372499β=(1μ)(α+β)=(1.27)5.083333=3.710831

您可以在 R 中检查此答案:

> mean(rbeta(10000000, 1.372499, 3.710831))
[1] 0.2700334
> var(rbeta(10000000, 1.372499, 3.710831))
[1] 0.03241907

我想将此添加为对出色答案的评论,但它运行时间很长,并且使用答案格式看起来会更好。

需要记住的是,并非所有都是可能的。很明显,但不太清楚的限制。(μ,σ2)μ[0,1]σ2

使用与大卫相同的推理,我们可以表达

σ2(α,μ)=μ2(1μ)α+μ

这相对于,最大的可以是:ασ2μ

limα0σ2(α,μ)=μ(1μ)

这只是一个上定论,因为有效的集是开放的(即,对于 Beta,我们必须有);此限制本身在处最大化。αα>0μ=12

注意与相应伯努利 RV 的关系。的 Beta 分布,因为它被迫取 0 和 1 之间的所有值,因此必须比具有相同均值的伯努利 RV(其所有质量都在末端)更分散(即具有更低的方差)区间)。事实上,将发送到 0 并修复相当于将 PDF 的质量越来越接近 0 和 1,即越来越接近一个伯努利分布,这就是为什么方差的上项恰好是对应的伯努利方差。μαβ=1μμα

综上所述,这是 Beta 的一组有效均值和方差:

在此处输入图像描述

(事实上​​,这在Beta 的 Wikipedia 页面上有说明)