NONMEM VI 1.0 Bug Alert #3
NONMEM Code.
3. Background. When the NWPRI utility is used, the values output in
the variance-covariance,
correlation and inverse covariance matrices of the estimates, other than
those for THETA's alone,
may not correspond to the printed labels and the output may be
incomplete.
Workaround: Ignore the values other than those for THETA's alone.
When NWPRI is used, covariances of the estimates that involve OMEGA's
are not meaningful.
Fix: To limit output to those values for THETA's alone make the
following two changes, one for
the default output and one for compressed output.
1) In PRNT1 locate the statement
430 CONTINUE
and immediately after it add
C ADDED 2/2007
KK=NV
2) In PRNT5 locate the statement
20 CONTINUE
and immediately after it add
C ADDED 2/2007
IF (K1.NE.M1.OR.K2.NE.M2.OR.K3.NE.M3) GO TO 100
Be sure to recompile NONMEM after making any code changes.
___________________________________________________
NMTRAN Code
3. Background: When the control stream includes $NONPAR OMIT
then the $COVARIANCE record is, in effect, ignored.
Workaround: Remove "$NONPAR OMIT" from any control stream that
also includes $COVAR.
Fix:
In GENERA
Locate the first instance of "GO TO 600".
Replace with "GO TO 550".
Original code:
WRITE (LINE,FORMS(7)) 'NONP',2
CALL PUT(LINE)
GO TO 600
Changed code:
WRITE (LINE,FORMS(7)) 'NONP',2
CALL PUT(LINE)
C CHANGED 1/2007
GO TO 550
4. Background: The following abbreviated code does not work correctly:
IF (ICALL.EQ.0) THEN
DO WHILE (DATA)
... transgeneration statements ...
ENDDO
ENDIF
NM-TRAN produces a spurious error message:
370 "DO WHILE" MAY NOT BE USED WHEN ICALL IS 2.
Workaround:
Change the abbreviated code from
IF (ICALL.EQ.0) THEN
to
IF (ICALL.EQ.0.AND.ICALL.NE.2) THEN
Fix: Two changes are needed to NM-TRAN's subroutine GETITM.
Locate:
IF (NDODAT.NE.0.AND..NOT.T(3).AND..NOT.T(1)) THEN
Change to:
C Modified 2/2007
IF (NDODAT.NE.0.AND..NOT.T(3).AND..NOT.T(1).AND..NOT.T(0)) THEN
Locate:
IF (.NOT.T(3).AND..NOT.T(1)) THEN
Change to:
C Modified 2/2007
IF (.NOT.T(3).AND..NOT.T(1).AND..NOT.T(0)) THEN
Be sure to recompile NM-TRAN after making any code changes.
_______________________________________________________________________________________________________________________________________