SIGMA estimates
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