在 pymc3 中定义一个 4 参数 Beta 先验答案

【问题标题】：Define a 4-parameters Beta prior in pymc3在 pymc3 中定义一个 4 参数 Beta 先验
【发布时间】：2021-11-24 23:11:58
【问题描述】：

我想使用 pymc3 进行贝叶斯分析。我的一个参数具有 Beta 分布，a=28.78, b=0.98, loc=-0.22, scale= 0.32。有谁知道如何在 pymc3 模型中定义 4 参数 Beta 分布？比如：

with pm.Model() as model_g:
    n=pm.Beta(‘n’, 28.78, 0.98, -0.22, 0.32)

【问题讨论】：

尝试使用 docs.scipy.org/doc/scipy/reference/generated/… 中的定义并仅将 a 和 b 传递给 pymc，我认为这是唯一受支持的

标签： python bayesian pymc3 beta-distribution

【解决方案1】：

这是一个确定性变换：

import pymc3 as pm

theta_a = 28.78
theta_b = 0.98
theta_loc = -0.22
theta_scale = 0.32

with pm.Model() as model_g:
  n_raw = pm.Beta('n_raw', theta_a, theta_b)
  n = pm.Deterministic('n', theta_scale * n_raw + theta_loc)

【讨论】：