NONMEM code for Mixture Model
Dear NMusers:
I am working on dataset with high variability on clearance and patients in
different disease states seems have different distribution of clearance. So I
want to try the mixture model but I don't know how to do the coding. I listed
the code of my base model as below:
$PROB 1hr IV INFUSION SINGLE DOSE WITHOUT COVARIATES
$DATA data01.CSV IGNORE=C
$INPUT ID TIME DV AMT RATE MDV
$SUBROUTINES ADVAN9 TRANS1 TOL=5
$MODEL NPAR=3, NCOMP=1, COMP=(CENTRAL,DEFOBS)
$PK
V1 = THETA(1)*EXP(ETA(1))
VM = THETA(2)*EXP(ETA(2))
KM = THETA(3)
S1 = V1
$ERROR
Y=F+F*EPS(1)+EPS(2)
IPRED=F
IRES=DV-IPRED
IF(AMT.NE.0)W=1
IF(AMT.EQ.0)W=F
IWRES=IRES/W
$DES
C1 = A(1)/V1
DADT(1) = - C1*VM/(KM+C1)
$THETA (0, 4.47) (0, 155) (0, 1380)
$OMEGA BLOCK (2)
0.5
0.3 0.9
$SIGMA (0.01) (0.1)
$EST POSTHOC METHOD=1 MAXEVAL=9990 PRINT=5
$COV
$TABLE ID TIME DV AMT RATE V1 VM KM IWRES IPRED NOPRINT FILE=TAB4
ONEHEADER
$SCAT (RES WRES) VS TIME BY ID
Any suggestion will be highly appreciated!
Best regards,
Huali