RE: output SE of population PK parameters in Nonmem
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: <
>
> [email protected]
>
> > 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. > >--------------------------------------------------------------------------- >---