PRED capped

From: Markus Joerger Date: February 23, 2014 technical Source: mail-archive.com
dear community, I use a previously validated model of capecitabine in a new data-set. While the model converges successfully, PRED for all metabolites/compartments are 'capped' for PRED (at 1.6 for compartment 2, 2.4 for compartment 3, 3.0 for compartment 4 and 1.2 for compartment 5). At the same time, IPRED's are not 'capped', and GoF look nice. In a previous data-set using the same model, this phenomenon did not occur. Any suggestion is very wellcome, best regards Markus Control stream (capecitabine = compartment 2, DFCR = compartment 3, DFUR = compartment 4, 5FU = compartment 5): $PROB RUN# 001 $INPUT C ID TIME AMT DV CMT MDV EVID STUDY NR SEX AGE LIVERM WT BSA PS BILI Krea ASAT ALAT T85C A496G G1601A A1627G T1896C EXON14S A2846T C24T G1249A A4145G C4488T GSTT1 GSTP1 GSTM1 TSTR TS3CG TYMSPHENO TS3UTR MTHFR ERCC1 XRCC1 GSTP1B GASTRECT $DATA cap5.CSV IGNORE=C $SUBROUTINES ADVAN6 TOL=6 $MODEL NCOMPARTMENTS=5 COMP=(DOS) COMP=(CAP) COMP=(DFCR) COMP=(DFUR) COMP=(5FU) $PK TVKA=THETA(1) KA=TVKA*EXP(ETA(1)) ;absorption rate constant.1 TVV2=THETA(2) V2=TVV2*EXP(ETA(2)) ;CAP distribution volume.2 TVV3=THETA(3) V3=TVV3*EXP(ETA(3)) ;DFCR distribution volume.3 TVV4=THETA(4) V4=TVV4*EXP(ETA(4)) ;DFUR distribution volume.4 TVV5=THETA(5) V5=TVV5*EXP(ETA(5)) ;5FU distribution volume.5 TVQ1=THETA(6) Q1=TVQ1*EXP(ETA(6)) ;CAP > DFCR.6 TVQ2=THETA(7) Q2=TVQ2*EXP(ETA(7)) ;DFCR > DFUR.7 TVQ3=THETA(8) Q3=TVQ3*EXP(ETA(8)) ;DFUR > 5FU.8 TVCL1=THETA(9) CL1=TVCL1*EXP(ETA(9)) ;CL1-CAP.9 TVCL2=THETA(10) CL2=TVCL2*EXP(ETA(10)) ;CL2-5FU.10 ALAG1=THETA(11)*EXP(ETA(11)) ;ALAG1.11 F1=0.95 S2=V2/1000 S3=V3/1000 S4=V4/1000 S5=V5/1000 K23=Q1/V2 K34=Q2/V3 K45=Q3/V4 K20=CL1/V2 K50=CL2/V5 $ERROR CALLFL=0 IPRE=-3 FLG2=0 FLG3=0 FLG4=0 FLG5=0 IF(CMT.EQ.2)FLG2=1 IF(CMT.EQ.3)FLG3=1 IF(CMT.EQ.4)FLG4=1 IF(CMT.EQ.5)FLG5=1 IF(F.GT.0) IPRE=LOG(F) Y=LOG(F)+EPS(1)*FLG2+EPS(2)*FLG3+EPS(3)*FLG4+EPS(4)*FLG5 W=1 IRES=DV-IPRE IWRES=IRES/W $DES DADT(1) = -KA*A(1) DADT(2) = KA*A(1)-K23*A(2)-K20*A(2) DADT(3) = K23*A(2)-K34*A(3) DADT(4) = K34*A(3)-K45*A(4) DADT(5) = K45*A(4)-K50*A(5) $THETA (0,16) ;KA.1 (0,580) ;V-CAP.2 1 FIX ;V-DFCR.3 1 FIX ;V-DFUR.4 1 FIX ;V-5FU.5 (0,24) ;CAP > DFCR.6 (0,16) ;DFCR > DFUR.7 (0,8) ;DFUR > 5FU.8 (0,180) ;CL1-CAP.9 (0,75) ;CL2-5FU.10 (0,0.2) ;ALAG1.11 $OMEGA 0.14 ;KA.1 0 FIX ;V-CAP.2 0 FIX ;V-DFCR.3 0 FIX ;V-DFUR.4 0.14 ;V-5FU.5 0.14 ;CAP > DFCR.6 0.14 ;DFCR > DFUR.7 0.14 ;DFUR > 5FU.8 0.14 ;CL1-CAP.9 0.14 ;CL2-5FU.10 0 FIX ;ALAG1.11 $SIGMA 0.14 ;CMP2 0.14 ;CMP3 0.14 ;CMP4 0.14 ;CMP5 $ESTIMATION PRINT=5 MAXEVAL=5000 METHOD=1 INTER NOABORT POSTHOC SIGDIG=3 MSFO=MSF020B $COVARIANCE MATRIX=R $TABLE ID TIME AMT CMT IPRE IWRES CL1 CL2 V2 V3 V4 V5 NOPRINT ONEHEADER FILE=cap8.TAB -- Markus Joerger MD-PhD ClinPharm Medical Oncology&Clinical Pharmacology Department of Internal Medicine Kantonsspital St. Gallen Switzerland [email protected] Phone: +41-765591070 Fax: +41-714946325
Feb 23, 2014 Markus Joerger PRED capped
Feb 23, 2014 Bill Denney Re: PRED capped