PK/PD model with a categorical variable
Dear NMusers,
I tried to fit a PK/PD model with a categorical variable (0/1).
The PK part was estimated in a previous run (2 steps process with fixed theta,
omega and sigma for PK).
The hypothesis for the categorical variable is a value of 1 with high
concentrations and 0 otherwise.
Something is wrong in the $ERROR as I observed that the F values used for
computing EXPP were not realistic (far to high).
Any advice will be much appreciate.
Best regards
Nicolas
$PROB PK/PD
$INPUT ID TIME DAT1=DROP DV AMT EVID MDV TYPE MTT PERIOD DOSE GROUPE RATE
$DATA ..\data.csv
IGNORE=#
$ESTIM MAXEVAL=9000 SIGDIG=2 PRINT=10 NOABORT METHOD=CONDITIONAL INTERACTION
LAPLACIAN
$COV
$THETA (440 FIXED) ;POPV1 ;L
$THETA (0.852 FIXED) ;POPK10 ;1/H
$THETA (2.13 FIXED) ;POPK12 ;1/H
$THETA (0.069 FIXED) ;POPK21 ;1/H
$THETA (3.98 FIXED) ;POPK13 ;1/H
$THETA (1.67 FIXED) ;POPK31 ;1/H
$THETA (2.41 FIXED) ;COV_GR
$THETA (0.677 FIXED) ;COV_DO
$THETA (.1) ;int
$THETA (.1) ;slp
$OMEGA 0 FIXED ; ETA_V1
$OMEGA 1.12 FIXED ; ETA_K10
$OMEGA 0 FIXED ; ETA_K12
$OMEGA 0 FIXED ; ETA_K21
$OMEGA 0 FIXED ; ETA_K13
$OMEGA 0 FIXED ; ETA_K31
$OMEGA 0.717 FIXED ; ETA_F1
$SIGMA 0.126 FIXED ; ERR1
$SUBR ADVAN11 TRANS1
$PK
D1=0.1667
GRP=0
IF(GROUPE.EQ.2) GRP=1
GRP_DO=0
IF(DOSE.GT.15) GRP_DO=1
TVF1=THETA(7)**GRP * THETA(8)**GRP_DO
F1=TVF1*EXP(ETA(7))
V=THETA(1)*V1*EXP(ETA(1))
K=THETA(2)*EXP(ETA(2))
K12=THETA(3)*EXP(ETA(3))
K21=THETA(4)*EXP(ETA(4))
K13=THETA(5)*EXP(ETA(5))
K31=THETA(6)*EXP(ETA(6))
S1=V
CL=K*V
INT=THETA(9)
SLP=THETA(10)
$ERROR
DEL=0
IF (F.EQ.0) DEL=0.0001
W=F+DEL
EXPP=INT+F*SLP
IF(EXPP.GT.30.0) EXPP=30.0
IF (TYPE.EQ.0) THEN
F_FLAG=0
Y=F + W*ERR(1)
ENDIF
IF (TYPE.EQ.1) THEN
F_FLAG=1
AA=DEXP(EXPP)
BB=1+AA
Y=DV*AA/BB+(1-DV)/BB
ENDIF
$TABLE ID TIME AMT DOSE GROUPE DV CWRES PRED EVID MDV TYPE
FILE=r
#ID TIME DATE DV AMT EVID MDV TYPE PERIOD DOSETH GROUPE RATE 1 10:04:00
07/02/2014 0 89.9984099 1 1 0 1 30 2 -2 1 10:04:01 07/02/2014 0.00248052 0 0 0
0 1 30 2 0 1 10:04:01 07/02/2014 0 0 0 0 1 1 30 2 0 1 10:30:00 07/02/2014
0.05406265 0 0 0 0 1 30 2 0 1 10:30:00 07/02/2014 1 0 0 0 1 1 30 2 0 1 10:45:00
07/02/2014 0.04515821 0 0 0 0 1 30 2 0 1 10:45:00 07/02/2014 0 0 0 0 1 1 30 2 0
1 11:50:00 07/02/2014 0.03370965 0 0 0 0 1 30 2 0 1 11:50:00 07/02/2014 0 0 0 0
1 1 30 2 0 1 12:25:00 07/02/2014 0.01590078 0 0 0 0 1 30 2 0 1 12:25:00
07/02/2014 0 0 0 0 1 1 30 2 0 1 14:27:00 07/02/2014 0.00477023 0 0 0 0 1 30 2 0
1 14:27:00 07/02/2014 0 0 0 0 1 1 30 2 0 1 16:31:00 07/02/2014 0.00314835 0 0 0
0 1 30 2 0 1 16:31:00 07/02/2014 0 0 0 0 1 1 30 2 0 1 18:31:00 07/02/2014
0.00298935 0 0 0 0 1 30 2 0 1 20:36:00 07/02/2014 0.00219431 0 0 0 0 1 30 2 0 1
22:30:00 07/02/2014 0.0020353 0 0 0 0 1 30 2 0 1 10:25:00 08/02/2014 0.00162188
0 0 0 0 1 30 2 0 1 10:25:00 08/02/2014 0 0 0 0 1 1 30 2 0 1 09:49:00 21/02/2014
0 26.3952934 4 1 0 2 10 2 -2 1 09:49:01 21/02/2014 0.00159008 0 0 0 0 2 10 2 0
1 09:49:01 21/02/2014 0 0 0 0 1 2 10 2 0