RE: problem with acquiring post hoc estimates
From: Leonid Gibiansky lgibiansky@emmes.com
Subject: RE: [NMusers] problem with acquiring post hoc estimates
Date: Tue, May 18, 2004 1:41 pm
Kristel
I do not know why, but with your control stream and data, NONMEM was using
F in compartment 1 for EVID =2 (only!) even when CMT was equal to 2. You
may use values that correspond to EVID=0 (3d and 4th lines of your data
file) or modify the code as below and use C1, C2 or IPRED/PRED (that are
computed explicitly using A1 and A2 rather than F).
Leonid
$ERROR
C1 = A(1)/S1
C2 = A(2)/S2
IPRED=A(1)/S1
IF(CMT.EQ.2) IPRED=A(2)/S2
IF(CMT.EQ.1)THEN
Y=IPRED+IPRED*THETA(7)*ERR(1)+THETA(8)*ERR(2)
W=SQRT(IPRED*IPRED*THETA(7)**2+THETA(8)**2)
ELSE
Y=IPRED+IPRED*THETA(9)*ERR(3)+THETA(10)*ERR(4)
W=SQRT(IPRED*IPRED*THETA(9)**2+THETA(10)**2)
ENDIF
IRES=DV-IPRED
IWRES=IRES/W
$TABLE .... CMT C1 C2 IPRED ..
_______________________________________________________