$PROBLEM PROJECT singledose i.v infusion ;DATE 6-2-04 PROGRAMMER:XXXX
;UNITS: Time=hour, Concentration=ng/ml
;Dose=100mg, Clearance=L/hr, Volume = L
$INPUT ID TIME CONC=DV AMT RATE MDV
$DATA cs3_ivinfest.csv IGNORE=C 

;RATE specifies the infusion rate, for a dose of 100mg, rate of 16.7mg/hr
;means infusion for 6hrs
$SUBROUTINE ADVAN1 TRANS2 ;A one compartment model from PREDPP library
$PK
FLG1=0
FLG2=0
FLG3=0
FLG4=0
FLG5=0
FLG6=0
FLG7=0
FLG8=0
FLG9=0
FLG10=0
FLG11=0
FLG12=0
FLG13=0
FLG14=0
FLG15=0
FLG16=0
FLG17=0
FLG18=0
FLG19=0
FLG20=0

IF (ID .EQ. 1) FLG1=1
IF (ID .EQ. 2) FLG2=1
IF (ID .EQ. 3) FLG3=1
IF (ID .EQ. 4) FLG4=1
IF (ID .EQ. 5) FLG5=1
IF (ID .EQ. 6) FLG6=1
IF (ID .EQ. 7) FLG7=1
IF (ID .EQ. 8) FLG8=1
IF (ID .EQ. 9) FLG9=1
IF (ID .EQ. 10) FLG10=1
IF (ID .EQ. 11) FLG11=1
IF (ID .EQ. 12) FLG12=1
IF (ID .EQ. 13) FLG13=1
IF (ID .EQ. 14) FLG14=1
IF (ID .EQ. 15) FLG15=1
IF (ID .EQ. 16) FLG16=1
IF (ID .EQ. 17) FLG17=1
IF (ID .EQ. 18) FLG18=1
IF (ID .EQ. 19) FLG19=1
IF (ID .EQ. 20) FLG20=1




CL = THETA(1)*FLG1+THETA(2)*FLG2+THETA(1)*FLG1+THETA(1)*FLG1+THETA(1)*FLG1+THETA(1)*FLG1+THETA(1)*FLG1+THETA(1)*FLG1 ;Clearance in L/hr
V = THETA(2) ;Volume of distribution in L
S1 =V/1000 ;scaling factor to match concentration in ng/ml
$ERROR
IPRED=F
Y=F+F*ERR(1)+ERR(2) ;Combined residual error model

$THETA (5,25) ;POPCL
$THETA (10,86) ;POPV
$OMEGA 0.02
$OMEGA 10

;$OMEGA 0.09 ;BSVCL ;30% BSV for Clearance
;$OMEGA 0.09 ;BSVV ;30% BSV for volume
;$SIGMA 0.02 ;ERRCV ;Proportional error = 15%
;$SIGMA 10 ;ERRSD ;Additive error =3ng/ml
$ESTIMATION METHOD=0 MAXEVAL=9999 PRINT=10 NOABORT 
$COVARIANCE SPECIAL PRINT=E
$TABLE ID TIME DV IPRED AMT CL V NOPRINT ONEHEADER FILE= cs3_ivinfestind.fit
INCLUDE CS3_IVINFESTIND2.txt 19