From: Sam Liao <*sliao*>

Date: Fri, 21 Nov 2008 09:38:54 -0500

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
*

*>
*

*>
*

*>
*

Received on Fri Nov 21 2008 - 09:38:54 EST

Date: Fri, 21 Nov 2008 09:38:54 -0500

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

Received on Fri Nov 21 2008 - 09:38:54 EST