Effect sizes complement p-values: they answer \u201chow big is the effect?\u201d rather than \u201cis there an effect at all?\u201d. MORIE exposes the standard families used in carceral, public-health, and sociolegal research.
library(morie)
set.seed(7)
group_a <- rnorm(60, mean = 0.0)
group_b <- rnorm(60, mean = 0.6)
d <- cohens_d(group_a, group_b)
d
#> [1] -0.5043227
Cohen\u2019s d expresses the difference between two means in pooled standard-deviation units. Conventional benchmarks: 0.2 small, 0.5 medium, 0.8 large.
tab <- matrix(c(20, 10, 30,
15, 25, 35), nrow = 2, byrow = TRUE)
v <- cramers_v(tab)
v
#> [1] 0.209657
Cramer\u2019s V scales the chi-square statistic to a 0\u20131 association measure for contingency tables. It is particularly useful for the provincial-vs-federal Mandela-rate cross-comparisons in MRM.
Omega-squared is a less-biased effect-size estimator than eta-squared for a one-way ANOVA design.
omega_squared(f_stat = 5.2, df_between = 2, df_within = 87, n = 90)
#> [1] 0.08536585
For a binomial proportion, MORIE exposes Wilson, Clopper\u2013Pearson exact, and Wald CIs. Wilson is the default and is what we recommend in published papers (it has better small-sample coverage than Wald).
proportion_ci(35, 100) # Wilson, 95% CI
#> $p_hat
#> [1] 0.35
#>
#> $ci_lower
#> [1] 0.2636425
#>
#> $ci_upper
#> [1] 0.4474556
proportion_ci(35, 100, method = "exact") # Clopper-Pearson
#> $p_hat
#> [1] 0.35
#>
#> $ci_lower
#> [1] 0.2572938
#>
#> $ci_upper
#> [1] 0.4518494
proportion_ci(35, 100, method = "wald") # Wald
#> $p_hat
#> [1] 0.35
#>
#> $ci_lower
#> [1] 0.2565157
#>
#> $ci_upper
#> [1] 0.4434843
e_value(rr = 2.0)
#> $e_value
#> [1] 3.414214
#>
#> $e_value_ci
#> [1] NA
The E-value is the minimum strength of association that an unmeasured confounder would need on both treatment and outcome to fully explain away an observed risk ratio.
cohens_d and
cramers_v, see the survey-weighted
vignette.causal-inference
vignette.