NONMEM AR(1) residual bug?

From: Immanuel Freedman Date: December 23, 2003 technical Source: cognigencorp.com
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 _______________________________________________________
Dec 17, 2003 Immanuel Freedman AR(1) residuals
Dec 17, 2003 Marc Gastonguay RE: AR(1) residuals
Dec 19, 2003 Marc Gastonguay RE: AR(1) residuals
Dec 23, 2003 Immanuel Freedman NONMEM AR(1) residual bug?