从az.plot_ppc观察到的曲线与实际数据的kde图不匹配

编程问答 2022-06-08

问题描述

我正在建模我认为是伽马分布的数据，因此我编写了如下的PyMC3模型。

import pymc3 as pm
import arviz as az

with pm.Model() as model2:
    alpha = pm.Exponential('alpha',5)
    beta = pm.Exponential('beta',5)
    
    p_factor = pm.Gamma('p_factor',alpha = alpha,beta = beta,observed = d.productivity_factor)
    
    prior_checks = pm.sample_prior_predictive(samples=100,random_seed= seed)
    
    trace2 = pm.sample(2000,tune=2000)

然后我进行后验预测检查

with model2:
    ppc2 = pm.sample_posterior_predictive(trace2,var_names=["alpha","beta",'p_factor'],random_seed=seed)
    
idata = az.from_pymc3(trace2,posterior_predictive=ppc2)
az.plot_ppc(idata)

，但是观察到的曲线与我的预期不符。我在做错什么吗？

KDE plot of observed data

PPC Plot

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

arvizbayesianpymc3python