Laplace approximation for mixed models in Julia (TMB-like functionality)

The quap in StatisticalRethinking.jl appears to just fit a normal distribution to a pre-sampled posterior rather than marginalizing out variables from an arbitrary likelihood function.

1 Like