NONMEM code for Mixture Model

From: Huali Wu Date: March 12, 2009 technical Source: mail-archive.com
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
Mar 12, 2009 Huali Wu NONMEM code for Mixture Model
Mar 12, 2009 Nidal . Alhuniti Re: NONMEM code for Mixture Model
Mar 13, 2009 Joachim Grevel Re: NONMEM code for Mixture Model
Mar 13, 2009 Joachim . Grevel Re: NONMEM code for Mixture Model
Mar 13, 2009 Huali Wu Re: NONMEM code for Mixture Model
Mar 16, 2009 Joachim Grevel Re: NONMEM code for Mixture Model
Mar 16, 2009 Joachim . Grevel Re: NONMEM code for Mixture Model