Hi, Everybody,
I have used Andrew Hooker's code to calculate CWRES with NONMEM 6.1
before (the NONMEM verbatim code and the R code for calculation). Now
our NONMEM has been upgraded to 6.2, and I am having problems with the
calculation of CWRES using the same code. In my NONMEM code, for the
error structure I used the following code for the combined error model
$SIGMA
1 FIXED
...
$PK
...
PERR=THETA(5)
AERR=THETA(6)
...
$ERROR
Y=F+(F*PERR + AERR)*EPS(1)
with NONMEM 6.1 and 6.2 using the exact same control file (with the
verbatim code for generating the matrices used for CWRES calculation
included), the outputs of the matrices are different. When using NONMEM
6.1, the output sigma and sesigma matrices used for calculation of CWRES
are non zero, and the NONMEM 6.2 outputs for these matrices are zero
matrices. The zero matrices seems to be one of the problems in
calculation of CWRES.
I saw in one thread that Tom Ludden mentioned that NONMEM 6.2 can output
the matrices used for calculation of CWRES directly in table format
without using the verbatim code. Has anybody used that function in 6.2
and would be willing to share the code.
Any help would be much appreciated.
Regards,
Susan
Received on Thu Nov 20 2008 - 22:18:12 EST
