We discussed this in Stan gathering, and @charlesm93 suggested sharing here. 1D normal prior and Poisson model. Integrating out the latent theta works with integrate_1d(), but doesn’t work with laplace_marginal(). The likely reason is a loss of numerical accuracy in some part. The Stan code includ…

Embedded Laplace numerical problem

WardBrian June 12, 2025, 9:15pm 2

Thanks Aki, I was able to replicate on my machine. Also copying @stevebronder

Topic		Replies	Views
Algebraic sovler problems (Laplace approximation in Stan) Modeling	65	7070	April 13, 2018
Integrated Laplace approximation roadmap Developers	49	5492	January 13, 2025
Problem with laplace coding using brms package Modeling brms	1	130	March 4, 2025
Defining a custom lpdf for the Normal-Laplace convoluted density in Stan General	15	1117	January 24, 2020
CmdStan & Stan 2.37 release candidate General	20	1493	October 7, 2025