R2_bayes function

Bayesian R-squared for flexreg Objects

Bayesian version of R-squared for flexible regression models for bounded continuous and discrete responses.

R2_bayes(model)

Arguments

• model: an object (or a list of objects) of class flexreg, usually the result of flexreg or flexreg_binom functions.

Returns

A list with the same length as the number of objects of class flexreg passed in the model argument. Each element of the list contains a vector of Bayesian R-squared values with the same length as the Markov Chain(s) after warmup.

Details

The function provides a Bayesian version of the R-squared measure, defined as the variance of the predicted values divided by itself plus the expected variance of the errors.

Examples

data("Reading")
FB <- flexreg(accuracy.adj ~ iq, data = Reading, type = "FB", n.iter=1000)
hist(R2_bayes(FB)[[1]])

References

Gelman, A., Goodrich, B., Gabry, J., Vehtari, A. (2019). R-squared for Bayesian Regression Models, The American Statistician, 73:3, 307--309. doi: 10.1080/00031305.2018.1549100