我正在更新一些计算,其中我使用pymc2到pymc3,当我在模型上有一些离散的随机变量时,我在采样器行为方面遇到了一些问题.例如,请考虑使用pymc2的以下模型:

import pymc as pm

N = 100

data = 10

p = pm.Beta('p', alpha=1.0, beta=1.0)

q = pm.Beta('q', alpha=1.0, beta=1.0)

A = pm.Binomial('A', N, p)

X = pm.Binomial('x', A, q, observed=True, value=data)

它并不代表任何东西,它只是一个模型,其中一个未观察到的变量是离散的.当我用pymc2对这个模型进行采样时,我得到以下结果:

mcmc = pm.MCMC(model)

mcmc.sample(iter=100000, burn=50000, thin=100)

plot(mcmc)

但是当我尝试使用PYMC3时,我得到了这个:

with pm.Model() as model:

N = 100

p = pm.Beta('p', alpha=1.0, beta=1.0)

q = pm.Beta('q', alpha=1.0, beta=1.0)

A = pm.Binomial('A', N, p)

X = pm.Binomial('x', A, q, observed=10)

with model:

start = pm.find_MAP()

with model:

step = pm.NUTS()

trace = pm.sample(3000, step, start)

pm.traceplot(trace)

看起来变量A根本没有被采样.我没有读过很多关于pymc3中使用的采样方法,但我注意到它似乎特别针对连续模型.这是否意味着它排除了模型上的离散未观察变量,或者是否有某种方法可以做我想做的事情?

python离散变量_python – 当涉及离散变量时,pymc3与pymc2的困难相关推荐

  1. python绘制折线图中文图例不显示_python使用matplotlib绘图时图例显示问题的解决...

    前言 matplotlib是基于Python语言的开源项目,旨在为Python提供一个数据绘图包.在使用Python matplotlib库绘制数据图时,需要使用图例标注数据类别,但是传参时,会出现图 ...

  2. python编程语言_python编程语言的简要历史

    python编程语言 The Python programming language is a general-purpose programming language that has carved ...

  3. Python可视化matplotlib自定义:运行时参数修改、自定义样式、默认样式+plt.style.use()

    Python可视化matplotlib自定义:运行时参数修改.自定义样式.默认样式+plt.style.use() 目录 Python可视化matplotlib自定义:运行时参数修改.自定义样式.默认 ...

  4. python画代码-Python教程_Python画Mandelbrot集 代码

    Python教程_Python画Mandelbrot集 代码 作者:Comet 来源: 课课家 www.kokojia.com点击数:278发布时间:2015-06-19 11:17:19 曼德勃罗集 ...

  5. python面试常见问题-Python面试常见问题,涉及Python各个方面

    原标题:Python面试常见问题,涉及Python各个方面 爬虫面试常见问题 一.项目问题: 1.你写爬虫的时候都遇到过什么反爬虫措施,你是怎样解决的 2.用的什么框架.为什么选择这个框架 二.框架问 ...

  6. python字符集_PYTHON 中的字符集

    Python中的字符编码是个老生常谈的话题,今天来梳理一下相关知识,希望给其他人些许帮助. Python2的 默认编码 是ASCII,不能识别中文字符,需要显式指定字符编码:Python3的 默认编码 ...

  7. python将EXCEL数据导入数据库时日期型数据变成数字并加.0的问题一行代码解决方案方案

    [问题描述]:python将EXCEL数据导入数据库时日期变成文本型数据并显示为数字格式 [解决方案] 数据源: codes: #!/usr/bin/python3 -- coding: utf-8 ...

  8. 在python中 函数赋值给变量时,需要注意的几个事项

    变量赋值是我们在日常开发中经常会遇到的一个问题,本文主要给大家介绍的是关于python将函数赋值给变量时需要注意的一些问题,分享出来供大家参考学习,下面话不多说了,来一起看看详细的介绍: 见过两种函数 ...

  9. python中将函数赋值给变量时需要注意的一些问题

    变量赋值是我们在日常开发中经常会遇到的一个问题,本文主要给大家介绍的是关于python将函数赋值给变量时需要注意的一些问题,分享出来供大家参考学习,下面话不多说了,来一起看看详细的介绍: 见过两种函数 ...

最新文章

  1. 空车上路,Waymo拿下加州首个“真”无人驾驶许可证
  2. leetcode算法第8题
  3. 排列与组合的一些定理(二)
  4. Tungsten Fabric SDN — Device Manager
  5. 第五章:管理数据库存储结构
  6. MATLAB格式化输出控制
  7. python用递归方式实现最大公约数_关于python最大递归深度 - 998
  8. Angular应用里的Template Reference变量
  9. 为长篇文档快速创建目录
  10. les物流执行系统_【精益运营】细化运营 精益求精 将“精准物流”落到实处...
  11. 云主机挂载硬盘 - 开机自动挂载
  12. asp.net 报表页面模板_Stimulsoft ASP.NET MVC报表教程:在设计器中保存报表模板
  13. Linux学习日志--共享内存
  14. 量化研究: Julia还是Matlab?
  15. cocos creator 文字描边效果
  16. paypal无法提现?最新解决办法(实战教程)!
  17. 神经网络训练集和测试集,神经网络验证集作用
  18. 《自卑与超越》的优秀读后感作文4000字
  19. 网传腾讯大规模裁员测试工程师,腾讯相关人员否认
  20. 屏蔽第三方SDK权限

热门文章

  1. PgSQL · 应用案例 · 逻辑订阅给业务架构带来了什么?
  2. elixir 关键字列表
  3. 下拉刷新和UITableView的section headerView冲突的原因分析与解决方案
  4. 微软宣布免费 Web 版 Office 2010发布日期
  5. 某公司PIX 520防火墙系统和NAT的实施
  6. 数据结构与算法题整理
  7. 使用线程池应该注意的问题
  8. hdu-2032杨辉三角
  9. LightOj_1030 Discovering Gold
  10. [linux驱动]proc学习笔记(一)