;Model Desc: base (minimal) model run 1 PROPORTIONAL residual error
;Project Name: phenocproject2nocooling
;Project ID: NO PROJECT DESCRIPTION
;WT AS COVARIATE WITH cL AND v WITH BASE MODEL AS SUGGESTED BY AB and follows chapter 11

$PROB RUN# pnc203 (PHENOBARBITAL POPULATION PK MODEL FOR non COOLED)
$INPUT C PT=DROP ID DATE=DROP DOSE=DROP TIME RT=DROP AMT DV EVID WT COOL=DROP APGT APGR=DROP  ALT AST PNA HIE 
$DATA 201.csv IGNORE=C 

$SUBROUTINES ADVAN1 
$PK
   
  TVCL = THETA(1) + THETA (3) * WT
 CL=TVCL*(1+ETA(1))
 TVVD=THETA(2) + THETA (4) * WT
 V = TVVD*(1+ETA(2))
 K=CL/V
 T=0.693/K
 S1=V

$ERROR ;PROPORTIONALERROR MODEL
   Y=F*(1+EPS(1)) 

$THETA
   (0, 0.0105) ;[CL, L/hr]
   (0, 1.05)   ;[V, L]
   (0, 0.0105) ; [CL WITH WT AS COVARIATE L/HR]
   (0, 1.05)   ;[V WITH WT AS COVARIATE L/HR]
  

$OMEGA
 .25 ;[P]
 .25 ;[P]

$SIGMA
 0.04;[A] ADDITIVECOMPONENT OF RESIDUAL VARIABILITY
 
$EST MAXEVAL=9999 PRINT=20 NOABORT POSTHOC MSF=pnc203.MSF


$TABLE ID  TIME  CL V WT K T PRED WRES  NOPRINT ONEHEADER FILE=pnc203.TAB
$SCATTERPLOT PRED VS DV UNIT
$SCATTERPLOT RES VS (PRED, WT, APGT)
$SCATTERPLOT WRES VS (PRED, WT, APGT)