96d0eebd5f75592deb850b8af84af2af9a8ae677,pymc3/examples/custom_dists.py,,,#,101

Before Change

    plot_MCMC_trace(ax[0], xdata, ydata, trace, True, colors=colors)
    plot_MCMC_model(ax[1], xdata, ydata, trace)

pymc3_trace = [trace["intercept"],
               trace["slope"],
               trace["sigma"]]

After Change

    return -tt.log(tt.abs_(value))

with pm.Model() as model:
    alpha = pm.Normal("intercept", mu=0, sd=100)
    // Create custom densities
    beta = pm.DensityDist("slope", loglike1, testval=0)
    sigma = pm.DensityDist("sigma", loglike2, testval=1)
    // Create likelihood
    like = pm.Normal("y_est", mu=alpha + beta *
                        xdata, sd=sigma, observed=ydata)

    trace = pm.sample(2000, njobs=2)

Instances