Hi everyone -
I am trying to obtain Bayes Factors for model comparison (more than two models), and I used cmdstan_stan
to fit all models. I’ve tried functions from different packages such as bayestestR
and brms
but failed because the CmdStanR fit object is incompatible with those functions according to the error message.
I also tried to convert the CmdStanR fit objects into RStan fit objects. Though the converting process is successful, the Bayes Factor using the bayesfactor()
function takes forever and I never get the output.
I wonder how I should proceed, or at least try. Any answer/direction/advice/idea would be much, much appreciative. Thanks much.
Matt