[NMusers] Confidence intervals

From: Venkat Pokiri <pokiri75_at_gmail.com>
Date: Thu, 1 May 2008 12:54:51 -0600

hi all:
I am relatively new to Nonmem and trying to incorporate the confidence
interevals in the control stream. I amgetting the following errror;
AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.

 195 $TABLE OR $SCAT REFERS TO NON-RESERVED ELEMENT OF COM
ARRAY.
fsubs did not get created by NM-TRAN.
Could some one advice me in this regard. I have alos attached the control
stream below.

Thanks in advance,
 Vpokiri

$PROB RUN# 108
$INPUT C ID TIME DV AMT ADDL II MDV
$DATA 108.csv IGNORE=C
$SUBROUTINE ADVAN2 TRANS2
$PK
  TVCL=THETA(1)
  CL=TVCL*EXP(ETA(1))
  TVV=THETA(2)
  V=TVV*EXP(ETA(2))
  TVKA=THETA(3)
  KA=TVKA*EXP(ETA(3))
  S2=V

$ERROR
"FIRST
" COMMON /rocm6/THETAF(40),OMEGAF(30,30),SIGMAF(30,30)

   IPRED = F
   IRES = DV-IPRED
   Y = IPRED+IPRED*EPS(1)+EPS(2)

" LAST
" IF(COMACT.EQ.1) THEN
" VARCP = G(1,1)**2*OMEGAF(1,1)
" VARCP = G(2,1)**2*OMEGAF(2,2)+VARCP
" VARCP = G(3,1)**2*OMEGAF(3,3)+VARCP
" VARCP = G(4,1)**2*OMEGAF(4,4)+VARCP
" VARCP = G(1,1)*G(2,1)*OMEGAF(1,2)+VARCP
" VARCP = G(1,1)*G(3,1)*OMEGAF(1,3)+VARCP
" VARCP = G(1,1)*G(4,1)*OMEGAF(1,4)+VARCP
" VARCP = G(2,1)*G(3,1)*OMEGAF(2,3)+VARCP
" VARCP = G(2,1)*G(4,1)*OMEGAF(2,4)+VARCP
" VARCP = G(3,1)*G(4,1)*OMEGAF(3,4)+VARCP
" SDCP = SQRT(VARCP)
" COM(1)= IPRED+1.96*SDCP
" COM(2)= IPRED-1.96*SDCP
" ENDIF

$THETA
  (0,24.6) ;[CL]
  (0,450) ;[V]
  (0,1.32) ;[KA]

$OMEGA
  0.141 ;[P] INTERINDIVIDUAL VARIABILITY IN CL
  0.233 ;[P] INTERINDIVIDUAL VARIABILITY IN V
  0.547 ;[P] INTERINDIVIDUAL VARIABILITY IN KA

$SIGMA
  1.20 ;[P] PROPORTIONAL COMPONENT OF RESIDUAL VARIABILITY
  0.243 ;[A] ADDITIVE COMPONENT OF RESIDUAL VARIABILITY

$ESTIMATION METHOD=1 PRINT=1 MAXEVAL=9999 NOABORT SIGDIGITS=3 POSTHOC
MSF=108.MSF
$SIMULATION (9283745)
$COV
$TABLE ID TIME AMT FILE=SIMDATA8 NOPRINT
;$TABLE ID NOPRINT CL V KA ONEHEADER FILE=patab108
$TABLE ID NOPRINT ONEHEADER FILE=cotab108
;$TABLE ID NOPRINT ONEHEADER FILE=catab108
$TABLE ID TIME IPRED NOPRINT ONEHEADER FILE=sdtab108
;$TABLE ID TIME DV IPRED IRES IWRES NOPRINT ONEHEADER FILE=108.txt
$TABLE ID AMT CL V KA NOPRINT ONEHEADER FILE=sdtab108.txt
$TABLE ID TIME IPRED COM(1)=UPP COM(2)=DOWN NOPRINT ONEHEADER FILE=citab108

Received on Thu May 01 2008 - 14:54:51 EDT

This archive was generated by hypermail 2.2.0 : Fri May 02 2008 - 11:22:07 EDT