Tables and Graphs for Mixed Models for Repeated Measures (MMRM) • tern.mmrm

GitHub commit activity GitHub contributors GitHub last commit GitHub pull requests GitHub repo size GitHub language count

Overview

tern.mmrm provides an interface for mixed model repeated measures (MMRM) within the tern framework to produce commonly used tables (using rtables and graphs. It builds on the R-package mmrm for the actual MMRM calculations.

When to use this package

If you would like to use the tern framework for tabulation and graphs, then this package is ideal for your MMRM fits. However if you use another reporting framework then it will be better to directly use mmrm and perform the tabulation and plots differently.

Main Features

Fitting of MMRM models to continuous longitudinal data collected over several time points (called visits) and optionally treatment arms.
Tabulation and plots of least square means per visit and treatment arm.
Tabulation of model diagnostics (e.g. BIC).
Diagnostic graphs (e.g. Q-Q plot of residuals).
Tabulation and plots of the covariance matrix estimate.
Subgroup specific refitting of the MMRM model and resulting forest plot.

Installation

tern.mmrm is available on CRAN and you can install the latest released version with:

install.packages("tern.mmrm")

or you can install the latest development version directly from GitHub by running the following:

# install.packages("pak")
pak::pak("insightsengineering/tern.mmrm")

Getting started

You can get started by trying out the example:

library(tern.mmrm)
fit <- fit_mmrm(
  vars = list(
    response = "FEV1",
    covariates = c("RACE", "SEX"),
    id = "USUBJID",
    arm = "ARMCD",
    visit = "AVISIT"
  ),
  data = mmrm_test_data,
  cor_struct = "unstructured",
  weights_emmeans = "equal",
  averages_emmeans = list(
    "VIS1+2" = c("VIS1", "VIS2")
  )
)

This specifies an MMRM with the FEV1 outcome and the RACE and SEX covariates for subjects identified by USUBJID and treatment arm ARMCD observed over time points identified by AVISIT in the mmrm_test_data data set. An unstructured covariance matrix is assumed. Note that only restricted maximum likelihood (REML) can be used. The least square means assume equal weights for factor combinations, and in addition to the single visit specific results an average across visits VIS1 and VIS2 will be computed.

Details

For a more detailed introduction to all of the features of this package, look at the introduction vignette:

vignette("introduction")

tern.mmrm