2016-12-06 5 views
0

Я следовал инструкциям о том, как запустить байесовский «т-тест» с использованием априорных по умолчанию в пакете BayesFactor в R.значение Байеса фактор в R пакете BayesFactor

Некоторые из возвращаемых значений являются астрономическими.

Вот пример сравнения с огромным коэффициентом Байеса:

#install.packages('BayesFactor') 
library(BayesFactor) 

condition1 <- c(0.94, 0.9, 0.96, 0.74, 1, 0.98, 0.86, 0.92, 0.918367346938776, 
       0.96, 0.4, 0.816326530612245, 0.8, 0.836734693877551, 0.56, 0.66, 
       0.605263157894737, 0.836734693877551, 0.84, 0.9, 0.92, 0.714285714285714, 
       0.82, 0.5, 0.565217391304348, 0.8, 0.62) 

condition2 <- c(0.34, 0.16, 0.23, 0.19, 0.71, 0.36, 0.02, 0.83, 0.11, 0.06, 
       0.27, 0.347368421052632, 0.21, 0.13953488372093, 0.11340206185567, 
       0.14, 0.142857142857143, 0.257731958762887, 0.15, 0.29, 0.67, 
       0.0515463917525773, 0.272727272727273, 0.0895522388059701, 0.0204081632653061, 
       0.13, 0.0612244897959184) 

bf = ttestBF(x = condition1, condition2, paired = TRUE) 
bf 

Это возвращает:

Bayes factor analysis 
-------------- 
[1] Alt., r=0.707 : 144035108289 ±0% 

Against denominator: 
    Null, mu = 0 
--- 
Bayes factor type: BFoneSample, JZS 

По большей части сравнения в диапазоне от менее 1 до нескольких сотен. Но я обеспокоен тем, что это значение (144035108289!) Свидетельствует о чем-то ошибочном с моей стороны.

FYI: значение p в критериях нулевой гипотезы по тем же данным, что указано выше = 4.649279e-14.

Любые заверения или сведения об этом возвращенном BF будут высоко оценены.

+1

Этот вопрос может быть лучше подходит для сообщества [Cross Validated] (http://stats.stackexchange.com), поскольку он является статистическим, а не вопросом программирования. – nya

ответ

0

Я вычислил BF, используя ручной ввод Т-значение и размер выборки, как это, используя один и тот же пакет:

exp(ttest.tstat(t=14.63, n1=27, rscale = 0.707)[['bf']]) 

Это дает тот же BF. Похоже, что это во многом связано с относительно большим размером выборки (27). Возвращенный BF, похоже, находится вверх и вверх.

 Смежные вопросы

  • Нет связанных вопросов^_^