Description Usage Arguments Details Author(s) Examples
Use final parameter estimates of a model to simulate data, calculate statistics for visual predictive check, and collect the results.
Important Note: in the current version of phxvpc.sim, the model.file (.mdl) must follow a specific format.
See the Details section: Instruction on model file.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 |
path |
System directory for location of the model run folder |
vpcpath |
System directory where vpc simulation results will be stored. If NULL(default), a subfolder vpc n will be created under the model run folder (path), in which n is a number. If vpc n subfolder already exists, n will automatically increase by 1 until no identical existing folder can be found. |
ivar |
Independent variable (x-axis). Default is "t", which is the independent variable when the structure model is set to any of Phoenix build-in PK models. Other options include "TAD", "PRED", or any user specified input. If user specified input is used, this variable must be in data set and defined in cols.file for model fit. |
nsim |
An integer value specifying the number of simulation replicates. Default is 200. |
pstrat |
A character vector that provides a maximum of 3 stratification variables. The vpc will stratify the data on unique values of the specified variable, and perform separate analyses on each set. Variables must be in data set and defined in cols.file for model fit. (Refer to Phoenix NLME manual for additional details regarding cols.file) |
setseed |
An integer value that provides a fixed seed for random resampling. If omitted, a random seed will be assigned automatically. |
pred.corr |
Optional. A character string specifying the options of how prediction correction is performed on dependent variable values before computing vpc results. Options are "proportional" or "additive", which can be shorten, with the minimal length of "prop" and "add", respectively. Not case sensitive. Default is NULL. |
var.corr |
A logical value. If TRUE, perform variability correction on dependent variable values before computing VPC results. Default is FALSE. Note: this option is only functional when the argument pred.corr is used. |
pi |
Default is c(0.025,0.5, 0.975). A vector of values that describe the prediction percentile that should be calculated. The prediction percentiles are displayed by the Q.pred.line option in |
pi.ci |
Default is c(0.025, 0.975). A vector of two values that describe the confidence interval of the prediction percentile that should be calculated. These values are used to specify the boundaries of the CI.Q.pred option in |
bin.option |
A character string that provide the binning options for simulation. Default is NULL. Alternative options are "K-means", "centers", and "boundaries", which can be shorten, with the minimal length of "k", "cent", and "bound", respectively. Not case sensitive. |
bin.bound |
A numeric vector that provides binning boundaries when bin.option "boundaries" is used. This option will be ignored when other bin.option is selected. Default is NULL. |
bin.center |
A numeric vector that provides centers for all bins when bin.option "centers" is used. This option will be ignored when other bin.option is selected. Default is NULL. |
clean |
A logical value. Default is FALSE. If TRUE, the NLME executable file in the vpcpath folder will be deleted after model execution, the results cannot be updated/modified. If FALSE, it makes possible to apply other estimation/simulation functions in the same folder. |
hold |
A logical value that determines command window behavior. Currently, it is not implemented. |
modsp.file |
An output file from |
out.file |
An output file from |
phxnlme
must be executed before using phxvpc.sim
.
Instruction on model file Model file (.mdl) must follow a certain format for phxvpc.sim to use final model estimates for simulation. First, the following blocks (fixed, ranef, stparm, error, observe) must be in the order of fixed->ranef(optional)->stparm->error->observe. Second, if the ranef block only contains diagnal elements, each element must be in a separate row. Examples are as following:
##correct format of model file:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | test(){
covariate(FEMALE)
covariate(DOSE) ##no restriction for covariate
fixef(
tvE0 = c(, 20,)
tvEMAX = c(, 120,)
ED50MALE = c(, 15,)
ED50FACTOR = c(, 1,)
)
ranef(
## each element must be in a separate row
diag(nE0, nEMAX, nED50) = c(1,
1,
1)
)
stparm(
E0 = tvE0 * exp(nE0)
EMAX = tvEMAX * exp(nEMAX)
ED50 = ED50MALE * exp(nED50) * ED50FACTOR^FEMALE
)
E = E0 + EMAX * DOSE / (DOSE + ED50)
error (EPS1 = 10)
observe(EOBS = E + EPS1)
}
|
Shuang Liang
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 | ## Note: .mdl file must be in the format as specified in the details section.
## Alternatively, use simmodel for VPC simulation.
## Note that pheno2.mdl needs to be modified to follow the specified format prior to
## running examples below
if(!is.null(checkphxnlme(testchk=TRUE))){
path="C:/Program Files (x86)/Pharsight/Phoenix/application/Examples/NLME Command Line/Model 5"
model.file="pheno2.mdl"
cols.file="colspheno2.txt"
data="pheno2.csv"
## Run model fit
phxnlme(path=path,model.file=model.file,cols.file=cols.file,data=data)
## VPC simulation
#phxvpc.sim(path)
}
## Change confidence interval of prediction percentiles
#phxvpc.sim(path, pi.ci=c(0.05, 0.95))
## Bin by boundaries
#phxvpc.sim(path, bin.option="bound", bin.bound=c(0, 0.5, 4, 8, 12))
## Note: lyon04.mdl needs to be modified to specified format prior to running example below
## For models not using build-in PK structure model
## Run model fit
if(!is.null(checkphxnlme(testchk=TRUE))){
path="C:/Program Files (x86)/Pharsight/Phoenix/application/Examples/NLME Command Line/Model 1"
model.file="lyon04.mdl"
cols.file="COLS04.txt"
data="EMAX02.csv"
phxnlme(path=path,model.file=model.file,cols.file=cols.file,data=data)
#phxvpc.sim(path, ivar="DOSE")
## Stratified VPC, 1 covariate
#phxvpc.sim(path, ivar="DOSE", pstrat="FEMALE")
## Stratified VPC, 3 covariates (covariates must be included in bot data and cols.file)
#phxvpc.sim(path, pstrat=c("SEX", "AGE", "DOSE"))
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.