Autoinduction

From: "Sreenivasa Rao Vanapalli" <svanapal@blue.weeg.uiowa.edu> Subject: Autoinduction Date: Mon, 14 May 2001 17:33:55 -0500 Dear NMUsers I'm trying to fit multiple oral dose data with autoinduction. The data fits well to a two compartment model with first order absorption. But when I'm trying to fit the data to a two compartment with hypothetical enzyme compartment, I'm getting error message 'floating-point error:divide by zero' The control file I'm using is, $PROBLEM MOXIDECTIN POOLED DATA $INPUT ID DOSE=AMT TIME CP=DV WT CMT SEX EVID $DATA MOXIMULTIEDIT $SUBROUTINES ADVAN6 TOL=5 $MODEL NCOMPARTMENTS=4 COMP=(DEPOT,INITIALOFF,DEFDOSE) COMP=(CENTRAL,DEFOBS,NOOFF) COMP=(ENZ) COMP=(PERIPH) $PK TVCLIN=THETA(1) CLIN=TVCLIN*EXP(ETA(1)) KA=THETA(2)*EXP(ETA(2)) Q=THETA(3)*EXP(ETA(3)) TVVD=THETA(4) V2=TVVD*EXP(ETA(4)) V3=THETA(5)*EXP(ETA(5)) KENZ=THETA(6)*EXP(ETA(6)) IC50=THETA(7)*EXP(ETA(7)) $DES DADT(1)=-KA*A(1) DADT(2)=KA*A(1)-CLIN*DV DADT(3)=KENZ-KENZ*A(3)*(1-(CP/(CP+IC50))) $ERROR DEL=0 IF (F.EQ.0) DEL=1 W=F+DEL Y=F*(1+ERR(1))+ERR(2) IPRED=F IRES=DV-IPRED IWRES=IRES/W CL=CLIN/A(3) $THETA (0,.3,10) ;CLIN (0,.8,100) ;ABSOR (0,1.5,10) ;Q (0,13,100) ;V2 (0,100,500);V3 (0,.001,100) ;KENZ (0,500,1000) ;IC50 $OMEGA .002 .0001 .0001 .0001 .001 .001 .001 $SIGMA .05 .01 $ESTIMATION MAXEVAL=9999 POSTHOC $COVARIANCE $TABLE ID WT TIME V2 CL KA $SCAT CP VS TIME $SCAT PRED VS TIME $SCAT (RES WRES) VS TIME $SCAT WRES VS CP UNIT $SCAT PRED VS CP UNIT Any help would be appreciated. Regards Sreenivasa Rao Vanapalli, Ph.D, Janssen Postdoctoral Research Scholar, S411 PHAR, College of Pharmacy, University of Iowa, Iowa City, IA-52242 319-353-5157 (Office) 319-337-2687 (Home)