 | Rik, I don't know anything about the non-parametric estimates: 'expected value of ETA' and
'Covariance matrix of ETAI. But, think the correlation matrix is CORRM in COMMON /CM12/ COVM(LPAR3),COVINM(LPAR3),STHTA(LTH),SEN(LVR,LVR), 1 CORRM(LPAR3),NVLS,VLS(LPAR) and you can get at it something like this: C DIMENSTION OF CORRLEATION MATRIX, ONE FOR EACH THETA AND FULL MATRIX FOR EACH OMEGA AND SIGMA C EVEN IF NOT USED WRITE(44,'(A10)') 'R^2 MATRIX' C NUMPAR IS TOTAL SIZE OF THE RANDOM EFFECTS MATRIX, DEFAULT IS 70X70 NUMPAR = LTH + LVR CORR = 0 C P IS THE POSITION IN CORRM, IT IS A 1XNUMPAR ARRAY P = 0 &nbs!
p; DO 66 I = 1,NUMPAR DO 65 N = 1,I P = P + 1 C CORRELATION MATRIX IS IN CORRM, LOWER TRIANGULAR ONLY WRITE(44,70) CORRM(P) ABSVAL = ABS(CORRM(P)) C IF (OFF DIAGONAL, I.NE.N) IS > 0.95, FAIL CORRELATION IF((ABSVAL.GT.0.95).AND.(I.NE.N).AND.(ABSVAL.LE.1.00)) THEN CORR = 1 END IF 65 CONTINUE C new line WRITE(44,'(A2)') ' ' 66 CONTINUE 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: "Rik Schoemaker" <[EMAIL PROTECTED]>
Date: Wed, February 06, 2008 8:30 am
To: <nmusers@globomaxnm.com>
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@globomaxnm.com [mailto:owner-nmusers@globomaxnm.com] On
Behalf Of GIRARD PASCAL
Sent: 01 February 2008 18:08
To: Jin, Bo; nmusers@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@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.
>
>---------------------------------------------------------------------------
>---
|
