贝叶斯因子边界计算
使用R语言中的“BayesFactor”的包可以计算贝叶斯因子边界,该包的安装方式为:
install.packages("BayesFactor")
计算贝叶斯因子边界的例子代码如下:
library(BayesFactor)
# 生成两个区分度不同的正态分布随机数
set.seed(123)
x1 <- rnorm(25, mean = 0, sd = 1)
x2 <- rnorm(25, mean = 0.5, sd = 1)
# 运用“比值积分对数”方法计算Bayes Factor
bf <- ttestBF(x1, x2, nullInterval = c(-Inf, 0))
# 计算贝叶斯因子边界
BayesFactor::bfInterval(bf, whichInterval = "ROPE")
代码中使用了正态分布生成随机数,并运用“比值积分对数”方法计算Bayes Factor,最后再调用“bfInterval”函数计算贝叶斯因子边界。