RE: [NMusers] different COMRES statements

Hi Frederik,

 

It is possible to compute the CWRES using the $PRIOR functionality in
NONMEM. However, it's a little tricky. The problem is this code:

 

$INFN

     IF (ICALL.EQ.3) THEN

         OPEN(50,FILE=CWTAB1.EST)

         WRITE(50,*) 'ETAS'

         DO WHILE(DATA)


               IF (NEWIND.LE.1) WRITE (50,*) ETA


             ENDDO


             WRITE(50,*) 'THETAS'

             WRITE(50,*) THETA

             WRITE(50,*) 'OMEGAS'

             WRITE(50,*) OMEGA(BLOCK)

             WRITE(50,*) 'SIGMAS'

             WRITE(50,*) SIGMA(BLOCK)

     ENDIF

 

Which writes out all the thetas and omegas and sigmas including the prior
stuff (which we don't need for the CWRES calculation). What I do is just
remove, by hand, the extra values in the "cwtab1.est" file to have only the
values for the etas, thetas, omegas and sigmas for the data and not the
prior (only the "first" thetas, omegas and sigmas). After that, calculating
the CWRES using XPOSE 4 should not be a problem.

 

For the COMRES values: If you need x COMRES values for the CWRES
calculation and y COMRES values for some other portion of your NONMEM model
then you just need to use the statement

 

$ABB COMRES= z

 

Where z = x+y.

 

Hope this helps!

 

-Andy

 

 

Andrew Hooker, Ph.D.

Assistant Professor of Pharmacometrics

Div. of Pharmacokinetics and Drug Therapy

Dept. of Pharmaceutical Biosciences

Uppsala University

Box 591

751 24 Uppsala

Sweden

Tel: +46 18 471 4355

www.farmbio.uu.se/research.php?avd=5

 

  _____

From: owner-nmusers_at_globomaxnm.com [mailto:owner-nmusers_at_globomaxnm.com] On
Behalf Of Frederik Gorter de Vries
Sent: Monday, October 15, 2007 1:40 PM
To: nmusers_at_globomaxnm.com
Subject: [NMusers] different COMRES statements

 

Dear NMuser,

 

I have a problem with implementing the CWRES in a model that's using $PRIOR.
Both need a $ABB COMRES= x statement, but with different numbers of x. Which
number should I take here? NONMEM seems to have no problem when using the
one from CWRES, but the CWRES cannot be calculated anymore in R using Xpose.
The following message appears:

 

> cwres <- compute.cwres(1, printToOutFile = FALSE,onlyNonZero = FALSE)

Error in `[.data.frame`(x, r, vars, drop = drop) :

        undefined columns selected

 

When I use the COMRES number from the $PRIOR NMTRAN gives an error regarding
the table with the wrong number of COM ARRAY .

 

I hope someone have seen this problem before and can tell me more about the
COMRES statement.

 

Thank you in advance and best regards,

 

Frederik Gorter de Vries

 

 

_______________________________________

Frederik Gorter de Vries

LACDR / Pharmacology

Gorleaus Laboratories

Einsteinweg 55 2333 CC Leiden

The Netherlands

Tel : ++31 71 527 6142

Email: frederik_at_gorterdevries.nl

 

 



Received on Mon Oct 15 2007 - 09:50:17 EDT