import pymc as pm, theano.tensor as T
with Model() as success_model:
    loss_lambda_rate = pm.Flat('loss_lambda_rate')
error = Poisson('error', mu=totalCounts * T.exp(loss_lambda_rate), observed
    =successCounts)
pymc3: несколько наблюдаемых значений