getDesignInverseNormal(informationRates = c(0.2, 0.6, 1)) |>
getSimulationMultiArmRates(
doseLevels = c(2.1, 3, 5.25, 9), # New!
typeOfShape = "sigmoidEmax", # Define dose-response relationship
gED50 = 2, # Dose with half of the maximum effect
slope = 0.5, # Slope of the dose-response curve
activeArms = 4,
plannedSubjects = c(10, 30, 50),
piControl = 0.1,
piMaxVector = seq(0.35, 0.45, 0.1),
intersectionTest = "Sidak",
conditionalPower = 0.8,
minNumberOfSubjectsPerStage = c(10, 4, 4),
maxNumberOfSubjectsPerStage = c(10, 100, 100),
maxNumberOfIterations = 10
)