Bayes Factor with Robustness Region (normal)

Bayes Factor with Robustness Region (normal)

bfrr(sample_mean = 0, sample_se = 0.1, sample_df = 99,
  model = "normal", ..., criterion = 3, rr_interval = NA,
  precision = 0.05)

Arguments

sample_mean	the observed sample mean
sample_se	the observed sample standard error
sample_df	the observed sample degrees of freedom
model	the model under which to calculate likelihood (H0, normal or uniform
...	model parameters (mean, sd and tail for normal, lower and upper for uniform)
criterion	the cutoff Bayes Factor for concluding evidence in favour of H0 or H1
rr_interval	the parameter intervals within which to test for robustness
precision	the step size for calculating the robustness region (default 0.05)

Value

list

Examples

rr <- bfrr(0.5, 1/sqrt(30), 29, "normal", mean = 0, sd = .7, tail = 1)
summary(rr)
#> The likelihood of your data under the theoretical distribution N(0, 0.7) is 0.16. The likelihood of your data under the null distribution T(29) is 0.01. The Bayes Factor is 12.53; this test finds evidence for H1 with a criterion of 3. The region of theoretical model parameters that give the same conclusion is `HN(0, [0.1, 1.4])`. 
plot(rr)