print_model_code#

pharmpy.modeling.print_model_code(model)[source]#

Print the model code of the underlying model language

Parameters:: model (Model) – Pharmpy model

Examples

>>> from pharmpy.modeling import print_model_code, load_example_model
>>> model = load_example_model("pheno")
>>> print_model_code(model)
$PROBLEM PHENOBARB SIMPLE MODEL
$DATA pheno.dta IGNORE=@
$INPUT ID TIME AMT WGT APGR DV FA1 FA2
$SUBROUTINE ADVAN1 TRANS2
$ABBREV REPLACE ETA_CL=ETA(1)
$ABBREV REPLACE ETA_VC=ETA(2)

$PK
TVCL = THETA(1)*WGT
TVV = THETA(2)*WGT
IF(APGR.LT.5) TVV = TVV*(1 + THETA(3))
CL = TVCL*EXP(ETA_CL)
VC = TVV*EXP(ETA_VC)
V = VC
S1 = VC

$ERROR
Y = F + F*EPS(1)

$THETA  (0,0.00469307) ; POP_CL
$THETA  (0,1.00916) ; POP_VC
$THETA  (-.99,.1) ; COVAPGR

$OMEGA  0.0309626 ; IIV_CL
$OMEGA  0.031128 ; IIV_VC

$SIGMA  0.0130865  ; SIGMA

$ESTIMATION METHOD=1 INTERACTION MAXEVALS=99999
$COVARIANCE UNCONDITIONAL PRINT=E
$TABLE ID TIME DV CIPREDI PRED RES CWRES NOAPPEND NOPRINT ONEHEADER FILE=pheno.tab