$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 RECS=ID

;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
CL = THETA(1) ;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