我是PyMC3的新手。在学习PyMC3时,我从一个博客(可能是为作者运行的)中重新创建了一个示例程序,并遇到以下错误:AttributeError:无法腌制本地对象'Op.make_py_thunk..rval'。我完全被卡住,任何指导将不胜感激。相关代码如下:
log_dose = np.array([-.86, -.3, -.05, .73])
log_dose_shared = shared(log_dose)
n = 5 * np.ones(4, dtype = int)
n_shared = shared(n)
deaths = np.array([0, 1, 3, 5])
with Model() as bioassay_model:
# Logit model parameters
alpha = Normal('alpha', 0, sd = 100)
beta = Normal('beta', 0, sd = 100)
# Calculate probabilities of death
theta = invlogit(alpha + beta * log_dose_shared)
# Data likelihood
obs_death = Binomial('obs_death', n = n_shared, p = theta, observed = deaths)
with bioassay_model:
# Obtain starting values via MAP
start = find_MAP(model = bioassay_model)
# Instantiate sampler
step = pm.Metropolis()
# Draw 2000 posterior samples
bioassay_trace = sample(50000, step = step, start = start)
logp = -13.034,|| grad || = 0.00043389:100%|████████████████████████████████████████████| 14/14 [00:00 <00:00,398.61it / s]
多过程采样(4个工作中的4个链)
CompoundStep
大都会:β
大都会:α
_RemoteTraceback追溯(最近一次通话最后一次)
_RemoteTraceback:
追溯(最近一次通话):
_feed中的文件“ C:\ Users \ bikim \ AppData \ Local \ conda \ conda \ envs \ pymc3p36 \ lib \ site-packages \ joblib \ externals \ loky \ backend \ queues.py”,第151行
obj,reducers = reducers)
在转储中的文件“ C:\ Users \ bikim \ AppData \ Local \ conda \ conda \ envs \ pymc3p36 \ lib \ site-packages \ joblib \ externals \ loky \ backend \ reduction.py”,第145行
p.dump(obj)
AttributeError:无法腌制本地对象'Op.make_py_thunk..rval'
为了完整起见,这是我的环境:
我使用conda安装了PyMC3
3.6.5 | Anaconda自定义(64位)| (默认值,2018年3月29日,13:32:41)[MSC v.1900 64位(AMD64)]
在PyMC3版本3.4.1上运行
Theano版本:1.0.2
Windows 10