analyzeSP23Design: Analyses the results of running a design. If a trial history,...

In sp23design: Design and Simulation of Seamless Phase II-III Clinical Trials

Analyses the results of running a design. If a trial history, such as the result of the function executeSP23Design or a history of the actual conduct of a single trial is provided, it returns the analysis results.

Description

Produces analysis results from the run of a single trial or a number of simulations.

Usage

analyzeSP23Design(sp23Design, trialHistory = NULL, data = NULL,
  col=c("red", "red", "brown", "brown"), lty=c(1,2,1,2))

Arguments

sp23Design	The design object typically produced by calling generateSP23Design.
trialHistory	Typically the result of calling executeSP23Design which is a single data frame, or, the result produced by calling exploreSP23Design which is a list of data frames, one for each simulation.
data	This is only used when the argument trialHistory is a single data frame, in which case, it should be the data generated in the clinical trial.
col	Colors used for the survival plots
lty	Line types for the survival plots

Details

If trialHistory is a single data frame, the naive response estimates and a survival plot are produced. Otherwise, the counts of the number of times the various hypothesis are rejected and other details are returned.

Value

A list of two items named responseSummary and designSummary If trialHistory is a single data frame, the naive response estimates are returned in responseSummary and a survival plot is produced. Otherwise, the counts of the number of times the various hypothesis are rejected and other details are returned in designSummary.

Author(s)

Mei-Chiung Shih, Balasubramanian Narasimhan, Pei He

References

Lai, Tze Leung and Lavori, Philip W. and Shih, Mei-Chiung. Sequential Design of Phase II-III Cancer Trials, Statistics in Medicine, Volume 31, issue 18, p.1944-1960, 2012.

Examples

## Not run: 
  catn <- function(...) cat(..., "\n")

  trialParameters <- list(minimumNumberOfEvents = 20,
                          minimumIncreaseInV = 0.2,
                          numberRecruitedEachYear = c(80, 120, 160, 160),
                          followupTime = 3,
                          adminCensoringTime = 7,
                          interimLookTime = c(1, 2, 3, 5, 7),
                          type1ErrorForResponse = 0.05,
                          type2ErrorForResponse = 0.01,
                          glrBoundarySidedness = "one", # one sided or two-sided
                          type1Error = 0.05,
                          type2Error = 0.10,
                          epsType1 = 1/3,
                          epsType2 = 1/3)

  trueParameters <- list(p0 = 0.3,
                         p1 = 0.3,
                         pdiffHyp=0.3,
                         theta = list(
                             alpha = 0,
                             beta = 0,
                             gamma = 0),
                         baselineLambda = 0.35,
                         etaHyp = 0.25)
  rngSeed <- 9872831
  sp23Design <- generateSP23Design(trueParameters, trialParameters)
  print(sp23Design)
  trialHistory <- exploreSP23Design(sp23Design, numberOfSimulations=25, rngSeed=rngSeed)
  result <- analyzeSP23Design(sp23Design, trialHistory)$designSummary
  catn("numberOfTimesH0RIsRejectedAtFirstLook",result[["numberOfTimesH0RIsRejectedAtFirstLook"]])
  catn("numberOfTimesH0RIsRejected", result[["numberOfTimesH0RIsRejected"]])
  catn("numberOfTimesStoppedForFutility", result[["numberOfTimesStoppedForFutility"]])
  catn("numberOfTimesH0SIsAccepted", result[["numberOfTimesH0SIsAccepted"]])
  catn("numberOfTimesH0SIsRejected", result[["numberOfTimesH0SIsRejected"]])
  catn("numberOfTimesFutilityDecidedAtLastLook",result[["numberOfTimesFutilityDecidedAtLastLook"]])
  catn("numberOfTimesTrialEndedAtLook", result[["numberOfTimesTrialEndedAtLook"]])
  catn("avgExitTime", result[["avgExitTime"]])

## End(Not run)

sp23design documentation built on April 19, 2022, 5:08 p.m.