dfpk (version 1.0)

scenarios: Generate New Scenarios

Description

This function can be used to generate several scenarios.

Usage

scenarios(param_pk, omega2, omega_a, sigma2, preal, limit_tox, time_sampling, N)

Arguments

param_pk

Subject's pharmacokinetics(PK) parameters.

omega2

The standard deviation of the log-normal distribution that the PK parameters clearance(CL) and volume(V) come from; it's equal to 70% or 30% in different scenarios.

omega_a

The standard deviation of the log-normal distribution which the sensitivity parameter $\alpha$ comes from.

sigma2

The additive or proportional error.

preal

The true probabilities of toxicity.

limit_tox

The threashold on toxicity.

time_sampling

Time sampling.

N

The total number of patients.

References

Ursino, M., et al, (2016) Dose-finding methods using pharmacokinetics in small populations (under review).

See Also

sim

Examples

#####################################
######### Scenario Inputs  ##########
#####################################

limit_tox <- 10.96                  
N <- 42                               
# number of patients
param_pk <- c(2,10,100)   # PK parameters ka,CL,V
omega2 <- 0.7             # Inter-individual
omega_a <- 0              # IIV   
preal <- c(0.001,.05,0.1,.2,0.35,0.45)      
# true prob of toxicity
time_sampling <- seq(0,24,length.out=48)    # time sampling
sigma2 <- c(rep(0.2,length(time_sampling)),rep(0.2,length(time_sampling)))           
# sigma2: Additive or proportional error

#############################################
############ Test the scenario ##############
#############################################

gen.scen <- scenarios(param_pk,omega2,omega_a,sigma2,preal,limit_tox,time_sampling, N) 
gen.scen