BayesianWorkflow.jl

Welcome!

using BayesianWorkflow

@model function f()
    x ~ Normal()
    y ~ Normal(x)
end

bm = BWModel(f(), (; y = 1.0))

prior_predictive_check(bm, 1000)

post = posterior(bm, 1000)

posterior_predictive_check(bm, post)