RE: output SE of population PK parameters in Nonmem
Hi Bo,
One simple thing you could do would be to insert the following code into
your control stream right before the $ERROR block. This will output the
point estimates and the SE into a ASCII dataset to read into other
applications.
$INFN
IF (ICALL.EQ.3) THEN
OPEN(50,FILE=outparms.txt')
WRITE (50,*) THETA
WRITE (50,*) SETHET
WRITE (50,*) OMEGA(DIAG)
WRITE (50,*) SEOMEG(DIAG)
WRITE (50,*) SIGMA(DIAG)
WRITE (50,*) SESIGM(DIAG)
ENDIF
Best wishes,
Scott Van Wart
Quoted reply history
-----Original Message-----
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Jin, Bo
Sent: Friday, February 01, 2008 11:48 AM
To: [email protected]
Subject: [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.
------------------------------------------------------------------------
------