Background

Log normal

f = brm(
  bf(GPT~1), data = d, family = lognormal(),
  chains = 4, cores = 4
)
## Compiling the C++ model
## Start sampling
pp_check(f, type="dens_overlay", nsamples = 100)