Skip to contents

Median causal effect via 1:1 nearest-neighbour PS matching

Usage

mrm_median_causal_effect(data, treatment_col, outcome_col, covariates)

Arguments

data

data.frame.

treatment_col

Binary 0/1 column.

outcome_col

Outcome column name.

covariates

Character vector of covariate columns.

Value

Named list with median_y1, median_y0, median_treatment_effect, n_matched, interpretation.

Examples

set.seed(2026)
n <- 200L
x <- rnorm(n)
D <- rbinom(n, 1, plogis(0.5 * x))
y <- 0.7 * D + 0.3 * x + rnorm(n, 0, 0.5)
df <- data.frame(D = D, y = y, age = x)
res <- mrm_median_causal_effect(df,
  treatment_col = "D",
  outcome_col = "y",
  covariates = "age"
)
res$median_treatment_effect
#> [1] 0.9529
res$n_matched
#> [1] 99