Re: [NMusers] 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
>
> -----Original Message-----
> From: owner-nmusers_at_globomaxnm.com [mailto:owner-nmusers_at_globomaxnm.com] On
> Behalf Of GIRARD PASCAL
> Sent: 01 February 2008 18:08
> To: Jin, Bo; nmusers_at_globomaxnm.com
> 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
> pascal.girard_at_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_at_globomaxnm.com [mailto:owner-nmusers_at_globomaxnm.com] De
>> la part de Jin, Bo
>> Envoyé : vendredi 1 février 2008 17:48
>> À : nmusers_at_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.
>>
>> ---------------------------------------------------------------------------
>> ---
>
>
>
Received on Wed Feb 06 2008 - 09:53:24 EST