MCPModSimulation: MCPMod-based simulation of dose-finding trial designs

In MCPModPack: Simulation-Based Design and Analysis of Dose-Finding Trials

MCPMod-based simulation of dose-finding trial designs

Description

This function implements the simulation-based analysis of dose-finding clinical trials with normally distributed, binary and count endpoints using the MCPMod methodology. For more information, see the technical manual in the package's doc folder.

Usage

MCPModSimulation(endpoint_type, models, alpha, direction, 
                 model_selection, Delta, theta, sim_models, sim_parameters)

Arguments

endpoint_type	Character value defining the primary endpoint's type. Possible values: "Normal": Normally distributed primary endpoint. "Binary": Binary primary endpoint. "Count": Count-type primary endpoint.
models	List of candidate dose-response models with initial values of the non-linear model parameters. The package supports the following dose-response models: linear, quadratic, exponential, emax, logistic and sigemax. No initial value is required for the linear model, a single initial value is required for the quadratic, exponential and Emax models, and two initial values are required for the logistic and sigEmax models.
alpha	Numeric value defining the one-sided significance level (default value is 0.025).
direction	Character value defining the direction of the dose-response relationship. Possible values: "Increasing": A larger value of the treatment difference corresponds to a beneficial treatment effect. "Decreasing": A smaller value of the treatment difference indicates a beneficial treatment effect.
model_selection	Character value defining the criterion for selecting the best dose-response model. Possible values: "AIC": Akaike information criterion (AIC). "maxT": Most significant test statistic. "aveAIC": Weighted AIC-based criterion.
Delta	Numeric value specifying the treatment effect for identifying the target dose. The treatment effect is defined relative to the placebo effect. A positive value is required if direction = "Increasing" and a negative value is required otherwise.
theta	Numeric vector defining the overdispersion parameter in each trial arm (required only with count-type primary endpoints).
sim_models	List defining the assumed dose-response model and initial values of the non-linear parameters in the simulated trial. The package supports the following dose-response models: linear, quadratic, exponential, Emax, logistic and sigEmax. No initial value is required for the linear model, a single initial value is required for the quadratic, exponential and Emax models, and two initial values are required for the logistic and sigEmax models. The following components are required: "placebo_effect": Numeric value defining the placebo effect in the assumed dose-response model. "max_effect": Numeric vector defining the effects at the maximum dose in the assumed dose-response model. These effects are defined relative to the placebo effect. Positive values are required if direction = "Increasing" and negative values are required otherwise. "sd": Numeric vector defining the standard deviations of the response variable in each trial arm (required for normally distributed endpoints).
sim_parameters	List defining the design and simulation parameters in the simulated trial. The following components are required: "n": Integer vector defining the number of patients in each trial arm. "doses": Numeric vector defining the dose levels in each trial arm. "dropout_rate": Numeric value defining the dropout rate in the simulated trial (between 0 and 1). "nsims": Integer value defining the number of simulation runs. "go_threshold": Numeric value specifying the threshold for computing go probabilities, i.e., probabilities that the maximum effect for the best dose-response model corresponding to a significant contrast exceeds a pre-defined value. The threshold is defined relative to the placebo effect. A positive value is required if direction = "Increasing" and a negative value is required otherwise.

Value

The function returns an object of class ‘⁠MCPModSimulationResults⁠’. This object is a list with the following components:

input_parameters	a list containing the user-specified parameters, e.g, endpoint type, model selection criteria, etc.
selected_models	a logical vector defining the candidate dose-response models.
sim_results	a list containing the simulation results based on the assumed dose-response model, e.g., power, target dose estimates, etc.

A detailed summary of the simulation results can be generated using the SimulationReport function.

Author(s)

Alex Dmitrienko <[e-mail]>

See Also

MCPModAnalysis

Examples

  
# Simulation-based evaluation of dose-finding trials with a count endpoint

# Endpoint type
endpoint_type = "Count"

# Select the candidate dose-response models and initial values 
# of the non-linear model parameters (linear, quadratic, exponential, 
# emax, logistic and sigemax)
models = list(linear = NA, 
              quadratic = -0.5, 
              exponential = 0.3, 
              emax = 0.3, 
              logistic = c(0.5, 0.1), 
              sigemax = c(0.5, 5))

# One-sided Type I error rate
alpha = 0.025

# Direction of the dose-response relationship
direction = "increasing"

# Model selection criterion
model_selection = "AIC"

# The treatment effect for identifying the target dose 
# (this effect is defined relative to the placebo effect)
Delta = 2

# Vector of overdispersion parameters
theta = c(2, 2, 2, 2, 2)

# Select the assumed dose-response model and values of the non-linear model parameters
sim_models = list(emax = 1, 
                  placebo_effect = 3, 
                  max_effect = seq(from = 0, to = 3, by = 1))

# Simulation parameters 
# (go threshold is defined relative to the placebo effect)
sim_parameters = list(n = c(50, 50, 50, 50, 50),
                      doses = c(0, 0.05, 0.2, 0.6, 1),
                      dropout_rate = 0.05,
                      nsims = 1000,
                      go_threshold = 2)

# Perform an MCPMod-based simulation
results = MCPModSimulation(endpoint_type = endpoint_type, 
                           models = models, 
                           alpha = alpha, 
                           direction = direction, 
                           model_selection = model_selection, 
                           Delta = Delta,
                           theta = theta,
                           sim_models = sim_models,
                           sim_parameters = sim_parameters)

# Simple summary of the MCPMod simulation results
results

# Detailed summary of the MCPMod simulation results (remove tempfile)
SimulationReport(results, 
  "MCPMod simulation summary (Count endpoint)", 
  tempfile("MCPMod simulation summary (Count endpoint).docx", fileext=".docx"))

  

MCPModPack documentation built on May 31, 2023, 5:23 p.m.