AR(1) residuals
From: drfreedman@drfreedmaninc.com
Subject: [NMusers] AR(1) residuals
Date: 12/17/2003 6:42 PM
Colleagues
Viewing the presentation
NONMEM Alphabet Soup: L2, AR(1), and INFN
by Marc Gastonguay posted at
; http://www.globomaxservice.com/presentations/presentations.h
tml , I tried the following code (all drug-specific)
information removed. The code including the "FIRST statement
appears to be inserted into FSUBS.FOR in an executable
section not permitting specification statements.
Can you suggest any changes? Can this approach to serial
correlation work with the ADBVAN subroutines?
Thank-you in advance
Dr. Immanuel Freedman
; run001.mod
$DATA input.csv IGNORE=C
;
$SUBROUTINE ADVAN2 TRANS2
;
$PK
;PK MODEL
F1=value; Bioavailability
CL=THETA(1)*EXP(ETA(1)); Clearance (mL/kg/d)
V=THETA(2)*EXP(ETA(2)); Central volume (mL/kg)
KA=THETA(3)*EXP(ETA(3)); Absorption rate (/d)
S2=V; Scale factor
;
$ERROR
" FIRST
" COMMON /NMPRD5/ C(50,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,I) = EXP(-THETA(4)*(TIME-T(J)))
" ENDIF
CP=A(2)/V;
Y=CP*EXP(ERR(1))+ERR(2);
;
$THETA
;values
$OMEGA
;values
$SIGMA
;values
$ESTIMATION NOABORT MAXEVAL=9999 PRINT=10 METHOD=0 POSTHOC
MSFO=run001.msf
$COVARIANCE
;
$TABLE NOPRINT ONEHEADER FILE=sdtab001
C ID TIME AMT DV CMT TYPE ROUT
F1 CL V KA;
;