微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何计算2D样本的MCMC Geweke准则收敛性?

如何解决如何计算2D样本的MCMC Geweke准则收敛性?

所以我正在使用pymc3来为MCMC链计算Geweke和Gelmen-Rubin诊断,例如:

chain = sampler(data,samples=5000,mu_init=9,sigma=sigma,proposal_width=2,mu_prior_mu=mu_prior,mu_prior_sd = sigma_prior)

gw_plot = pm3.geweke(chain)

这在1D模式中很好,但是现在我想在2D模式中做到这一点,但是我的链长5000 x 2,我不确定是否应该输出一个或两个数字?即我是分别进行每个维度还是希望可以得到一个摘要诊断?但是目前,我得到2 x 5000 x 2的形状作为gw_plot,我不确定该怎么做?谢谢

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。