Я хочу вычислить байесовский коэффициент для теста пропорции одной выборки, где нулевая гипотеза состоит в том, что -
Пропорции для уровней для номинальной переменной не отличаются от теоретических пропорций
Я знаю, как использовать proportionBF из пакета BayesFactor, чтобы сделать это, когда есть две пропорции с нулевой гипотезой, согласно которой оба уровня одинаково вероятны (p = 0.5).Это эквивалентно stats::binom.test.
# frequentist test (p-value > 0.05)
broom::tidy(stats::binom.test(x = 21, n = 44, p = 0.5))
#> # A tibble: 1 x 8
#> estimate statistic p.value parameter conf.low conf.high method
#> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <chr>
#> 1 0.477 21 0.880 44 0.325 0.633 Exact~
#> # ... with 1 more variable: alternative <chr>
# quantifying evidence in favor of null hypothesis
library(BayesFactor)
1/ proportionBF(y = 21, N = 44, p = 0.5, rscale = 0.707)
#> Bayes factor analysis
#> --------------
#> [1] Null, p=0.5 : 3.724132 ±0.02%
#>
#> Against denominator:
#> Alternative, p0 = 0.5, r = 0.707, p =/= p0
#> ---
#> Bayes factor type: BFproportion, logistic
Но я не знаю, как выполнить этот анализ, когда существует более двух пропорций, что эквивалентно stats::chisq.test (тест на соответствие пригодности).
# frequentist test (p-value > 0.05)
table(mtcars$cyl)
#>
#> 4 6 8
#> 11 7 14
broom::tidy(stats::chisq.test(x = table(mtcars$cyl)))
#> # A tibble: 1 x 4
#> statistic p.value parameter method
#> <dbl> <dbl> <dbl> <chr>
#> 1 2.31 0.315 2 Chi-squared test for given probabilities
Как вычислить коэффициент Байеса для этого теста с использованием пакета BayesFactor?
(PS Если это невозможно, мне также будет интересно узнать о любом другом пакете, который можетто же самое.)