________________________________
De: Vanesa Escudero Ortiz
Enviado el: 14 June 2011 10:4 AM
Para: '[email protected]'
Asunto: [NMusers] Question about correlations with NWPRI
Hi all,
I have problems running the following model trying to implement correlations in
NWPRI subroutine. The model runs fine (minimization successful) when
correlations are not implement. However, when correlations are implemented an
error message appeared. Model a priori does not consider the correlation
between CL and V2.
The error message is:
INPUTS SPECIFIED TO ROUTINE NWPRI ARE INAPPROPRIATE
Model file
$PROB MODELO
$INPUT ID NHC CYCL TIME AMT DOSE DUR CONC DV CNR RATE EVID BSA AGE SEX
$DATA data.csv IGNORE=#
$SUBROUTINE ADVAN5
$PRIOR NWPRI NTHETA=6 NETA=9 NEPS=1 NTHP=6 NETP=9 NPEXP=1
$MODEL
NCOMPARTMENTS=3
COMP(CENTRAL)
COMP(PERIPH1)
COMP(PERIPH2)
$PK
OCC3=0
OCC4=0
IF(CYCL.EQ.1) OCC3=1
IF(CYCL.EQ.2) OCC4=1
IOV1 = OCC3*ETA(4)+OCC4*ETA(5)
IOV2 = OCC3*ETA(6)+OCC4*ETA(7)
IOV3 = OCC3*ETA(8)+OCC4*ETA(9)
CL = THETA(1)*EXP(ETA(1)+IOV1)
V1 = THETA(2)
Q2 = THETA(3)*EXP(ETA(3)+IOV2)
V2 = THETA(4)*EXP(ETA(2)+IOV3)
Q3 = THETA(5)
V3 = THETA(6)
S1 = V1
K10 = CL/V1
K12 = Q2/V1
K21 = Q2/V2
K13 = Q3/V1
K31 = Q3/V3
$ERROR
IPRED = LOG(0.000001)
IF (F.GT.0) IPRED = LOG(F)
IRES = DV - IPRED
Y= IPRED + EPS(1)
$THETA (0 57.6)
$THETA (0 21.2)
$THETA (0 84.7)
$THETA (0 2170)
$THETA (0 34.7)
$THETA (0 101)
$OMEGA BLOCK(2)
0.045
0.0067 0.001 ;CL_V2
$OMEGA 0.011 ;IIV V2
$OMEGA BLOCK(1) 0.0292 ;IOV CL
$OMEGA BLOCK(1) SAME
$OMEGA BLOCK(1) 0.001 ;IOV Q2
$OMEGA BLOCK(1) SAME
$OMEGA BLOCK(1) 0.001 ;IOV V2
$OMEGA BLOCK(1) SAME
$THETA (57.6) FIX
$THETA (21.2) FIX
$THETA (84.7) FIX
$THETA (2170) FIX
$THETA (34.7) FIX
$THETA (101) FIX
$OMEGA BLOCK(6) 2.828
0 1.719
0 0 8.216
0 0 0 126.077
0 0 0 0 3.78
0 0 0 0 0 12.12 FIX
$OMEGA BLOCK(2)
0.0402
0 0.0053 FIX ;CL_V2
$OMEGA 0.0112 FIX ;IIV V2
$OMEGA BLOCK(1) 0.032 FIX ;IOV CL
$OMEGA BLOCK(1) SAME
$OMEGA BLOCK(1) 0.258 FIX ;IOV Q2
$OMEGA BLOCK(1) SAME
$OMEGA BLOCK(1) 0.068 FIX ;IOV V2
$OMEGA BLOCK(1) SAME
$THETA 10 FIX
$THETA 10 FIX
$THETA 10 FIX
$THETA 10 FIX
$THETA 10 FIX
$THETA 10 FIX
$SIGMA 0.05
$ESTIMATION METHOD=1 INTERACTION MAXEVALS=9990 PRINT=10 POSTHOC
$COVARIANCE PRINT=E MATRIX=R
$TABLE ID TIME AMT CYCL DOSE IPRED IRES NOPRINT ONEHEADER
FILE=TAB_MOD_C_CORRE.TXT
$TABLE ID CL V1 V2 Q2 Q3 V3 ETA1 ETA2 ETA3
ETA4 ETA5 ETA6 ETA7 ETA8 ETA9
NOPRINT ONEHEADER FILE=PAR_MOD_C_CORRE.TXT
Any help would be much appreciated
Best regards
_________________________________________
Vanesa Escudero Ortíz
Personalized Pharmacotherapy Unit
Platform of Oncology
USP Hospital San Jaime
Ptda. de la Loma, s/n
03180 - Torrevieja (Alicante), Spain
Tel.: +34 96 692 13 13 ext. 3197
Fax.: +34 96 692 27 06
Email: [email protected]
Web: www.uspsanjaime.com
_________________________________________