# List of all parameters
parameters = list()
# Endpoint type
parameters$endpoint_type = "Normal"
# Direction of beneficial effect
parameters$direction = "Lower"
# Number of dose-placebo comparisons
parameters$n_comparisons = 2
# Number of endpoints
parameters$n_endpoints = 2
# Number of enrolled patients (placebo, Dose L and Dose H)
parameters$sample_size = c(120, 120, 120)
# Patient dropout rate
parameters$dropout_rate = 0.10
# Endpoint information (rows corresponds to endpoints and
# columns corresponds to dose-placebo comparisons)
parameters$control_mean = c(0, 0)
parameters$treatment_mean = matrix(c(-0.30, -0.35,
-0.35, -0.40), 2, 2, byrow = TRUE)
parameters$control_sd = c(1, 1)
parameters$treatment_sd = matrix(c(1, 1,
1, 1), 2, 2, byrow = TRUE)
# Endpoint correlation matrix
parameters$endpoint_correlation = matrix(c(1, 0.5,
0.5, 1), 2, 2)
# Component procedure to be used in the gatekeeping procedure
parameters$mult_test = "Hochberg"
# Mixture method used in the gatekeeping procedure
parameters$mult_method = "Modified"
# Truncation parameters in the gatekeeping procedure
parameters$mult_test_gamma = c(0.8, 1)
# Overall one-sided Type I error rate
parameters$alpha = 0.025
# Number of simulations
parameters$nsims = 100
# Number of cores for parallel calculations
parameters$ncores = 1
# Run simulations to perform power calculations
results = MultAdj(parameters)
# Generate a simulation report (remove tempfile)
GenerateReport(results,
tempfile("MultAdj Multiple endpoints and dose-placebo comparisons.docx",
fileext=".docx"))
Run the code above in your browser using DataLab