Re: output SE of population PK parameters in Nonmem
Hi Rik,
$COV PRINT=E
should give you correlation matrix of the
population parameter estimates
As far as I know, there is no way to get SE (and hence covariance matrix) for individual parameter estimates (ETA values).
Leonid
--------------------------------------
Leonid Gibiansky, Ph.D.
President, QuantPharm LLC
web: www.quantpharm.com
e-mail: LGibiansky at quantpharm.com
tel: (301) 767 5566
Rik Schoemaker wrote:
> 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.
> >
> > ---------------------------------------------------------------------------
> > ---