BayesFactor Package R:两种不同的输出
BayesFactor Package R: Two different Output
我正在使用 BayesFactor
包,但我得到 两个不同的输出 对于 相同的数据 使用 两个非常相似的代码。我想知道哪个是正确的?
if(!require(BayesFactor)){install.packages('BayesFactor')}
require(BayesFactor)
##################################################
exp(ttest.tstat(t= 2 , n1=40, n2=40, nullInterval =c(0, Inf), rscale = sqrt(2)/2,
complement = FALSE, simple = FALSE)$bf) ### !CHECK THIS OUTPUT! ###
exp(ttest.tstat(t= 2 , n1=40, n2=40, nullInterval =c(0, Inf), rscale = sqrt(2)/2,
complement = FALSE, simple = TRUE)) ### !CHECK THIS OUTPUT! ###
在?ttest.tstat
的文档中有提到
If simple is TRUE, returns the Bayes factor (against the null). If
FALSE, the function returns a vector of length 3 containing the
computed log(e) Bayes factor, along with a proportional error estimate
on the Bayes factor and the method used to compute it.
本质上,我们不需要在simple=TRUE
时取exp
ttest.tstat(t= 2 , n1=40, n2=40, nullInterval =c(0, Inf),
rscale = sqrt(2)/2, complement = FALSE, simple = TRUE)
# B10
# 2.502954
而 simple=FALSE
exp(ttest.tstat(t= 2 , n1=40, n2=40, nullInterval =c(0, Inf),
rscale = sqrt(2)/2, complement = FALSE, simple = FALSE)$bf)
#[1] 2.502954
我正在使用 BayesFactor
包,但我得到 两个不同的输出 对于 相同的数据 使用 两个非常相似的代码。我想知道哪个是正确的?
if(!require(BayesFactor)){install.packages('BayesFactor')}
require(BayesFactor)
##################################################
exp(ttest.tstat(t= 2 , n1=40, n2=40, nullInterval =c(0, Inf), rscale = sqrt(2)/2,
complement = FALSE, simple = FALSE)$bf) ### !CHECK THIS OUTPUT! ###
exp(ttest.tstat(t= 2 , n1=40, n2=40, nullInterval =c(0, Inf), rscale = sqrt(2)/2,
complement = FALSE, simple = TRUE)) ### !CHECK THIS OUTPUT! ###
在?ttest.tstat
If simple is TRUE, returns the Bayes factor (against the null). If FALSE, the function returns a vector of length 3 containing the computed log(e) Bayes factor, along with a proportional error estimate on the Bayes factor and the method used to compute it.
本质上,我们不需要在simple=TRUE
exp
ttest.tstat(t= 2 , n1=40, n2=40, nullInterval =c(0, Inf),
rscale = sqrt(2)/2, complement = FALSE, simple = TRUE)
# B10
# 2.502954
而 simple=FALSE
exp(ttest.tstat(t= 2 , n1=40, n2=40, nullInterval =c(0, Inf),
rscale = sqrt(2)/2, complement = FALSE, simple = FALSE)$bf)
#[1] 2.502954