Hello, I do not think I have seen this solved but new to Stan and needed a confirmation: For a normal hierarchical case say Y ~ N(theta, sigma^2) with theta unknown and sigma^2 unknown. We know that the sample mean or Y-bar ~ N(theta, sigma^2 /n) where n is the sample size and (n-1)*S^2 / sigm…

I’ve re-explored this and can confirm that you can use sufficient statistics for the target += command. Below script 1 uses the normal likelihood and script 2 uses the normal sufficient statistics. Hope this can clear any doubts and help others better understand the target += command. By running bo…

Thank you Mike! This helped. Do you know why I can not use target += as coded now? The model runs. Could also just not be a good model, but convergence is the issue and I think it’s because of the target +=.

Sufficient Statistics for Normal Case

Modeling

mike-lawrence June 24, 2021, 1:09pm 2

See this post. And maybe this one too.

Topic		Replies	Views
Normal likelihood via sufficient statistics Modeling	5	265	June 17, 2024
PSA: using stantargets for SBC(-esque) checks, & big speedups for big gaussian likelihoods Publicity techniques , specification , performance	1	1008	May 30, 2021
Fitting a simple normal model in Stan conditional on sufficient statistics & Jeffreys prior? Modeling techniques	10	611	April 30, 2024
PSA: Sufficient-Stats for Gaussian Likelihoods should speed up most Hierarchical models with IID observation subgroups Modeling techniques , specification , psa	3	1167	July 18, 2023
Normal_sufficient signature Developers	7	1074	March 30, 2017

Sufficient Statistics for Normal Case

Related topics