RE: 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
_____
Quoted reply history
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Frederik Gorter de Vries
Sent: Monday, October 15, 2007 1:40 PM
To: [email protected]
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: [EMAIL PROTECTED]