Wed Apr  9 13:28:36 CEST 2008
$PROBLEM PSMValidation

$INPUT ID TIME EVID DV CMT AMT

$DATA NONMEMinputdata.csv  IGNORE "

$SUBROUTINE ADVAN6 TOL 6 DP

$MODEL
COMP=(CENTRAL)
COMP=(PERIPH)


$THETA (0 0.6) ; W
$THETA (0 5 50) ; V1

$OMEGA 0.5     ; V1

$SIGMA 1 FIX

$PK
K1 = 0.005
K2 = 0.005
K3 = 0.002  

V1 = THETA(2) * EXP(ETA(1))

$DES
DADT(1)= -(K1+K3)*A(1) + K2*A(2)
DADT(2)= K1*A(1)-K2*A(2) 

$ERROR
IPRED=A(1)/V1
IRES=DV-IPRED
W=THETA(1)
IWRES=IRES/W
Y=IPRED+W*EPS(1)

$ESTIMATION MAXEVALS=9999 METHOD=1 INTERACTION NOABORT SIGDIGITS=3 PRINT=5

$COVARIANCE

$TABLE ID TIME EVID IPRED IRES IWRES W ETA(1) 
       NOPRINT ONEHEADER FILE=TABLE.PAR

 
 
1NONLINEAR MIXED EFFECTS MODEL PROGRAM (NONMEM)    DOUBLE PRECISION NONMEM    VERSION VI LEVEL 1.2  
 DEVELOPED AND PROGRAMMED BY STUART BEAL AND LEWIS SHEINER
 
 PROBLEM NO.:         1
 PSMValidation                                                           
0DATA CHECKOUT RUN:              NO 
 DATA SET LOCATED ON UNIT NO.:    2
 THIS UNIT TO BE REWOUND:        NO 
 NO. OF DATA RECS IN DATA SET:  440
 NO. OF DATA ITEMS IN DATA SET:   7
 ID DATA ITEM IS DATA ITEM NO.:   1
 DEP VARIABLE IS DATA ITEM NO.:   4
 MDV DATA ITEM IS DATA ITEM NO.:  7
0INDICES PASSED TO SUBROUTINE PRED:
  3  2  6  0  0  0  5  0  0
  0  0
0LABELS FOR DATA ITEMS:
   ID    TIME    EVID      DV     CMT     AMT     MDV
0(NONBLANK) LABELS FOR PRED-DEFINED ITEMS:
 IPRE    IRES       W    IWRE
0FORMAT FOR DATA:
 (6E9.0,1F2.0)                                                                   
 
 TOT. NO. OF OBS RECS:     410
 TOT. NO. OF INDIVIDUALS:   10
0LENGTH OF THETA:  2
0OMEGA HAS SIMPLE DIAGONAL FORM WITH DIMENSION:  1
0SIGMA HAS SIMPLE DIAGONAL FORM WITH DIMENSION:  1
0INITIAL ESTIMATE OF THETA:
 LOWER BOUND    INITIAL EST    UPPER BOUND
  0.0000E+00     0.6000E+00     0.1000E+07
  0.0000E+00     0.5000E+01     0.5000E+02
0INITIAL ESTIMATE OF OMEGA:
 0.5000E+00
0INITIAL ESTIMATE OF SIGMA:
 0.1000E+01
0SIGMA CONSTRAINED TO BE THIS INITIAL ESTIMATE
0ESTIMATION STEP OMITTED:           NO 
 CONDITIONAL ESTIMATES USED:       YES 
 CENTERED ETA:                      NO 
 EPS-ETA INTERACTION:              YES 
 LAPLACIAN OBJ. FUNC.:              NO 
 NO. OF FUNCT. EVALS. ALLOWED:    9999
 NO. OF SIG. FIGURES REQUIRED:       3
 INTERMEDIATE PRINTOUT:            YES 
 ESTIMATE OUTPUT TO MSF:            NO 
 ABORT WITH PRED EXIT CODE 1:       NO 
 IND. OBJ. FUNC. VALUES SORTED:     NO 
0COVARIANCE STEP OMITTED:    NO 
 EIGENVLS. PRINTED:          NO 
 SPECIAL COMPUTATION:        NO 
 COMPRESSED FORMAT:          NO 
0TABLES STEP OMITTED:    NO 
 NO. OF TABLES:           1
0-- TABLE  1 --
 PRINTED:                NO 
 HEADER:                YES 
 FILE TO BE FORWARDED:   NO 
0USER-CHOSEN ITEMS 
 IN THE ORDER THEY WILL APPEAR IN THE TABLE:
   ID    TIME    EVID    IPRE    IRES    IWRE       W    ETA1
1DOUBLE PRECISION PREDPP VERSION V LEVEL 1.0  
 
 GENERAL NONLINEAR KINETICS MODEL (ADVAN6)
0MODEL SUBROUTINE USER-SUPPLIED - ID NO. 9999
0MAXIMUM NO. OF BASIC PK PARAMETERS:   3
0COMPARTMENT ATTRIBUTES 
 COMPT. NO.   FUNCTION   INITIAL    ON/OFF      DOSE      DEFAULT    DEFAULT
                         STATUS     ALLOWED    ALLOWED    FOR DOSE   FOR OBS.
    1         CENTRAL      ON         YES        YES        YES        YES
    2         PERIPH       ON         YES        YES        NO         NO 
    3         OUTPUT       OFF        YES        NO         NO         NO 
0NRD VALUE FROM SUBROUTINE TOL:   6
1
 ADDITIONAL PK PARAMETERS - ASSIGNMENT OF ROWS IN GG
 COMPT. NO.                             INDICES
              SCALE      BIOAVAIL.   ZERO-ORDER  ZERO-ORDER  ABSORB
                         FRACTION    RATE        DURATION    LAG  
    1           *           *           *           *           *
    2           *           *           *           *           *
    3           *           -           -           -           -
             - PARAMETER IS NOT ALLOWED FOR THIS MODEL
             * PARAMETER IS NOT SUPPLIED BY PK SUBROUTINE;
               WILL DEFAULT TO ONE IF APPLICABLE
0DATA ITEM INDICES USED BY PRED ARE:
   EVENT ID DATA ITEM IS DATA ITEM NO.:      3
   TIME DATA ITEM IS DATA ITEM NO.:          2
   DOSE AMOUNT DATA ITEM IS DATA ITEM NO.:   6
   COMPT. NO. DATA ITEM IS DATA ITEM NO.:    5
 
0PK SUBROUTINE CALLED WITH EVERY EVENT RECORD.
 PK SUBROUTINE NOT CALLED AT NONEVENT (ADDITIONAL OR LAGGED) DOSE TIMES.
0ERROR SUBROUTINE CALLED WITH EVERY EVENT RECORD.
0ERROR SUBROUTINE INDICATES THAT DERIVATIVES OF COMPARTMENT AMOUNTS ARE USED.
0DES SUBROUTINE USES COMPACT STORAGE MODE.
1
 MONITORING OF SEARCH:

0ITERATION NO.:    0    OBJECTIVE VALUE:  0.90026E+05    NO. OF FUNC. EVALS.: 4
 CUMULATIVE NO. OF FUNC. EVALS.:    4
 PARAMETER:  0.1000E+00  0.1000E+00  0.1000E+00
 GRADIENT:  -0.1798E+06  0.1036E+05  0.4645E+01
0ITERATION NO.:    5    OBJECTIVE VALUE:  0.23065E+04    NO. OF FUNC. EVALS.: 7
 CUMULATIVE NO. OF FUNC. EVALS.:   40
 PARAMETER:  0.2809E+01  0.1045E+00 -0.3769E-01
 GRADIENT:  -0.3693E+00  0.4974E+00 -0.2237E+00
0ITERATION NO.:   10    OBJECTIVE VALUE:  0.23065E+04    NO. OF FUNC. EVALS.: 7
 CUMULATIVE NO. OF FUNC. EVALS.:   76
 PARAMETER:  0.2810E+01  0.9302E-01 -0.3236E-01
 GRADIENT:   0.2042E-01  0.6421E-02 -0.3603E-02
0ITERATION NO.:   13    OBJECTIVE VALUE:  0.23065E+04    NO. OF FUNC. EVALS.: 4
 CUMULATIVE NO. OF FUNC. EVALS.:   94
 PARAMETER:  0.2810E+01  0.9283E-01 -0.3225E-01
 GRADIENT:  -0.3577E-02 -0.1625E-02  0.8766E-03
0MINIMIZATION SUCCESSFUL
 NO. OF FUNCTION EVALUATIONS USED:   94
 NO. OF SIG. DIGITS IN FINAL EST.:  3.4

 ETABAR IS THE ARITHMETIC MEAN OF THE ETA-ESTIMATES,
 AND THE P-VALUE IS GIVEN FOR THE NULL HYPOTHESIS THAT THE TRUE MEAN IS 0.

 ETABAR:  -0.13E-03
 SE:       0.20E+00

 P VAL.:   0.10E+01
1
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 ************************************************************************************************************************
 ********************                                                                                ********************
 ********************                           MINIMUM VALUE OF OBJECTIVE FUNCTION                  ********************
 ********************                                                                                ********************
 ************************************************************************************************************************
 





 **************************************************       2306.517     **************************************************
1
 ************************************************************************************************************************
 ********************                                                                                ********************
 ********************                                  FINAL PARAMETER ESTIMATE                      ********************
 ********************                                                                                ********************
 ************************************************************************************************************************
 


 THETA - VECTOR OF FIXED EFFECTS PARAMETERS   *********


            TH 1      TH 2
 
         9.01E+00  4.97E+00
 


 OMEGA - COV MATRIX FOR RANDOM EFFECTS - ETAS  ********


            ETA1
 
 ETA1
+        3.84E-01
 


 SIGMA - COV MATRIX FOR RANDOM EFFECTS - EPSILONS  ****


            EPS1
 
 EPS1
+        1.00E+00
 
1
 ************************************************************************************************************************
 ********************                                                                                ********************
 ********************                             STANDARD ERROR OF ESTIMATE                         ********************
 ********************                                                                                ********************
 ************************************************************************************************************************
 


 THETA - VECTOR OF FIXED EFFECTS PARAMETERS   *********


            TH 1      TH 2
 
         8.48E-01  9.75E-01
 


 OMEGA - COV MATRIX FOR RANDOM EFFECTS - ETAS  ********


            ETA1
 
 ETA1
+        1.12E-01
 


 SIGMA - COV MATRIX FOR RANDOM EFFECTS - EPSILONS  ****


            EPS1
 
 EPS1
+       .........
 
1
 ************************************************************************************************************************
 ********************                                                                                ********************
 ********************                          COVARIANCE MATRIX OF ESTIMATE                         ********************
 ********************                                                                                ********************
 ************************************************************************************************************************
 

            TH 1      TH 2      OM11      SG11
 
 TH 1
+        7.18E-01
 
 TH 2
+       -6.53E-01  9.50E-01
 
 OM11
+        2.42E-02  2.22E-02  1.25E-02
 
 SG11
+       ......... ......... ......... .........
 
1
 ************************************************************************************************************************
 ********************                                                                                ********************
 ********************                         CORRELATION MATRIX OF ESTIMATE                         ********************
 ********************                                                                                ********************
 ************************************************************************************************************************
 

            TH 1      TH 2      OM11      SG11
 
 TH 1
+        1.00E+00
 
 TH 2
+       -7.91E-01  1.00E+00
 
 OM11
+        2.56E-01  2.04E-01  1.00E+00
 
 SG11
+       ......... ......... ......... .........
 
1
 ************************************************************************************************************************
 ********************                                                                                ********************
 ********************                  INVERSE COVARIANCE MATRIX OF ESTIMATE                         ********************
 ********************                                                                                ********************
 ************************************************************************************************************************
 

            TH 1      TH 2      OM11      SG11
 
 TH 1
+        7.22E+00
 
 TH 2
+        5.52E+00  5.32E+00
 
 OM11
+       -2.39E+01 -2.02E+01  1.63E+02
 
 SG11
+       ......... ......... ......... .........
 
Wed Apr  9 13:28:38 CEST 2008
