r - 贝叶斯因子作为 ratioBF 的非逻辑输出

Question

我正在使用 R 中的命令 ratioBF() 估计比例的贝叶斯因子。我想复制维基百科https://en.wikipedia.org/wiki/Bayes_factor中的示例。ratioBF 命令的输出在 log 中，所以我将 log 的输出转换为试图与 wikipedia 的示例一致的结果，因为我犯了一个错误。我的代码是：

library(BayesFactor)
bf = proportionBF(115, 200, p = 1/2)
exp(bf@bayesFactor$bf)

任何人都可以帮助我了解哪个是我的错误。

标签： rbayesianhypothesis-test

Answer 1

根据 R 文档（您可以通过键入 ?proportionBF 找到它）。

默认选择是 λ~logistic(λ_0,r)

其中 lambda_0=logit(p_0) 和 lambda=logit(p)。r 用作先验比例参数。

但从 Wiki 来看，q 的先验分布来自均匀分布。