Неправильная экспоненциальная выборка в PyMC ⇐ Python

[Anonymous]

Очевидно, я делаю что-то не так... Пожалуйста, взгляните на следующую программу. Он работает хорошо, но дает мне параметр лямбда для экспоненциального распределения, который далек от параметра, который я использовал для генерации случайных наблюдений:

Код: Выделить всё

import numpy as np
import arviz as az
import pymc as pm

lambda_param = 0.25
random_size = 1000
x = np.random.exponential(lambda_param, random_size)

basic_model = pm.Model()

with basic_model:
_lam_ = pm.HalfNormal("lambda", sigma = 1)
Y_obs = pm.Exponential("Y_obs", lam = _lam_, observed = x)
start = pm.find_MAP(model = basic_model)
idata = pm.sample(1000, start = start)

summary = az.summary(idata, round_to = 6)
summary

После моего последнего запуска программы я обнаружил, что среднее значение лямбда превышает 4..., где лямбда=0,25 как Я использовал его.
Буду очень признателен, если укажу пальцем на мои ошибки в программировании.

Подробнее здесь: https://stackoverflow.com/questions/754 ... ng-in-pymc