From: "michel.tod" <michel.tod@avc.ap-hop-paris.fr>
Subject: SIGMA estimates
Date: Mon, 26 Feb 2001 10:40:28 +0100
Dear NONMEM users,
In a $PRED block of statements, I need to read the current SIGMA estimates. In the users guide, I found that THETA and OMEGA estimates can be read in some read-only COMMONs. However, I did not found mention of SIGMA estimates.
Any suggestion or comment would be welcome...
Michel Tod
Dpt of Pharmacotoxicology
Avicenne Hospital
Bobigny, France
SIGMA estimates
From: "Pascal Girard" <pgirard@pharsight.com>
Subject: Re: SIGMA estimates
Date: Mon, 26 Feb 2001 11:13:01 +0100
Hi Michel,
You'l find final THETA, OMEGA and SIGMA estimates in ROCM6 and their S.E. in ROCM7.
Best,
Pascal Girard
From: Mats Karlsson <Mats.Karlsson@biof.uu.se>
Subject: Re: SIGMA estimates
Date: Mon, 26 Feb 2001 12:29:05 +0100
Dear Michel,
It is possible you could get it somehow, but one easy way of solving your problem is to FIX sigma to 1 and multiply it by a THETA:
Y = F + EPS(1)*THETA(.)
...
$THETA (0,.1)
$SIGMA 1 FIX
This may be slightly slower (Sturat has told me), but makes the estimate available (among other advantages).
Best regards,
Mats
Mats Karlsson, PhD
Professor of Pharmacometrics
Div. of Pharmacokinetics and Drug Therapy
Dept. of Pharmaceutical Biosciences
Faculty of Pharmacy
Uppsala University
Box 591
SE-751 24 Uppsala
Sweden
phone +46 18 471 4105
fax +46 18 471 4003
mats.karlsson@farmbio.uu.se
From: R.Port@DKFZ-Heidelberg.de
Subject: SIGMA estimates
Date: Mon, 26 Feb 2001 12:48:25 +0100
Hi, Dr. Tod:
Here is a control file where THETA, OMEGA, and SIGMA estimates are pulled from ROCM6 and are later used for the output table:
$PROBLEM urinary excretion
$INPUT ID TRTM AEXC=DV ; AEXC := amount excreted (mg)
$DATA test.data1 IGNORE=#
$SUBROUTINES
$PRED
"FIRST
" COMMON /ROCM6/ THETAF(40),OMEGAF(30,30),SIGMAF(30,30)
" DOUBLE PRECISION THETAF,OMEGAF,SIGMAF
" TTH1 = THETAF(1)
" OOM11 = OMEGAF(1,1)
" SSG11 = SIGMAF(1,1)
Y = THETA(1) + ETA(1) + EPS(1)
TH1 = TTH1
OM11 = OOM11
SG11 = SSG11
$THETA 30
$OMEGA 25
$SIGMA 25
;$MSFI test.msf1
$ESTIM PRINT=1 MAXEVALS=9999 ;POSTHOC ;MSF=test.msf1
$TABLE ID TRTM TH1 OM11 SG11 ETA1 NOPRINT ONEHEADER FILE=anmtable
With best regards, Ruedi Port
-----------------------------------------------------------------------------
Dr. R.E. Port, German Cancer Research Center, D-0200
P.O. Box 10 19 49, D-69009 Heidelberg, Germany
phone: +49-6221 42-3385
-3347
fax: -3382
e-mail: r.port@dkfz.de