Two case studies:

Splines in Stan, by Milad Kharratzadeh.

Spatial Models in Stan: Intrinsic Auto-Regressive Models for Areal Data, by Mitzi Morris.

This is great. Thanks, Mitzi! Thanks, Milad!

Posted by Andrew on 10 October 2017, 3:23 pm

Two case studies:

Splines in Stan, by Milad Kharratzadeh.

Spatial Models in Stan: Intrinsic Auto-Regressive Models for Areal Data, by Mitzi Morris.

This is great. Thanks, Mitzi! Thanks, Milad!

## Recent Comments

- Chris on It’s not matching
*or*regression, it’s matching*and*regression. - Andre on Are female scientists worse mentors? This study pretends to know
- Ryan Bernstein on Further formalization of the “multiverse” idea in statistical modeling
- Max on Greek statistician is in trouble for . . . telling the truth!
- gec on Further formalization of the “multiverse” idea in statistical modeling
- Martha (Smith) on Greek statistician is in trouble for . . . telling the truth!
- Ryan Bernstein on Further formalization of the “multiverse” idea in statistical modeling
- Bob Carpenter on Further formalization of the “multiverse” idea in statistical modeling
- gec on Greek statistician is in trouble for . . . telling the truth!
- Anoneuoid on Further formalization of the “multiverse” idea in statistical modeling
- Further formalization of the “multiverse” idea in statistical modeling « Statistical Modeling, Causal Inference, and Social Science on It’s not matching
*or*regression, it’s matching*and*regression. - Andrew on Are female scientists worse mentors? This study pretends to know
- anon on Follow-up on yesterday’s posts: some maps are less misleading than others.
- Dan Weeks on Are female scientists worse mentors? This study pretends to know
- Keith O'Rourke on Greek statistician is in trouble for . . . telling the truth!
- Jukka on Today in spam
- Eric de Souza on Greek statistician is in trouble for . . . telling the truth!
- Xi'an on Greek statistician is in trouble for . . . telling the truth!
- Anonymous on Are female scientists worse mentors? This study pretends to know
- Brendan on Are female scientists worse mentors? This study pretends to know

- Chris on It’s not matching
## Categories

As a spatial statistician, this puts me one step closer to using Stan in my MS-level Bayes course, so kudos to you! (The bigger step for is me deciding whether or not I want to teach Gibbs/Metropolis and then go to software that doesn’t use Gibbs/Metropolis, but that’s my problem not yours…)

One minor comment: In the ICAR section of Mitzi’s write-up, the joint distribution is not a normal distribution (due to D-W being singular, as you point out a few lines later), and the pairwise difference formulation is off slightly — in particular, the “1/tau^{n/2}” piece outside of the exponent should have “(n-1)/2” if you only have one “spatial island” or, more generally, “(n-I)/2” if you have I spatial islands. Similarly, if you have I spatial islands, you need to have I sum-to-zero constraints (e.g., one for the contiguous US, one for Alaska, and one for Hawaii). Jim Hodges has a 2003 paper related to this in Biometrics (“On the precision of the conditionally autoregressive prior in spatial models”).

thank you for noticing that! and thanks for pointer to the Hodges paper which is admirably clear – available at: https://pdfs.semanticscholar.org/db52/43bcee9b8148695167460a51e513a27df368.pdf

There is a neat proposal on how to treat disconnected graphs in Freni-Sterrantino et al. (2017) on ArXiv: “A note on intrinsic Conditional Autoregressive models for disconnected graphs” https://arxiv.org/pdf/1705.04854.pdf

On the splines example, the math near the end shows a_1 ~ N(0,tau) and a_i ~ N(a_i-1,tau) . I dont think this is right, as it seems to assume the same variance should be used for the starting value of the spline as well as the first derivative. Also, if we are interested in smoothing, I would be tempted to smooth the second derivative rather than the first (ie curvature rather than slope), but maybe it is an empirical question as to which is more appropriate. Finally, it would be useful to link this to the established literature on p-splines (ie penalized B-splines), as it seems that is what is going on here. These are all just suggestions (and may well be wrong) — thanks for the great case study!

PS — sorry to put the comment here, there doesnt appear to be a comment or contact link on the case studies themselves. Maybe this would be useful?