Questions Regarding the calculation of CWRES with NONMEM version 6.2

1 messages 1 people Latest: Nov 21, 2008
Dear Susan: We have no problem to run CWRES in xpose4 with NONMEM 6.2. But we did not try the same code with NONMEM 6.1. We made some changes to the control stream for nm 6.1 before running nm 6.2 as follows: 1) remove the "LAST block 2) change the $TABLE for FILE=cwtabx.deriv to take out all the "COM(x)=" Best regards, Sam Liao Pharmax Research > 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 > > >