Rstanarm(optimize = TRUE) diagnostics

[bgoodri]

@avehtari added a few to 2.19

New releases of rstanarm and rstantools R packages

Pareto diagnostic and Pareto smoothed importance sampling https://arxiv.org/abs/1507.02646 is used to diagnose and improve inference when using stan_glm with algorithm=‘optimizing’. Demonstration of timings with n=100,000 and p=100, and Gaussian and logistic regression at https://avehtari.github.io/RAOS-Examples/BigData/bigdata.html Using PSIS and importance resampling means also that we can use PSIS-LOO also when algorithm=‘optimizing’ (I recommend to increase the number of draws from the def…