using BFpack

BFpack.R

library(tidyverse)
library(BFpack)

n <- 12
kat_data <- tibble(condition = c(rep("reliable", n), 
                                 rep("unreliable", n)),
                   harder = c(rbernoulli(n, p = .8), 
                              rbernoulli(n, p = .3)))

kat_data %>% group_by(condition) %>% summarise(mean_harder = mean(harder))

mod <- glm(harder ~ condition, family = "binomial", data = kat_data)
summary(mod)

BF(mod)
BF(mod, hypothesis = "condition < 0")