FEV1 data comparison between Bayesian and frequentist MMRMs
2026-20-10
About
This vignette provides an example comparison of a Bayesian MMRM fit,
obtained by brms.mmrm::brm_model(), and a frequentist MMRM
fit, obtained by mmrm::mmrm(). An overview of parameter
estimates and differences by type of MMRM is given in the summary (Tables 4 and 5) at the end.
1 Prerequisites
This comparison workflow requires the following packages.
> packages <- c(
+ "dplyr",
+ "tidyr",
+ "ggplot2",
+ "gt",
+ "gtsummary",
+ "purrr",
+ "parallel",
+ "brms.mmrm",
+ "mmrm",
+ "posterior"
+ )
> invisible(lapply(packages, library, character.only = TRUE))We set a seed for the random number generator to ensure statistical reproducibility.
> set.seed(123L)2 Data
2.1 Pre-processing
This analysis exercise uses the fev_dat dataset
contained in the mmrm-package:
> data(fev_data, package = "mmrm")It is an artificial (simulated) dataset of a clinical trial
investigating the effect of an active treatment on FEV1
(forced expired volume in one second), compared to placebo.
FEV1 is a measure of how quickly the lungs can be emptied
and low levels may indicate chronic obstructive pulmonary disease
(COPD).
The dataset is a tibble with 800 rows and the following
notable variables:
USUBJID(subject ID)AVISIT(visit number, factor)VISITN(visit number, numeric)ARMCD(treatment,TRTorPBO)RACE(3-category race)SEX(female or male)FEV1_BL(FEV1 at baseline, %)FEV1(FEV1 at study visits)WEIGHT(weighting variable)
The primary endpoint for the analysis is change from baseline in
FEV1, which we derive below and denote
FEV1_CHG.
> fev_data <- fev_data |>
+ mutate("FEV1_CHG" = FEV1 - FEV1_BL)The rest of the pre-processing steps create factors for the study arm
and visit and apply the usual checking and standardization steps of
brms.mmrm::brm_data().
> fev_data <- brm_data(
+ data = fev_data,
+ outcome = "FEV1_CHG",
+ group = "ARMCD",
+ time = "AVISIT",
+ patient = "USUBJID",
+ baseline = "FEV1_BL",
+ reference_group = "PBO",
+ covariates = c("RACE", "SEX")
+ ) |>
+ brm_data_chronologize(order = "VISITN")The following table shows the first rows of the dataset.
> head(fev_data) |>
+ gt() |>
+ tab_caption(caption = md("Table 1. First rows of the pre-processed `fev_dat` dataset."))| USUBJID | AVISIT | ARMCD | RACE | SEX | FEV1_BL | FEV1 | WEIGHT | VISITN | VISITN2 | FEV1_CHG |
|---|---|---|---|---|---|---|---|---|---|---|
| PT2 | VIS1 | PBO | Asian | Male | 45.02477 | NA | 0.4651848 | 1 | 0.3295078 | NA |
| PT2 | VIS2 | PBO | Asian | Male | 45.02477 | 31.45522 | 0.2330974 | 2 | -0.8204684 | -13.569552 |
| PT2 | VIS3 | PBO | Asian | Male | 45.02477 | 36.87889 | 0.3600763 | 3 | 0.4874291 | -8.145878 |
| PT2 | VIS4 | PBO | Asian | Male | 45.02477 | 48.80809 | 0.5073795 | 4 | 0.7383247 | 3.783324 |
| PT3 | VIS1 | PBO | Black or African American | Female | 43.50070 | NA | 0.6821642 | 1 | 0.5757814 | NA |
| PT3 | VIS2 | PBO | Black or African American | Female | 43.50070 | 35.98699 | 0.8917896 | 2 | -0.3053884 | -7.513705 |
2.2 Descriptive statistics
Table of baseline characteristics:
> fev_data |>
+ select(ARMCD, USUBJID, SEX, RACE, FEV1_BL) |>
+ distinct() |>
+ select(-USUBJID) |>
+ tbl_summary(
+ by = c(ARMCD),
+ statistic = list(
+ all_continuous() ~ "{mean} ({sd})",
+ all_categorical() ~ "{n} / {N} ({p}%)"
+ )
+ ) |>
+ modify_caption("Table 2. Baseline characteristics.")| Characteristic | PBO N = 1051 |
TRT N = 951 |
|---|---|---|
| SEX | ||
| Male | 50 / 105 (48%) | 44 / 95 (46%) |
| Female | 55 / 105 (52%) | 51 / 95 (54%) |
| RACE | ||
| Asian | 38 / 105 (36%) | 32 / 95 (34%) |
| Black or African American | 46 / 105 (44%) | 29 / 95 (31%) |
| White | 21 / 105 (20%) | 34 / 95 (36%) |
| FEV1_BL | 40 (9) | 40 (9) |
| 1 n / N (%); Mean (SD) | ||
Table of change from baseline in FEV1 over 52 weeks:
> fev_data |>
+ pull(AVISIT) |>
+ unique() |>
+ sort() |>
+ purrr::map(
+ .f = ~ fev_data |>
+ filter(AVISIT %in% .x) |>
+ tbl_summary(
+ by = ARMCD,
+ include = FEV1_CHG,
+ type = FEV1_CHG ~ "continuous2",
+ statistic = FEV1_CHG ~ c(
+ "{mean} ({sd})",
+ "{median} ({p25}, {p75})",
+ "{min}, {max}"
+ ),
+ label = list(FEV1_CHG = paste("Visit ", .x))
+ )
+ ) |>
+ tbl_stack(quiet = TRUE) |>
+ modify_caption("Table 3. Change from baseline.")| Characteristic | PBO N = 105 |
TRT N = 95 |
|---|---|---|
| Visit VIS1 | ||
| Mean (SD) | -8 (9) | -2 (10) |
| Median (Q1, Q3) | -9 (-16, 0) | -4 (-9, 7) |
| Min, Max | -26, 12 | -24, 20 |
| Unknown | 37 | 29 |
| Visit VIS2 | ||
| Mean (SD) | -3 (8) | 2 (9) |
| Median (Q1, Q3) | -3 (-10, 1) | 2 (-4, 8) |
| Min, Max | -20, 15 | -22, 23 |
| Unknown | 36 | 24 |
| Visit VIS3 | ||
| Mean (SD) | 2 (8) | 5 (9) |
| Median (Q1, Q3) | 2 (-2, 8) | 6 (0, 11) |
| Min, Max | -15, 20 | -19, 30 |
| Unknown | 34 | 37 |
| Visit VIS4 | ||
| Mean (SD) | 8 (12) | 13 (13) |
| Median (Q1, Q3) | 6 (1, 20) | 12 (5, 23) |
| Min, Max | -20, 39 | -14, 47 |
| Unknown | 38 | 28 |
The following figure shows the primary endpoint over the four study visits in the data.
> fev_data |>
+ group_by(ARMCD) |>
+ ggplot(aes(x = AVISIT, y = FEV1_CHG, fill = factor(ARMCD))) +
+ geom_hline(yintercept = 0, col = "grey", linewidth = 1.2) +
+ geom_boxplot(na.rm = TRUE) +
+ labs(
+ x = "Visit",
+ y = "Change from baseline in FEV1",
+ fill = "Treatment"
+ ) +
+ scale_fill_manual(values = c("darkgoldenrod2", "coral2")) +
+ theme_bw()
Figure 1. Change from baseline in FEV1 over 4 visit time points.
3 Fitting MMRMs
3.1 Bayesian model
The formula for the Bayesian model includes additive effects for baseline, study visit, race, sex, and study-arm-by-visit interaction.
> b_mmrm_formula <- brm_formula(
+ data = fev_data,
+ intercept = TRUE,
+ baseline = TRUE,
+ group = FALSE,
+ time = TRUE,
+ baseline_time = FALSE,
+ group_time = TRUE,
+ correlation = "unstructured"
+ )
> print(b_mmrm_formula)
#> FEV1_CHG ~ FEV1_BL + ARMCD:AVISIT + AVISIT + RACE + SEX + unstr(time = AVISIT, gr = USUBJID)
#> sigma ~ 0 + AVISITWe fit the model using brms.mmrm::brm_model(). To ensure
a good basis of comparison with the frequentist model, we put an
extremely diffuse prior on the intercept. The parameters already have
diffuse flexible priors by default.
> b_mmrm_fit <- brm_model(
+ data = filter(fev_data, !is.na(FEV1_CHG)),
+ formula = b_mmrm_formula,
+ prior = brms::prior(class = "Intercept", prior = "student_t(3, 0, 1000)"),
+ iter = 10000,
+ warmup = 2000,
+ chains = 4,
+ cores = 4,
+ seed = 1,
+ refresh = 0
+ )
#> Running MCMC with 4 parallel chains...
#>
#> Chain 2 finished in 22.6 seconds.
#> Chain 4 finished in 23.0 seconds.
#> Chain 1 finished in 23.4 seconds.
#> Chain 3 finished in 23.3 seconds.
#>
#> All 4 chains finished successfully.
#> Mean chain execution time: 23.1 seconds.
#> Total execution time: 23.5 seconds.Here is a posterior summary of model parameters, including fixed effects and pairwise correlation among visits within patients.
> summary(b_mmrm_fit)
#> Family: gaussian
#> Links: mu = identity; sigma = log
#> Formula: FEV1_CHG ~ FEV1_BL + ARMCD:AVISIT + AVISIT + RACE + SEX + unstr(time = AVISIT, gr = USUBJID)
#> sigma ~ 0 + AVISIT
#> Data: modeled_data (Number of observations: 537)
#> Draws: 4 chains, each with iter = 10000; warmup = 2000; thin = 1;
#> total post-warmup draws = 32000
#>
#> Correlation Structures:
#> Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
#> cortime(VIS1,VIS2) 0.36 0.08 0.19 0.52 1.00 55170 25674
#> cortime(VIS1,VIS3) 0.14 0.10 -0.06 0.33 1.00 57265 24050
#> cortime(VIS2,VIS3) 0.04 0.10 -0.15 0.23 1.00 50963 24352
#> cortime(VIS1,VIS4) 0.17 0.11 -0.06 0.38 1.00 57250 24491
#> cortime(VIS2,VIS4) 0.11 0.09 -0.06 0.28 1.00 57480 25753
#> cortime(VIS3,VIS4) 0.01 0.10 -0.18 0.21 1.00 49870 24634
#>
#> Regression Coefficients:
#> Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
#> Intercept 24.36 1.41 21.57 27.12 1.00 48546 25877
#> FEV1_BL -0.84 0.03 -0.90 -0.79 1.00 64473 23059
#> AVISIT2 4.79 0.81 3.22 6.38 1.00 37872 26950
#> AVISIT3 10.37 0.82 8.77 12.00 1.00 34320 25281
#> AVISIT4 15.19 1.32 12.59 17.76 1.00 37729 25306
#> RACEBlackorAfricanAmerican 1.41 0.59 0.25 2.57 1.00 52830 25732
#> RACEWhite 5.45 0.62 4.24 6.65 1.00 51233 25125
#> SEXFemale 0.34 0.51 -0.65 1.33 1.00 57472 26156
#> AVISITVIS1:ARMCDTRT 3.98 1.07 1.87 6.09 1.00 35035 26023
#> AVISITVIS2:ARMCDTRT 3.94 0.82 2.32 5.53 1.00 53609 23736
#> AVISITVIS3:ARMCDTRT 2.98 0.68 1.65 4.31 1.00 55724 24704
#> AVISITVIS4:ARMCDTRT 4.41 1.70 1.11 7.73 1.00 45487 25110
#> sigma_AVISITVIS1 1.83 0.06 1.71 1.95 1.00 57928 25818
#> sigma_AVISITVIS2 1.59 0.06 1.47 1.71 1.00 53397 23759
#> sigma_AVISITVIS3 1.33 0.06 1.20 1.46 1.00 55513 23227
#> sigma_AVISITVIS4 2.28 0.06 2.16 2.41 1.00 61076 23707
#>
#> Draws were sampled using sample(hmc). For each parameter, Bulk_ESS
#> and Tail_ESS are effective sample size measures, and Rhat is the potential
#> scale reduction factor on split chains (at convergence, Rhat = 1).3.2 Frequentist model
The formula for the frequentist model is the same, except for the different syntax for specifying the covariance structure of the MMRM. We fit the model below.
> f_mmrm_fit <- mmrm::mmrm(
+ formula = FEV1_CHG ~ FEV1_BL + ARMCD:AVISIT + AVISIT + RACE + SEX +
+ us(AVISIT | USUBJID),
+ data = mutate(
+ fev_data,
+ AVISIT = factor(as.character(AVISIT), ordered = FALSE)
+ )
+ )The parameter summaries of the frequentist model are below.
> summary(f_mmrm_fit)
#> mmrm fit
#>
#> Formula: FEV1_CHG ~ FEV1_BL + ARMCD:AVISIT + AVISIT + RACE + SEX + us(AVISIT | USUBJID)
#> Data: mutate(fev_data, AVISIT = factor(as.character(AVISIT), ordered = FALSE)) (used 537
#> observations from 197 subjects with maximum 4 timepoints)
#> Covariance: unstructured (10 variance parameters)
#> Method: Satterthwaite
#> Vcov Method: Asymptotic
#> Inference: REML
#>
#> Model selection criteria:
#> AIC BIC logLik deviance
#> 3381.4 3414.2 -1680.7 3361.4
#>
#> Coefficients:
#> Estimate Std. Error df t value Pr(>|t|)
#> (Intercept) 24.35372 1.40754 257.97000 17.302 < 2e-16 ***
#> FEV1_BL -0.84022 0.02777 190.27000 -30.251 < 2e-16 ***
#> AVISITVIS2 4.79036 0.79848 144.82000 5.999 1.51e-08 ***
#> AVISITVIS3 10.36601 0.81318 157.08000 12.748 < 2e-16 ***
#> AVISITVIS4 15.19231 1.30857 139.25000 11.610 < 2e-16 ***
#> RACEBlack or African American 1.41921 0.57874 169.56000 2.452 0.015211 *
#> RACEWhite 5.45679 0.61626 157.54000 8.855 1.65e-15 ***
#> SEXFemale 0.33812 0.49273 166.43000 0.686 0.493529
#> AVISITVIS1:ARMCDTRT 3.98329 1.04540 142.32000 3.810 0.000206 ***
#> AVISITVIS2:ARMCDTRT 3.93076 0.81351 142.26000 4.832 3.46e-06 ***
#> AVISITVIS3:ARMCDTRT 2.98372 0.66567 129.61000 4.482 1.61e-05 ***
#> AVISITVIS4:ARMCDTRT 4.40400 1.66049 132.88000 2.652 0.008970 **
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#>
#> Covariance estimate:
#> VIS1 VIS2 VIS3 VIS4
#> VIS1 37.8301 11.3255 3.4796 10.6844
#> VIS2 11.3255 23.5476 0.7760 5.5103
#> VIS3 3.4796 0.7760 13.8037 0.5683
#> VIS4 10.6844 5.5103 0.5683 92.96254 Comparison
This section compares the Bayesian posterior parameter estimates from
brms.mmrm to the frequentist parameter estimates of the
mmrm package.
4.1 Extract estimates from Bayesian model
We extract and standardize the Bayesian estimates.
> b_mmrm_draws <- b_mmrm_fit |>
+ as_draws_df()
> visit_levels <- sort(unique(as.character(fev_data$AVISIT)))
> for (level in visit_levels) {
+ name <- paste0("b_sigma_AVISIT", level)
+ b_mmrm_draws[[name]] <- exp(b_mmrm_draws[[name]])
+ }
> b_mmrm_summary <- b_mmrm_draws |>
+ summarize_draws() |>
+ select(variable, mean, sd) |>
+ filter(!(variable %in% c("Intercept", "lprior", "lp__"))) |>
+ rename(bayes_estimate = mean, bayes_se = sd) |>
+ mutate(
+ variable = variable |>
+ tolower() |>
+ gsub(pattern = "b_", replacement = "") |>
+ gsub(pattern = "b_sigma_AVISIT", replacement = "sigma_") |>
+ gsub(pattern = "cortime", replacement = "correlation") |>
+ gsub(pattern = "__", replacement = "_") |>
+ gsub(pattern = "avisitvis", replacement = "avisit")
+ )4.2 Extract estimates from frequentist model
We extract and standardize the frequentist estimates.
> f_mmrm_fixed <- summary(f_mmrm_fit)$coefficients |>
+ as_tibble(rownames = "variable") |>
+ mutate(variable = tolower(variable)) |>
+ mutate(variable = gsub("(", "", variable, fixed = TRUE)) |>
+ mutate(variable = gsub(")", "", variable, fixed = TRUE)) |>
+ mutate(variable = gsub("avisitvis", "avisit", variable)) |>
+ rename(freq_estimate = Estimate, freq_se = `Std. Error`) |>
+ select(variable, freq_estimate, freq_se)> f_mmrm_variance <- tibble(
+ variable = paste0("sigma_AVISIT", visit_levels) |>
+ tolower() |>
+ gsub(pattern = "avisitvis", replacement = "avisit"),
+ freq_estimate = sqrt(diag(f_mmrm_fit$cov))
+ )> f_diagonal_factor <- diag(1 / sqrt(diag(f_mmrm_fit$cov)))
> f_corr_matrix <- f_diagonal_factor %*% f_mmrm_fit$cov %*% f_diagonal_factor
> colnames(f_corr_matrix) <- visit_levels> f_mmrm_correlation <- f_corr_matrix |>
+ as.data.frame() |>
+ as_tibble() |>
+ mutate(x1 = visit_levels) |>
+ pivot_longer(
+ cols = -any_of("x1"),
+ names_to = "x2",
+ values_to = "freq_estimate"
+ ) |>
+ filter(
+ as.numeric(gsub("[^0-9]", "", x1)) < as.numeric(gsub("[^0-9]", "", x2))
+ ) |>
+ mutate(variable = sprintf("correlation_%s_%s", x1, x2)) |>
+ select(variable, freq_estimate)> f_mmrm_summary <- bind_rows(
+ f_mmrm_fixed,
+ f_mmrm_variance,
+ f_mmrm_correlation
+ ) |>
+ mutate(variable = gsub("\\s+", "", variable) |> tolower())4.3 Summary
The first table below summarizes the parameter estimates from each model and the differences between estimates (Bayesian minus frequentist). The second table shows the standard errors of these estimates and differences between standard errors. In each table, the “Relative” column shows the relative difference (the difference divided by the frequentist quantity).
Because of the different statistical paradigms and estimation
procedures, especially regarding the covariance parameters, it would not
be realistic to expect the Bayesian and frequentist approaches to yield
virtually identical results. Nevertheless, the absolute and relative
differences in the table below show strong agreement between
brms.mmrm and mmrm.
> b_f_comparison <- full_join(
+ x = b_mmrm_summary,
+ y = f_mmrm_summary,
+ by = "variable"
+ ) |>
+ mutate(
+ diff_estimate = bayes_estimate - freq_estimate,
+ diff_relative_estimate = diff_estimate / freq_estimate,
+ diff_se = bayes_se - freq_se,
+ diff_relative_se = diff_se / freq_se
+ ) |>
+ select(variable, ends_with("estimate"), ends_with("se"))> table_estimates <- b_f_comparison |>
+ select(variable, ends_with("estimate"))
> gt(table_estimates) |>
+ fmt_number(decimals = 4) |>
+ tab_caption(
+ caption = md(
+ paste(
+ "Table 4. Comparison of parameter estimates between",
+ "Bayesian and frequentist MMRMs."
+ )
+ )
+ ) |>
+ cols_label(
+ variable = "Variable",
+ bayes_estimate = "Bayesian",
+ freq_estimate = "Frequentist",
+ diff_estimate = "Difference",
+ diff_relative_estimate = "Relative"
+ )| Variable | Bayesian | Frequentist | Difference | Relative |
|---|---|---|---|---|
| intercept | 24.3554 | 24.3537 | 0.0017 | 0.0001 |
| fev1_bl | −0.8402 | −0.8402 | 0.0000 | −0.0001 |
| avisit2 | 4.7922 | 4.7904 | 0.0018 | 0.0004 |
| avisit3 | 10.3706 | 10.3660 | 0.0046 | 0.0004 |
| avisit4 | 15.1916 | 15.1923 | −0.0007 | 0.0000 |
| raceblackorafricanamerican | 1.4121 | 1.4192 | −0.0071 | −0.0050 |
| racewhite | 5.4472 | 5.4568 | −0.0096 | −0.0018 |
| sexfemale | 0.3426 | 0.3381 | 0.0045 | 0.0132 |
| avisit1:armcdtrt | 3.9763 | 3.9833 | −0.0069 | −0.0017 |
| avisit2:armcdtrt | 3.9363 | 3.9308 | 0.0055 | 0.0014 |
| avisit3:armcdtrt | 2.9781 | 2.9837 | −0.0056 | −0.0019 |
| avisit4:armcdtrt | 4.4076 | 4.4040 | 0.0036 | 0.0008 |
| sigma_avisit1 | 6.2339 | 6.1506 | 0.0833 | 0.0135 |
| sigma_avisit2 | 4.9149 | 4.8526 | 0.0623 | 0.0128 |
| sigma_avisit3 | 3.7760 | 3.7153 | 0.0606 | 0.0163 |
| sigma_avisit4 | 9.8025 | 9.6417 | 0.1608 | 0.0167 |
| correlation_vis1_vis2 | 0.3605 | 0.3795 | −0.0189 | −0.0499 |
| correlation_vis1_vis3 | 0.1415 | 0.1523 | −0.0108 | −0.0709 |
| correlation_vis2_vis3 | 0.0388 | 0.0430 | −0.0042 | −0.0975 |
| correlation_vis1_vis4 | 0.1679 | 0.1802 | −0.0123 | −0.0683 |
| correlation_vis2_vis4 | 0.1110 | 0.1178 | −0.0068 | −0.0576 |
| correlation_vis3_vis4 | 0.0135 | 0.0159 | −0.0023 | −0.1472 |
> table_se <- b_f_comparison |>
+ select(variable, ends_with("se")) |>
+ filter(!is.na(freq_se))
> gt(table_se) |>
+ fmt_number(decimals = 4) |>
+ tab_caption(
+ caption = md(
+ paste(
+ "Table 5. Comparison of parameter standard errors between",
+ "Bayesian and frequentist MMRMs."
+ )
+ )
+ ) |>
+ cols_label(
+ variable = "Variable",
+ bayes_se = "Bayesian",
+ freq_se = "Frequentist",
+ diff_se = "Difference",
+ diff_relative_se = "Relative"
+ )| Variable | Bayesian | Frequentist | Difference | Relative |
|---|---|---|---|---|
| intercept | 1.4129 | 1.4075 | 0.0053 | 0.0038 |
| fev1_bl | 0.0281 | 0.0278 | 0.0003 | 0.0104 |
| avisit2 | 0.8069 | 0.7985 | 0.0085 | 0.0106 |
| avisit3 | 0.8225 | 0.8132 | 0.0093 | 0.0115 |
| avisit4 | 1.3223 | 1.3086 | 0.0137 | 0.0105 |
| raceblackorafricanamerican | 0.5886 | 0.5787 | 0.0099 | 0.0171 |
| racewhite | 0.6210 | 0.6163 | 0.0047 | 0.0077 |
| sexfemale | 0.5065 | 0.4927 | 0.0138 | 0.0280 |
| avisit1:armcdtrt | 1.0656 | 1.0454 | 0.0202 | 0.0193 |
| avisit2:armcdtrt | 0.8238 | 0.8135 | 0.0103 | 0.0126 |
| avisit3:armcdtrt | 0.6765 | 0.6657 | 0.0108 | 0.0162 |
| avisit4:armcdtrt | 1.7000 | 1.6605 | 0.0395 | 0.0238 |
5 Session info
> sessionInfo()
#> R version 4.5.2 (2025-10-31)
#> Platform: aarch64-apple-darwin20
#> Running under: macOS Tahoe 26.7.1
#>
#> Matrix products: default
#> BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib
#> LAPACK: /Library/Frameworks/R.framework/Versions/4.5-arm64/Resources/lib/libRlapack.dylib; LAPACK version 3.12.1
#>
#> locale:
#> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
#>
#> time zone: America/Indiana/Indianapolis
#> tzcode source: internal
#>
#> attached base packages:
#> [1] parallel stats graphics grDevices utils datasets methods base
#>
#> other attached packages:
#> [1] posterior_1.7.0 mmrm_0.3.19 purrr_1.2.2 gtsummary_2.6.1 gt_1.4.0
#> [6] ggplot2_4.0.3 tidyr_1.3.2 knitr_1.52 dplyr_1.2.1 brms.mmrm_1.1.1.9001
#> [11] brms_2.23.0 Rcpp_1.1.2
#>
#> loaded via a namespace (and not attached):
#> [1] tidyselect_1.2.1 svUnit_1.0.8 farver_2.1.2 loo_2.10.1 tidybayes_3.0.7
#> [6] S7_0.2.2 fastmap_1.2.0 tensorA_0.36.2.1 promises_1.5.0 digest_0.6.39
#> [11] estimability_2.0.0 lifecycle_1.0.5 StanHeaders_2.39.1 processx_3.9.0 magrittr_2.0.5
#> [16] compiler_4.5.2 sass_0.4.10 rlang_1.3.0 tools_4.5.2 utf8_1.2.6
#> [21] data.table_1.18.6.1 labeling_0.4.3 bridgesampling_1.2-1 pkgbuild_1.4.8 curl_8.1.0
#> [26] plyr_1.8.9 xml2_1.6.0 RColorBrewer_1.1-3 cmdstanr_0.9.0 abind_1.4-8
#> [31] withr_3.0.3 grid_4.5.2 stats4_4.5.2 xtable_1.8-8 colorspace_2.1-4
#> [36] inline_0.3.21 emmeans_2.0.4 scales_1.4.0 gtools_3.9.5 MASS_7.3-66
#> [41] ggridges_0.5.7 cli_3.6.6 mvtnorm_1.4-2 generics_0.1.4 otel_0.2.0
#> [46] RcppParallel_6.2.1 binom_1.1-2 reshape2_1.4.5 commonmark_2.0.0 rstan_2.32.7
#> [51] stringr_1.6.0 bayesplot_1.15.0 matrixStats_1.5.0 vctrs_0.7.3 V8_8.2.0
#> [56] Matrix_1.7-6 jsonlite_2.0.0 litedown_0.11 arrayhelpers_1.1-2 ggdist_3.3.3
#> [61] glue_1.8.1 mcptools_1.0.3 codetools_0.2-20 ps_1.9.3 distributional_0.9.0
#> [66] stringi_1.8.9 gtable_0.3.6 later_1.4.8 QuickJSR_1.11.0 tibble_3.3.1
#> [71] pillar_1.11.1 htmltools_0.5.9 Brobdingnag_1.3-1 TMB_1.9.25 R6_2.6.1
#> [76] Rdpack_2.6.6 evaluate_1.0.5 lattice_0.23-1 markdown_2.0 cards_0.9.0
#> [81] rbibutils_2.4.1 backports_1.5.1 nanonext_1.10.3.9000 renv_1.3.1 trialr_0.1.6
#> [86] rstantools_2.7.1 coda_0.19-4.1 gridExtra_2.3.1 nlme_3.1-171 checkmate_2.3.4
#> [91] xfun_0.62 fs_2.1.0 zoo_1.9-1 pkgconfig_2.0.3