RE: output SE of population PK parameters in Nonmem
Dear all,
It's a neat trick and works very well, but can anyone tell me the keyword
(like for instance SETHET below) for exporting the correlation matrix of the
estimates and for the non-parametric estimates: 'expected value of ETA' and
'Covariance matrix of ETA'....
Thanks in advance,
Rik
Quoted reply history
-----Original Message-----
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of GIRARD PASCAL
Sent: 01 February 2008 18:08
To: Jin, Bo; [email protected]
Subject: RE : [NMusers] output SE of population PK parameters in Nonmem
Hi Jin,
See in NONMEM help:
___________________________________________________________________
| |
| FINALIZATION EXAMPLE ($PRED AND $INFN) |
|_________________________________________________________________|
This example contains abbreviated code which can be inserted in a
$PRED or $INFN block. This code outputs final parameter estimates,
standard errors, minimum value of the objective function, and condi-
tional estimates of etas to various user files. The return codes from
Estimation and Covariance steps (zero for normal termination) are also
output.
IF (ICALL.EQ.3) THEN
DO WHILE(DATA)
IF (NEWIND.LE.1) WRITE (50,*) ETA
ENDDO
WRITE (51,*) OBJECT
WRITE (52,*) THETA
WRITE (53,*) SETHET
WRITE (54,*) OMEGA(BLOCK)
WRITE (55,*) SEOMEG(BLOCK)
WRITE (56,*) SIGMA(BLOCK)
WRITE (57,*) SESIGM(BLOCK)
WRITE (58,*) IERE,IERC
ENDIF
Alternatively, you take advantage of the fact that all those estimmates are
stored in Vectors and matrices of various COMMON:
COMMON /ROCM6/ THETAF(40),OMEGAF(30,30),SIGMAF(30,30)
COMMON /ROCM7/ SETH(40),SEOM(30,30),SESIG(30,30)
COMMON /ROCM8/ OBJECT
COMMON /ROCM9/ IERE,IERC
COMMON /CM12/ COV(2850)
and write your own INFN subroutine (see INFN help).
Best regards,
Pascal Girard, PhD
EA 3738, CTO
Fac Medecine Lyon-Sud, BP12
69921 OULLINS Cedex, France
[EMAIL PROTECTED]
Tel +33 (0)4 26 23 59 54 / Fax +33 (0)4 26 23 59 76
Master Recherche Lyon 1 Santé et Populations, Spécialité PhIT
http://master-sante-pop.univ-lyon1.fr/
>-----Message d'origine-----
>De : [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] De
>la part de Jin, Bo
>Envoyé : vendredi 1 février 2008 17:48
>À : [email protected]
>Objet : [NMusers] output SE of population PK parameters in Nonmem
>
>Dear all: I have a question. Is there any way to output the SE or SD of
>a population PK parameter estimate in Nonmem (say typical value of Cl)
>into some files which can be read by SAS or S-Plus etc.? (e.g. some
>excel or TAB files?)
>
>thanks,
>
>- Bo
>
>_____________________________________
>
>BO JIN
>Clinical Pharmacology Statistics
>Merck Research Labs
>Phone: 267-305-7876
>
>
>
>---------------------------------------------------------------------------
>---
>Notice: This e-mail message, together with any attachments, contains
>information of Merck & Co., Inc. (One Merck Drive, Whitehouse Station,
>New Jersey, USA 08889), and/or its affiliates (which may be known
>outside the United States as Merck Frosst, Merck Sharp & Dohme or MSD
>and in Japan, as Banyu - direct contact information for affiliates is
>available at http://www.merck.com/contact/contacts.html) that may be
>confidential, proprietary copyrighted and/or legally privileged. It is
>intended solely for the use of the individual or entity named on this
>message. If you are not the intended recipient, and have received this
>message in error, please notify us immediately by reply e-mail and then
>delete it from your system.
>
>---------------------------------------------------------------------------
>---