Hi all,

I have run a bayesian linear mixed model using brm(), but when I try to use bayes_factor, I get a .Call error message. Specifically,

I created a brm fit object using default prior:

fit1<-brm(y~x+(1|gp1)+(1|g2)), data=data, save_all_pars=TRUE)

fit0<-brm(y~(1|gp1)+(1|gp2),data=data, save_all_pars=TRUE)

These all run well. But when I used:

bayes_factor(fit1,fit0)

I got two error messages:

Error in .Call(“dmvtCpp”, X_ = X, mu_ = mu, sigma_ = sigma, df = -1, log_ = log, :

“dmvtCpp” not resolved from current namespace (mvnfast)

Error: Bridgesampling failed. Did you set ‘save_all_pars’ to TRUE when fitting your model?

For the second one, I did set save_all_pars=TRUE so I guess it’s just an error that’s tied to the first one?

I tried googling the problem on .call error, but can’t find any good solutions. Any help is much appreciated!