RE: Problem running NWPRI

From: [email protected] [mailto:[email protected]] On Behalf Of Brogren, Jacob Sent: den 19 augusti 2010 15:15 To: [email protected] Subject: [NMusers] Problem running NWPRI Hi all, I have problems running the following model trying to implement NWPRI in NONMEM. The model runs fine (minimization successful and $COV run) when not using PRIOR. Any help is much appreciated. (NB. The prior modes are perturbed.) The error message is: ======================================== “Bareword found where operator expected at (eval 119) line 1, near "0INITIAL" (Missing operator before INITIAL?) Bareword found where operator expected at (eval 124) line 1, near "0SIGMA" (Missing operator before SIGMA?) Bareword found where operator expected at (eval 131) line 1, near "0ESTIMATION" (Missing operator before ESTIMATION?)” ========================================== Model file ================================== $PROB Test NWPRI $INPUT STUD ID PID=DROP EID=DROP DAT2=DROP TIME TAD=DROP DAY FLAG AMT CP=DROP DV MDV UVOL CMT EVID DOSE RACE=DROP SEX AGE WT HT ALAT=DROP ASAT=DROP ALB CLCR $DATA ../data/data.csv IGNORE=@ $SUBS ADVAN4 TRANS4 PRIOR=nwprior.f $PK IF(SEX.EQ.1) LBW = (1.1*WT) - 128*(WT**2/HT**2) IF(SEX.EQ.2) LBW = (1.07*WT)- 148*(WT**2/HT**2) TVLG = THETA(1) TVKA = THETA(2)*EXP(-THETA(3)*(DOSE/100)) TVV2 = THETA(4)*(LBW/58)**1 TVQ = THETA(5)*(LBW/58)**0.75 TVV3 = THETA(6)*(LBW/58)**1 TVAD = THETA(9) TVKCL = THETA(10) TVCL = (THETA(7)*(CLCR/90)**THETA(8) + TVAD*EXP(-TVKCL*TIME))*(LBW/58)**0.75 ALAG1= TVLG KA = TVKA*EXP(ETA(1)) CL = TVCL*EXP(ETA(2)) V2 = TVV2*EXP(ETA(3)) Q = TVQ V3 = TVV3 S2 = V2 S3 = V3 $ERROR DEL = 0 IF (F.EQ.0) DEL=1 IPRED=LOG(F +DEL) W = SQRT(THETA(11)**2) + DEL IRES = DV-IPRED IWRES= IRES/W Y = IPRED + ERR(1)*W ;--------- Initial estimates -------------- $THETA (0.01,0.35,3) ; 1. ALAG1 (0.01,3.5,15) ; 2. KA (0.01,0.55,3) ; 3. DOSE on KA (1,150,300) ; 4. V2 (0.1,15,50) ; 5. Q (0.1,40,200) ; 6. V3 (0.01,30,80) ; 7. CL at CLCR=90 mL/min (0.01,0.6,3) ; 8. Power on CLCR (0.01,6,20) ; 9. Addition to CL at first dose (0.01,0.05,3) ; 10. Rate constant of additional CL decay (0.01,0.2,2) ; 11. prop error $OMEGA 0.09 ; 1. KA $OMEGA BLOCK(2) 0.09 ; 2. CL 0.003 0.09 ; 3. V2 ;--------- Prior modes for THETA -------------- $THETA (0.3 FIX) ; 1. ALAG1 (lag time) (3 FIX) ; 2. KA (0.5 FIX) ; 3. DOSE on KA (175 FIX) ; 4. V2 (10 FIX) ; 5. Q (50 FIX) ; 6. V3 (35 FIX) ; 7. CL at CLCR = 90 mL/min (0.5 FIX) ; 8. Power on CLCR (5 FIX) ; 9. Addition to CL at first dose (0.1 FIX) ; 10. Rate constant of additional CL decay (0.2 FIX) ; 11. Proportional error ;--------- COV matrix for THETA --------------- $OMEGA BLOCK(11) 0.0001 0.0021 0.2628 0.0001 0.0233 0.0047 0.007 0.5331 0.0058 13.9852 -0.0009 -0.0707 0.0058 -2.061 0.7897 -0.0015 -0.126 0.0158 -4.7538 1.7982 4.3763 0.0007 0.0669 0.0088 1.7858 0.0403 0.1896 0.7052 0.0001 0.0009 0.0003 -0.0095 -0.0015 -0.0017 0.0015 0.0009 -0.0014 -0.1102 -0.0057 -0.1326 -0.0891 -0.4479 -0.37 -0.005 1.5127 0.0001 0.0005 0.0001 -0.0049 -0.0018 -0.0047 -0.001 0.0001 0.0064 0.0002 0.0001 0.0012 0.0003 0.0262 -0.0033 -0.0055 0.0047 0.0001 -0.005 0.0001 0.0004 FIX ;--------- Prior mode for OMEGA --------------- $OMEGA 0.31 FIX ; 1. KA $OMEGA BLOCK(2) 0.032 ; 2. CL 0.033 0.034 FIX ; 3. V2 ;--------- D.F. for OMEGA --------------------- $THETA 150 FIX $SIGMA 1 FIX $EST PRINT=1 MAXEVAL=9999 METHOD=1 INTER POSTHOC $COV PRINT=E ========================================== Here is the nwprior.f file: ========================================== SUBROUTINE PRIOR (ICALL,CNT) DOUBLE PRECISION CNT NTHETA=11 NETA=3 NEPS=1 NTHP=11 NETP=3 NEPP=0 NPEXP=1 PLEV=0 CALL NWPRI (NTHETA,NETA,NEPS,NTHP,NETP,NEPP,NPEXP,ITYP,PLEV, 1 NSAM,ISS,CNT) RETURN END ============================= I have counted the number of leading spaces in the nwprior.f file to 6. The model is run with NONMEM v 6.2.1 using execute in PsN. Best Regards Jacob ------------------------------------------------- Jacob Brogren MSc Pharm Senior Research Scientist Pharmacometrics Clinical Pharmacology and DMPK AstraZeneca R&D Södertälje SE-151 85 Södertälje Sweden Tel. +46 8 553 236 33 Fax +46 8 553 274 91 [email protected] www.astrazeneca.com ________________________________ Confidentiality Notice: This message is private and may contain confidential and proprietary information. If you have received this message in error, please notify us and remove it from your system and note that you must not copy, distribute or take any action in reliance on it. Any unauthorized use or disclosure of the contents of this message is not permitted and may be unlawful. -------------------------------------------------------------------------- Confidentiality Notice: This message is private and may contain confidential and proprietary information. If you have received this message in error, please notify us and remove it from your system and note that you must not copy, distribute or take any action in reliance on it. Any unauthorized use or disclosure of the contents of this message is not permitted and may be unlawful.