-
Notifications
You must be signed in to change notification settings - Fork 5
/
Copy pathOutput_simulated_merop_PK_run3.lst
542 lines (426 loc) · 21.8 KB
/
Output_simulated_merop_PK_run3.lst
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
Mon Apr 15 17:23:02 CEST 2019
;; 1. Based on: 3
;; 2. Description: covariate model + WT on V2 + CPRED in pred3 ($TABLE)
;; x1. Author: user
$PROBLEM MEROPENEM IV INFUSION 3COMP DESCRIPTION PLASMA AND ELF
; CONCENTRATIONS FINAL COV MODEL
; ------------dataset------------
$INPUT ID WT GFR DV TIME MDV CMT AMT RATE SS II EVID
; WT [kg], GFRC [0: poor hepatic function, 1: good hepatic function],
; GFR [mL/min], DV [mg/L], TIME [h], CMT [1: plasma conc., 2: ELF conc],
; AMT (DOSE) [g], RATE (K) [g/h], II (T) [h], GRP=AMT/RATE (delta) [h]
; CL [L/h], V1 V2 V3 [L], Q2 Q3 [L/h], S1 S2 [L]
$DATA simulated_dataset.csv IGNORE=@
; ------------model------------
$SUBROUTINE ADVAN11 TRANS4
$PK
; parmacokinetic parameters
TVCL=THETA(1)*((GFR/65)**THETA(2)) ; CENTRAL
TVV1=THETA(3)*((WT/75)**THETA(4))
TVQ2=THETA(5) ; ELF
TVV2=THETA(6)*((WT/75)**THETA(9))
TVQ3=THETA(7) ; PERIPHERAL
TVV3=THETA(8)
; interindividual variance model
CL=TVCL*EXP(ETA(1))
V1=TVV1*EXP(ETA(2))
Q2=TVQ2
V2=TVV2*EXP(ETA(3))
Q3=TVQ3*EXP(ETA(4))
V3=TVV3
; scaling factor
S1=V1/1000 ; dose [g] and conc. [mg/L]
S2=V2/THETA(10)
$ERROR
; calcultate de result (i.e. model prediction)
H1=0
H2=0
IF(CMT.EQ.1) H1=1
IF(CMT.EQ.2) H2=1
Y=H1*(F*(1+EPS(1))+EPS(2))+H2*(F*(1+EPS(3))) ; +EPS(4)
W=F
IPRED=F ; prediction individuelle
IRES=DV-IPRED ; (individual-specific residual)
IWRES=IRES/W ; (individual-specific weighted residual)
; Initial conditions
$THETA (0,7.28847922347366) ; CL
$THETA 0.741447648249536 ; GFRCL
$THETA (0,14.1662561197661) ; V1
$THETA (0,0.889886352071436); WTV1
$THETA (0,7.35442828009044) ; Q2; ELF
$THETA (0,4.11640593262913) ; V2
$THETA (0,7.73152957383366) ; Q3; PERIPHERAL
$THETA (0,10.033350553562) ; V3
$THETA (0,0.944423036580702); WTV2
$THETA 251.758138000927 ; V2/S2
$OMEGA BLOCK(4)
0.130531253539076 ; IIVCL
0 0.149440414536403 ; IIVV1
0 0 1.71967661710853 ; IIVV2
0 0 0 0.171838259641997 ; IIVQ3
$SIGMA 0.0263659853631903 ; epsPROP1
0.219976467635949 ; epsADD1
0.368481652473912 ; epsPROP2
$ESTIMATION METHOD=1 INTER MAXEVAL=9999 SIGDIGITS=3 POSTHOC PRINT=5 ; PRINT=1 ;
; precision des estimation? standard error of estimates & matrice de correlation
$COVARIANCE PRINT=E
$TABLE ID WT GFR TIME MDV CMT RATE PRED CPRED RES WRES IPRED IRES
IWRES CWRES EVID NOPRINT FILE=pred4new ; population and individual predictions and residuals
$TABLE ID WT GFR CL V1 Q2 V2 Q3 V3 ETA(1) ETA(2) ETA(3) ETA(4)
TVCL TVV1 TVQ2 TVV2 TVQ3 TVV3 NOPRINT NOAPPEND FIRSTONLY
FILE=param4new ; individual PK parameters (bayesiens,posthoc)
NM-TRAN MESSAGES
WARNINGS AND ERRORS (IF ANY) FOR PROBLEM 1
(WARNING 2) NM-TRAN INFERS THAT THE DATA ARE POPULATION.
License Registered to: Universite de Namur
Expiration Date: 14 MAR 2020
Current Date: 15 APR 2019
Days until program expires : 334
1NONLINEAR MIXED EFFECTS MODEL PROGRAM (NONMEM) VERSION 7.4.3
ORIGINALLY DEVELOPED BY STUART BEAL, LEWIS SHEINER, AND ALISON BOECKMANN
CURRENT DEVELOPERS ARE ROBERT BAUER, ICON DEVELOPMENT SOLUTIONS,
AND ALISON BOECKMANN. IMPLEMENTATION, EFFICIENCY, AND STANDARDIZATION
PERFORMED BY NOUS INFOSYSTEMS.
PROBLEM NO.: 1
MEROPENEM IV INFUSION 3COMP DESCRIPTION PLASMA AND ELF
0DATA CHECKOUT RUN: NO
DATA SET LOCATED ON UNIT NO.: 2
THIS UNIT TO BE REWOUND: NO
NO. OF DATA RECS IN DATA SET: 898
NO. OF DATA ITEMS IN DATA SET: 12
ID DATA ITEM IS DATA ITEM NO.: 1
DEP VARIABLE IS DATA ITEM NO.: 4
MDV DATA ITEM IS DATA ITEM NO.: 6
0INDICES PASSED TO SUBROUTINE PRED:
12 5 8 9 10 11 7 0 0 0 0
0LABELS FOR DATA ITEMS:
ID WT GFR DV TIME MDV CMT AMT RATE SS II EVID
0(NONBLANK) LABELS FOR PRED-DEFINED ITEMS:
TVCL TVV1 TVQ2 TVV2 TVQ3 TVV3 CL V1 Q2 V2 Q3 V3 IPRED IRES IWRES
0FORMAT FOR DATA:
(E3.0,E5.0,E6.0,E8.0,E4.0,7E2.0)
TOT. NO. OF OBS RECS: 433
TOT. NO. OF INDIVIDUALS: 60
0LENGTH OF THETA: 10
0DEFAULT THETA BOUNDARY TEST OMITTED: NO
0OMEGA HAS BLOCK FORM:
1
1 1
1 1 1
1 1 1 1
0DEFAULT OMEGA BOUNDARY TEST OMITTED: NO
0SIGMA HAS SIMPLE DIAGONAL FORM WITH DIMENSION: 3
0DEFAULT SIGMA BOUNDARY TEST OMITTED: NO
0INITIAL ESTIMATE OF THETA:
LOWER BOUND INITIAL EST UPPER BOUND
0.0000E+00 0.7288E+01 0.1000E+07
-0.1000E+07 0.7414E+00 0.1000E+07
0.0000E+00 0.1417E+02 0.1000E+07
0.0000E+00 0.8899E+00 0.1000E+07
0.0000E+00 0.7354E+01 0.1000E+07
0.0000E+00 0.4116E+01 0.1000E+07
0.0000E+00 0.7732E+01 0.1000E+07
0.0000E+00 0.1003E+02 0.1000E+07
0.0000E+00 0.9444E+00 0.1000E+07
-0.1000E+07 0.2518E+03 0.1000E+07
0INITIAL ESTIMATE OF OMEGA:
BLOCK SET NO. BLOCK FIXED
1 NO
0.1305E+00
0.0000E+00 0.1494E+00
0.0000E+00 0.0000E+00 0.1720E+01
0.0000E+00 0.0000E+00 0.0000E+00 0.1718E+00
0INITIAL ESTIMATE OF SIGMA:
0.2637E-01
0.0000E+00 0.2200E+00
0.0000E+00 0.0000E+00 0.3685E+00
0COVARIANCE STEP OMITTED: NO
EIGENVLS. PRINTED: YES
SPECIAL COMPUTATION: NO
COMPRESSED FORMAT: NO
GRADIENT METHOD USED: NOSLOW
SIGDIGITS ETAHAT (SIGLO): -1
SIGDIGITS GRADIENTS (SIGL): -1
EXCLUDE COV FOR FOCE (NOFCOV): NO
TURN OFF Cholesky Transposition of R Matrix (CHOLROFF): NO
KNUTHSUMOFF: -1
RESUME COV ANALYSIS (RESUME): NO
SIR SAMPLE SIZE (SIRSAMPLE): -1
NON-LINEARLY TRANSFORM THETAS DURING COV (THBND): 1
PRECONDTIONING CYCLES (PRECOND): 0
PRECONDTIONING TYPES (PRECONDS): TOS
FORCED PRECONDTIONING CYCLES (PFCOND):0
PRECONDTIONING TYPE (PRETYPE): 0
FORCED POS. DEFINITE SETTING: (FPOSDEF):0
0TABLES STEP OMITTED: NO
NO. OF TABLES: 2
SEED NUMBER (SEED): 11456
RANMETHOD: 3U
MC SAMPLES (ESAMPLE): 300
WRES SQUARE ROOT TYPE (WRESCHOL): EIGENVALUE
0-- TABLE 1 --
0RECORDS ONLY: ALL
04 COLUMNS APPENDED: YES
PRINTED: NO
HEADER: YES
FILE TO BE FORWARDED: NO
FORMAT: S1PE11.4
LFORMAT:
RFORMAT:
FIXED_EFFECT_ETAS:
0USER-CHOSEN ITEMS:
ID WT GFR TIME MDV CMT RATE CPRED IPRED IRES IWRES CWRES EVID
0-- TABLE 2 --
0RECORDS ONLY: FIRSTONLY
04 COLUMNS APPENDED: NO
PRINTED: NO
HEADER: YES
FILE TO BE FORWARDED: NO
FORMAT: S1PE11.4
LFORMAT:
RFORMAT:
FIXED_EFFECT_ETAS:
0USER-CHOSEN ITEMS:
ID WT GFR CL V1 Q2 V2 Q3 V3 ETA1 ETA2 ETA3 ETA4 TVCL TVV1 TVQ2 TVV2 TVQ3 TVV3
1DOUBLE PRECISION PREDPP VERSION 7.4.3
THREE COMPARTMENT MODEL (ADVAN11)
0MAXIMUM NO. OF BASIC PK PARAMETERS: 6
0BASIC PK PARAMETERS (AFTER TRANSLATION):
BASIC PK PARAMETER NO. 1: ELIMINATION RATE (K)
BASIC PK PARAMETER NO. 2: CENTRAL-TO-PERIPH1 RATE (K12)
BASIC PK PARAMETER NO. 3: PERIPH1-TO-CENTRAL RATE (K21)
BASIC PK PARAMETER NO. 4: CENTRAL-TO-PERIPH2 RATE (K13)
BASIC PK PARAMETER NO. 5: PERIPH2-TO-CENTRAL RATE (K31)
TRANSLATOR WILL CONVERT PARAMETERS
CL, V1, Q2, V2, Q3, V3 TO K, K12, K21, K13, K31 (TRANS4)
0COMPARTMENT ATTRIBUTES
COMPT. NO. FUNCTION INITIAL ON/OFF DOSE DEFAULT DEFAULT
STATUS ALLOWED ALLOWED FOR DOSE FOR OBS.
1 CENTRAL ON NO YES YES YES
2 PERIPH1 ON NO YES NO NO
3 PERIPH2 ON NO YES NO NO
4 OUTPUT OFF YES NO NO NO
1
ADDITIONAL PK PARAMETERS - ASSIGNMENT OF ROWS IN GG
COMPT. NO. INDICES
SCALE BIOAVAIL. ZERO-ORDER ZERO-ORDER ABSORB
FRACTION RATE DURATION LAG
1 7 * * * *
2 8 * * * *
3 * * * * *
4 * - - - -
- 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.: 12
TIME DATA ITEM IS DATA ITEM NO.: 5
DOSE AMOUNT DATA ITEM IS DATA ITEM NO.: 8
DOSE RATE DATA ITEM IS DATA ITEM NO.: 9
STEADY STATE DATA ITEM IS DATA ITEM NO.: 10
INTERVAL DATA ITEM IS DATA ITEM NO.: 11
COMPT. NO. DATA ITEM IS DATA ITEM NO.: 7
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.
1
#TBLN: 1
#METH: First Order Conditional Estimation with Interaction
ESTIMATION STEP OMITTED: NO
ANALYSIS TYPE: POPULATION
NUMBER OF SADDLE POINT RESET ITERATIONS: 0
GRADIENT METHOD USED: NOSLOW
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
IND. OBJ. FUNC. VALUES SORTED: NO
NUMERICAL DERIVATIVE
FILE REQUEST (NUMDER): NONE
MAP (ETAHAT) ESTIMATION METHOD (OPTMAP): 0
ETA HESSIAN EVALUATION METHOD (ETADER): 0
INITIAL ETA FOR MAP ESTIMATION (MCETA): 0
SIGDIGITS FOR MAP ESTIMATION (SIGLO): 100
GRADIENT SIGDIGITS OF
FIXED EFFECTS PARAMETERS (SIGL): 100
NOPRIOR SETTING (NOPRIOR): OFF
NOCOV SETTING (NOCOV): OFF
DERCONT SETTING (DERCONT): OFF
FINAL ETA RE-EVALUATION (FNLETA): ON
EXCLUDE NON-INFLUENTIAL (NON-INFL.) ETAS
IN SHRINKAGE (ETASTYPE): NO
NON-INFL. ETA CORRECTION (NONINFETA): OFF
RAW OUTPUT FILE (FILE): psn.ext
EXCLUDE TITLE (NOTITLE): NO
EXCLUDE COLUMN LABELS (NOLABEL): NO
FORMAT FOR ADDITIONAL FILES (FORMAT): S1PE12.5
PARAMETER ORDER FOR OUTPUTS (ORDER): TSOL
WISHART PRIOR DF INTERPRETATION (WISHTYPE):0
KNUTHSUMOFF: 0
INCLUDE LNTWOPI: NO
INCLUDE CONSTANT TERM TO PRIOR (PRIORC): NO
INCLUDE CONSTANT TERM TO OMEGA (ETA) (OLNTWOPI):NO
ADDITIONAL CONVERGENCE TEST (CTYPE=4)?: NO
EM OR BAYESIAN METHOD USED: NONE
THE FOLLOWING LABELS ARE EQUIVALENT
PRED=PREDI
RES=RESI
WRES=WRESI
IWRS=IWRESI
IPRD=IPREDI
IRS=IRESI
MONITORING OF SEARCH:
0ITERATION NO.: 0 OBJECTIVE VALUE: 1521.92722102785 NO. OF FUNC. EVALS.: 15
CUMULATIVE NO. OF FUNC. EVALS.: 15
NPARAMETR: 7.2885E+00 7.4145E-01 1.4166E+01 8.8989E-01 7.3544E+00 4.1164E+00 7.7315E+00 1.0033E+01 9.4442E-01 2.5176E+02
1.3053E-01 1.4944E-01 1.7197E+00 1.7184E-01 2.6366E-02 2.1998E-01 3.6848E-01
PARAMETER: 1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01
1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01 1.0000E-01
GRADIENT: -1.1068E+02 -4.2828E+01 -3.9451E+01 -9.9723E+00 1.1085E+01 -1.5085E+00 -1.5360E+01 5.8721E+01 -3.1456E+00 -5.1710E+01
8.1158E+00 1.9005E+01 1.9216E+01 6.3121E+00 -8.4011E+00 5.0305E+00 -5.3393E-01
0ITERATION NO.: 5 OBJECTIVE VALUE: 1502.67971788807 NO. OF FUNC. EVALS.: 82
CUMULATIVE NO. OF FUNC. EVALS.: 97
NPARAMETR: 8.6094E+00 7.1339E-01 1.5277E+01 9.6899E-01 6.7561E+00 4.2798E+00 7.9967E+00 8.0926E+00 9.7133E-01 2.4939E+02
1.0847E-01 1.1404E-01 1.3470E+00 1.5251E-01 2.6669E-02 2.0300E-01 3.5486E-01
PARAMETER: 2.6656E-01 9.6216E-02 1.7549E-01 1.8516E-01 1.5138E-02 1.3894E-01 1.3373E-01 -1.1497E-01 1.2810E-01 9.9058E-02
7.4510E-03 -3.5165E-02 -2.2132E-02 4.0345E-02 1.0571E-01 5.9855E-02 8.1161E-02
GRADIENT: 5.0302E+01 -7.2331E+01 -5.2819E+00 -1.0223E+01 4.0387E+00 -9.0166E+00 6.2115E+00 -1.5462E+01 -3.4401E+00 -6.7748E+01
1.3826E+00 7.2962E+00 8.1247E+00 5.7442E+00 -1.1119E+00 2.3071E-02 -2.8154E-01
0ITERATION NO.: 10 OBJECTIVE VALUE: 1493.45477299108 NO. OF FUNC. EVALS.: 80
CUMULATIVE NO. OF FUNC. EVALS.: 177
NPARAMETR: 8.5618E+00 7.1304E-01 1.5206E+01 1.3740E+00 9.0886E+00 6.6413E+00 5.6481E+00 7.5094E+00 1.7464E+00 2.6388E+02
1.0406E-01 1.0389E-01 3.4671E-01 4.6296E-02 2.6521E-02 3.0404E-01 3.0036E-01
PARAMETER: 2.6102E-01 9.6168E-02 1.7085E-01 5.3439E-01 3.1172E-01 5.7832E-01 -2.1399E-01 -1.8976E-01 7.1475E-01 1.0481E-01
-1.3330E-02 -8.1791E-02 -7.0070E-01 -5.5575E-01 1.0293E-01 2.6181E-01 -2.2013E-03
GRADIENT: 4.3530E+01 -9.5716E+01 1.1884E+01 8.1816E+00 4.1203E+00 1.5949E+00 2.8705E+00 7.3020E+00 -2.8925E+00 -1.4428E+02
-4.8437E+00 1.8031E+00 -4.1494E+00 8.5714E-01 -6.7771E+00 1.4025E+01 -8.7364E+00
0ITERATION NO.: 15 OBJECTIVE VALUE: 1490.61413521149 NO. OF FUNC. EVALS.: 81
CUMULATIVE NO. OF FUNC. EVALS.: 258
NPARAMETR: 8.4555E+00 7.4305E-01 1.5165E+01 1.1970E+00 8.8290E+00 6.4683E+00 5.7279E+00 7.5589E+00 1.8333E+00 2.7119E+02
1.0839E-01 1.0890E-01 4.6400E-01 9.0101E-03 2.6847E-02 2.5821E-01 3.2888E-01
PARAMETER: 2.4852E-01 1.0022E-01 1.6814E-01 3.9647E-01 2.8274E-01 5.5193E-01 -1.9997E-01 -1.8319E-01 7.6329E-01 1.0772E-01
7.0777E-03 -5.8216E-02 -5.5500E-01 -1.3741E+00 1.0903E-01 1.8012E-01 4.3147E-02
GRADIENT: 3.0413E+01 2.8694E+01 2.1496E+00 -1.0740E+00 7.7790E-01 2.8661E+00 -8.2499E-01 5.9800E+00 -2.5045E+00 5.4485E+01
1.4961E+00 5.4445E+00 2.2327E+00 2.1267E-01 -3.0160E+00 3.5853E+00 3.5476E+00
0ITERATION NO.: 20 OBJECTIVE VALUE: 1489.66348132977 NO. OF FUNC. EVALS.: 81
CUMULATIVE NO. OF FUNC. EVALS.: 339
NPARAMETR: 8.2078E+00 7.4020E-01 1.5180E+01 1.1930E+00 8.2202E+00 5.9170E+00 6.0691E+00 7.6891E+00 2.1998E+00 2.6921E+02
1.0690E-01 9.6701E-02 4.5891E-01 8.4018E-03 2.7400E-02 2.4580E-01 3.1564E-01
PARAMETER: 2.1879E-01 9.9832E-02 1.6910E-01 3.9313E-01 2.1129E-01 4.6285E-01 -1.4210E-01 -1.6611E-01 9.4554E-01 1.0693E-01
1.2416E-04 -1.1764E-01 -5.6052E-01 -1.4091E+00 1.1924E-01 1.5549E-01 2.2612E-02
GRADIENT: -2.9645E-01 3.4623E+00 -3.0699E-01 -3.3537E-01 -1.8480E-01 -1.6735E-01 -3.2983E-01 -6.1281E-02 1.9294E-02 8.0821E+00
3.4757E-01 7.1304E-02 1.5247E-01 2.1254E-01 6.4120E-01 -2.2148E-04 2.5363E-01
0ITERATION NO.: 25 OBJECTIVE VALUE: 1489.55913928762 NO. OF FUNC. EVALS.: 107
CUMULATIVE NO. OF FUNC. EVALS.: 446
NPARAMETR: 8.2140E+00 7.3989E-01 1.5162E+01 1.2032E+00 8.2228E+00 5.9268E+00 6.0864E+00 7.6999E+00 2.2100E+00 2.6867E+02
1.0649E-01 9.6767E-02 4.5987E-01 3.6901E-04 2.7382E-02 2.4651E-01 3.1541E-01
PARAMETER: 2.1955E-01 9.9789E-02 1.6793E-01 4.0161E-01 2.1161E-01 4.6450E-01 -1.3925E-01 -1.6470E-01 9.5019E-01 1.0672E-01
-1.7957E-03 -1.1730E-01 -5.5947E-01 -2.9717E+00 1.1891E-01 1.5694E-01 2.2240E-02
GRADIENT: -4.0163E+00 -4.6498E+01 -8.9421E-01 7.6512E-02 -1.4165E-01 2.7651E-02 -2.3505E-01 8.1752E-02 1.3079E-01 -1.5283E+00
-8.4631E-02 -2.3000E-02 -1.0997E-02 9.1996E-03 7.7080E-02 6.7505E-02 -7.4985E-03
0ITERATION NO.: 30 OBJECTIVE VALUE: 1489.51334345201 NO. OF FUNC. EVALS.: 191
CUMULATIVE NO. OF FUNC. EVALS.: 637
NPARAMETR: 8.2385E+00 7.5069E-01 1.5183E+01 1.2026E+00 8.2259E+00 5.9232E+00 6.0968E+00 7.7038E+00 2.2007E+00 2.6881E+02
1.0652E-01 9.6888E-02 4.6209E-01 1.7184E-05 2.7380E-02 2.4636E-01 3.1523E-01
PARAMETER: 2.2252E-01 1.0125E-01 1.6933E-01 4.0117E-01 2.1198E-01 4.6389E-01 -1.3754E-01 -1.6420E-01 9.4595E-01 1.0677E-01
-1.6324E-03 -1.1667E-01 -5.5707E-01 -4.5829E+00 1.1887E-01 1.5665E-01 2.1953E-02
GRADIENT: -1.1332E-01 -2.8278E-01 3.1412E-02 1.0000E-02 7.5806E-03 -5.9031E-03 -9.9965E-03 2.4286E-02 7.8972E-03 -2.8697E-01
-1.6040E-02 -2.0025E-04 1.2800E-03 0.0000E+00 -3.6162E-02 1.4409E-02 -4.4356E-03
0ITERATION NO.: 33 OBJECTIVE VALUE: 1489.51332594737 NO. OF FUNC. EVALS.: 102
CUMULATIVE NO. OF FUNC. EVALS.: 739
NPARAMETR: 8.2394E+00 7.5074E-01 1.5183E+01 1.2025E+00 8.2285E+00 5.9264E+00 6.0947E+00 7.7014E+00 2.1993E+00 2.6883E+02
1.0654E-01 9.6883E-02 4.6178E-01 1.7184E-05 2.7383E-02 2.4630E-01 3.1521E-01
PARAMETER: 2.2263E-01 1.0125E-01 1.6928E-01 4.0105E-01 2.1230E-01 4.6444E-01 -1.3788E-01 -1.6451E-01 9.4532E-01 1.0678E-01
-1.5553E-03 -1.1670E-01 -5.5740E-01 -4.5601E+00 1.1892E-01 1.5651E-01 2.1925E-02
GRADIENT: 7.0950E-04 8.9136E-03 8.1387E-04 -6.0649E-04 2.5205E-05 -1.4630E-04 -8.0484E-04 1.8746E-03 -4.8299E-05 -6.6194E-03
1.8365E-03 -5.7066E-04 2.7411E-04 0.0000E+00 -7.9842E-04 -9.8131E-05 -1.6826E-04
#TERM:
0MINIMIZATION SUCCESSFUL
NO. OF FUNCTION EVALUATIONS USED: 739
NO. OF SIG. DIGITS IN FINAL EST.: 3.6
0PARAMETER ESTIMATE IS NEAR ITS BOUNDARY
THIS MUST BE ADDRESSED BEFORE THE COVARIANCE STEP CAN BE IMPLEMENTED
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: 3.7198E-04 -6.7102E-03 1.5416E-02 -1.7335E-06
SE: 4.1181E-02 3.1359E-02 5.2205E-02 1.9963E-06
N: 60 60 60 60
P VAL.: 9.9279E-01 8.3056E-01 7.6776E-01 3.8519E-01
ETASHRINKSD(%) 1.4459E+00 2.1301E+01 3.9991E+01 9.9624E+01
ETASHRINKVR(%) 2.8708E+00 3.8065E+01 6.3989E+01 9.9999E+01
EBVSHRINKSD(%) 2.1758E+00 2.1566E+01 3.8644E+01 9.9588E+01
EBVSHRINKVR(%) 4.3043E+00 3.8481E+01 6.2354E+01 9.9998E+01
EPSSHRINKSD(%) 1.4402E+01 1.4402E+01 4.9416E+00
EPSSHRINKVR(%) 2.6730E+01 2.6730E+01 9.6390E+00
TOTAL DATA POINTS NORMALLY DISTRIBUTED (N): 433
N*LOG(2PI) CONSTANT TO OBJECTIVE FUNCTION: 795.80076975524651
OBJECTIVE FUNCTION VALUE WITHOUT CONSTANT: 1489.5133259473719
OBJECTIVE FUNCTION VALUE WITH CONSTANT: 2285.3140957026185
REPORTED OBJECTIVE FUNCTION DOES NOT CONTAIN CONSTANT
TOTAL EFFECTIVE ETAS (NIND*NETA): 240
#TERE:
Elapsed estimation time in seconds: 19.04
Elapsed postprocess time in seconds: 0.09
1
************************************************************************************************************************
******************** ********************
******************** FIRST ORDER CONDITIONAL ESTIMATION WITH INTERACTION ********************
#OBJT:************** MINIMUM VALUE OF OBJECTIVE FUNCTION ********************
******************** ********************
************************************************************************************************************************
#OBJV:******************************************** 1489.513 **************************************************
1
************************************************************************************************************************
******************** ********************
******************** FIRST ORDER CONDITIONAL ESTIMATION WITH INTERACTION ********************
******************** FINAL PARAMETER ESTIMATE ********************
******************** ********************
************************************************************************************************************************
THETA - VECTOR OF FIXED EFFECTS PARAMETERS *********
TH 1 TH 2 TH 3 TH 4 TH 5 TH 6 TH 7 TH 8 TH 9 TH10
8.24E+00 7.51E-01 1.52E+01 1.20E+00 8.23E+00 5.93E+00 6.09E+00 7.70E+00 2.20E+00 2.69E+02
OMEGA - COV MATRIX FOR RANDOM EFFECTS - ETAS ********
ETA1 ETA2 ETA3 ETA4
ETA1
+ 1.07E-01
ETA2
+ 0.00E+00 9.69E-02
ETA3
+ 0.00E+00 0.00E+00 4.62E-01
ETA4
+ 0.00E+00 0.00E+00 0.00E+00 1.72E-05
SIGMA - COV MATRIX FOR RANDOM EFFECTS - EPSILONS ****
EPS1 EPS2 EPS3
EPS1
+ 2.74E-02
EPS2
+ 0.00E+00 2.46E-01
EPS3
+ 0.00E+00 0.00E+00 3.15E-01
1
OMEGA - CORR MATRIX FOR RANDOM EFFECTS - ETAS *******
ETA1 ETA2 ETA3 ETA4
ETA1
+ 3.26E-01
ETA2
+ 0.00E+00 3.11E-01
ETA3
+ 0.00E+00 0.00E+00 6.80E-01
ETA4
+ 0.00E+00 0.00E+00 0.00E+00 4.15E-03
SIGMA - CORR MATRIX FOR RANDOM EFFECTS - EPSILONS ***
EPS1 EPS2 EPS3
EPS1
+ 1.65E-01
EPS2
+ 0.00E+00 4.96E-01
EPS3
+ 0.00E+00 0.00E+00 5.61E-01
Elapsed finaloutput time in seconds: 0.04
#CPUT: Total CPU Time in Seconds, 19.192
Stop Time:
Mon Apr 15 17:23:21 CEST 2019