$THETA
Gives initial estimates and bounds for THETAs.
|
|
Discussion
This record gives initial estimates and bounds for elements of the THETA matrix.
THETAs are numbered in the order in which they are defined. When a $MSFI record appears in the problem
specification, the $THETA record should not used.
Options
Each value defines a theta and gives its initial estimate and bounds. A value has one of four forms:
init [FIXED]
Init is the initial estimate. If FIXED is used, the final parameter estimate is to be constrained to equal the initial parameter estimate.
([low,] init [,up] [FIXED])
Low and up are lower and upper bounds respectively. They are the boundaries for the minimization search. Commas are optional. If an upper bound is used, a lower bound must also be used (e.g., -INF; see below). The lower and upper bounds (or if just the lower bound is used, then just this bound) may be omitted, in which case this form differs from the one described above only in so far as with this form, the initial estimate and the FIXED attribute are enclosed in parentheses. When FIXED is used, and a bound appears the bound must equal the initial estimate. FIXED is implied when all three values are equal. The lower bound can be -INF (negative infinity), and the upper bound can be INF (infinity). They are the defaults for lower and upper bounds. 1
([low,] init [,up]) [FIXED]
This is just like the form described above except that if the FIXED attribute is used, the attribute occurs outside the parentheses, and then if moreover, a bound appears, the bound need not equal the initial estimate.
(low,,up)
The commas are required. Because no initial estimate is given, a search for an initial estimate (the Initial Estimates Step) is undertaken by NONMEM. With NONMEM 7.4, when initial THETAs are to be estimated, evaluations can now be done for FOCE and Laplace, not just for FO.
(value)xn
Any initial value or group of initial values may be enclosed in parentheses and followed by "xn", which means to replicate the values within parentheses n times ("repeated value"). The values within the parenthesis may have any of the above forms. For example, the following two are equivalent:
|
|
UNINT (NM75)
UNINT is used during the Optimal Design Step to identify a theta as uninteresting. UNINT may be used anywhere that FIXED may be used.
label=value [FIXED] (NM75)
This is a compact method of defining an element of theta, specifying
its initial estimate, and specifying a symbolic index for this
element of THETA. The symbolic index may be used for THETA in
abbreviated code, and will also identify this element of THETA in
the NONMEM output. (Only the first 9 characters of the label will
appear). If new $THETA records change the ordering, the
abbreviated code does not have to be changed. For example, suppose
the third element of THETA that is defined happens to be
|
|
Once defined, the $PK or $PRED code can use these definitions:
|
|
The NONMEM report will describe the relationship,
|
|
Alternatively one can specify the label using
|
|
control record but the NONMEM report will not identify the
relationship and will not generate the symbolic label THETA(CL) in
the report. For dynamic, or implicit mapping of labels, such as for
various occasions, these still need to be done via the $ABBR REPLACE
record.
$THETA records must be placed ahead of any records that use the
symbolic label. Therefore these records must be placed before
$THETAI, $THETAR, and abbreviated code.
NAMES (label …)value … (NM75)
A compact way of defining one or more THETAs with labels and initial values. For example
|
|
NUMBERPOINTS=n
During NONMEM's search for an initial estimate, a number of points will be examined. This number will be automatically determined by NONMEM, or it can be specified with this option. May also be coded NUMBERPTS, NUMPOINTS, NUMPTS.
ABORT | NOABORT | NOABORTFIRST
ABORT: during the Initial Estimates step, NONMEM is to abort when PRED sets the error return code to 1. (The PRED error return code n is set by the statement "EXIT n [k]" in abbreviated code, or by the statement IERPRD=n in user-supplied code, or by PREDPP when it detects an error.) This is the default.
NOABORT: during the Initial Estimates step, NONMEM is simply to ignore values of the theta vector that result in PRED error return code 1. (Ordinarily the first value of the theta vector is never ignored.) These will not be feasible values for an initial estimate.
NOABORTFIRST: Same as NOABORT option, but also applies to the first value of the theta vector that is tried. It cannot be shortened; all characters must be coded. May be used with the NOABORT option, in which case the stronger condition (NOABORTFIRST) takes precedence.
They are communicated to NONMEM as numeric values -1000000 and 1000000. Therefore one can also use these two numbers to indicate negative and positive infinity here.