Phase12Compare (version 1.5)

SimSpsoContour: Simulates trial replications from the SPSO model with desirability optimality function.

Description

Simulates replications from a Bayesian adaptive phase 12 clinical trial design using the SPSO model and a optimality function based on the desirability tradeoff contour.

Usage

SimSpsoContour(
  NSims,
  PE,
  PT,
  corET,
  Nmax,
  cohort,
  Contour,
  CutE,
  CutT,
  AcceptE,
  AcceptT,
  HypermeansE,
  HypermeansT,
  Hypervars
)

Arguments

NSims

Number of trial simulations to run.

PE

True efficacy probability for each dose.

PT

True toxicity probaiblity for each dose.

corET

Correlation parameter between efficacy and toxicity probability.

Nmax

Maximum trial sample size.

cohort

Patient cohort size.

Contour

Contour vector for desirability function. Contains in order: (pi_1,E,pi_2,T,pi_3,E,pi_4,T).

CutE

Cutoff for efficacy probability acceptability. Dose-efficacy probabilities must be larger than this for patient assignment.

CutT

Cutoff for toxicity probability acceptability. Dose-toxicity probabilities must be smaller than this for patient assignment.

AcceptE

Posterior probability threshold for efficacy acceptability.

AcceptT

Posterior probability threshold for toxicity acceptability.

HypermeansE

Hypermeans for dose-specific efficacy parameters.

HypermeansT

Hypermeans for dose-specific toxcity parameters.

Hypervars

Hypervariances needed for the SPSO model. Contains, in order (sigma_0^2, sigma_mu^2,tau).

Value

A list with the first entry corresponding to a matrix with: (1) True toxicity probabilities at each dose, (2) True efficacy probabilities at each dose, (3) True Desirability of each dose, (4) Optimal dose selection probability, (5) Average sample size of patients treated at each dose. The second entry of the list contains a vector with rows corresponding to (1) the true binary bivariate correlation between efficacy and toxicity, (2) Stopping probability of the trial, (3) Average number of efficacy events, (4) Average number of toxicity events, and (5) Delta.

Examples

# NOT RUN {
library(mvtnorm)
##True toxicity probability
PT=c(.05,.10,.15,.20,.30)
##True Efficacy Probability
PE=c(.2,.4,.6,.65,.7)
#True Correlation
corET=.5
##Number of simulations
NSims=1 ##Increase this when using
##Hypermeans and hypervariances
HypermeansE=c(-1.189, -0.357,  0.360,  0.546,  0.743)
HypermeansT=c(-2.325, -1.811, -1.464, -1.189, -0.740)
Hypervars=c(1,16,1)
##Trial Parameters
##Cohort Size, N^F and N_ET
cohort=3
##Starting Dose
DoseStart=1
Nmax=30 ##Max Sample Size
#Acceptability Criterion
CutE=.3
CutT=.4
##Limits on acceptability
AcceptE=.1
AcceptT=.1
##Contour vector
Contour = c(.35, .75,.7,.4)
RESULTS=SimSpsoContour(NSims, PE, PT, corET, Nmax, cohort, Contour,
CutE, CutT, AcceptE, AcceptT,HypermeansE, HypermeansT, Hypervars)
RESULTS
# }