step = pm.CompoundStep([pm.Metropolis([total_lambda_mu]), pm.Metropolis([
    total_lambda_tau]), pm.Metropolis([total_lambda]), pm.Metropolis([
    loss_lambda_factor])])
pymc3: несколько наблюдаемых значений