 | Jin, Another options it the Fortran app nmsee ( ftp://ftp.globomaxnm.com/Public/nonmem/nmsee/) or a version of nmsee I recently put together (cleverly called nmsee2) at http://www.nextlevelsolns.com/downloads.html. These are both post-processors that read the NONMEM output file (and so you only get 3 digits for the parameters). Either of these can be put near the end of the NMFE?.bat file to run automatically after NONMEM. Mark Mark Sale MD
Next Level Solutions, LLC
www.NextLevelSolns.com
919-846-9185
-------- Original Message --------
Subject: RE : [NMusers] output SE of population PK parameters in Nonmem
From: "GIRARD PASCAL" <[EMAIL PROTECTED]>
Date: Fri, February 01, 2008 12:07 pm
To: "Jin, Bo" <[EMAIL PROTECTED]>, <[email protected]>
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
pascal.girard@adm.univ-lyon1.fr
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 : owner-nmusers@globomaxnm.com [mailto:owner-nmusers@globomaxnm.com] De
>la part de Jin, Bo
>Envoyé : vendredi 1 février 2008 17:48
>À : nmusers@globomaxnm.com
>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.
>
>---------------------------------------------------------------------------
>---
|
