RE: Ln(DV)
From:"Bachman, William"
Subject: RE: [NMusers] Ln(DV)
Date:Thu, 13 Mar 2003 14:15:39 -050
Have you constrained Q to be positive?
example code for constraining PK parameters to be positive:
$SUBROUTINES ADVAN4 TRANS4
$PK
TVCL=THETA(1)
CL=TVCL*EXP(ETA(1))
TVV2=THETA(2)
V2=TVV2*EXP(ETA(2))
TVQ=THETA(3)
Q=TVQ*EXP(ETA(3))
TVv3=THETA(4)
V3=TVv3*EXP(ETA(4))
TVKA=THETA(5)
KA=TVKA*EXP(ETA(5))
S2=V2
;TLAG1=THETA(6)
;ALAG1=TLAG1*EXP(ETA(6))
$ERROR
IPRED=F
IF(IPRED.EQ.0) IPRED=EXP(-9)
Y=LOG(IPRED)+ERR(1)
$THETA
(0, 4.6) ;[CL]
(0, 6.5) ;[V2]
(0, 3.9) ;[Q]
(0, 3.0) ;[V3]
(0, 0.25);[KA]
;(0,0.2) ;[ALAG1]