NONMEM AR(1) residual bug?
From: drfreedman@drfreedmaninc.com
Subject: [NMusers] NONMEM AR(1) residual bug?
Date: 12/23/2003 8:59 PM
Colleagues
As a follow-up to my previous message on this topic, the
residual serial correlation control stream below yields the
following FSUBS.FOR fragment which does not compile with the
Compaq Fortran 6.6a compiler.
Is this a bug in NONMEM?
Thank-you in advance
Immanuel Freedman, Ph.D., MIEEE
(619) 884-1347
Control Stream:
$ERROR
" FIRST
" COMMON /NMPRD5/ C(250,29)
" DIMENSION T(50)
" IF (NEWIND.NE.2) I=0
" IF (MDV.EQ.0) THEN
" I=I+1
" T(I)=TIME
" DO 5 J=1,I
" 5 C(J,1)=EXP(-THETA(11)*(TIME-T(J)))
" ENDIF
FSUBS.FOR fragment:
COMMON /NMPRD5/ C(250,29)
DIMENSION T(50)
IF (NEWIND.NE.2) I=0
IF (EVTREC(N,14).EQ.0) THEN
I=I+1
T(I)=EVTREC(N,03)
DO 5 J=1,I
5 EVTREC(N,01)(J,1)=EXP(-THETA(11)*(EVTREC(N,03)-T
(J)))
ENDIF
_______________________________________________________