2017-02-17 43 views
1

我使用的是BayesFactor包,但我用兩個非常相似的代碼得到的相同的數據兩種不同的輸出。我想知道哪一個是正確的?BayesFactor包裝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

據的?ttest.tstat

如果簡單爲TRUE,則返回貝葉斯因子(針對零)的文檔中提到的。如果爲 FALSE,則函數返回長度爲3的向量,其中包含貝葉斯因子以及貝葉斯因子上的比例誤差估計值 以及用於計算它的方法。

從本質上說,我們並不需要採取expsimple=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 
相關問題