使用有序贝叶斯概率 - 如何计算贝叶斯因子？

Question

我试图为我运行的三个模型找到bayesfactor（下面的代码）。

std::union_members_count<>

但是，当我使用包

时

library(Zelig)
DV1.1 <- zelig(as.factor(DV1) ~ IV1, model = "oprobit.bayes", mcmc=500000, data = DataCL)
DV1.2 <- zelig(as.factor(DV1) ~ IV1 + IV2, model = "oprobit.bayes", mcmc=500000, data = DataCL)
DV1.3 <- zelig(as.factor(DV1) ~ IV1 + IV2 + IV3 , model = "oprobit.bayes", mcmc=500000, data = DataCL)

我收到错误代码：

library(BayesFactor)

BayesFactorDV1 <- BayesFactor (DV1.1, DV1.2, DV1.3)

有人能帮我吗？谢谢

我的数据结构如下所示：

 Error in BayesFactor(DV1.1, DV1.2, DV1.3) : argument not of class mcmc