;; 1. Based on: final
;; 2. Description: base_1-comp_add, linear elim
;; x1. Author: user
;; 3. Label:
$PROBLEM    PK
$INPUT      ID DAT2=DROP DAY TIME AMT RATE DV MDV EVID SEX AGE WT
TYPE=DROP SECTION=DROP ALT=DROP AST=DROP AKP=DROP ALB=DROP
AG=DROP TBA=DROP UREA TBIL DBIL=DROP URIC CRP=DROP CREA
CLCR=DROP EGFR=DROP TP=DROP GGT=DROP WBC RBC=DROP HGB PLT
NEUT1 NEUT2 TRIG=DROP CHOL=DROP HDLC=DROP LDLC=DROP
APOA=DROP APOB=DROP LFP=DROP QJSSXS=DROP FSM=DROP
CVVH=DROP SX=DROP
$DATA      elderlynoCVVH.csv IGNORE=#
$SUBROUTINE ADVAN1 TRANS2
$PK
;;; CLURIC-DEFINITION START
CLURIC = EXP(THETA(3)*(URIC - 262.52))
;;; CLURIC-DEFINITION END

;;; CL-RELATION START
CLCOV=CLURIC
;;; CL-RELATION END



TVCL = THETA(1) 

TVCL = CLCOV*TVCL
CL=TVCL* EXP(ETA(1))
TVV  = THETA(2)
V=TVV* EXP(ETA(2))
S1= V

$ERROR

IPRED = F
Y= F+EPS(1)



$THETA
(0, 4.22) ; CL
(0, 45.8) ; V
(-0.0046, -0.00116,0.0046) ; CLURIC1

$OMEGA 
0.155 ; IIV/BSV CL
0 FIX  ;  IIV/BSV V

$SIGMA 
4.63

$ESTIMATION METHOD=1 INTER MAXEVAL=9999 NOABORT SIG=3 PRINT=1 POSTHOC
MCETA=1
; Xpose
$COVARIANCE
$TABLE ID TIME DV MDV EVID IPRED CWRES CIWRES ETA1 ONEHEADER
NOPRINT
FILE=sdtabf02
$TABLE CL V ONEHEADER NOPRINT FIRSTONLY
FILE=patabf02
;$TABLE ID DAY TIME AMT RATE DV MDV EVID CL V ETA1 ETA2 SEX AGE WT TYPE SECTION ALT AST AKP ALB AG TBA UREA TBIL DBIL URIC CRP CREA CLCR EGFR TP GGT WBC RBC HGB PLT NEUT1 NEUT2 TRIG CHOL HDLC LDLC APOA APOB INR PT FIB DEJT LFP QJSSXS FSM CVVH SX ONEHEADER NOPRINT FILE=.fit
$ETAS       FILE=CLTBA1.phi