Kfold returning NaN

[bgoodri]

Models with a standardized outcome don’t make sense in a Bayesian context and standardized predictors tend to confuse rstanarm. The best thing is to estimate models with the raw variables and then

Multiple outcomes in rstanarm

The lmer results are suspect because the correlations are 1 or NaN. The stan_lmer ones might be okay but you get a very small intercept (I think) due to standardizing your variables, which isn’t recommended. Just do it with unstandardized variables, and if you have to interpret things in terms of standard deviations. Do posterior_predict once with the actual data and once with a one standard deviation increase for all observations on the variable of interest. The distribution of the difference…