fit_maihda: Fit MAIHDA Model

In MAIHDA: Multilevel Analysis of Individual Heterogeneity and Discriminatory Accuracy

Fit MAIHDA Model

Description

Fits a multilevel model for MAIHDA (Multilevel Analysis of Individual Heterogeneity and Discriminatory Accuracy) using either lme4 or brms.

Usage

fit_maihda(formula, data, engine = "lme4", family = "gaussian", ...)

Arguments

formula	A formula specifying the model. Should include random effect for stratum (e.g., outcome ~ fixed_vars + (1 | stratum)).
data	A data frame containing the variables in the formula.
engine	Character string specifying which engine to use: "lme4" (default) or "brms".
family	Character string or family object specifying the model family. Common options: "gaussian", "binomial", "poisson". Default is "gaussian".
...	Additional arguments passed to lmer/glmer (lme4) or brm (brms).

Value

A maihda_model object containing:

model	The fitted model object (lme4 or brms)
engine	The engine used ("lme4" or "brms")
formula	The model formula
data	The data used for fitting
family	The family used
strata_info	The strata information from make_strata() if available, NULL otherwise

Examples

# Create strata
strata_result <- make_strata(maihda_sim_data, vars = c("gender", "race", "education"))

# Fit model with lme4
model <- fit_maihda(health_outcome ~ age + (1 | stratum),
                    data = strata_result$data,
                    engine = "lme4")

# Fit model with brms (if brms is available)
# model_brms <- fit_maihda(health_outcome ~ age + (1 | stratum),
#                          data = strata_result$data,
#                          engine = "brms")

MAIHDA documentation built on April 5, 2026, 5:06 p.m.