Anyway to calculate the bayes factor of the model conveniently?

[Krozard]

I am try to tuning the emission function for my data and have some candidate, is there a way to calculate the bayes factor of the pyhsmm outputs? Or in other words calculate the marginal likelihood for a fitmodel?