Case label: Pachepsky. Additive model with evaporation. All SWRC parameters except alpha fixed. alpha and saturated K log-fitted.   
Parameterization = ADV.
Nr of fixed and fitting parameters: 9. Seed1 = 98877665. LogK = T.
The maximum relative difference in the objective function between consecutive shuffles is  1.000E-06.
The maximum allowed value or the root mean square error of the best fit is:  1.000E-02.
The maximum number of convergence criteria any parameter may fail is 7.
Maximum number of evaluations of the objective function is 3000. Few complexes (T/F): T.
The best map points (T) or random points (F) fill the first complexes: T.

There are 26 observations. K is given as a function of the matric potential (T/F): T.
Temperature:   20.0     degrees Celcius.

There will be 3 parameter optimization runs performed.

The numbers of the flexible fitting parameters:   4  6  7  8  9

The cumulative probability density function of the Root Mean Square Error in the parameter space (  9.999000E+03 is the error value)
          RMSE      cdf
_______________________
  1.841647E-01   0.0000
  2.245686E-01   0.0001
  3.741769E-01   0.0010
  4.930027E-01   0.0100
  8.496346E-01   0.0500
  1.024536E+00   0.1000
  1.222132E+00   0.1500
  1.425398E+00   0.2000
  1.521475E+00   0.2500
  1.705450E+00   0.3000
  1.991854E+00   0.3500
  2.372777E+00   0.4000
  2.721537E+00   0.4500
  2.820863E+00   0.5000
  2.953368E+00   0.5500
  3.245820E+00   0.6000
  3.964017E+00   0.6500
  4.207346E+00   0.7000
  4.216956E+00   0.7500
  4.611498E+00   0.8000
  5.054093E+00   0.8500
  5.328468E+00   0.9000
  5.777869E+00   0.9500
  6.869744E+00   0.9900
  6.872544E+00   0.9990
  6.873540E+00   0.9999
  6.874783E+00   1.0000

Average RMSE:   3.036129E+00. Standard deviation of the RMSE:   1.649791E+00.

Histogram of the mapped RMSE-values. Parameter 0 refers to the RMSE.

The histogram of parameter nr. 0.

Bin height:             0.0000    0.2608E-01    0.7779E-01    0.2427E+00    0.2445E+00    0.3484E+00    0.2419E+00    0.1642E+00    0.1224E+00    0.1495E+00    0.3334E+00    0.2346E+00    0.1176E+00    0.4573E-01    0.6123E-01    0.3796E+00
Upper bin limit:    0.1842E+00    0.4530E+00    0.7219E+00    0.9907E+00    0.1260E+01    0.1528E+01    0.1797E+01    0.2066E+01    0.2335E+01    0.2604E+01    0.2873E+01    0.3142E+01    0.3410E+01    0.3679E+01    0.3948E+01    0.4217E+01
Out of 59049 data points, 14763 are beyond the last bin.

Convergence criteria:
1. in the best fits from the last set of shuffles, the range of the parameter does not exceed both the absolute and the relative user-specified tolerance
2. in the best fits from the last set of shuffles, the range of the objective function does not exceed its absolute user-specified tolerance
3. the parameter range in the final complexes does not exceed the maximum internally set permissible value
4. the volume of the space enveloping the final complexes does not exceed the internally set maximum permissible value
5. the parameter range in the most successful complex (minus 1 member) does not exceed the internally set maximum permissible value
6. the volume of the space enveloping the most successful complex (minus 1 member) does not exceed the internally set maximum permissible value
7. the parameter does not exceed both the absolute and the relative user-specified tolerance in the final complexes
8. the parameter does not exceed both the absolute and the relative user-specified tolerance in the most successful complex (minus 1 member)
9. the change of the objective function between consecutive shuffles does not exceed the user-specified tolerance
10. the RMSE of the best fit does not exceed the user-specified maximum


The optimization of run 1 converged after 724 evaluations of the objective function and 20 shuffles.
440 reflection points, 248 contraction points, and 36 random points were calculated.
The RMSE for run 1 is   1.615669E-01. The RMSE normalized by the maximum observed conductivity is   2.076588E-02.

A parameter did (T) or did not (F) achieve convergence on a particular criterion.

Parameter nr  Criterion 1  2  3  4  5  6  7  8  9 10
____________________________________________________
           1            T  F  T  T  T  T  T  T  T  F
           2            T  F  T  T  T  T  T  T  T  F
           3            T  F  T  T  T  T  T  T  T  F
           4            T  F  T  T  T  T  T  T  T  F
           5            T  F  T  T  T  T  T  T  T  F
           6            T  F  F  T  T  T  T  T  T  F
           7            T  F  T  T  T  T  T  T  T  F
           8            F  F  F  T  T  T  F  T  T  F
           9            F  F  F  T  F  T  F  F  T  F

Fitted parameter values of run 1.

Parameter name   Parameter value
________________________________
theta-s           4.28882800E-01
air-entry value  -2.02020833E+00
h oven-dryness   -6.30938045E+06
alpha             7.67000712E-04
n                 1.47972963E+00
Ks-cap            6.43138217E+00
Ks-ad             2.14158141E-03
gamma             6.90342028E+00
tau               1.45258838E-01

The correlations between parameters.
N.B. For those parameters whose logarithmic values were fitted, the correlations apply to the logarithms of their absolute values.

Nr of data points from which to draw the subset for estimating the correlation and, if desired, the covariance matrix: 724. Nr of data points in that subset: 500.

Correlation matrix (9.9990E+03 is an error value)

Variable nr               1            2            3            4            5            6            7            8            9
       1         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       2         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       3         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       4         9.9990E+03   9.9990E+03   9.9990E+03   1.0000E+00   9.9990E+03   7.5163E-01  -3.5374E-02  -8.0767E-01   5.7916E-01
       5         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       6         9.9990E+03   9.9990E+03   9.9990E+03   7.5163E-01   9.9990E+03   1.0000E+00  -6.8173E-02  -5.1861E-01   4.0059E-01
       7         9.9990E+03   9.9990E+03   9.9990E+03  -3.5374E-02   9.9990E+03  -6.8173E-02   1.0000E+00   1.2792E-01  -6.2777E-02
       8         9.9990E+03   9.9990E+03   9.9990E+03  -8.0767E-01   9.9990E+03  -5.1861E-01   1.2792E-01   1.0000E+00  -6.9423E-01
       9         9.9990E+03   9.9990E+03   9.9990E+03   5.7916E-01   9.9990E+03   4.0059E-01  -6.2777E-02  -6.9423E-01   1.0000E+00

The optimization of run 2 converged after 725 evaluations of the objective function and 21 shuffles.
462 reflection points, 235 contraction points, and 28 random points were calculated.
The RMSE for run 2 is   1.615654E-01. The RMSE normalized by the maximum observed conductivity is   2.076570E-02.

A parameter did (T) or did not (F) achieve convergence on a particular criterion.

Parameter nr  Criterion 1  2  3  4  5  6  7  8  9 10
____________________________________________________
           1            T  F  T  T  T  T  T  T  T  F
           2            T  F  T  T  T  T  T  T  T  F
           3            T  F  T  T  T  T  T  T  T  F
           4            T  F  T  T  T  T  T  T  T  F
           5            T  F  T  T  T  T  T  T  T  F
           6            T  F  F  T  T  T  T  T  T  F
           7            T  F  T  T  T  T  T  T  T  F
           8            F  F  T  T  T  T  T  T  T  F
           9            F  F  F  T  F  T  F  F  T  F

Fitted parameter values of run 2.

Parameter name   Parameter value
________________________________
theta-s           4.28882800E-01
air-entry value  -2.02020833E+00
h oven-dryness   -6.30938045E+06
alpha             7.67000712E-04
n                 1.47972963E+00
Ks-cap            6.44343778E+00
Ks-ad             2.13667821E-03
gamma             6.90773369E+00
tau               1.15513146E-01

The correlations between parameters.
N.B. For those parameters whose logarithmic values were fitted, the correlations apply to the logarithms of their absolute values.

Nr of data points from which to draw the subset for estimating the correlation and, if desired, the covariance matrix: 725. Nr of data points in that subset: 500.

Correlation matrix (9.9990E+03 is an error value)

Variable nr               1            2            3            4            5            6            7            8            9
       1         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       2         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       3         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       4         9.9990E+03   9.9990E+03   9.9990E+03   1.0000E+00   9.9990E+03   7.7502E-01  -2.3032E-01  -8.2092E-01   4.0163E-01
       5         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       6         9.9990E+03   9.9990E+03   9.9990E+03   7.7502E-01   9.9990E+03   1.0000E+00  -1.5951E-01  -4.9387E-01   2.4158E-01
       7         9.9990E+03   9.9990E+03   9.9990E+03  -2.3032E-01   9.9990E+03  -1.5951E-01   1.0000E+00   8.0726E-02  -6.7303E-02
       8         9.9990E+03   9.9990E+03   9.9990E+03  -8.2092E-01   9.9990E+03  -4.9387E-01   8.0726E-02   1.0000E+00  -5.2835E-01
       9         9.9990E+03   9.9990E+03   9.9990E+03   4.0163E-01   9.9990E+03   2.4158E-01  -6.7303E-02  -5.2835E-01   1.0000E+00

The optimization of run 3 converged after 984 evaluations of the objective function and 31 shuffles.
682 reflection points, 280 contraction points, and 22 random points were calculated.
The RMSE for run 3 is   1.615663E-01. The RMSE normalized by the maximum observed conductivity is   2.076580E-02.

A parameter did (T) or did not (F) achieve convergence on a particular criterion.

Parameter nr  Criterion 1  2  3  4  5  6  7  8  9 10
____________________________________________________
           1            T  F  T  T  T  T  T  T  T  F
           2            T  F  T  T  T  T  T  T  T  F
           3            T  F  T  T  T  T  T  T  T  F
           4            F  F  T  T  T  T  T  T  T  F
           5            T  F  T  T  T  T  T  T  T  F
           6            T  F  F  T  F  T  T  T  T  F
           7            T  F  T  T  T  T  T  T  T  F
           8            F  F  F  T  F  T  F  F  T  F
           9            F  F  F  T  F  T  F  F  T  F

Fitted parameter values of run 3.

Parameter name   Parameter value
________________________________
theta-s           4.28882800E-01
air-entry value  -2.02020833E+00
h oven-dryness   -6.30938045E+06
alpha             7.67000712E-04
n                 1.47972963E+00
Ks-cap            6.43315657E+00
Ks-ad             2.13691651E-03
gamma             6.91073062E+00
tau               5.52856147E-02

The correlations between parameters.
N.B. For those parameters whose logarithmic values were fitted, the correlations apply to the logarithms of their absolute values.

Nr of data points from which to draw the subset for estimating the correlation and, if desired, the covariance matrix: 984. Nr of data points in that subset: 500.

Correlation matrix (9.9990E+03 is an error value)

Variable nr               1            2            3            4            5            6            7            8            9
       1         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       2         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       3         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       4         9.9990E+03   9.9990E+03   9.9990E+03   1.0000E+00   9.9990E+03   5.6728E-01  -2.8991E-02  -5.1131E-01   6.6286E-01
       5         9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03   9.9990E+03
       6         9.9990E+03   9.9990E+03   9.9990E+03   5.6728E-01   9.9990E+03   1.0000E+00  -9.9022E-02  -2.1515E-01   4.4903E-01
       7         9.9990E+03   9.9990E+03   9.9990E+03  -2.8991E-02   9.9990E+03  -9.9022E-02   1.0000E+00   6.1523E-02   2.0768E-01
       8         9.9990E+03   9.9990E+03   9.9990E+03  -5.1131E-01   9.9990E+03  -2.1515E-01   6.1523E-02   1.0000E+00  -4.2345E-01
       9         9.9990E+03   9.9990E+03   9.9990E+03   6.6286E-01   9.9990E+03   4.4903E-01   2.0768E-01  -4.2345E-01   1.0000E+00

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

The best fit is that of run 2, with the RMSE =   1.615654E-01, and the AICc =  -5.236689E+01.

Fitted parameter values of run 2.

Parameter name   Parameter value
________________________________
theta-s           4.28882800E-01
air-entry value  -2.02020833E+00
h oven-dryness   -6.30938045E+06
alpha             7.67000712E-04
n                 1.47972963E+00
Ks-cap            6.44343778E+00
Ks-ad             2.13667821E-03
gamma             6.90773369E+00
tau               1.15513146E-01

tau            : The (non-logarithmic) fitted values differ by more than 10% of the mean of the smallest and largest value.
