Learn R Programming

Reference Based Multiple Imputation (`rbmi`)

Overview

The rbmi package is used for the imputation of missing data in clinical trials with continuous multivariate normal longitudinal outcomes. It supports imputation under a missing at random (MAR) assumption, reference-based imputation methods, and delta adjustments (as required for sensitivity analysis such as tipping point analyses). The package implements both Bayesian and approximate Bayesian multiple imputation combined with Rubin's rules for inference, and frequentist conditional mean imputation combined with (jackknife or bootstrap) resampling.

Installation

The package can be installed directly from CRAN via:

install.packages("rbmi")

Note that the usage of Bayesian multiple imputation requires the installation of the suggested package rstan.

install.packages("rstan")

Usage

The package is designed around its 4 core functions:

draws() - Fits multiple imputation models
impute() - Imputes multiple datasets
analyse() - Analyses multiple datasets
pool() - Pools multiple results into a single statistic

The basic usage of these core functions is described in the quickstart vignette:

vignette(topic = "quickstart", package = "rbmi")

Validation

For clarification on the current validation status of rbmi please see the FAQ vignette.

Support

For any help with regards to using the package or if you find a bug please create a GitHub issue

Copy Link

Version

Install

install.packages('rbmi')

Monthly Downloads

676

Version

1.4.1

License

Apache License (>= 2)

Issues

Pull Requests

Stars

Forks

Repository

https://github.com/insightsengineering/rbmi

Homepage

https://insightsengineering.github.io/rbmi/

Maintainer

Isaac Gravestock

Last Published

March 3rd, 2025

Functions in rbmi (1.4.1)

Creates a "MMRM" ready dataset

Create an imputation object

Convert object to dataframe

Construct an analysis object

Convert indicator to index

Diagnostics of the MCMC based on ESS

check_hmc_diagn

Diagnostics of the MCMC based on HMC-related measures.

Create vector of Stratas

Creates a draws object

as_cropped_char

as_cropped_char

Diagnostics of the MCMC

Compute covariance matrix for some reference-based methods (JR, CIR)

Convert character variables to factor

assert_variables_exist

Assert that all variables exist within a dataset

as_simple_formula

Creates a simple formula object from a string

Expand data.frame into a design matrix

as_mmrm_formula

Create MMRM formula

Fit the base imputation model and get parameter estimates

Ensure rstan exists

Evaluate a call to mmrm

Expand and fill in missing data.frame rows

Get imputation strategies

Extract the Effective Sample Size (ESS) from a stanfit object

extract_covariates

Extract Variables from string vector

Fit the base imputation model using a Bayesian approach

get_jackknife_stack

Creates a stack object populated with jackknife samples

extract_data_nmar_as_na

Set to NA outcome values that would be MNAR if they were missing (i.e. which occur after an ICE handled using a reference-based imputation strategy)

Fit a MMRM model

get_pattern_groups

Determine patients missingness group

get_mmrm_sample

Fit MMRM and returns parameter estimates

Control the computational details of the imputation methods

convert_to_imputation_list_df

Convert list of imputation_list_single() objects to an imputation_list_df() object (i.e. a list of imputation_df() objects's)

Extract parameters from a MMRM model

extract_imputed_dfs

Extract imputed datasets

get_pattern_groups_unique

Get Pattern Summary

get_pool_components

Expected Pool Components

get_session_hash

Get session hash

format_method_descriptions

Format method descriptions

Create a valid imputation_df object

imputation_list_df

List of imputations_df

Is value absent

Is character or factor

Does object have a class ?

Fit the base imputation model on bootstrap samples

imputation_single

Create a valid imputation_single object

Is single character

get_delta_template

Get delta utility variables

imputation_list_single

A collection of imputation_singles() grouped by a single subjid ID

is_in_rbmi_development

Is package in development mode?

Calculate delta from a lagged scale coefficient

Extract draws from a stanfit object

Create a delta data.frame template

Invert and derive indexes

extract_imputed_df

Extract imputed dataset

get_bootstrap_stack

Creates a stack object populated with bootstrapped samples

generate_data_single

Generate data for a single group

longDataConstructor

R6 Class for Storing / Accessing & Sampling Longitudinal Data

get_conditional_parameters

Derive conditional multivariate normal parameters

Calculate design vector for the lsmeans

make_rbmi_cluster

Create a rbmi ready cluster

Least Square Means

Von Hippel and Bartlett pooling of BMLMI method

Get Compiled Stan Object

get_visit_distribution_parameters

Derive visit distribution parameters

impute_internal

Create imputed datasets

Sample outcome value

get_example_data

Simulate a realistic example dataset

Set the multiple imputation methodology

impute_data_individual

Impute data for a single subject

Create imputed datasets

Calculate parametric confidence intervals

random_effects_expr

Construct random effects formula

Sample random values from the multivariate normal distribution

Create and validate a sample_list object

pval_percentile

P-value of percentile bootstrap

Pool analysis results obtained from the imputed datasets

pool_bootstrap_normal

Bootstrap Pooling via normal approximation

pool_bootstrap_percentile

Bootstrap Pooling via Percentiles

Combine estimates using Rubin's rules

Parallelise Lapply

print.imputation

Print imputation object

R6 Class for printing current sampling progress

Sample Patient Ids

is_num_char_fact

Is character, factor or numeric

rbmi: Reference Based Multiple Imputation

Print analysis object

Last Observation Carried Forward

Simulate intercurrent event

simulate_dropout

Simulate drop-out

Internal Pool Methods

Capture all Output

prepare_stan_data

Prepare input data to run the Stan model

recursive_reduce

recursive_reduce

Does a string contain a substring

split_imputations

Split a flat list of imputation_single() into multiple imputation_df()'s by ID

simulate_test_data

Create simulated datasets

scalerConstructor

R6 Class for scaling (and un-scaling) design matrices

Create object of sample_single class

Sort data.frame

Transform array into list of arrays

validate.analysis

Validate analysis objects

Validate draws object

validate.sample_single

Validate sample_single object

validate.is_mar

Validate is_mar for a given subject

validate_datalong

Validate a longdata object

Validate inputs for vars

validate_strategies

Validate user specified strategies

Print draws object

validate_analyse_pars

Validate analysis results

validate.stan_data

Validate a stan_data object

Generic validation method

transpose_samples

Transpose samples

validate.simul_pars

Validate a simul_pars object

transpose_results

Transpose results object

Set key variables

transpose_imputations

Transpose imputations

remove_if_all_missing

Remove subjects from dataset if they have no observed values

Barnard and Rubin degrees of freedom adjustment

Set simulation parameters of a study group.

validate.sample_list

Validate sample_list object

validate.references

Validate user supplied references

adjust_trajectories

Adjust trajectories due to the intercurrent event (ICE)

QR decomposition

R6 Class for a FIFO stack

Analysis of Covariance

Implements an Analysis of Covariance (ANCOVA)

adjust_trajectories_single

Adjust trajectory of a subject's outcome due to the intercurrent event (ICE)

Analyse Multiple Imputed Datasets

antidepressant_data

Antidepressant trial data

Applies delta adjustment