From: "Olofsen, E. (ANST)"
Subject: [NMusers] correlation function for the epsilons
Date:Wed, 14 May 2003 22:35:57 +0200
Dear group,
I wanted to specify a correlation function for the epsilons as discussed in:
http://www.cognigencorp.com/nonmem/nm/98jul101998.html
Furthermore, I wanted to test this with single-subject data. At the end of
the above link it is noted that the single-subject data should be entered
as population data, but with OMEGA constrained to zero. However, when the
data file and $INPUT contain the ID item, the covariance step does not give
correct results (the estimation step does). If I remove `ID', NONMEM gives
a warning about the SIGMA matrix, it treats every observation as coming from
a different individual, and I wonder if the above technique will work.
Does anyone know how I should proceed?
Thanks in advance,
Erik Olofsen
Leiden University Medical Center
PS: An example (t.d contains e.g. [1 0; 1 1; 1 2; 1 3]):
$PROBLEM indeed
$DATA t.d
$INPUT ID DV
$PRED
A = THETA(1)+ETA(1)
Y = A+ERR(1)
$THETA
1
$OMEGA
0 FIX
$SIGMA
0.1
$ESTIMATION MAX=9999 SIG=4 PRINT=1
$COVR
_______________________________________________________