2017-02-17 9 views
1

私はBayesFactorパッケージを使用していますが、私は2つの非常に類似したコードを使用して同じデータため二つの異なる出力を取得します。 どちらが正しいのですか?ベイズ因子パッケージR:二つの異なる出力

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! ### 

答えて

0

単純がTRUEである場合、それは?ttest.tstat

の文書に記載されているが、(ヌルに対する)ベイズ因子を返します。 がFALSEの場合、 計算されたlog(e)ベイズ係数と、ベイズ係数の比例誤差推定 とそれを計算するために使用された方法を含む長さ3のベクトルを返します。

要するに、我々は行う必要はありませんexpときsimple=TRUE

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 
関連する問題