Confidence intervals
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