Here’s the plan for Stan core development that Bob presented at Stancon last week (that is, back at the end of August, 2018): Part I. Rear-View Mirror Stan 2.18 Released Multi-core Processing has Landed! Multi-Process Parallelism Map Function New Built-in Functions Manuals to HTML Improved Effective Sample Size Foreach Loops Data-qualified Arguments Bug Fixes and […]

# Category: Stan

## Stan This Month

So much is going on with Stan that it can be hard to keep track, so we (the Stan project) are starting a monthly update and newsletter. If you want to be included in the monthly mailing list, just type in your email here. Charles Margossian is the editor of Stan This Month and indeed […]

## Computational Bayesian Statistics [book review]

This Cambridge University Press book by M. Antónia Amaral Turkman, Carlos Daniel Paulino, and Peter Müller is an enlarged translation of a set of lecture notes in Portuguese. (Warning: I have known Peter Müller from his PhD years in Purdue University and cannot pretend to perfect objectivity. For one thing, Peter once brought me frozen-solid […]

## Principal Stratification on a Latent Variable (fitting a multilevel model using Stan)

Adam Sales points to this article with John Pane on principal stratification on a latent variable, and writes: Besides the fact that the paper uses Stan, and it’s about principal stratification, which you just blogged about, I thought you might like it because of its central methodological contribution. We had been trying to use computer […]

## Autodiff! (for the C++ jockeys in the audience)

Here’s a cool thread from Bob Carpenter on the Stan forums, “A new continuation-based autodiff by refactoring,” on the inner workings of Stan. Enjoy.

P.S. Just for laffs, see this post from 2010 which is the earliest mention of auto…

## Transforming parameters in a simple time-series model; debugging the Jacobian

So. This one is pretty simple. But the general idea could be useful to some of you. So here goes. We were fitting a model with an autocorrelation parameter, rho, which was constrained to be between 0 and 1. The model looks like this: eta_t ~ normal(rho*eta_{t-1}, sigma_res), for t = 2, 3, … T […]

## NYC Meetup Thursday: Under the hood: Stan’s library, language, and algorithms

I (Bob, not Andrew!) will be doing a meetup talk next Thursday in New York City. Here’s the link with registration and location and time details (summary: pizza unboxing at 6:30 pm in SoHo): Bayesian Data Analysis Meetup: Under the hood: Stan’s library, language, and algorithms After summarizing what Stan does, this talk will focus […]

The post NYC Meetup Thursday: Under the hood: Stan’s library, language, and algorithms appeared first on Statistical Modeling, Causal Inference, and Social Science.

## Reproducibility and Stan

Aki prepared these slides which cover a series of topics, starting with notebooks, open code, and reproducibility of code in R and Stan; then simulation-based calibration of algorithms; then model averaging and prediction. Lots to think about here: t…

## Reproducibility and Stan

Aki prepared these slides which cover a series of topics, starting with notebooks, open code, and reproducibility of code in R and Stan; then simulation-based calibration of algorithms; then model averaging and prediction. Lots to think about here: t…

## Prior distributions for covariance matrices

Someone sent me a question regarding the inverse-Wishart prior distribution for covariance matrix, as it is the default in some software he was using. Inverse-Wishart does not make sense for prior distribution; it has problems because the shape and scale are tangled. See this paper, “Visualizing Distributions of Covariance Matrices,” by Tomoki Tokuda, Ben Goodrich, […]

The post Prior distributions for covariance matrices appeared first on Statistical Modeling, Causal Inference, and Social Science.

## Prior distributions for covariance matrices

Someone sent me a question regarding the inverse-Wishart prior distribution for covariance matrix, as it is the default in some software he was using. Inverse-Wishart does not make sense for prior distribution; it has problems because the shape and scale are tangled. See this paper, “Visualizing Distributions of Covariance Matrices,” by Tomoki Tokuda, Ben Goodrich, […]

The post Prior distributions for covariance matrices appeared first on Statistical Modeling, Causal Inference, and Social Science.

## Multilevel models for multiple comparisons! Varying treatment effects!

Mark White writes: I have a question regarding using multilevel models for multiple comparisons, per your 2012 paper and many blog posts. I am in a situation where I do randomized experiments, and I have a lot of additional demographic information about people, as well. For the moment, let us just assume that all of […]

The post Multilevel models for multiple comparisons! Varying treatment effects! appeared first on Statistical Modeling, Causal Inference, and Social Science.

## Fitting the Besag, York, and Mollie spatial autoregression model with discrete data

Rudy Banerjee writes: I am trying to use the Besag, York & Mollie 1991 (BYM) model to study the sociology of crime and space/time plays a vital role. Since many of the variables and parameters are discrete in nature is it possible to develop a BYM that uses an Integer Auto-regressive (INAR) process instead of […]

The post Fitting the Besag, York, and Mollie spatial autoregression model with discrete data appeared first on Statistical Modeling, Causal Inference, and Social Science.

## Stan development in RStudio

Check this out! RStudio now has special features for Stan: – Improved, context-aware autocompletion for Stan files and chunks – A document outline, which allows for easy navigation between Stan code blocks – Inline diagnostics, which help to find issues while you develop your Stan model – The ability to interrupt Stan parallel workers launched […]

The post Stan development in RStudio appeared first on Statistical Modeling, Causal Inference, and Social Science.

## Stan on the web! (thanks to RStudio)

This is big news. Thanks to RStudio, you can now run Stan effortlessly on the web. So you can get started on Stan without any investment in set-up time, no need to install C++ on your computer, etc. As Ben Goodrich writes, “RStudio Cloud is particularly useful for Stan tutorials where a lot of time […]

The post Stan on the web! (thanks to RStudio) appeared first on Statistical Modeling, Causal Inference, and Social Science.

## Cool postdoc position in Arizona on forestry forecasting using tree ring models!

Margaret Evans sends in this cool job ad: Two-Year Post Doctoral Fellowship in Forest Ecological Forecasting, Data Assimilation A post-doctoral fellowship is available in the Laboratory of Tree-Ring Research (University of Arizona) to work on an NSF Macrosystems Biology-funded project assimilating together tree-ring and forest inventory data to analyze patterns and drivers of forest productivity […]

The post Cool postdoc position in Arizona on forestry forecasting using tree ring models! appeared first on Statistical Modeling, Causal Inference, and Social Science.

## The hot hand—in darts!

Roland Langrock writes: Since on your blog you’ve regularly been discussing hot hand literature – which we closely followed – I’m writing to share with you a new working paper we wrote on a potential hot hand pattern in professional darts. We use state-space models in which a continuous-valued latent “hotness” variable, modeled as an […]

The post The hot hand—in darts! appeared first on Statistical Modeling, Causal Inference, and Social Science.

## “Dynamically Rescaled Hamiltonian Monte Carlo for Bayesian Hierarchical Models”

Aki points us to this paper by Tore Selland Kleppe, which begins: Dynamically rescaled Hamiltonian Monte Carlo (DRHMC) is introduced as a computationally fast and easily implemented method for performing full Bayesian analysis in hierarchical statistical models. The method relies on introducing a modified parameterisation so that the re-parameterised target distribution has close to constant […]

The post “Dynamically Rescaled Hamiltonian Monte Carlo for Bayesian Hierarchical Models” appeared first on Statistical Modeling, Causal Inference, and Social Science.

## Against Winner-Take-All Attribution

This is the anti-Wolfram. I did not design or write the Stan language. I’m a user of Stan. Lots of people designed and wrote Stan, most notably Bob Carpenter (designed the language and implemented lots of the algorithms), Matt Hoffman (came up with the Nuts algorithm), and Daniel Lee (put together lots of the internals […]

The post Against Winner-Take-All Attribution appeared first on Statistical Modeling, Causal Inference, and Social Science.

## StanCon 2018 Helsinki tutorial videos online

StanCon 2018 Helsinki tutorial videos are now online at Stan YouTube channel List of tutorials at StanCon 2018 Helsinki Basics of Bayesian inference and Stan, parts 1 + 2, Jonah Gabry & Lauren Kennedy Hierarchical models, parts 1 + 2, Ben Goodrich Stan C++ development: Adding a new function to Stan, parts 1 + 2, […]

The post StanCon 2018 Helsinki tutorial videos online appeared first on Statistical Modeling, Causal Inference, and Social Science.