我正在使用R中的comand ratioBF()来估计比例的BayesFactor。我想复制维基百科https://en.wikipedia.org/wiki/Bayes_factor中的示例。比例BF命令的输出在日志中,因此我在输入错误时试图转换日志的输出,使其结果与维基百科的示例一致。我的代码是:
library(BayesFactor)
bf = proportionBF(115, 200, p = 1/2)
exp(bf@bayesFactor$bf)
任何人都可以帮助我了解哪个是我的错误对象。
答案 0 :(得分:0)
根据R文档(您可以通过键入?proportionBF来找到它)。
默认选择是λ〜logistic(λ_0,r)
其中lambda_0 = logit(p_0)和lambda = logit(p)。 r作为优先尺度参数。
但是从Wiki来说,q的先前分布来自均匀分布。