如何在PyMC3中定义四参数Beta先验分布

PyMC3原生的pm.Beta默认仅输出[0,1]区间的标准Beta分布，你期望的直接按位置传四个参数的写法不生效，有两种可行的实现方案： #

方案1：高版本PyMC3直接指定内置参数 #

PyMC3 3.9及以上版本已经为连续分布内置了loc（位置移位）和scale（缩放）参数，你只需要显式指定参数名传入即可，避免位置传参匹配错误：

with pm.Model() as model_g:
    n = pm.Beta('n', alpha=28.78, beta=0.98, loc=-0.22, scale=0.32)

方案2：全版本兼容的手动变换写法 #

如果你的PyMC3版本较低，或者需要更可控的变换逻辑，可以先定义标准Beta分布，再通过线性变换得到目标四参数Beta分布：

with pm.Model() as model_g:
    # 定义[0,1]区间的标准Beta分布
    std_n = pm.Beta('std_n', alpha=28.78, beta=0.98)
    # 线性变换得到四参数Beta，取值范围为[loc, loc+scale] = [-0.22, 0.1]
    n = pm.Deterministic('n', -0.22 + 0.32 * std_n)

注意：四参数Beta的有效取值范围是loc到loc+scale，请确认该范围匹配你要定义的参数的实际取值逻辑。

内容的提问来源于stack exchange，提问作者H2H