In debruine/bfrr: Bayes Factors and Robustness Regions

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
set.seed(8675309)

library(bfrr)

a <- rnorm(30, 0, 1)
b <- rnorm(30, 0.4, 1)

tt <- t.test(a, b)

sample_mean <- tt$estimate[[1]] - tt$estimate[[2]]
sample_se <- tt$stderr
sample_df <- tt$parameter[[1]]

rr <- bfrr(sample_mean, sample_se, sample_df, 
     model = "normal", mean = 0, sd = 0.5, tail = 1)

summary(rr)

plot(rr)

debruine/bfrr documentation built on March 7, 2020, 5:47 p.m.