TODO: * after the posterior mean and CI of the regression coefficients, also provide the posterior mean (SD) of the heterogeneity distribution of within-subjects variables
bug: add library(rstan) before rstan::stan_model and rstan::sampling, otherwise Error in cpp_object_initializer(.self, .refClassDef, …) : could not find function “cpp_object_initializer”
saving workspace problem