Re: how to implement AR residual error model in Nonmem with FOCE INT routine
The verbatim code provided below was written for NONMEM 6 or earlier. You'll find more detail here:
http://gaps.cpb.ouhsc.edu/nm/98jul1098.html
This code will not work with NONMEM 7, because COMMONS are no longer part of the NONMEM system. The NONMEM 7 help file ARnmprd5.for: " CORRELATION ACROSS L2 RECORDS " seems to address this type of auto- correlated residual model, but I have not yet tried this in version 7. Perhaps others have some experience.
Hope this helps.
Marc
Marc R. Gastonguay, Ph.D. < [email protected] >
President & CEO, Metrum Research Group LLC < metrumrg.com >
Scientific Director, Metrum Institute < metruminstitute.org >
2 Tunxis Rd, Suite 112, Tariffville, CT 06081 Direct: +1.860.670.0744 Main: +1.860.735.7043 Fax: +1.860.760.6014
Quoted reply history
On Jan 6, 2010, at 11:29 AM, [email protected] wrote:
> The following code causes errors. I use ADVAN13 to run a PD model.
>
> $ERROR
> " FIRST
> " COMMON /NMPRD5/ C(50,29)
> " DIMENSION T(300)
> " IF (NEWIND.NE.2) I=0
> " IF (CMT==4) THEN
> " I=I+1
> " T(I)=TIME
> " DO 5 J=1,I
> " 5 C(J,I) = EXP(-THETA(9)*(TIME-T(J)))
> " ENDIF
>
> Errors:
>
> FIRST
> 1
> Error: Unclassifiable statement at (1)
> fsubs.f90:365.25:
> COMMON /NMPRD5/ C(50,29)
> 1
> Error: Unexpected COMMON statement at (1)
> fsubs.f90:366.21:
> DIMENSION T(300)
> 1
> Error: Unexpected attribute declaration statement at (1)
> fsubs.f90:370.22:
> T(I)=EVTREC(NVNT,03)
> 1
> Error: Unexpected STATEMENT FUNCTION statement at (1)
> fsubs.f90:372.50:
> 5 C(J,I) = EXP(-THETA(9)*(EVTREC(NVNT,03)-T(J)))
> 1
> Error: Unexpected STATEMENT FUNCTION statement at (1)
> fsubs.f90:373.4:
> ENDIF
> 1
> Error: Expecting END DO statement at (1)
> fsubs.f90:492.6:
> END
> 1
> Error: END DO statement expected at (1)
> fsubs.f90:493.6:
> SUBROUTINE TOL(NRD)
> 1
> Error: Unclassifiable statement at (1)
> fsubs.f90:494.39:
> USE SIZES, ONLY: DPSIZE,ISIZE
> 1
> Error: Unexpected USE statement at (1)
> fsubs.f90:495.25:
> DIMENSION :: NRD(*)
> 1
> Error: Unexpected attribute declaration statement at (1)
> fsubs.f90:496.32:
> INTEGER(KIND=ISIZE) :: NRD
> 1
> Error: Unexpected data declaration statement at (1)
> fsubs.f90:497.6:
> NRD(1)=7
> 1
> Error: Unclassifiable statement at (1)
> fsubs.f90:500.6:
> END
> 1
> Error: END DO statement expected at (1)
> fsubs.f90:501.6:
> SUBROUTINE DES (A,P,T,DADT,IR,DA,DP,DT)
> 1
> Error: Unclassifiable statement at (1)
> fsubs.f90:502.17:
> USE NMPRD4P
> 1
> Error: Unexpected USE statement at (1)
> fsubs.f90:503.39:
> USE SIZES, ONLY: DPSIZE,ISIZE
> 1
> Error: Unexpected USE statement at (1)
> fsubs.f90:504.57:
> USE NMPRD_INT, ONLY: MSEC=>ISECDER,MFIRST=>IFRSTDER
> 1
> Error: Unexpected USE statement at (1)
> fsubs.f90:505.33:
> USE PROCM_REAL,ONLY: EVTREC
> 1
> Error: Unexpected USE statement at (1)
> fsubs.f90:506.39:
> USE PROCM_INT, ONLY: NVNT=>NEVENT
> 1
> Error: Unexpected USE statement at (1)
> fsubs.f90:507.52:
> USE PRMOD_INT, ONLY: ICALL=>ICALLD,IDEFD,IDEFA
> 1
> Error: Unexpected USE statement at (1)
> fsubs.f90:508.78:
> IMPLICIT REAL(KIND=DPSIZE) (A-Z)
> 1
> Error: Letter A already has an IMPLICIT type at (1)
> fsubs.f90:509.10:
> SAVE
> 1
> Error: Unexpected attribute declaration statement at (1)
> fsubs.f90:510.78:
> INTEGER(KIND=ISIZE) :: IR
> 1
> Error: Unexpected data declaration statement at (1)
> fsubs.f90:511.20:
> DIMENSION :: A(*),P(*),DADT(*),DA(IR,*),DP(IR,*),DT(*)
> 1
> Error: Cannot change attributes of USE-associated symbol amnt at (1)
> fsubs.f90:514.6:
> IDEFD(1)= 0
> 1
> Error: Unclassifiable statement at (1)
> Fatal Error: Error count reached limit of 25.
> gfortran: Internal error: Aborted (program f951)
> Please submit a full bug report.
> See http://gcc.gnu.org/bugs.html for instructions.
> gfortran: fsubs.o: No such file or directory
> No nonmem execution.
>
> ----- Original Message -----
> From: [email protected]
> Date: Wednesday, January 6, 2010 10:50 am
>
> Subject: Re: [NMusers] how to implement AR residual error model in Nonmem with FOCE INT routine
>
> To: [email protected]
>
> > Hello Nonmem Users,
> >
> > Is there a more detailed explanation of the code in
> >
> >
> >
> > http://cognigencorp.com/nonmem/nm/97dec172003.html somwhere? I
> > try to implement a second-order autoregressive process. If
> > C(J,1) is a correlation coefficient, can I use more complex
> > equations for C(j,1)? Will this code work if we use ADVAN13 (PD
> > model here)?
> >
> > Thanks!
> > Pavel
> >
> >
> >
> >
> > ----- Original Message -----
> > From: "Gastonguay, Marc"
> > Date: Tuesday, December 29, 2009 10:09 am
> > Subject: Re: [NMusers] how to implement AR residual error model
> > in Nonmem with FOCE INT routine
> > To: Ethan Wu
> > Cc: [email protected]
> >
> > > Hello Ethan,
> > > This archived NMusers discussion thread might be useful:
> > >
> > > http://cognigencorp.com/nonmem/nm/97dec172003.html
> > >
> > > Marc
> > >
> > > Marc R. Gastonguay, Ph.D.
> > > President & CEO, Metrum Research Group LLC < metrumrg.com >
> > > Scientific Director, Metrum Institute < metruminstitute.org >
> > > 2 Tunxis Rd, Suite 112, Tariffville, CT 06081 Direct:
> > > +1.860.670.0744 Main: +1.860.735.7043 Fax: +1.860.760.6014
> > >
> > >
> > >
> > > On Dec 29, 2009, at 9:14 AM, Ethan Wu wrote:
> > >
> > > > Dear all,
> > > > I am fitting a PK/PD model using FOCE INT. I do see a high
> > >
> > > > correlation between WRES
> > > > I wonder could someone help me with implementation of AR
> > > residual
> > > > error model in NONMEM.
> > > > Many thanks! Happy new year to all.
> > > >
> > > >
> > >
> > >
> >