PopED Results 

        2018-09-10 08:27:20

==============================================================================
Model description : PopED model 

Model Sizes : 
Number of individual model parameters                  g[j]    : Ng    = 5
Number of population model fixed parameters            bpop[j] : Nbpop = 4
Number of population model random effects parameters   b[j]    : Nb    = 3

Typical Population Parameters:
bpop[1]:  0.15 
bpop[2]:     8 
bpop[3]:     1 
bpop[4]:     1 

Between Subject Variability matrix D (variance units) 
0.07 0.00 0.00
0.00 0.02 0.00
0.00 0.00 0.60

Diagonal Elements of D [sqrt(param)]:
D[1,1]:  0.07 [0.2646] 
D[2,2]:  0.02 [0.1414] 
D[3,3]:   0.6 [0.7746] 

Residual Unexplained Variability matrix SIGMA (variance units) : 
0.01 0.00
0.00 0.25

Diagonal Elements of SIGMA [sqrt(param)]:
SIGMA[1,1]:  0.01 [  0.1] 
SIGMA[2,2]:  0.25 [  0.5] 

==============================================================================
Experiment description (design and design space)

Number of individuals: 32
Number of groups (individuals with same design): 1
Number of individuals per group:
     Group 1: 32
Number of samples per group:
 Number of discrete experimental variables: 0
Number of model covariates: 1

Initial Sampling Schedule
Group 1:    0.5      1      2      6     24     36     72    120

Minimum allowed sampling values
Group 1:   0.01   0.01   0.01   0.01   0.01   0.01   0.01   0.01

Maximum allowed sampling values
Group 1:    120    120    120    120    120    120    120    120

Covariates  (min, max):
Group 1: 70 (0.01, 100)

===============================================================================
Initial design evaluation

Initial OFV = 55.3964

Efficiency criterion [usually defined as OFV^(1/npar)]  = 1016.94

Initial design
expected relative standard error
(%RSE, rounded to nearest integer)
    Parameter   Values   RSE_0
      bpop[1]     0.15       5
      bpop[2]        8       3
      bpop[3]        1      14
       D[1,1]     0.07      30
       D[2,2]     0.02      37
       D[3,3]      0.6      27
   SIGMA[1,1]     0.01      32
   SIGMA[2,2]     0.25      26

==============================================================================
Criterion Specification

OFV calculation for FIM: 4 
  1=Determinant of FIM,
  4=log determinant of FIM,
  6=determinant of interesting part of FIM (Ds)

Approximation method: 0
  0=FO, 
  1=FOCE, 
  2=FOCEI, 
  3=FOI

Fisher Information Matrix type: 1
  0=Full FIM,
  1=Reduced FIM,
  2=weighted models,
  3=Loc models,
  4=reduced FIM with derivative of SD of sigma as pfim,
  5=FULL FIM parameterized with A,B,C matrices & derivative of variance,
  6=Calculate one model switch at a time, good for large matrices,
  7=Reduced FIM parameterized with A,B,C matrices & derivative of variance

Design family: 1
  D-family design (1) or 
  ED-family design (0) 
  (with or without parameter uncertainty)

==============================================================================
Optimization of design parameters

* Optimize Sampling Schedule
* Optimize Covariates

*******************************
Initial Value
 OFV(mf) = 55.3964
*******************************

RS - It. : 2   OFV : 56.1972

*******************************
RS Results
 OFV(mf) = 56.1972

Optimized Sampling Schedule
Group 1:   0.01 0.3033  2.658  6.079  7.281  46.91  87.02    120

Optimized Covariates:
Group 1: 69.948

*********************************

SG - It. : 1  OFV : 56.32   Diff. :     1
SG - It. : 2  OFV :  56.4   Diff. : 0.00131

SG - Iteration 2 --------- FINAL -------------------------
Normalized gradient: Grad_xt(OFV)/OFV
-0.00197814128670107
 -8.40339235184533e-05
 0.00861116807949714
 -0.000252583348418647
 0.00277151280773806
 6.63220120881358e-05
 1.18153736961139e-05
 -1.24891901433085e-05
xt opt:
0.01
 5.78882218016777
 0.593007944062438
 6.99089494498356
 2.94810085067416
 47.1956564701637
 87.3131069173753
 119.710318514651
Normalized gradient: Grad_a(OFV)/OFV
6.673559e-04
aopt:
7.018940e+01
OFV(mf)    : 56.3967
diff       : 0.00131013
*************************************************************
*****************************
            Line Search

group 1 -- xt[6] changed from  47.1957 to  60.005
     OFV(MF) changed from 56.3967 to 56.4447 
group 1 -- xt[6] changed from  60.005 to  120
     OFV(MF) changed from 56.4447 to 56.4594 
group 1 -- xt[1] changed from  0.01 to  60.005
     OFV(MF) changed from 56.4594 to 56.4722 
group 1 -- xt[1] changed from  60.005 to  120
     OFV(MF) changed from 56.4722 to 56.5348 

Best value for OFV(MF) = 56.5348

Best value for xt:
Group 1:  0.593  2.948  5.789  6.991  87.31  119.7    120    120

group 1 -- a[1] changed from  70.1894 to  100
     OFV(MF) changed from 56.5348 to 57.1558 
Best value for OFV(MF) = 57.1558

Best value for a: 
Group 1: 100


Line search run time: 0.193 seconds
***************************

===============================================================================
FINAL RESULTS
Optimized Sampling Schedule
Group 1:  0.593  2.948  5.789  6.991  87.31  119.7    120    120

Optimized Covariates:
Group 1: 100

 FIM: 
 1.847048e+04 -2.457849e+01 -9.535716e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
-2.457849e+01  2.008446e+01 -2.055966e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
-9.535716e+00 -2.055966e+00  4.980363e+01  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
 0.000000e+00  0.000000e+00  0.000000e+00  2.698618e+03  1.359229e+01  3.196753e-02  3.698485e+02  5.960048e+01
 0.000000e+00  0.000000e+00  0.000000e+00  1.359229e+01  2.581668e+04  4.226994e+00  1.004061e+04  8.990176e+01
 0.000000e+00  0.000000e+00  0.000000e+00  3.196753e-02  4.226994e+00  3.875628e+01  9.343952e+01  2.507505e+00
 0.000000e+00  0.000000e+00  0.000000e+00  3.698485e+02  1.004061e+04  9.343952e+01  2.313818e+05  3.511576e+03
 0.000000e+00  0.000000e+00  0.000000e+00  5.960048e+01  8.990176e+01  2.507505e+00  3.511576e+03  6.442719e+02


Inverse(FIM):
 5.423737e-05  6.772253e-05  1.318031e-05  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
 6.772253e-05  5.008559e-02  2.080572e-03  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
 1.318031e-05  2.080572e-03  2.016727e-02  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
 0.000000e+00  0.000000e+00  0.000000e+00  3.713202e-04 -4.759407e-08  2.080685e-06 -7.736619e-08 -3.392997e-05
 0.000000e+00  0.000000e+00  0.000000e+00 -4.759407e-08  3.940985e-05 -2.928598e-07 -1.773269e-06  4.171412e-06
 0.000000e+00  0.000000e+00  0.000000e+00  2.080685e-06 -2.928598e-07  2.582877e-02 -9.695160e-06 -4.783411e-05
 0.000000e+00  0.000000e+00  0.000000e+00 -7.736619e-08 -1.773269e-06 -9.695160e-06  4.795055e-06 -2.584290e-05
 0.000000e+00  0.000000e+00  0.000000e+00 -3.392997e-05  4.171412e-06 -4.783411e-05 -2.584290e-05  1.695738e-03

OFV = 57.1558

Efficiency criterion [usually defined as det(FIM)^(1/npar)]  = 1267.08

Efficiency: 
  ((exp(ofv_final) / exp(ofv_init))^(1/n_parameters)) = 1.246

Expected relative standard error
(%RSE, rounded to nearest integer):
    Parameter   Values   RSE_0   RSE
      bpop[1]     0.15       5     5
      bpop[2]        8       3     3
      bpop[3]        1      14    14
       D[1,1]     0.07       0     0
       D[2,2]     0.02      37    31
       D[3,3]      0.6       0     0
   SIGMA[1,1]     0.01      32    22
   SIGMA[2,2]     0.25      26    16

Total running time: 0.979 seconds
