Re: Questions Regarding the calculation of CWRES with NONMEM version 6.2

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