NONMEM Users Network Archive

Hosted by Cognigen

Questions Regarding the calculation of CWRES with NONMEM version 6.2

From: Susan Zhou <susan.zhou>
Date: Thu, 20 Nov 2008 22:18:12 -0500

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

The NONMEM Users Network is maintained by ICON plc. Requests to subscribe to the network should be sent to: nmusers-request@iconplc.com.

Once subscribed, you may contribute to the discussion by emailing: nmusers@globomaxnm.com.