Fit a non-population dynamic model

dynmodel(
  system,
  model,
  inits,
  data,
  fixPars = NULL,
  nlmixrObject = NULL,
  control = list(),
  ...
)

Arguments

system

RxODE object. See RxODE for more details.

model

Error model.

inits

Initial values of system parameters.

data

Dataset to estimate. Needs to be RxODE compatible in EVIDs.

fixPars

Fixed system parameters. Default is NULL.

nlmixrObject

nlmixr object. See nlmixr for more details. Default is NULL.

control

Control options for dynmodel dynmodelControl .

...

Other parameters (ignored)

Value

A dynmodel fit object

Author

Wenping Wang, Mason McComb and Matt Fidler

Examples


# \donttest{

# dynmodel example --------------------------------------------------------
ode <- "
      kel = CL/V;
      d/dt(X) = -kel*X;
      C=X/V;
      PRED = C
      "
ode_system <- RxODE(model = ode)
#>  
model_error_structure <- cp ~ C + add(0.01) + prop(0.01)
inits <- c(CL = 1, V = 10)
control <- dynmodelControl(method = "Nelder-Mead")
fit <-
  try(dynmodel(
    system = ode_system,
    model = model_error_structure,
    data = Bolus_1CPT,
    inits = inits,
    control = control
  ))
#> |-----+---------------+-----------+-----------+-----------+-----------|
#> |    #| Objective Fun |        CL |         V |       add |      prop |
#> |    1|     1487008.6 |     1.000 |     10.00 |   0.01000 |   0.01000 |
#> |    U|     1487008.6 |     1.000 |     10.00 |   0.01000 |   0.01000 |
#> |    2|     1487008.6 |    0.8000 |     10.00 |   0.01000 |   0.01000 |
#> |    U|     1487008.6 |    0.8000 |     10.00 |   0.01000 |   0.01000 |
#> |    3|     1487008.6 |     1.000 |     8.000 |   0.01000 |   0.01000 |
#> |    U|     1487008.6 |     1.000 |     8.000 |   0.01000 |   0.01000 |
#> |    4|     1493743.8 |     1.000 |     10.00 |  0.008000 |   0.01000 |
#> |    U|     1493743.8 |     1.000 |     10.00 |  0.008000 |   0.01000 |
#> |    5|     2302237.1 |     1.000 |     10.00 |   0.01000 |  0.008000 |
#> |    U|     2302237.1 |     1.000 |     10.00 |   0.01000 |  0.008000 |
#> |    6|     1043965.3 |    0.9000 |     9.000 |  0.009000 |   0.01200 |
#> |    U|     1043965.3 |    0.9000 |     9.000 |  0.009000 |   0.01200 |
#> |    7|     774746.27 |    0.8500 |     8.500 |  0.008500 |   0.01400 |
#> |    U|     774746.27 |    0.8500 |     8.500 |  0.008500 |   0.01400 |
#> |    8|     1039469.2 |    0.8250 |     8.250 |   0.01125 |   0.01200 |
#> |    U|     1039469.2 |    0.8250 |     8.250 |   0.01125 |   0.01200 |
#> |    9|     892071.79 |     1.037 |     7.375 |  0.009875 |   0.01300 |
#> |    U|     892071.79 |     1.037 |     7.375 |  0.009875 |   0.01300 |
#> |   10|     722274.73 |    0.8562 |     9.062 |  0.009813 |   0.01450 |
#> |    U|     722274.73 |    0.8562 |     9.062 |  0.009813 |   0.01450 |
#> |   11|     547640.04 |    0.7844 |     9.594 |  0.009719 |   0.01675 |
#> |    U|     547640.04 |    0.7844 |     9.594 |  0.009719 |   0.01675 |
#> |   12|     483911.36 |    0.7484 |     6.859 |  0.009672 |   0.01787 |
#> |    U|     483911.36 |    0.7484 |     6.859 |  0.009672 |   0.01787 |
#> |   13|     333073.34 |    0.6227 |     5.289 |  0.009508 |   0.02181 |
#> |    U|     333073.34 |    0.6227 |     5.289 |  0.009508 |   0.02181 |
#> |   14|     365491.93 |    0.8223 |     7.129 |  0.007551 |   0.02078 |
#> |    U|     365491.93 |    0.8223 |     7.129 |  0.007551 |   0.02078 |
#> |   15|     287159.92 |    0.5021 |     7.881 |  0.007764 |   0.02367 |
#> |    U|     287159.92 |    0.5021 |     7.881 |  0.007764 |   0.02367 |
#> |   16|     199785.52 |    0.2345 |     8.134 |  0.006708 |   0.02901 |
#> |    U|     199785.52 |    0.2345 |     8.134 |  0.006708 |   0.02901 |
#> |   17|     186206.44 |    0.3819 |     6.573 |  0.008243 |   0.03018 |
#> |    U|     186206.44 |    0.3819 |     6.573 |  0.008243 |   0.03018 |
#> |   18|     125691.57 |    0.1478 |     5.609 |  0.008114 |   0.03826 |
#> |    U|     125691.57 |    0.1478 |     5.609 |  0.008114 |   0.03826 |
#> |   19|     126339.01 |    0.1292 |     3.487 |  0.006222 |   0.03818 |
#> |    U|     126339.01 |    0.1292 |     3.487 |  0.006222 |   0.03818 |
#> |   20|     105729.15 |   -0.2552 |     4.130 |  0.007725 |   0.04285 |
#> |    U|     105729.15 |   -0.2552 |     4.130 |  0.007725 |   0.04285 |
#> |   21|     77192.666 |   -0.7939 |     2.631 |  0.007812 |   0.05389 |
#> |    U|     77192.666 |   -0.7939 |     2.631 |  0.007812 |   0.05389 |
#> |   22|     70930.129 |   -0.7638 |     4.641 |  0.004920 |   0.05786 |
#> |    U|     70930.129 |   -0.7638 |     4.641 |  0.004920 |   0.05786 |
#> |   23|     53932.881 |    -1.457 |     4.317 |  0.002626 |   0.07588 |
#> |    U|     53932.881 |    -1.457 |     4.317 |  0.002626 |   0.07588 |
#> |   24|     54975.876 |    -1.221 |   -0.1115 |  0.005679 |   0.07410 |
#> |    U|     54975.876 |    -1.221 |   -0.1115 |  0.005679 |   0.07410 |
#> |   25|     50213.136 |    -1.792 |     2.736 |  0.005894 |   0.08288 |
#> |    U|     50213.136 |    -1.792 |     2.736 |  0.005894 |   0.08288 |
#> |   26|     43493.772 |    -2.752 |     2.361 |  0.005730 |    0.1052 |
#> |    U|     43493.772 |    -2.752 |     2.361 |  0.005730 |    0.1052 |
#> |   27|     41723.672 |    -3.260 |    -1.010 |  0.002809 |    0.1163 |
#> |    U|     41723.672 |    -3.260 |    -1.010 |  0.002809 |    0.1163 |
#> |   28|     38792.534 |    -4.964 |    -4.319 | 0.0001571 |    0.1553 |
#> |    U|     38792.534 |    -4.964 |    -4.319 | 0.0001571 |    0.1553 |
#> |   29|     38945.709 |    -4.403 |    -1.507 |-0.0007160 |    0.1514 |
#> |    U|     38945.709 |    -4.403 |    -1.507 |-0.0007160 |    0.1514 |
#> |   30|     38374.346 |    -5.567 |    0.5376 | -0.001780 |    0.1698 |
#> |    U|     38374.346 |    -5.567 |    0.5376 | -0.001780 |    0.1698 |
#> |   31|     38027.291 |    -7.739 |    0.8621 | -0.005510 |    0.2176 |
#> |    U|     38027.291 |    -7.739 |    0.8621 | -0.005510 |    0.2176 |
#> |   32|     38128.019 |    -8.472 |    -5.619 | -0.002795 |    0.2389 |
#> |    U|     38128.019 |    -8.472 |    -5.619 | -0.002795 |    0.2389 |
#> |   33|     38466.546 |    -10.04 |    -7.653 |  -0.01016 |    0.2764 |
#> |    U|     38466.546 |    -10.04 |    -7.653 |  -0.01016 |    0.2764 |
#> |   34|     38649.017 |    -11.20 |    -6.857 | -0.008439 |    0.2927 |
#> |    U|     38649.017 |    -11.20 |    -6.857 | -0.008439 |    0.2927 |
#> |   35|     39454.074 |    -13.76 |    -5.314 |  -0.01361 |    0.3575 |
#> |    U|     39454.074 |    -13.76 |    -5.314 |  -0.01361 |    0.3575 |
#> |   36|     38020.121 |    -7.163 |    -4.568 | -0.003285 |    0.2059 |
#> |    U|     38020.121 |    -7.163 |    -4.568 | -0.003285 |    0.2059 |
#> |   37|     38245.755 |    -5.503 |    -1.632 | -0.002437 |    0.1766 |
#> |    U|     38245.755 |    -5.503 |    -1.632 | -0.002437 |    0.1766 |
#> |   38|     39336.269 |    -4.402 |     2.175 |  0.003149 |    0.1431 |
#> |    U|     39336.269 |    -4.402 |     2.175 |  0.003149 |    0.1431 |
#> |   39|     38157.426 |    -8.628 |    -5.196 | -0.006834 |    0.2431 |
#> |    U|     38157.426 |    -8.628 |    -5.196 | -0.006834 |    0.2431 |
#> |   40|     38462.561 |    -10.50 |    -5.629 | -0.006775 |    0.2761 |
#> |    U|     38462.561 |    -10.50 |    -5.629 | -0.006775 |    0.2761 |
#> |   41|     38028.571 |    -6.752 |    -2.631 | -0.003521 |    0.2015 |
#> |    U|     38028.571 |    -6.752 |    -2.631 | -0.003521 |    0.2015 |
#> |   42|     38100.333 |    -6.435 |   -0.7819 |-0.0007213 |    0.1889 |
#> |    U|     38100.333 |    -6.435 |   -0.7819 |-0.0007213 |    0.1889 |
#> |   43|     38409.707 |    -5.573 |     2.060 | -0.003723 |    0.1680 |
#> |    U|     38409.707 |    -5.573 |     2.060 | -0.003723 |    0.1680 |
#> |   44|     38037.842 |    -7.747 |    -3.699 | -0.003027 |    0.2212 |
#> |    U|     38037.842 |    -7.747 |    -3.699 | -0.003027 |    0.2212 |
#> |   45|     38098.156 |    -8.266 |    -4.236 | -0.006950 |    0.2342 |
#> |    U|     38098.156 |    -8.266 |    -4.236 | -0.006950 |    0.2342 |
#> |   46|     38043.181 |    -7.808 |    -3.373 | -0.005393 |    0.2229 |
#> |    U|     38043.181 |    -7.808 |    -3.373 | -0.005393 |    0.2229 |
#> |   47|     38033.280 |    -6.893 |    -1.645 | -0.002279 |    0.2002 |
#> |    U|     38033.280 |    -6.893 |    -1.645 | -0.002279 |    0.2002 |
#> |   48|     38077.121 |    -6.526 |   -0.2918 | -0.004270 |    0.1914 |
#> |    U|     38077.121 |    -6.526 |   -0.2918 | -0.004270 |    0.1914 |
#> |   49|     38020.744 |    -7.442 |    -2.848 | -0.003338 |    0.2137 |
#> |    U|     38020.744 |    -7.442 |    -2.848 | -0.003338 |    0.2137 |
#> |   50|     38031.184 |    -7.656 |    -2.947 | -0.005548 |    0.2192 |
#> |    U|     38031.184 |    -7.656 |    -2.947 | -0.005548 |    0.2192 |
#> |   51|     38021.231 |    -7.465 |    -2.621 | -0.004731 |    0.2144 |
#> |    U|     38021.231 |    -7.465 |    -2.621 | -0.004731 |    0.2144 |
#> |   52|     38048.786 |    -8.153 |    -1.957 | -0.004910 |    0.2243 |
#> |    U|     38048.786 |    -8.153 |    -1.957 | -0.004910 |    0.2243 |
#> |   53|     38018.540 |    -7.102 |    -2.462 | -0.003869 |    0.2072 |
#> |    U|     38018.540 |    -7.102 |    -2.462 | -0.003869 |    0.2072 |
#> |   54|     38025.751 |    -6.847 |    -7.112 | -0.002101 |    0.2030 |
#> |    U|     38025.751 |    -6.847 |    -7.112 | -0.002101 |    0.2030 |
#> |   55|     38019.452 |    -7.070 |    -5.118 | -0.002953 |    0.2066 |
#> |    U|     38019.452 |    -7.070 |    -5.118 | -0.002953 |    0.2066 |
#> |   56|     38027.420 |    -6.924 |    -4.877 | -0.001991 |    0.2023 |
#> |    U|     38027.420 |    -6.924 |    -4.877 | -0.001991 |    0.2023 |
#> |   57|     38018.257 |    -7.330 |    -3.185 | -0.004046 |    0.2114 |
#> |    U|     38018.257 |    -7.330 |    -3.185 | -0.004046 |    0.2114 |
#> |   58|     38027.633 |    -6.891 |    -4.820 | -0.003738 |    0.2018 |
#> |    U|     38027.633 |    -6.891 |    -4.820 | -0.003738 |    0.2018 |
#> |   59|     38018.175 |    -7.304 |    -3.341 | -0.003438 |    0.2108 |
#> |    U|     38018.175 |    -7.304 |    -3.341 | -0.003438 |    0.2108 |
#> |   60|     38018.862 |    -7.240 |    -2.485 | -0.003868 |    0.2121 |
#> |    U|     38018.862 |    -7.240 |    -2.485 | -0.003868 |    0.2121 |
#> |   61|     38020.804 |    -7.418 |   -0.6183 | -0.004657 |    0.2141 |
#> |    U|     38020.804 |    -7.418 |   -0.6183 | -0.004657 |    0.2141 |
#> |   62|     38018.089 |    -7.157 |    -3.993 | -0.003379 |    0.2085 |
#> |    U|     38018.089 |    -7.157 |    -3.993 | -0.003379 |    0.2085 |
#> |   63|     38019.049 |    -7.207 |    -4.006 | -0.003498 |    0.2068 |
#> |    U|     38019.049 |    -7.207 |    -4.006 | -0.003498 |    0.2068 |
#> |   64|     38018.059 |    -7.231 |    -2.865 | -0.003776 |    0.2108 |
#> |    U|     38018.059 |    -7.231 |    -2.865 | -0.003776 |    0.2108 |
#> |   65|     38020.426 |    -7.409 |    -4.230 | -0.003451 |    0.2135 |
#> |    U|     38020.426 |    -7.409 |    -4.230 | -0.003451 |    0.2135 |
#> |   66|     38017.847 |    -7.179 |    -2.904 | -0.003764 |    0.2088 |
#> |    U|     38017.847 |    -7.179 |    -2.904 | -0.003764 |    0.2088 |
#> |   67|     38018.385 |    -7.106 |    -3.366 | -0.003133 |    0.2080 |
#> |    U|     38018.385 |    -7.106 |    -3.366 | -0.003133 |    0.2080 |
#> |   68|     38017.946 |    -7.274 |    -3.231 | -0.003818 |    0.2106 |
#> |    U|     38017.946 |    -7.274 |    -3.231 | -0.003818 |    0.2106 |
#> |   69|     38017.830 |    -7.116 |    -3.156 | -0.003930 |    0.2086 |
#> |    U|     38017.830 |    -7.116 |    -3.156 | -0.003930 |    0.2086 |
#> |   70|     38018.225 |    -7.022 |    -3.064 | -0.004177 |    0.2075 |
#> |    U|     38018.225 |    -7.022 |    -3.064 | -0.004177 |    0.2075 |
#> |   71|     38017.893 |    -7.243 |    -2.085 | -0.004265 |    0.2108 |
#> |    U|     38017.893 |    -7.243 |    -2.085 | -0.004265 |    0.2108 |
#> |   72|     38017.753 |    -7.175 |    -2.823 | -0.004113 |    0.2086 |
#> |    U|     38017.753 |    -7.175 |    -2.823 | -0.004113 |    0.2086 |
#> |   73|     38018.191 |    -7.146 |    -2.801 | -0.004281 |    0.2075 |
#> |    U|     38018.191 |    -7.146 |    -2.801 | -0.004281 |    0.2075 |
#> |   74|     38017.999 |    -7.083 |    -2.253 | -0.004218 |    0.2078 |
#> |    U|     38017.999 |    -7.083 |    -2.253 | -0.004218 |    0.2078 |
#> |   75|     38017.743 |    -7.226 |    -2.986 | -0.003918 |    0.2099 |
#> |    U|     38017.743 |    -7.226 |    -2.986 | -0.003918 |    0.2099 |
#> |   76|     38018.782 |    -7.105 |    -3.850 | -0.003598 |    0.2071 |
#> |    U|     38018.782 |    -7.105 |    -3.850 | -0.003598 |    0.2071 |
#> |   77|     38017.676 |    -7.209 |    -2.526 | -0.004098 |    0.2099 |
#> |    U|     38017.676 |    -7.209 |    -2.526 | -0.004098 |    0.2099 |
#> |   78|     38017.585 |    -7.184 |    -2.841 | -0.004265 |    0.2097 |
#> |    U|     38017.585 |    -7.184 |    -2.841 | -0.004265 |    0.2097 |
#> |   79|     38017.557 |    -7.186 |    -2.810 | -0.004516 |    0.2101 |
#> |    U|     38017.557 |    -7.186 |    -2.810 | -0.004516 |    0.2101 |
#> |   80|     38017.769 |    -7.281 |    -2.416 | -0.004392 |    0.2107 |
#> |    U|     38017.769 |    -7.281 |    -2.416 | -0.004392 |    0.2107 |
#> |   81|     38017.651 |    -7.240 |    -2.601 | -0.004276 |    0.2101 |
#> |    U|     38017.651 |    -7.240 |    -2.601 | -0.004276 |    0.2101 |
#> |   82|     38018.198 |    -7.256 |    -2.639 | -0.004291 |    0.2114 |
#> |    U|     38018.198 |    -7.256 |    -2.639 | -0.004291 |    0.2114 |
#> |   83|     38017.621 |    -7.195 |    -2.777 | -0.004157 |    0.2093 |
#> |    U|     38017.621 |    -7.195 |    -2.777 | -0.004157 |    0.2093 |
#> |   84|     38017.477 |    -7.189 |    -2.371 | -0.004606 |    0.2098 |
#> |    U|     38017.477 |    -7.189 |    -2.371 | -0.004606 |    0.2098 |
#> |   85|     38017.353 |    -7.171 |    -2.063 | -0.004950 |    0.2098 |
#> |    U|     38017.353 |    -7.171 |    -2.063 | -0.004950 |    0.2098 |
#> |   86|     38017.381 |    -7.188 |    -2.599 | -0.004852 |    0.2098 |
#> |    U|     38017.381 |    -7.188 |    -2.599 | -0.004852 |    0.2098 |
#> |   87|     38017.318 |    -7.130 |    -2.523 | -0.004961 |    0.2094 |
#> |    U|     38017.318 |    -7.130 |    -2.523 | -0.004961 |    0.2094 |
#> |   88|     38017.223 |    -7.074 |    -2.484 | -0.005304 |    0.2090 |
#> |    U|     38017.223 |    -7.074 |    -2.484 | -0.005304 |    0.2090 |
#> |   89|     38017.165 |    -7.114 |    -2.201 | -0.005653 |    0.2101 |
#> |    U|     38017.165 |    -7.114 |    -2.201 | -0.005653 |    0.2101 |
#> |   90|     38017.060 |    -7.074 |    -1.914 | -0.006401 |    0.2105 |
#> |    U|     38017.060 |    -7.074 |    -1.914 | -0.006401 |    0.2105 |
#> |   91|     38016.905 |    -7.067 |    -1.720 | -0.006238 |    0.2094 |
#> |    U|     38016.905 |    -7.067 |    -1.720 | -0.006238 |    0.2094 |
#> |   92|     38016.670 |    -7.007 |    -1.175 | -0.007098 |    0.2091 |
#> |    U|     38016.670 |    -7.007 |    -1.175 | -0.007098 |    0.2091 |
#> |   93|     38016.684 |    -6.976 |    -1.219 | -0.007025 |    0.2094 |
#> |    U|     38016.684 |    -6.976 |    -1.219 | -0.007025 |    0.2094 |
#> |   94|     38016.452 |    -6.895 |    -1.333 | -0.007964 |    0.2091 |
#> |    U|     38016.452 |    -6.895 |    -1.333 | -0.007964 |    0.2091 |
#> |   95|     38016.166 |    -6.757 |   -0.9678 | -0.009471 |    0.2088 |
#> |    U|     38016.166 |    -6.757 |   -0.9678 | -0.009471 |    0.2088 |
#> |   96|     38016.174 |    -6.833 |   -0.1534 | -0.009694 |    0.2099 |
#> |    U|     38016.174 |    -6.833 |   -0.1534 | -0.009694 |    0.2099 |
#> |   97|     38016.187 |    -6.712 |    0.1561 |  -0.01024 |    0.2081 |
#> |    U|     38016.187 |    -6.712 |    0.1561 |  -0.01024 |    0.2081 |
#> |   98|     38015.922 |    -6.679 |    0.1485 |  -0.01123 |    0.2085 |
#> |    U|     38015.922 |    -6.679 |    0.1485 |  -0.01123 |    0.2085 |
#> |   99|     38015.788 |    -6.531 |    0.8319 |  -0.01333 |    0.2081 |
#> |    U|     38015.788 |    -6.531 |    0.8319 |  -0.01333 |    0.2081 |
#> |  100|     38015.643 |    -6.410 |     1.109 |  -0.01427 |    0.2083 |
#> |    U|     38015.643 |    -6.410 |     1.109 |  -0.01427 |    0.2083 |
#> |  101|     38015.601 |    -6.111 |     2.250 |  -0.01786 |    0.2080 |
#> |    U|     38015.601 |    -6.111 |     2.250 |  -0.01786 |    0.2080 |
#> |  102|     38015.541 |    -6.404 |    0.8244 |  -0.01493 |    0.2093 |
#> |    U|     38015.541 |    -6.404 |    0.8244 |  -0.01493 |    0.2093 |
#> |  103|     38015.587 |    -6.249 |     1.159 |  -0.01728 |    0.2098 |
#> |    U|     38015.587 |    -6.249 |     1.159 |  -0.01728 |    0.2098 |
#> |  104|     38015.943 |    -6.069 |     1.623 |  -0.01810 |    0.2072 |
#> |    U|     38015.943 |    -6.069 |     1.623 |  -0.01810 |    0.2072 |
#> |  105|     38015.897 |    -5.800 |     3.733 |  -0.02264 |    0.2075 |
#> |    U|     38015.897 |    -5.800 |     3.733 |  -0.02264 |    0.2075 |
#> |  106|     38015.489 |    -6.354 |     2.197 |  -0.01628 |    0.2092 |
#> |    U|     38015.489 |    -6.354 |     2.197 |  -0.01628 |    0.2092 |
#> |  107|     38015.750 |    -6.497 |     2.484 |  -0.01537 |    0.2102 |
#> |    U|     38015.750 |    -6.497 |     2.484 |  -0.01537 |    0.2102 |
#> |  108|     38016.370 |    -6.900 |   -0.6808 | -0.008559 |    0.2098 |
#> |    U|     38016.370 |    -6.900 |   -0.6808 | -0.008559 |    0.2098 |
#> |  109|     38015.581 |    -6.075 |     2.629 |  -0.01912 |    0.2080 |
#> |    U|     38015.581 |    -6.075 |     2.629 |  -0.01912 |    0.2080 |
#> |  110|     38015.533 |    -5.941 |     3.118 |  -0.02076 |    0.2092 |
#> |    U|     38015.533 |    -5.941 |     3.118 |  -0.02076 |    0.2092 |
#> |  111|     38015.597 |    -6.276 |     2.134 |  -0.01769 |    0.2099 |
#> |    U|     38015.597 |    -6.276 |     2.134 |  -0.01769 |    0.2099 |
#> |  112|     38015.493 |    -6.235 |     2.163 |  -0.01773 |    0.2094 |
#> |    U|     38015.493 |    -6.235 |     2.163 |  -0.01773 |    0.2094 |
#> |  113|     38015.838 |    -6.392 |     1.522 |  -0.01573 |    0.2105 |
#> |    U|     38015.838 |    -6.392 |     1.522 |  -0.01573 |    0.2105 |
#> |  114|     38015.470 |    -6.154 |     2.352 |  -0.01827 |    0.2086 |
#> |    U|     38015.470 |    -6.154 |     2.352 |  -0.01827 |    0.2086 |
#> |  115|     38015.557 |    -5.938 |     4.091 |  -0.02159 |    0.2089 |
#> |    U|     38015.557 |    -5.938 |     4.091 |  -0.02159 |    0.2089 |
#> |  116|     38015.480 |    -6.287 |     1.641 |  -0.01660 |    0.2092 |
#> |    U|     38015.480 |    -6.287 |     1.641 |  -0.01660 |    0.2092 |
#> |  117|     38015.610 |    -6.574 |     1.058 |  -0.01368 |    0.2091 |
#> |    U|     38015.610 |    -6.574 |     1.058 |  -0.01368 |    0.2091 |
#> |  118|     38015.477 |    -6.099 |     2.603 |  -0.01899 |    0.2091 |
#> |    U|     38015.477 |    -6.099 |     2.603 |  -0.01899 |    0.2091 |
#> |  119|     38015.477 |    -6.099 |     2.603 |  -0.01899 |    0.2091 |
#> |    U|     38015.477 |    -6.099 |     2.603 |  -0.01899 |    0.2091 |
#> |  120|     38015.470 |    -6.154 |     2.352 |  -0.01827 |    0.2086 |
#> |    U|     38015.470 |    -6.154 |     2.352 |  -0.01827 |    0.2086 |
#> |  121|     38015.493 |    -6.235 |     2.163 |  -0.01773 |    0.2094 |
#> |    U|     38015.493 |    -6.235 |     2.163 |  -0.01773 |    0.2094 |
#> |  122|     38015.480 |    -6.287 |     1.641 |  -0.01660 |    0.2092 |
#> |    U|     38015.480 |    -6.287 |     1.641 |  -0.01660 |    0.2092 |
#> |  123|     38015.489 |    -6.354 |     2.197 |  -0.01628 |    0.2092 |
#> |    U|     38015.489 |    -6.354 |     2.197 |  -0.01628 |    0.2092 |
#> 
#> Error in solve.default(fit$hessian) : 
#>   Lapack routine dgesv: system is exactly singular: U[1,1] = 0

# nlmixr model example ----------------------------------------------------------
model_onecmt_bolus <- function() {
  ini({
    CL <- c(0, 5, 10) # Clearance (L/hr)
    V <- c(0, 50, 100) # Volume of Distribution
    prop.err <- c(0, 0.01, 1)
  })
  model({
    kel <- CL / V
    d / dt(X) <- -kel * X
    cp <- X / V
    cp ~ prop(prop.err)
  })
}

# note on some platforms this fit is not successful
fit <- try(nlmixr(object = model_onecmt_bolus, data = Bolus_1CPT, est = "dynmodel"))
#>  parameter labels from comments will be replaced by 'label()'
#> → creating full model...
#> → pruning branches (`if`/`else`)...
#>  done
#> → loading into symengine environment...
#>  done
#> → finding duplicate expressions in EBE model...
#> 
#> → optimizing duplicate expressions in EBE model...
#> 
#> → compiling EBE model...
#>  
#>  done
#>  
#>  
#> |-----+---------------+-----------+-----------+-----------|
#> |    #| Objective Fun |        CL |         V |      prop |
#> |    1| 4.0228158e+09 |   -0.8004 |     1.000 |    -1.000 |
#> |    U| 4.0228158e+09 |     5.000 |     50.00 |   0.01000 |
#> |    2|     10088613. |   -0.8004 |     1.000 |   -0.8100 |
#> |    U|     10088613. |     5.000 |     50.00 |    0.2000 |
#> |    3|     11311242. |   -0.6004 |     1.000 |   -0.8100 |
#> |    U|     11311242. |     5.040 |     50.00 |    0.2000 |
#> |    4|     10078717. |   -0.8004 |     1.200 |   -0.8100 |
#> |    U|     10078717. |     5.000 |     50.00 |    0.2000 |
#> |    5|     2550672.7 |   -0.8004 |     1.000 |   -0.6100 |
#> |    U|     2550672.7 |     5.000 |     50.00 |    0.4000 |
#> |    6|     8998173.4 |    -1.000 |     1.000 |   -0.8100 |
#> |    U|     8998173.4 |     4.960 |     50.00 |    0.2000 |
#> |    7|     10098520. |   -0.8004 |    0.8000 |   -0.8100 |
#> |    U|     10098520. |     5.000 |     50.00 |    0.2000 |
#> |    8| 1.8117103e+21 |   -0.8004 |     1.000 |    -1.010 |
#> |    U| 1.8117103e+21 |     5.000 |     50.00 | 1.490e-08 |
#> |    9|     4504202.9 |   -0.8004 |     1.000 |   -0.7100 |
#> |    U|     4504202.9 |     5.000 |     50.00 |    0.3000 |
#> |   10|     2114977.2 |   -0.8004 |     1.000 |   -0.5700 |
#> |    U|     2114977.2 |     5.000 |     50.00 |    0.4400 |
#> |   11|     2505577.9 |   -0.8306 |     1.000 |   -0.6098 |
#> |    U|     2505577.9 |     4.994 |     50.00 |    0.4002 |
#> |   12|     2098337.6 |   -0.7808 |    0.9993 |   -0.5657 |
#> |    U|     2098337.6 |     5.004 |     50.00 |    0.4443 |
#> |   13|     2256365.1 |   -0.7789 |     1.011 |   -0.5816 |
#> |    U|     2256365.1 |     5.004 |     50.00 |    0.4284 |
#> |   14|     1963459.6 |   -0.7796 |     1.012 |   -0.5503 |
#> |    U|     1963459.6 |     5.004 |     50.00 |    0.4597 |
#> |   15|     1824149.0 |   -0.7776 |     1.003 |   -0.5324 |
#> |    U|     1824149.0 |     5.005 |     50.00 |    0.4776 |
#> |   16|     1686944.4 |   -0.7752 |     1.004 |   -0.5125 |
#> |    U|     1686944.4 |     5.005 |     50.00 |    0.4975 |
#> |   17|     1604269.6 |   -0.7916 |     1.008 |   -0.5019 |
#> |    U|     1604269.6 |     5.002 |     50.00 |    0.5081 |
#> |   18|     1498903.6 |   -0.7904 |     1.000 |   -0.4837 |
#> |    U|     1498903.6 |     5.002 |     50.00 |    0.5263 |
#> |   19|     1418723.5 |   -0.7870 |     1.012 |   -0.4681 |
#> |    U|     1418723.5 |     5.003 |     50.00 |    0.5419 |
#> |   20|     1328386.5 |   -0.7958 |     1.017 |   -0.4508 |
#> |    U|     1328386.5 |     5.001 |     50.00 |    0.5592 |
#> |   21|     1245001.5 |   -0.7950 |     1.022 |   -0.4316 |
#> |    U|     1245001.5 |     5.001 |     50.00 |    0.5784 |
#> |   22|     1095740.8 |   -0.8003 |     1.031 |   -0.3929 |
#> |    U|     1095740.8 |     5.000 |     50.00 |    0.6171 |
#> |   23|     1071789.7 |   -0.7861 |     1.067 |   -0.3833 |
#> |    U|     1071789.7 |     5.003 |     50.00 |    0.6267 |
#> |   24|     1065384.0 |   -0.8038 |     1.077 |   -0.3845 |
#> |    U|     1065384.0 |     4.999 |     50.00 |    0.6255 |
#> |   25|     1006131.6 |   -0.8035 |     1.070 |   -0.3655 |
#> |    U|     1006131.6 |     4.999 |     50.00 |    0.6445 |
#> |   26|     953406.83 |   -0.8062 |     1.062 |   -0.3476 |
#> |    U|     953406.83 |     4.999 |     50.00 |    0.6624 |
#> |   27|     900141.19 |   -0.8106 |     1.064 |   -0.3282 |
#> |    U|     900141.19 |     4.998 |     50.00 |    0.6818 |
#> |   28|     855151.42 |   -0.8218 |     1.061 |   -0.3118 |
#> |    U|     855151.42 |     4.996 |     50.00 |    0.6982 |
#> |   29|     839879.22 |   -0.8263 |     1.080 |   -0.3061 |
#> |    U|     839879.22 |     4.995 |     50.00 |    0.7039 |
#> |   30|     879872.78 |   -0.8333 |     1.084 |   -0.3245 |
#> |    U|     879872.78 |     4.993 |     50.00 |    0.6855 |
#> |   31|     809362.43 |   -0.8149 |     1.078 |   -0.2898 |
#> |    U|     809362.43 |     4.997 |     50.00 |    0.7202 |
#> |   32|     768838.21 |   -0.8171 |     1.082 |   -0.2704 |
#> |    U|     768838.21 |     4.997 |     50.00 |    0.7396 |
#> |   33|     754779.11 |   -0.8151 |     1.101 |   -0.2627 |
#> |    U|     754779.11 |     4.997 |     50.00 |    0.7473 |
#> |   34|     742385.17 |   -0.8167 |     1.099 |   -0.2565 |
#> |    U|     742385.17 |     4.997 |     50.00 |    0.7535 |
#> |   35|     722947.92 |   -0.8212 |     1.107 |   -0.2468 |
#> |    U|     722947.92 |     4.996 |     50.00 |    0.7632 |
#> |   36|     691210.45 |   -0.8308 |     1.126 |   -0.2305 |
#> |    U|     691210.45 |     4.994 |     50.00 |    0.7795 |
#> |   37|     653200.62 |   -0.8414 |     1.151 |   -0.2090 |
#> |    U|     653200.62 |     4.992 |     50.00 |    0.8010 |
#> |   38|     618137.03 |   -0.8538 |     1.179 |   -0.1879 |
#> |    U|     618137.03 |     4.989 |     50.00 |    0.8221 |
#> |   39|     563861.59 |   -0.8814 |     1.238 |   -0.1529 |
#> |    U|     563861.59 |     4.984 |     50.00 |    0.8571 |
#> |   40|     543896.98 |   -0.8916 |     1.260 |   -0.1385 |
#> |    U|     543896.98 |     4.982 |     50.01 |    0.8715 |
#> |   41|     491772.79 |   -0.9124 |     1.316 |  -0.09466 |
#> |    U|     491772.79 |     4.978 |     50.01 |    0.9153 |
#> |   42|     409686.72 |   -0.9575 |     1.429 |  -0.01000 |
#> |    U|     409686.72 |     4.969 |     50.01 |     1.000 |
#> |   43|     404075.37 |   -0.9840 |     1.471 |  -0.01000 |
#> |    U|     404075.37 |     4.963 |     50.01 |     1.000 |
#> |   44|     406175.79 |   -0.9741 |     1.455 |  -0.01000 |
#> |    U|     406175.79 |     4.965 |     50.01 |     1.000 |
#> |   45|     412238.46 |   -0.9701 |     1.480 |  -0.01724 |
#> |    U|     412238.46 |     4.966 |     50.01 |    0.9928 |
#> |   46|     407188.09 |   -0.9855 |     1.488 |  -0.01480 |
#> |    U|     407188.09 |     4.963 |     50.01 |    0.9952 |
#> |   47|     409356.74 |   -0.9873 |     1.468 |  -0.01824 |
#> |    U|     409356.74 |     4.963 |     50.01 |    0.9918 |
#> |   48|     402190.92 |   -0.9932 |     1.471 |  -0.01000 |
#> |    U|     402190.92 |     4.961 |     50.01 |     1.000 |
#> |   49|     398451.30 |    -1.011 |     1.471 |  -0.01000 |
#> |    U|     398451.30 |     4.958 |     50.01 |     1.000 |
#> |   50|     391088.50 |    -1.048 |     1.471 |  -0.01000 |
#> |    U|     391088.50 |     4.950 |     50.01 |     1.000 |
#> |   51|     376914.13 |    -1.120 |     1.480 |  -0.01000 |
#> |    U|     376914.13 |     4.936 |     50.01 |     1.000 |
#> |   52|     350997.04 |    -1.261 |     1.520 |  -0.01000 |
#> |    U|     350997.04 |     4.908 |     50.01 |     1.000 |
#> |   53|     303854.10 |    -1.551 |     1.558 |  -0.01000 |
#> |    U|     303854.10 |     4.850 |     50.01 |     1.000 |
#> |   54|     235434.59 |    -2.085 |     1.319 |  -0.01000 |
#> |    U|     235434.59 |     4.743 |     50.01 |     1.000 |
#> |   55|     181463.51 |    -2.669 |    0.3060 |  -0.01000 |
#> |    U|     181463.51 |     4.626 |     49.99 |     1.000 |
#> |   56|     140110.01 |    -3.299 |   -0.6788 |  -0.01000 |
#> |    U|     140110.01 |     4.500 |     49.97 |     1.000 |
#> |   57|     146973.88 |    -3.197 |    -3.015 |  -0.01000 |
#> |    U|     146973.88 |     4.521 |     49.92 |     1.000 |
#> |   58|     248822.36 |    -3.273 |   -0.6957 |   -0.3241 |
#> |    U|     248822.36 |     4.505 |     49.97 |    0.6859 |
#> |   59|     127628.14 |    -3.550 |    -1.821 |  -0.01000 |
#> |    U|     127628.14 |     4.450 |     49.94 |     1.000 |
#> |   60|     109368.45 |    -3.988 |    -2.905 |  -0.01000 |
#> |    U|     109368.45 |     4.362 |     49.92 |     1.000 |
#> |   61|     106026.09 |    -4.089 |    -4.070 |  -0.01000 |
#> |    U|     106026.09 |     4.342 |     49.90 |     1.000 |
#> |   62|     114080.69 |    -3.884 |    -5.221 |  -0.01000 |
#> |    U|     114080.69 |     4.383 |     49.88 |     1.000 |
#> |   63|     88485.074 |    -4.673 |    -4.070 |  -0.01000 |
#> |    U|     88485.074 |     4.225 |     49.90 |     1.000 |
#> |   64|     76050.713 |    -5.258 |    -4.067 |  -0.01000 |
#> |    U|     76050.713 |     4.108 |     49.90 |     1.000 |
#> |   65|     65082.367 |    -6.022 |    -3.181 |  -0.01000 |
#> |    U|     65082.367 |     3.956 |     49.92 |     1.000 |
#> |   66|     61057.453 |    -6.412 |   -0.8756 |  -0.01000 |
#> |    U|     61057.453 |     3.878 |     49.96 |     1.000 |
#> |   67|     205362.83 |    -6.828 |     1.306 |   -0.7418 |
#> |    U|     205362.83 |     3.795 |     50.01 |    0.2682 |
#> |   68|     99165.127 |    -6.465 |   -0.9734 |   -0.5003 |
#> |    U|     99165.127 |     3.867 |     49.96 |    0.5097 |
#> |   69|     55514.873 |    -7.228 |    -1.714 |  -0.01000 |
#> |    U|     55514.873 |     3.714 |     49.95 |     1.000 |
#> |   70|     55138.370 |    -7.334 |    -1.304 |  -0.02500 |
#> |    U|     55138.370 |     3.693 |     49.95 |    0.9850 |
#> |   71|     58127.346 |    -6.793 |    -1.524 |  -0.01000 |
#> |    U|     58127.346 |     3.801 |     49.95 |     1.000 |
#> |   72|     62383.088 |    -7.216 |    -1.278 |   -0.2910 |
#> |    U|     62383.088 |     3.717 |     49.95 |    0.7190 |
#> |   73|     54994.481 |    -7.329 |    -1.350 |  -0.01000 |
#> |    U|     54994.481 |     3.694 |     49.95 |     1.000 |
#> |   74|     55429.616 |    -7.243 |    -1.466 |  -0.01000 |
#> |    U|     55429.616 |     3.712 |     49.95 |     1.000 |
#> |   75|     55075.824 |    -7.381 |    -1.393 |  -0.03880 |
#> |    U|     55075.824 |     3.684 |     49.95 |    0.9712 |
#> |   76|     54649.162 |    -7.401 |    -1.340 |  -0.01000 |
#> |    U|     54649.162 |     3.680 |     49.95 |     1.000 |
#> |   77|     54319.458 |    -7.473 |    -1.349 |  -0.01000 |
#> |    U|     54319.458 |     3.665 |     49.95 |     1.000 |
#> |   78|     53751.144 |    -7.607 |    -1.289 |  -0.01000 |
#> |    U|     53751.144 |     3.639 |     49.95 |     1.000 |
#> |   79|     52859.970 |    -7.843 |    -1.117 |  -0.01000 |
#> |    U|     52859.970 |     3.591 |     49.96 |     1.000 |
#> |   80|     51385.727 |    -8.349 |   -0.8242 |  -0.01000 |
#> |    U|     51385.727 |     3.490 |     49.96 |     1.000 |
#> |   81|     94108.838 |    -9.212 |   -0.6474 |   -0.7793 |
#> |    U|     94108.838 |     3.318 |     49.97 |    0.2307 |
#> |   82|     51161.841 |    -8.447 |   -0.9933 |  -0.01000 |
#> |    U|     51161.841 |     3.471 |     49.96 |     1.000 |
#> |   83|     50277.644 |    -8.918 |   -0.6465 |  -0.01000 |
#> |    U|     50277.644 |     3.377 |     49.97 |     1.000 |
#> |   84|     49747.822 |    -9.327 |   -0.4212 |  -0.01000 |
#> |    U|     49747.822 |     3.295 |     49.97 |     1.000 |
#> |   85|     49478.559 |    -9.617 |   -0.2578 |  -0.01000 |
#> |    U|     49478.559 |     3.237 |     49.97 |     1.000 |
#> |   86|     49172.841 |    -10.14 |    0.1596 |  -0.01000 |
#> |    U|     49172.841 |     3.133 |     49.98 |     1.000 |
#> |   87|     49151.295 |    -10.19 |    0.1503 |  -0.01000 |
#> |    U|     49151.295 |     3.121 |     49.98 |     1.000 |
#> |   88|     49170.455 |    -10.14 |    0.1282 |  -0.01000 |
#> |    U|     49170.455 |     3.131 |     49.98 |     1.000 |
#> |   89|     48974.917 |    -10.20 |    0.1498 |  -0.06514 |
#> |    U|     48974.917 |     3.121 |     49.98 |    0.9449 |
#> |   90|     48872.476 |    -10.24 |    0.1698 |  -0.09206 |
#> |    U|     48872.476 |     3.112 |     49.98 |    0.9179 |
#> |   91|     48823.825 |    -10.29 |    0.1852 |  -0.09974 |
#> |    U|     48823.825 |     3.102 |     49.98 |    0.9103 |
#> |   92|     48823.818 |    -10.34 |    0.2075 |  -0.09251 |
#> |    U|     48823.818 |     3.092 |     49.98 |    0.9175 |
#> |   93|     48856.687 |    -10.33 |    0.2300 |  -0.08399 |
#> |    U|     48856.687 |     3.095 |     49.98 |    0.9260 |
#> |   94|     48739.019 |    -10.35 |    0.2134 |   -0.1190 |
#> |    U|     48739.019 |     3.091 |     49.98 |    0.8910 |
#> |   95|     48613.464 |    -10.38 |    0.2377 |   -0.1561 |
#> |    U|     48613.464 |     3.084 |     49.98 |    0.8539 |
#> |   96|     48427.837 |    -10.46 |    0.2859 |   -0.2074 |
#> |    U|     48427.837 |     3.067 |     49.99 |    0.8026 |
#> |   97|     48122.020 |    -10.65 |    0.3758 |   -0.2925 |
#> |    U|     48122.020 |     3.031 |     49.99 |    0.7175 |
#> |   98|     47863.190 |    -10.99 |    0.5747 |   -0.4789 |
#> |    U|     47863.190 |     2.961 |     49.99 |    0.5311 |
#> |   99|     48310.028 |    -11.42 |    0.5559 |   -0.6030 |
#> |    U|     48310.028 |     2.877 |     49.99 |    0.4070 |
#> |  100|     47919.935 |    -11.05 |    0.6582 |   -0.2853 |
#> |    U|     47919.935 |     2.950 |     49.99 |    0.7247 |
#> |  101|     47579.884 |    -11.20 |    0.6475 |   -0.4404 |
#> |    U|     47579.884 |     2.920 |     49.99 |    0.5696 |
#> |  102|     47758.451 |    -11.01 |    0.5707 |   -0.3647 |
#> |    U|     47758.451 |     2.959 |     49.99 |    0.6453 |
#> |  103|     47849.948 |    -11.16 |    0.7123 |   -0.5188 |
#> |    U|     47849.948 |     2.929 |     49.99 |    0.4912 |
#> |  104|     47514.933 |    -11.29 |    0.7101 |   -0.4470 |
#> |    U|     47514.933 |     2.902 |     49.99 |    0.5630 |
#> |  105|     47484.743 |    -11.39 |    0.7361 |   -0.4863 |
#> |    U|     47484.743 |     2.882 |     49.99 |    0.5237 |
#> |  106|     47426.720 |    -11.49 |    0.7634 |   -0.4974 |
#> |    U|     47426.720 |     2.861 |     50.00 |    0.5126 |
#> |  107|     47350.846 |    -11.71 |    0.8034 |   -0.5249 |
#> |    U|     47350.846 |     2.818 |     50.00 |    0.4851 |
#> |  108|     47489.604 |    -11.91 |    0.8761 |   -0.5789 |
#> |    U|     47489.604 |     2.778 |     50.00 |    0.4311 |
#> |  109|     47808.634 |    -11.70 |    0.7744 |   -0.5942 |
#> |    U|     47808.634 |     2.819 |     50.00 |    0.4158 |
#> |  110|     47284.956 |    -11.78 |    0.8383 |   -0.4433 |
#> |    U|     47284.956 |     2.805 |     50.00 |    0.5667 |
#> |  111|     47256.567 |    -11.79 |    0.9459 |   -0.4634 |
#> |    U|     47256.567 |     2.802 |     50.00 |    0.5466 |
#> |  112|     47241.631 |    -11.84 |    0.9207 |   -0.4641 |
#> |    U|     47241.631 |     2.792 |     50.00 |    0.5459 |
#> |  113|     47252.856 |    -11.79 |    0.8998 |   -0.4699 |
#> |    U|     47252.856 |     2.803 |     50.00 |    0.5401 |
#> |  114|     47265.850 |    -11.86 |    0.9216 |   -0.4433 |
#> |    U|     47265.850 |     2.789 |     50.00 |    0.5667 |
#> |  115|     47234.957 |    -11.85 |    0.9451 |   -0.4677 |
#> |    U|     47234.957 |     2.790 |     50.00 |    0.5423 |
#> |  116|     47230.279 |    -11.87 |    0.9617 |   -0.4663 |
#> |    U|     47230.279 |     2.786 |     50.00 |    0.5437 |
#> |  117|     47212.952 |    -11.92 |    0.9902 |   -0.4728 |
#> |    U|     47212.952 |     2.776 |     50.00 |    0.5372 |
#> |  118|     47183.913 |    -12.01 |     1.050 |   -0.4859 |
#> |    U|     47183.913 |     2.758 |     50.00 |    0.5241 |
#> |  119|     47168.795 |    -12.09 |     1.131 |   -0.4889 |
#> |    U|     47168.795 |     2.743 |     50.00 |    0.5211 |
#> |  120|     47149.655 |    -12.19 |     1.149 |   -0.5014 |
#> |    U|     47149.655 |     2.721 |     50.00 |    0.5086 |
#> |  121|     47143.706 |    -12.26 |     1.233 |   -0.5015 |
#> |    U|     47143.706 |     2.707 |     50.00 |    0.5085 |
#> |  122|     47160.872 |    -12.27 |     1.231 |   -0.4819 |
#> |    U|     47160.872 |     2.706 |     50.00 |    0.5281 |
#> |  123|     47141.587 |    -12.28 |     1.217 |   -0.5074 |
#> |    U|     47141.587 |     2.705 |     50.00 |    0.5026 |
#> |  124|     47140.716 |    -12.29 |     1.233 |   -0.5116 |
#> |    U|     47140.716 |     2.703 |     50.00 |    0.4984 |
#> |  125|     47141.773 |    -12.27 |     1.220 |   -0.5107 |
#> |    U|     47141.773 |     2.706 |     50.00 |    0.4993 |
#> |  126|     47141.386 |    -12.28 |     1.241 |   -0.5053 |
#> |    U|     47141.386 |     2.703 |     50.00 |    0.5047 |
#> |  127|     47140.116 |    -12.30 |     1.237 |   -0.5101 |
#> |    U|     47140.116 |     2.701 |     50.00 |    0.4999 |
#> |  128|     47139.330 |    -12.31 |     1.250 |   -0.5116 |
#> |    U|     47139.330 |     2.698 |     50.00 |    0.4984 |
#> |  129|     47138.325 |    -12.35 |     1.273 |   -0.5150 |
#> |    U|     47138.325 |     2.691 |     50.01 |    0.4950 |
#> |  130|     47139.086 |    -12.41 |     1.320 |   -0.5217 |
#> |    U|     47139.086 |     2.677 |     50.01 |    0.4883 |
#> |  131|     47141.927 |    -12.36 |     1.252 |   -0.5012 |
#> |    U|     47141.927 |     2.688 |     50.01 |    0.5088 |
#> |  132|     47137.817 |    -12.39 |     1.276 |   -0.5155 |
#> |    U|     47137.817 |     2.683 |     50.01 |    0.4945 |
#> |  133|     47143.653 |    -12.39 |     1.275 |   -0.5284 |
#> |    U|     47143.653 |     2.683 |     50.01 |    0.4816 |
#> |  134|     47138.022 |    -12.37 |     1.279 |   -0.5109 |
#> |    U|     47138.022 |     2.687 |     50.01 |    0.4991 |
#> |  135|     47137.916 |    -12.40 |     1.279 |   -0.5165 |
#> |    U|     47137.916 |     2.681 |     50.01 |    0.4935 |
#> |  136|     47137.784 |    -12.38 |     1.273 |   -0.5142 |
#> |    U|     47137.784 |     2.685 |     50.01 |    0.4958 |
#> |  137|     47137.866 |    -12.37 |     1.270 |   -0.5128 |
#> |    U|     47137.866 |     2.687 |     50.01 |    0.4972 |
#> |  138|     47138.027 |    -12.38 |     1.270 |   -0.5111 |
#> |    U|     47138.027 |     2.684 |     50.01 |    0.4989 |
#> |  139|     47137.778 |    -12.38 |     1.274 |   -0.5138 |
#> |    U|     47137.778 |     2.684 |     50.01 |    0.4962 |
#> |  140|     47137.785 |    -12.38 |     1.272 |   -0.5138 |
#> |    U|     47137.785 |     2.684 |     50.01 |    0.4962 |
#> |  141|     47137.777 |    -12.38 |     1.275 |   -0.5138 |
#> |    U|     47137.777 |     2.684 |     50.01 |    0.4962 |
#> |  142|     47137.781 |    -12.38 |     1.274 |   -0.5137 |
#> |    U|     47137.781 |     2.684 |     50.01 |    0.4963 |
#> |  143|     47137.776 |    -12.38 |     1.275 |   -0.5142 |
#> |    U|     47137.776 |     2.684 |     50.01 |    0.4958 |
#> |  144|     47137.773 |    -12.38 |     1.275 |   -0.5141 |
#> |    U|     47137.773 |     2.684 |     50.01 |    0.4959 |
#> |  145|     47137.771 |    -12.38 |     1.276 |   -0.5141 |
#> |    U|     47137.771 |     2.684 |     50.01 |    0.4959 |
#> |  146|     47137.768 |    -12.38 |     1.276 |   -0.5142 |
#> |    U|     47137.768 |     2.684 |     50.01 |    0.4958 |
#> |  147|     47137.763 |    -12.38 |     1.278 |   -0.5143 |
#> |    U|     47137.763 |     2.684 |     50.01 |    0.4957 |
#> |  148|     47137.760 |    -12.38 |     1.280 |   -0.5145 |
#> |    U|     47137.760 |     2.684 |     50.01 |    0.4955 |
#> |  149|     47137.758 |    -12.38 |     1.283 |   -0.5148 |
#> |    U|     47137.758 |     2.683 |     50.01 |    0.4952 |
#> |  150|     47137.750 |    -12.38 |     1.285 |   -0.5147 |
#> |    U|     47137.750 |     2.684 |     50.01 |    0.4953 |
#> |  151|     47137.739 |    -12.38 |     1.287 |   -0.5144 |
#> |    U|     47137.739 |     2.684 |     50.01 |    0.4956 |
#> |  152|     47137.725 |    -12.38 |     1.294 |   -0.5147 |
#> |    U|     47137.725 |     2.684 |     50.01 |    0.4953 |
#> |  153|     47137.708 |    -12.38 |     1.302 |   -0.5148 |
#> |    U|     47137.708 |     2.684 |     50.01 |    0.4952 |
#> |  154|     47137.698 |    -12.38 |     1.315 |   -0.5155 |
#> |    U|     47137.698 |     2.683 |     50.01 |    0.4945 |
#> |  155|     47137.651 |    -12.38 |     1.328 |   -0.5153 |
#> |    U|     47137.651 |     2.684 |     50.01 |    0.4947 |
#> |  156|     47137.619 |    -12.38 |     1.357 |   -0.5158 |
#> |    U|     47137.619 |     2.684 |     50.01 |    0.4942 |
#> |  157|     47137.604 |    -12.38 |     1.381 |   -0.5163 |
#> |    U|     47137.604 |     2.684 |     50.01 |    0.4937 |
#> |  158|     47137.628 |    -12.38 |     1.406 |   -0.5172 |
#> |    U|     47137.628 |     2.684 |     50.01 |    0.4928 |
#> |  159|     47139.602 |    -12.38 |     1.381 |   -0.5053 |
#> |    U|     47139.602 |     2.684 |     50.01 |    0.5047 |
#> |  160|     47137.393 |    -12.38 |     1.404 |   -0.5137 |
#> |    U|     47137.393 |     2.685 |     50.01 |    0.4963 |
#> |  161|     47137.455 |    -12.38 |     1.381 |   -0.5139 |
#> |    U|     47137.455 |     2.684 |     50.01 |    0.4961 |
#> |  162|     47140.343 |    -12.37 |     1.404 |   -0.5237 |
#> |    U|     47140.343 |     2.686 |     50.01 |    0.4863 |
#> |  163|     47137.384 |    -12.37 |     1.415 |   -0.5144 |
#> |    U|     47137.384 |     2.686 |     50.01 |    0.4956 |
#> |  164|     47137.405 |    -12.37 |     1.403 |   -0.5141 |
#> |    U|     47137.405 |     2.685 |     50.01 |    0.4959 |
#> |  165|     47137.851 |    -12.37 |     1.418 |   -0.5097 |
#> |    U|     47137.851 |     2.685 |     50.01 |    0.5003 |
#> |  166|     47137.349 |    -12.38 |     1.419 |   -0.5144 |
#> |    U|     47137.349 |     2.685 |     50.01 |    0.4956 |
#> |  167|     47137.340 |    -12.37 |     1.424 |   -0.5141 |
#> |    U|     47137.340 |     2.685 |     50.01 |    0.4959 |
#> |  168|     47137.317 |    -12.37 |     1.430 |   -0.5139 |
#> |    U|     47137.317 |     2.685 |     50.01 |    0.4961 |
#> |  169|     47137.276 |    -12.38 |     1.442 |   -0.5139 |
#> |    U|     47137.276 |     2.685 |     50.01 |    0.4961 |
#> |  170|     47137.216 |    -12.37 |     1.465 |   -0.5139 |
#> |    U|     47137.216 |     2.685 |     50.01 |    0.4961 |
#> |  171|     47137.140 |    -12.38 |     1.488 |   -0.5147 |
#> |    U|     47137.140 |     2.684 |     50.01 |    0.4953 |
#> |  172|     47137.087 |    -12.38 |     1.511 |   -0.5148 |
#> |    U|     47137.087 |     2.685 |     50.01 |    0.4952 |
#> |  173|     47137.075 |    -12.38 |     1.535 |   -0.5158 |
#> |    U|     47137.075 |     2.684 |     50.01 |    0.4942 |
#> |  174|     47136.926 |    -12.38 |     1.558 |   -0.5144 |
#> |    U|     47136.926 |     2.685 |     50.01 |    0.4956 |
#> |  175|     47136.781 |    -12.38 |     1.605 |   -0.5140 |
#> |    U|     47136.781 |     2.684 |     50.01 |    0.4960 |
#> |  176|     47136.501 |    -12.38 |     1.699 |   -0.5138 |
#> |    U|     47136.501 |     2.684 |     50.01 |    0.4962 |
#> |  177|     47135.937 |    -12.38 |     1.888 |   -0.5139 |
#> |    U|     47135.937 |     2.684 |     50.02 |    0.4961 |
#> |  178|     47134.843 |    -12.38 |     2.265 |   -0.5133 |
#> |    U|     47134.843 |     2.683 |     50.03 |    0.4967 |
#> |  179|     47132.845 |    -12.40 |     3.019 |   -0.5125 |
#> |    U|     47132.845 |     2.681 |     50.04 |    0.4975 |
#> |  180|     47129.411 |    -12.40 |     4.527 |   -0.5091 |
#> |    U|     47129.411 |     2.680 |     50.07 |    0.5009 |
#> |  181|     47129.647 |    -12.43 |     7.543 |   -0.4992 |
#> |    U|     47129.647 |     2.674 |     50.13 |    0.5108 |
#> |  182|     60348.283 |    -12.43 |     4.527 |   -0.8196 |
#> |    U|     60348.283 |     2.674 |     50.07 |    0.1904 |
#> |  183|     49108.982 |    -11.56 |     5.678 |  -0.01000 |
#> |    U|     49108.982 |     2.848 |     50.09 |     1.000 |
#> |  184|     47475.948 |    -12.57 |     4.290 |   -0.6120 |
#> |    U|     47475.948 |     2.646 |     50.07 |    0.3980 |
#> |  185|     49618.569 |    -12.86 |     4.849 |  -0.01000 |
#> |    U|     49618.569 |     2.588 |     50.08 |     1.000 |
#> |  186|     47180.880 |    -12.17 |     4.339 |   -0.4647 |
#> |    U|     47180.880 |     2.727 |     50.07 |    0.5453 |
#> |  187|     47493.469 |    -12.76 |     4.527 |   -0.4014 |
#> |    U|     47493.469 |     2.608 |     50.07 |    0.6086 |
#> |  188|     47130.518 |    -12.32 |     4.658 |   -0.5035 |
#> |    U|     47130.518 |     2.695 |     50.07 |    0.5065 |
#> |  189|     47643.587 |    -12.38 |     4.612 |   -0.3419 |
#> |    U|     47643.587 |     2.684 |     50.07 |    0.6681 |
#> |  190|     47129.048 |    -12.43 |     4.550 |   -0.5147 |
#> |    U|     47129.048 |     2.675 |     50.07 |    0.4953 |
#> |  191|     47219.819 |    -12.46 |     4.615 |   -0.4566 |
#> |    U|     47219.819 |     2.667 |     50.07 |    0.5534 |
#> |  192|     47129.305 |    -12.40 |     4.529 |   -0.5097 |
#> |    U|     47129.305 |     2.680 |     50.07 |    0.5003 |
#> |  193|     47183.681 |    -12.43 |     4.538 |   -0.4694 |
#> |    U|     47183.681 |     2.675 |     50.07 |    0.5406 |
#> |  194|     47130.096 |    -12.44 |     4.531 |   -0.5220 |
#> |    U|     47130.096 |     2.672 |     50.07 |    0.4880 |
#> |  195|     47135.192 |    -12.44 |     4.533 |   -0.5024 |
#> |    U|     47135.192 |     2.673 |     50.07 |    0.5076 |
#> |  196|     47131.728 |    -12.43 |     4.550 |   -0.5264 |
#> |    U|     47131.728 |     2.675 |     50.07 |    0.4836 |
#> |  197|     47129.031 |    -12.43 |     4.562 |   -0.5184 |
#> |    U|     47129.031 |     2.674 |     50.07 |    0.4916 |
#> |  198|     47128.834 |    -12.42 |     4.565 |   -0.5178 |
#> |    U|     47128.834 |     2.675 |     50.07 |    0.4922 |
#> |  199|     47128.901 |    -12.43 |     4.571 |   -0.5180 |
#> |    U|     47128.901 |     2.675 |     50.07 |    0.4920 |
#> |  200|     47128.854 |    -12.43 |     4.567 |   -0.5158 |
#> |    U|     47128.854 |     2.675 |     50.07 |    0.4942 |
#> |  201|     47128.812 |    -12.42 |     4.563 |   -0.5176 |
#> |    U|     47128.812 |     2.675 |     50.07 |    0.4924 |
#> |  202|     47128.776 |    -12.42 |     4.561 |   -0.5173 |
#> |    U|     47128.776 |     2.675 |     50.07 |    0.4927 |
#> |  203|     47128.826 |    -12.42 |     4.560 |   -0.5181 |
#> |    U|     47128.826 |     2.675 |     50.07 |    0.4919 |
#> |  204|     47128.776 |    -12.42 |     4.561 |   -0.5170 |
#> |    U|     47128.776 |     2.675 |     50.07 |    0.4930 |
#> |  205|     47128.760 |    -12.42 |     4.562 |   -0.5171 |
#> |    U|     47128.760 |     2.675 |     50.07 |    0.4929 |
#> |  206|     47128.725 |    -12.42 |     4.563 |   -0.5170 |
#> |    U|     47128.725 |     2.676 |     50.07 |    0.4930 |
#> |  207|     47128.693 |    -12.42 |     4.564 |   -0.5168 |
#> |    U|     47128.693 |     2.676 |     50.07 |    0.4932 |
#> |  208|     47128.686 |    -12.42 |     4.566 |   -0.5159 |
#> |    U|     47128.686 |     2.676 |     50.07 |    0.4941 |
#> |  209|     47128.663 |    -12.42 |     4.566 |   -0.5162 |
#> |    U|     47128.663 |     2.676 |     50.07 |    0.4938 |
#> |  210|     47128.639 |    -12.42 |     4.566 |   -0.5164 |
#> |    U|     47128.639 |     2.676 |     50.07 |    0.4936 |
#> |  211|     47128.598 |    -12.42 |     4.566 |   -0.5166 |
#> |    U|     47128.598 |     2.676 |     50.07 |    0.4934 |
#> |  212|     47128.530 |    -12.42 |     4.568 |   -0.5169 |
#> |    U|     47128.530 |     2.677 |     50.07 |    0.4931 |
#> |  213|     47128.463 |    -12.41 |     4.567 |   -0.5170 |
#> |    U|     47128.463 |     2.677 |     50.07 |    0.4930 |
#> |  214|     47128.409 |    -12.41 |     4.569 |   -0.5171 |
#> |    U|     47128.409 |     2.678 |     50.07 |    0.4929 |
#> |  215|     47128.364 |    -12.41 |     4.569 |   -0.5174 |
#> |    U|     47128.364 |     2.679 |     50.07 |    0.4926 |
#> |  216|     47128.218 |    -12.40 |     4.571 |   -0.5168 |
#> |    U|     47128.218 |     2.680 |     50.07 |    0.4932 |
#> |  217|     47128.023 |    -12.39 |     4.576 |   -0.5162 |
#> |    U|     47128.023 |     2.682 |     50.07 |    0.4938 |
#> |  218|     47127.918 |    -12.37 |     4.584 |   -0.5161 |
#> |    U|     47127.918 |     2.686 |     50.07 |    0.4939 |
#> |  219|     47127.851 |    -12.37 |     4.583 |   -0.5150 |
#> |    U|     47127.851 |     2.686 |     50.07 |    0.4950 |
#> |  220|     47127.850 |    -12.37 |     4.584 |   -0.5149 |
#> |    U|     47127.850 |     2.687 |     50.07 |    0.4951 |
#> |  221|     47127.853 |    -12.37 |     4.584 |   -0.5150 |
#> |    U|     47127.853 |     2.687 |     50.07 |    0.4950 |
#> |  222|     47127.858 |    -12.37 |     4.584 |   -0.5151 |
#> |    U|     47127.858 |     2.687 |     50.07 |    0.4949 |
#> |  223|     47127.846 |    -12.37 |     4.584 |   -0.5148 |
#> |    U|     47127.846 |     2.686 |     50.07 |    0.4952 |
#> |  224|     47127.842 |    -12.37 |     4.584 |   -0.5146 |
#> |    U|     47127.842 |     2.686 |     50.07 |    0.4954 |
#> |  225|     47127.841 |    -12.37 |     4.584 |   -0.5146 |
#> |    U|     47127.841 |     2.686 |     50.07 |    0.4954 |
#> |  226|     47127.842 |    -12.37 |     4.584 |   -0.5147 |
#> |    U|     47127.842 |     2.686 |     50.07 |    0.4953 |
#> |  227|     47127.843 |    -12.37 |     4.584 |   -0.5148 |
#> |    U|     47127.843 |     2.686 |     50.07 |    0.4952 |
#> |  228|     47127.841 |    -12.37 |     4.584 |   -0.5146 |
#> |    U|     47127.841 |     2.686 |     50.07 |    0.4954 |
#> |  229|     47127.841 |    -12.37 |     4.584 |   -0.5146 |
#> |    U|     47127.841 |     2.686 |     50.07 |    0.4954 |
#> |  230|     47127.840 |    -12.37 |     4.584 |   -0.5145 |
#> |    U|     47127.840 |     2.686 |     50.07 |    0.4955 |
#> |  231|     47127.839 |    -12.37 |     4.585 |   -0.5144 |
#> |    U|     47127.839 |     2.686 |     50.07 |    0.4956 |
#> |  232|     47127.838 |    -12.37 |     4.585 |   -0.5142 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4958 |
#> |  233|     47127.838 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4957 |
#> |  234|     47127.838 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4957 |
#> |  235|     47127.840 |    -12.37 |     4.585 |   -0.5141 |
#> |    U|     47127.840 |     2.686 |     50.07 |    0.4959 |
#> |  236|     47127.838 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4957 |
#> |  237|     47127.838 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4957 |
#> |  238|     47127.838 |    -12.37 |     4.585 |   -0.5142 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4958 |
#> |  239|     47127.838 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4957 |
#> |  240|     47127.838 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4957 |
#> |  241|     47127.838 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4957 |
#> |  242|     47127.838 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4957 |
#> |  243|     47127.838 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4957 |
#> |  244|     47127.838 |    -12.37 |     4.585 |   -0.5144 |
#> |    U|     47127.838 |     2.686 |     50.07 |    0.4956 |
#> |  245|     47127.837 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.837 |     2.686 |     50.07 |    0.4957 |
#> |  246|     47127.837 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.837 |     2.686 |     50.07 |    0.4957 |
#> |  247|     47127.837 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.837 |     2.686 |     50.07 |    0.4957 |
#> |  248|     47127.837 |    -12.37 |     4.585 |   -0.5143 |
#> |    U|     47127.837 |     2.686 |     50.07 |    0.4957 |
#> |  249|     47127.837 |    -12.37 |     4.585 |   -0.5144 |
#> |    U|     47127.837 |     2.686 |     50.07 |    0.4956 |
#> |  250|     47127.837 |    -12.37 |     4.585 |   -0.5144 |
#> |    U|     47127.837 |     2.686 |     50.07 |    0.4956 |
#> |  251|     47127.837 |    -12.37 |     4.585 |   -0.5145 |
#> |    U|     47127.837 |     2.686 |     50.07 |    0.4955 |
#> |  252|     47127.837 |    -12.37 |     4.585 |   -0.5144 |
#> |    U|     47127.837 |     2.686 |     50.07 |    0.4956 |
#> |  253|     47127.836 |    -12.37 |     4.585 |   -0.5144 |
#> |    U|     47127.836 |     2.686 |     50.07 |    0.4956 |
#> |  254|     47127.836 |    -12.37 |     4.585 |   -0.5144 |
#> |    U|     47127.836 |     2.686 |     50.07 |    0.4956 |
#> |  255|     47127.836 |    -12.37 |     4.586 |   -0.5143 |
#> |    U|     47127.836 |     2.686 |     50.07 |    0.4957 |
#> |  256|     47127.836 |    -12.37 |     4.586 |   -0.5143 |
#> |    U|     47127.836 |     2.686 |     50.07 |    0.4957 |
#> |  257|     47127.836 |    -12.37 |     4.586 |   -0.5142 |
#> |    U|     47127.836 |     2.686 |     50.07 |    0.4958 |
#> |  258|     47127.836 |    -12.37 |     4.586 |   -0.5142 |
#> |    U|     47127.836 |     2.686 |     50.07 |    0.4958 |
#> |  259|     47127.834 |    -12.37 |     4.586 |   -0.5143 |
#> |    U|     47127.834 |     2.686 |     50.07 |    0.4957 |
#> |  260|     47127.833 |    -12.37 |     4.586 |   -0.5144 |
#> |    U|     47127.833 |     2.686 |     50.07 |    0.4956 |
#> |  261|     47127.833 |    -12.37 |     4.586 |   -0.5145 |
#> |    U|     47127.833 |     2.686 |     50.07 |    0.4955 |
#> |  262|     47127.832 |    -12.37 |     4.587 |   -0.5145 |
#> |    U|     47127.832 |     2.686 |     50.07 |    0.4955 |
#> |  263|     47127.830 |    -12.37 |     4.588 |   -0.5145 |
#> |    U|     47127.830 |     2.686 |     50.07 |    0.4955 |
#> |  264|     47127.827 |    -12.37 |     4.589 |   -0.5144 |
#> |    U|     47127.827 |     2.685 |     50.07 |    0.4956 |
#> |  265|     47127.826 |    -12.38 |     4.592 |   -0.5145 |
#> |    U|     47127.826 |     2.685 |     50.07 |    0.4955 |
#> |  266|     47127.826 |    -12.37 |     4.591 |   -0.5144 |
#> |    U|     47127.826 |     2.685 |     50.07 |    0.4956 |
#> |  267|     47127.826 |    -12.38 |     4.592 |   -0.5150 |
#> |    U|     47127.826 |     2.685 |     50.07 |    0.4950 |
#> |  268|     47127.853 |    -12.38 |     4.592 |   -0.5157 |
#> |    U|     47127.853 |     2.685 |     50.07 |    0.4943 |
#> |  269|     47127.825 |    -12.38 |     4.593 |   -0.5147 |
#> |    U|     47127.825 |     2.685 |     50.07 |    0.4953 |
#> |  270|     47127.825 |    -12.38 |     4.592 |   -0.5148 |
#> |    U|     47127.825 |     2.685 |     50.07 |    0.4952 |
#> |  271|     47127.827 |    -12.38 |     4.593 |   -0.5144 |
#> |    U|     47127.827 |     2.685 |     50.07 |    0.4956 |
#> |  272|     47127.823 |    -12.38 |     4.593 |   -0.5147 |
#> |    U|     47127.823 |     2.685 |     50.07 |    0.4953 |
#> |  273|     47127.822 |    -12.38 |     4.594 |   -0.5147 |
#> |    U|     47127.822 |     2.685 |     50.07 |    0.4953 |
#> |  274|     47127.822 |    -12.38 |     4.594 |   -0.5147 |
#> |    U|     47127.822 |     2.685 |     50.07 |    0.4953 |
#> |  275|     47127.822 |    -12.38 |     4.595 |   -0.5147 |
#> |    U|     47127.822 |     2.685 |     50.07 |    0.4953 |
#> |  276|     47127.822 |    -12.38 |     4.595 |   -0.5147 |
#> |    U|     47127.822 |     2.685 |     50.07 |    0.4953 |
#> |  277|     47127.821 |    -12.38 |     4.596 |   -0.5147 |
#> |    U|     47127.821 |     2.685 |     50.07 |    0.4953 |
#> |  278|     47127.820 |    -12.38 |     4.597 |   -0.5147 |
#> |    U|     47127.820 |     2.685 |     50.07 |    0.4953 |
#> |  279|     47127.819 |    -12.38 |     4.598 |   -0.5147 |
#> |    U|     47127.819 |     2.684 |     50.07 |    0.4953 |
#> |  280|     47127.818 |    -12.38 |     4.600 |   -0.5147 |
#> |    U|     47127.818 |     2.684 |     50.07 |    0.4953 |
#> |  281|     47127.822 |    -12.38 |     4.601 |   -0.5145 |
#> |    U|     47127.822 |     2.684 |     50.07 |    0.4955 |
#> |  282|     47127.830 |    -12.38 |     4.600 |   -0.5141 |
#> |    U|     47127.830 |     2.684 |     50.07 |    0.4959 |
#> |  283|     47127.816 |    -12.38 |     4.601 |   -0.5147 |
#> |    U|     47127.816 |     2.684 |     50.07 |    0.4953 |
#> |  284|     47127.815 |    -12.38 |     4.602 |   -0.5148 |
#> |    U|     47127.815 |     2.684 |     50.07 |    0.4952 |
#> |  285|     47127.816 |    -12.38 |     4.605 |   -0.5148 |
#> |    U|     47127.816 |     2.684 |     50.07 |    0.4952 |
#> |  286|     47127.844 |    -12.38 |     4.602 |   -0.5138 |
#> |    U|     47127.844 |     2.684 |     50.07 |    0.4962 |
#> |  287|     47127.813 |    -12.38 |     4.604 |   -0.5148 |
#> |    U|     47127.813 |     2.684 |     50.07 |    0.4952 |
#> |  288|     47127.811 |    -12.38 |     4.607 |   -0.5148 |
#> |    U|     47127.811 |     2.684 |     50.07 |    0.4952 |
#> |  289|     47127.809 |    -12.38 |     4.613 |   -0.5149 |
#> |    U|     47127.809 |     2.683 |     50.07 |    0.4951 |
#> |  290|     47127.809 |    -12.39 |     4.617 |   -0.5150 |
#> |    U|     47127.809 |     2.683 |     50.07 |    0.4950 |
#> |  291|     47127.873 |    -12.38 |     4.613 |   -0.5135 |
#> |    U|     47127.873 |     2.683 |     50.07 |    0.4965 |
#> |  292|     47127.808 |    -12.38 |     4.613 |   -0.5149 |
#> |    U|     47127.808 |     2.683 |     50.07 |    0.4951 |
#> |  293|     47127.805 |    -12.38 |     4.615 |   -0.5149 |
#> |    U|     47127.805 |     2.683 |     50.07 |    0.4951 |
#> |  294|     47127.802 |    -12.38 |     4.619 |   -0.5149 |
#> |    U|     47127.802 |     2.683 |     50.07 |    0.4951 |
#> |  295|     47127.802 |    -12.39 |     4.626 |   -0.5150 |
#> |    U|     47127.802 |     2.683 |     50.07 |    0.4950 |
#> |  296|     47128.007 |    -12.39 |     4.626 |   -0.5126 |
#> |    U|     47128.007 |     2.682 |     50.07 |    0.4974 |
#> |  297|     47127.806 |    -12.39 |     4.628 |   -0.5149 |
#> |    U|     47127.806 |     2.683 |     50.07 |    0.4951 |
#> |  298|     47127.804 |    -12.39 |     4.627 |   -0.5143 |
#> |    U|     47127.804 |     2.683 |     50.07 |    0.4957 |
#> |  299|     47127.801 |    -12.39 |     4.628 |   -0.5146 |
#> |    U|     47127.801 |     2.683 |     50.07 |    0.4954 |
#> |  300|     47127.801 |    -12.39 |     4.626 |   -0.5146 |
#> |    U|     47127.801 |     2.683 |     50.07 |    0.4954 |
#> |  301|     47127.800 |    -12.39 |     4.625 |   -0.5146 |
#> |    U|     47127.800 |     2.683 |     50.07 |    0.4954 |
#> |  302|     47127.802 |    -12.39 |     4.624 |   -0.5148 |
#> |    U|     47127.802 |     2.683 |     50.07 |    0.4952 |
#> |  303|     47127.801 |    -12.39 |     4.625 |   -0.5145 |
#> |    U|     47127.801 |     2.683 |     50.07 |    0.4955 |
#> |  304|     47127.800 |    -12.39 |     4.625 |   -0.5148 |
#> |    U|     47127.800 |     2.683 |     50.07 |    0.4952 |
#> |  305|     47127.802 |    -12.39 |     4.625 |   -0.5147 |
#> |    U|     47127.802 |     2.683 |     50.07 |    0.4953 |
#> |  306|     47127.798 |    -12.39 |     4.625 |   -0.5148 |
#> |    U|     47127.798 |     2.683 |     50.07 |    0.4952 |
#> |  307|     47127.796 |    -12.39 |     4.625 |   -0.5149 |
#> |    U|     47127.796 |     2.683 |     50.07 |    0.4951 |
#> |  308|     47127.792 |    -12.39 |     4.625 |   -0.5150 |
#> |    U|     47127.792 |     2.683 |     50.07 |    0.4950 |
#> |  309|     47127.787 |    -12.39 |     4.626 |   -0.5150 |
#> |    U|     47127.787 |     2.683 |     50.07 |    0.4950 |
#> |  310|     47127.777 |    -12.38 |     4.626 |   -0.5149 |
#> |    U|     47127.777 |     2.683 |     50.07 |    0.4951 |
#> |  311|     47127.760 |    -12.38 |     4.626 |   -0.5146 |
#> |    U|     47127.760 |     2.684 |     50.07 |    0.4954 |
#> |  312|     47127.749 |    -12.38 |     4.627 |   -0.5139 |
#> |    U|     47127.749 |     2.685 |     50.07 |    0.4961 |
#> |  313|     47127.738 |    -12.38 |     4.623 |   -0.5143 |
#> |    U|     47127.738 |     2.685 |     50.07 |    0.4957 |
#> |  314|     47127.740 |    -12.37 |     4.623 |   -0.5140 |
#> |    U|     47127.740 |     2.685 |     50.07 |    0.4960 |
#> |  315|     47127.749 |    -12.38 |     4.623 |   -0.5155 |
#> |    U|     47127.749 |     2.685 |     50.07 |    0.4945 |
#> |  316|     47127.733 |    -12.37 |     4.622 |   -0.5147 |
#> |    U|     47127.733 |     2.685 |     50.07 |    0.4953 |
#> |  317|     47127.728 |    -12.37 |     4.623 |   -0.5146 |
#> |    U|     47127.728 |     2.685 |     50.07 |    0.4954 |
#> |  318|     47127.720 |    -12.37 |     4.625 |   -0.5143 |
#> |    U|     47127.720 |     2.686 |     50.07 |    0.4957 |
#> |  319|     47127.727 |    -12.37 |     4.628 |   -0.5136 |
#> |    U|     47127.727 |     2.687 |     50.07 |    0.4964 |
#> |  320|     47127.721 |    -12.37 |     4.624 |   -0.5145 |
#> |    U|     47127.721 |     2.686 |     50.07 |    0.4955 |
#> |  321|     47127.724 |    -12.37 |     4.624 |   -0.5141 |
#> |    U|     47127.724 |     2.686 |     50.07 |    0.4959 |
#> |  322|     47127.719 |    -12.37 |     4.626 |   -0.5142 |
#> |    U|     47127.719 |     2.686 |     50.07 |    0.4958 |
#> |  323|     47127.719 |    -12.37 |     4.627 |   -0.5140 |
#> |    U|     47127.719 |     2.686 |     50.07 |    0.4960 |
#> |  324|     47127.719 |    -12.37 |     4.627 |   -0.5141 |
#> |    U|     47127.719 |     2.686 |     50.07 |    0.4959 |
#> |  325|     47127.718 |    -12.37 |     4.627 |   -0.5141 |
#> |    U|     47127.718 |     2.686 |     50.07 |    0.4959 |
#> |  326|     47127.718 |    -12.37 |     4.627 |   -0.5140 |
#> |    U|     47127.718 |     2.686 |     50.07 |    0.4960 |
#> |  327|     47127.719 |    -12.37 |     4.628 |   -0.5139 |
#> |    U|     47127.719 |     2.686 |     50.07 |    0.4961 |
#> |  328|     47127.723 |    -12.37 |     4.627 |   -0.5139 |
#> |    U|     47127.723 |     2.686 |     50.07 |    0.4961 |
#> |  329|     47127.712 |    -12.37 |     4.627 |   -0.5144 |
#> |    U|     47127.712 |     2.686 |     50.07 |    0.4956 |
#> |  330|     47127.709 |    -12.37 |     4.628 |   -0.5145 |
#> |    U|     47127.709 |     2.686 |     50.07 |    0.4955 |
#> |  331|     47127.706 |    -12.37 |     4.629 |   -0.5145 |
#> |    U|     47127.706 |     2.686 |     50.07 |    0.4955 |
#> |  332|     47127.705 |    -12.37 |     4.630 |   -0.5145 |
#> |    U|     47127.705 |     2.685 |     50.07 |    0.4955 |
#> |  333|     47127.705 |    -12.37 |     4.630 |   -0.5145 |
#> |    U|     47127.705 |     2.685 |     50.07 |    0.4955 |
#> |  334|     47127.705 |    -12.37 |     4.631 |   -0.5145 |
#> |    U|     47127.705 |     2.685 |     50.07 |    0.4955 |
#> |  335|     47127.706 |    -12.37 |     4.630 |   -0.5148 |
#> |    U|     47127.706 |     2.685 |     50.07 |    0.4952 |
#> |  336|     47127.703 |    -12.37 |     4.631 |   -0.5146 |
#> |    U|     47127.703 |     2.685 |     50.07 |    0.4954 |
#> |  337|     47127.700 |    -12.37 |     4.631 |   -0.5146 |
#> |    U|     47127.700 |     2.685 |     50.07 |    0.4954 |
#> |  338|     47127.696 |    -12.37 |     4.633 |   -0.5145 |
#> |    U|     47127.696 |     2.686 |     50.07 |    0.4955 |
#> |  339|     47127.692 |    -12.37 |     4.634 |   -0.5144 |
#> |    U|     47127.692 |     2.686 |     50.07 |    0.4956 |
#> |  340|     47127.689 |    -12.37 |     4.635 |   -0.5143 |
#> |    U|     47127.689 |     2.686 |     50.07 |    0.4957 |
#> |  341|     47127.687 |    -12.37 |     4.635 |   -0.5143 |
#> |    U|     47127.687 |     2.686 |     50.07 |    0.4957 |
#> |  342|     47127.687 |    -12.37 |     4.636 |   -0.5141 |
#> |    U|     47127.687 |     2.686 |     50.07 |    0.4959 |
#> |  343|     47127.685 |    -12.37 |     4.636 |   -0.5143 |
#> |    U|     47127.685 |     2.686 |     50.07 |    0.4957 |
#> |  344|     47127.683 |    -12.37 |     4.637 |   -0.5144 |
#> |    U|     47127.683 |     2.686 |     50.07 |    0.4956 |
#> |  345|     47127.678 |    -12.37 |     4.638 |   -0.5144 |
#> |    U|     47127.678 |     2.686 |     50.07 |    0.4956 |
#> |  346|     47127.670 |    -12.37 |     4.641 |   -0.5146 |
#> |    U|     47127.670 |     2.686 |     50.07 |    0.4954 |
#> |  347|     47127.652 |    -12.37 |     4.647 |   -0.5145 |
#> |    U|     47127.652 |     2.686 |     50.07 |    0.4955 |
#> |  348|     47127.617 |    -12.37 |     4.659 |   -0.5146 |
#> |    U|     47127.617 |     2.686 |     50.07 |    0.4954 |
#> |  349|     47127.550 |    -12.37 |     4.682 |   -0.5143 |
#> |    U|     47127.550 |     2.687 |     50.07 |    0.4957 |
#> |  350|     47127.406 |    -12.37 |     4.729 |   -0.5146 |
#> |    U|     47127.406 |     2.686 |     50.07 |    0.4954 |
#> |  351|     47127.124 |    -12.37 |     4.823 |   -0.5143 |
#> |    U|     47127.124 |     2.687 |     50.08 |    0.4957 |
#> |  352|     47126.558 |    -12.37 |     5.012 |   -0.5145 |
#> |    U|     47126.558 |     2.686 |     50.08 |    0.4955 |
#> |  353|     47125.519 |    -12.36 |     5.389 |   -0.5156 |
#> |    U|     47125.519 |     2.689 |     50.09 |    0.4944 |
#> |  354|     47123.170 |    -12.37 |     6.143 |   -0.5147 |
#> |    U|     47123.170 |     2.687 |     50.10 |    0.4953 |
#> |  355|     47118.704 |    -12.35 |     7.651 |   -0.5144 |
#> |    U|     47118.704 |     2.691 |     50.13 |    0.4956 |
#> |  356|     47109.646 |    -12.35 |     10.67 |   -0.5156 |
#> |    U|     47109.646 |     2.690 |     50.19 |    0.4944 |
#> |  357|     47091.713 |    -12.32 |     16.70 |   -0.5161 |
#> |    U|     47091.713 |     2.696 |     50.31 |    0.4939 |
#> |  358|     47056.010 |    -12.29 |     28.77 |   -0.5185 |
#> |    U|     47056.010 |     2.702 |     50.56 |    0.4915 |
#> |  359|     46985.647 |    -12.19 |     52.90 |   -0.5213 |
#> |    U|     46985.647 |     2.722 |     51.04 |    0.4887 |
#> |  360|     46847.910 |    -12.04 |     101.2 |   -0.5258 |
#> |    U|     46847.910 |     2.752 |     52.00 |    0.4842 |
#> |  361|     46586.818 |    -11.73 |     197.7 |   -0.5339 |
#> |    U|     46586.818 |     2.814 |     53.93 |    0.4761 |
#> |  362|     46124.668 |    -11.02 |     390.8 |   -0.5455 |
#> |    U|     46124.668 |     2.956 |     57.80 |    0.4645 |
#> |  363|     45491.819 |    -9.827 |     776.9 |   -0.5936 |
#> |    U|     45491.819 |     3.195 |     65.52 |    0.4164 |
#> |  364|     45442.746 |    -9.803 |     776.9 |   -0.5274 |
#> |    U|     45442.746 |     3.199 |     65.52 |    0.4826 |
#> |  365|     45077.334 |    -8.358 |     1163. |   -0.5630 |
#> |    U|     45077.334 |     3.489 |     73.24 |    0.4470 |
#> |  366|     45442.746 |    -9.803 |     776.9 |   -0.5274 |
#> |    U|     45442.746 |     3.199 |     65.52 |    0.4826 |
#> 
#> → creating full model...
#> → pruning branches (`if`/`else`)...
#>  done
#> → loading into symengine environment...
#>  done
#> → finding duplicate expressions in EBE model...
#> 
#> → optimizing duplicate expressions in EBE model...
#> 
#> → compiling EBE model...
#>  
#>  done
#> Calculating residuals/tables
#> done

if (inherits(fit, "nlmixrDynmodel")) {
 as.dynmodel(fit)
}
#>             est       %cv
#> CL    3.1994401 0.5474521
#> V    65.5187671 0.7156383
#> prop  0.4825785 1.0149641
#> 
#>  -loglik      AIC      BIC 
#> 45442.75 90891.49 90912.03 
#> 

# method = "focei" is slightly more flexible and well tested

fit <- try(nlmixr(object = model_onecmt_bolus, data = Bolus_1CPT, est = "focei"))
#>  parameter labels from comments will be replaced by 'label()'
#> Key: U: Unscaled Parameters; X: Back-transformed parameters; G: Gill difference gradient approximation
#> F: Forward difference gradient approximation
#> C: Central difference gradient approximation
#> M: Mixed forward and central difference gradient approximation
#> Unscaled parameters for Omegas=chol(solve(omega));
#> Diagonals are transformed, as specified by foceiControl(diagXform=)
#> |-----+---------------+-----------+-----------+-----------|
#> |    #| Objective Fun |        CL |         V |  prop.err |
#> |    1| 8.0454138e+09 |   -0.8004 |     1.000 |    -1.000 |
#> |    U| 8.0454138e+09 |     5.000 |     50.00 |   0.01000 |
#> |    X| 8.0454138e+09 |     5.000 |     50.00 |   0.01000 |
#> |    G|    Gill Diff. | 1.154e+11 |-9.901e+10 |-8.044e+09 |
#> |    2|     25787707. |    -1.558 |     1.650 |   -0.9472 |
#> |    U|      25787707 |     1.210 |     82.51 |   0.01026 |
#> |    X|      25787707 |     1.210 |     82.51 |   0.01026 |
#> |    M|   Mixed Diff. |-5.748e+07 | 1.415e+06 |-2.509e+07 |
#> |    3|     12967373. |    -1.286 |     1.943 |   -0.6466 |
#> |    U|      12967373 |     2.572 |     97.13 |   0.01177 |
#> |    X|      12967373 |     2.572 |     97.13 |   0.01177 |
#> |    M|   Mixed Diff. |-1.429e+07 | 7.740e+06 |-1.100e+07 |
#> |    4|     7423175.1 |    -1.178 |     1.725 |   -0.2093 |
#> |    U|     7423175.1 |     3.111 |     86.27 |   0.01395 |
#> |    X|     7423175.1 |     3.111 |     86.27 |   0.01395 |
#> |    M|   Mixed Diff. |-1.180e+06 | 4.378e+06 |-5.298e+06 |
#> |    5|     5484755.4 |    -1.144 |     1.430 |    0.1923 |
#> |    U|     5484755.4 |     3.281 |     71.48 |   0.01596 |
#> |    X|     5484755.4 |     3.281 |     71.48 |   0.01596 |
#> |    M|   Mixed Diff. | 1.408e+07 |-2.824e+06 |-3.417e+06 |
#> |    6|     8000501.1 |    -1.609 |     1.111 |     1.018 |
#> |    U|     8000501.1 |    0.9563 |     55.53 |   0.02009 |
#> |    X|     8000501.1 |    0.9563 |     55.53 |   0.02009 |
#> |    7|     6450161.2 |    -1.418 |     1.370 |    0.4573 |
#> |    U|     6450161.2 |     1.914 |     68.48 |   0.01729 |
#> |    X|     6450161.2 |     1.914 |     68.48 |   0.01729 |
#> |    8|     5616884.6 |    -1.296 |     1.460 |    0.2292 |
#> |    U|     5616884.6 |     2.522 |     73.00 |   0.01615 |
#> |    X|     5616884.6 |     2.522 |     73.00 |   0.01615 |
#> |    9|     5075129.6 |    -1.216 |     1.444 |    0.2098 |
#> |    U|     5075129.6 |     2.921 |     72.20 |   0.01605 |
#> |    X|     5075129.6 |     2.921 |     72.20 |   0.01605 |
#> |    M|   Mixed Diff. |-2.409e+06 | 9.218e+05 |-3.142e+06 |
#> |   10|     4863902.1 |    -1.192 |     1.417 |    0.2758 |
#> |    U|     4863902.1 |     3.044 |     70.86 |   0.01638 |
#> |    X|     4863902.1 |     3.044 |     70.86 |   0.01638 |
#> |    M|   Mixed Diff. | 2.610e+06 |-2.612e+05 |-2.950e+06 |
#> |   11|     4566749.8 |    -1.247 |     1.373 |    0.4087 |
#> |    U|     4566749.8 |     2.765 |     68.64 |   0.01704 |
#> |    X|     4566749.8 |     2.765 |     68.64 |   0.01704 |
#> |    M|   Mixed Diff. |-4.993e+06 | 5.568e+05 |-2.659e+06 |
#> |   12|     4266337.6 |    -1.207 |     1.213 |    0.6609 |
#> |    U|     4266337.6 |     2.967 |     60.64 |   0.01830 |
#> |    X|     4266337.6 |     2.967 |     60.64 |   0.01830 |
#> |    M|   Mixed Diff. | 1.123e+07 |-5.104e+06 |-2.312e+06 |
#> |   13|     4128738.8 |    -1.335 |     1.077 |    0.8975 |
#> |    U|     4128738.8 |     2.324 |     53.85 |   0.01949 |
#> |    X|     4128738.8 |     2.324 |     53.85 |   0.01949 |
#> |    M|   Mixed Diff. |-8.655e+06 |-9.469e+05 |-2.099e+06 |
#> |   14|     3928525.9 |    -1.248 |     1.270 |    0.6829 |
#> |    U|     3928525.9 |     2.760 |     63.50 |   0.01841 |
#> |    X|     3928525.9 |     2.760 |     63.50 |   0.01841 |
#> |    M|   Mixed Diff. |-2.176e+06 |-7.457e+05 |-2.114e+06 |
#> |   15|     4854860.6 |    -1.180 |     1.139 |    0.9455 |
#> |    U|     4854860.6 |     3.103 |     56.95 |   0.01973 |
#> |    X|     4854860.6 |     3.103 |     56.95 |   0.01973 |
#> |   16|     3876327.0 |    -1.207 |     1.284 |    0.7232 |
#> |    U|       3876327 |     2.968 |     64.22 |   0.01862 |
#> |    X|       3876327 |     2.968 |     64.22 |   0.01862 |
#> |    M|   Mixed Diff. | 5.070e+06 |-2.312e+06 |-2.064e+06 |
#> |   17|     3754785.9 |    -1.245 |     1.267 |    0.7659 |
#> |    U|     3754785.9 |     2.778 |     63.34 |   0.01883 |
#> |    X|     3754785.9 |     2.778 |     63.34 |   0.01883 |
#> |    M|   Mixed Diff. |-1.408e+06 |-8.831e+05 |-1.975e+06 |
#> |   18|     3650705.3 |    -1.223 |     1.225 |    0.8756 |
#> |    U|     3650705.3 |     2.888 |     61.26 |   0.01938 |
#> |    X|     3650705.3 |     2.888 |     61.26 |   0.01938 |
#> |    M|   Mixed Diff. | 4.818e+06 |-2.812e+06 |-1.866e+06 |
#> |   19|     3602499.8 |    -1.302 |     1.089 |     1.055 |
#> |    U|     3602499.8 |     2.490 |     54.47 |   0.02028 |
#> |    X|     3602499.8 |     2.490 |     54.47 |   0.02028 |
#> |    M|   Mixed Diff. |-3.795e+06 |-1.661e+06 |-1.758e+06 |
#> |   20|     3236491.4 |    -1.231 |     1.316 |     1.032 |
#> |    U|     3236491.4 |     2.847 |     65.80 |   0.02016 |
#> |    X|     3236491.4 |     2.847 |     65.80 |   0.02016 |
#> |    M|   Mixed Diff. |-7.476e+05 |-4.760e+05 |-1.588e+06 |
#> |   21|     3375170.1 |    -1.133 |     1.379 |     1.240 |
#> |    U|     3375170.1 |     3.338 |     68.93 |   0.02120 |
#> |    X|     3375170.1 |     3.338 |     68.93 |   0.02120 |
#> |   22|     3150224.1 |    -1.194 |     1.340 |     1.111 |
#> |    U|     3150224.1 |     3.033 |     66.99 |   0.02055 |
#> |    X|     3150224.1 |     3.033 |     66.99 |   0.02055 |
#> |    M|   Mixed Diff. | 3.903e+06 |-1.324e+06 |-1.515e+06 |
#> |   23|     3023016.2 |    -1.215 |     1.413 |     1.159 |
#> |    U|     3023016.2 |     2.926 |     70.66 |   0.02079 |
#> |    X|     3023016.2 |     2.926 |     70.66 |   0.02079 |
#> |    M|   Mixed Diff. |-9.123e+05 | 2.721e+05 |-1.436e+06 |
#> |   24|     3000246.8 |    -1.145 |     1.532 |     1.276 |
#> |    U|     3000246.8 |     3.277 |     76.60 |   0.02138 |
#> |    X|     3000246.8 |     3.277 |     76.60 |   0.02138 |
#> |    M|   Mixed Diff. | 4.214e+06 | 1.159e+05 |-1.387e+06 |
#> |   25|     2706408.3 |    -1.206 |     1.411 |     1.396 |
#> |    U|     2706408.3 |     2.971 |     70.53 |   0.02198 |
#> |    X|     2706408.3 |     2.971 |     70.53 |   0.02198 |
#> |    M|   Mixed Diff. | 5.628e+04 | 7.816e+04 |-1.215e+06 |
#> |   26|     2322580.8 |    -1.223 |     1.387 |     1.756 |
#> |    U|     2322580.8 |     2.887 |     69.37 |   0.02378 |
#> |    X|     2322580.8 |     2.887 |     69.37 |   0.02378 |
#> |   27|     1649814.2 |    -1.271 |     1.321 |     2.789 |
#> |    U|     1649814.2 |     2.648 |     66.05 |   0.02894 |
#> |    X|     1649814.2 |     2.648 |     66.05 |   0.02894 |
#> |    M|   Mixed Diff. |-2.486e+06 | 1.283e+05 |-5.554e+05 |
#> |   28| 2.3670848e+29 |    -1.233 |   0.08498 |     3.440 |
#> |    U| 2.3670848e+29 |     2.839 |     4.249 |   0.03220 |
#> |    X| 2.3670848e+29 |     2.839 |     4.249 |   0.03220 |
#> |   29|     1713126.7 |    -1.199 |     1.221 |     2.855 |
#> |    U|     1713126.7 |     3.008 |     61.06 |   0.02927 |
#> |    X|     1713126.7 |     3.008 |     61.06 |   0.02927 |
#> |   30|     1583758.2 |    -1.217 |     1.318 |     2.801 |
#> |    U|     1583758.2 |     2.915 |     65.91 |   0.02900 |
#> |    X|     1583758.2 |     2.915 |     65.91 |   0.02900 |
#> |    M|   Mixed Diff. | 5.091e+05 |-4.111e+05 |-5.319e+05 |
#> |   31|     1565278.3 |    -1.230 |     1.365 |     2.827 |
#> |    U|     1565278.3 |     2.850 |     68.23 |   0.02913 |
#> |    X|     1565278.3 |     2.850 |     68.23 |   0.02913 |
#> |    M|   Mixed Diff. |-8.564e+05 | 3.677e+04 |-5.231e+05 |
#> |   32|     1525940.2 |    -1.191 |     1.436 |     2.899 |
#> |    U|     1525940.2 |     3.049 |     71.82 |   0.02950 |
#> |    X|     1525940.2 |     3.049 |     71.82 |   0.02950 |
#> |    M|   Mixed Diff. | 6.319e+05 | 2.047e+04 |-5.036e+05 |
#> |   33|     1470157.6 |    -1.210 |     1.401 |     3.001 |
#> |    U|     1470157.6 |     2.950 |     70.04 |   0.03000 |
#> |    X|     1470157.6 |     2.950 |     70.04 |   0.03000 |
#> |   34|     1384039.9 |    -1.208 |     1.288 |     3.277 |
#> |    U|     1384039.9 |     2.962 |     64.40 |   0.03139 |
#> |    X|     1384039.9 |     2.962 |     64.40 |   0.03139 |
#> |    M|   Mixed Diff. | 1.597e+06 |-7.525e+05 |-4.276e+05 |
#> |   35|     1209842.6 |    -1.234 |     1.291 |     3.683 |
#> |    U|     1209842.6 |     2.830 |     64.56 |   0.03341 |
#> |    X|     1209842.6 |     2.830 |     64.56 |   0.03341 |
#> |   36|     1031572.3 |    -1.249 |     1.284 |     4.290 |
#> |    U|     1031572.3 |     2.756 |     64.21 |   0.03645 |
#> |    X|     1031572.3 |     2.756 |     64.21 |   0.03645 |
#> |    M|   Mixed Diff. |-7.101e+05 |-1.302e+05 |-2.707e+05 |
#> |   37|     799114.69 |    -1.239 |     1.302 |     5.303 |
#> |    U|     799114.69 |     2.808 |     65.12 |   0.04152 |
#> |    X|     799114.69 |     2.808 |     65.12 |   0.04152 |
#> |    M|   Mixed Diff. |-3.514e+05 |-1.064e+05 |-1.813e+05 |
#> |   38|     526437.70 |    -1.226 |     1.332 |     7.361 |
#> |    U|      526437.7 |     2.872 |     66.60 |   0.05180 |
#> |    X|      526437.7 |     2.872 |     66.60 |   0.05180 |
#> |    M|   Mixed Diff. |-9.739e+04 |-6.044e+04 |-9.211e+04 |
#> |   39|     377147.49 |    -1.217 |     1.356 |     9.482 |
#> |    U|     377147.49 |     2.918 |     67.81 |   0.06241 |
#> |    X|     377147.49 |     2.918 |     67.81 |   0.06241 |
#> |    M|   Mixed Diff. |-1.745e+04 |-3.311e+04 |-5.215e+04 |
#> |   40|     270007.16 |    -1.208 |     1.382 |     12.24 |
#> |    U|     270007.16 |     2.963 |     69.11 |   0.07620 |
#> |    X|     270007.16 |     2.963 |     69.11 |   0.07620 |
#> |    M|   Mixed Diff. | 1.555e+04 |-1.490e+04 |-2.832e+04 |
#> |   41|     201117.97 |    -1.199 |     1.408 |     15.49 |
#> |    U|     201117.97 |     3.006 |     70.40 |   0.09247 |
#> |    X|     201117.97 |     3.006 |     70.40 |   0.09247 |
#> |    M|   Mixed Diff. | 2.471e+04 |    -4757. |-1.562e+04 |
#> |   42|     155198.10 |    -1.190 |     1.436 |     19.46 |
#> |    U|      155198.1 |     3.051 |     71.78 |    0.1123 |
#> |    X|      155198.1 |     3.051 |     71.78 |    0.1123 |
#> |    M|   Mixed Diff. | 2.435e+04 |     770.8 |    -8547. |
#> |   43|     127246.75 |    -1.182 |     1.463 |     23.79 |
#> |    U|     127246.75 |     3.094 |     73.14 |    0.1339 |
#> |    X|     127246.75 |     3.094 |     73.14 |    0.1339 |
#> |    M|   Mixed Diff. | 2.069e+04 |     3320. |    -4912. |
#> |   44|     106714.22 |    -1.171 |     1.496 |     29.55 |
#> |    U|     106714.22 |     3.148 |     74.82 |    0.1628 |
#> |    X|     106714.22 |     3.148 |     74.82 |    0.1628 |
#> |    M|   Mixed Diff. | 1.557e+04 |     4547. |    -2629. |
#> |   45|     94284.677 |    -1.159 |     1.532 |     36.06 |
#> |    U|     94284.677 |     3.206 |     76.61 |    0.1953 |
#> |    X|     94284.677 |     3.206 |     76.61 |    0.1953 |
#> |    M|   Mixed Diff. | 1.083e+04 |     4827. |    -1438. |
#> |   46|     86446.913 |    -1.146 |     1.573 |     43.70 |
#> |    U|     86446.913 |     3.271 |     78.64 |    0.2335 |
#> |    X|     86446.913 |     3.271 |     78.64 |    0.2335 |
#> |    F| Forward Diff. |     6660. |     4657. |    -772.9 |
#> |   47|     81819.577 |    -1.132 |     1.617 |     52.27 |
#> |    U|     81819.577 |     3.342 |     80.85 |    0.2764 |
#> |    X|     81819.577 |     3.342 |     80.85 |    0.2764 |
#> |    F| Forward Diff. |     3268. |     4288. |    -411.9 |
#> |   48|     79279.580 |    -1.117 |     1.663 |     61.29 |
#> |    U|      79279.58 |     3.417 |     83.13 |    0.3215 |
#> |    X|      79279.58 |     3.417 |     83.13 |    0.3215 |
#> |    F| Forward Diff. |     690.6 |     3876. |    -219.3 |
#> |   49|     77976.485 |    -1.102 |     1.708 |     70.32 |
#> |    U|     77976.485 |     3.490 |     85.38 |    0.3666 |
#> |    X|     77976.485 |     3.490 |     85.38 |    0.3666 |
#> |    F| Forward Diff. |    -1203. |     3496. |    -115.1 |
#> |   50|     77350.191 |    -1.088 |     1.752 |     79.34 |
#> |    U|     77350.191 |     3.564 |     87.59 |    0.4117 |
#> |    X|     77350.191 |     3.564 |     87.59 |    0.4117 |
#> |    F| Forward Diff. |    -2627. |     3159. |    -55.21 |
#> |   51|     77140.054 |    -1.076 |     1.786 |     86.42 |
#> |    U|     77140.054 |     3.621 |     89.31 |    0.4471 |
#> |    X|     77140.054 |     3.621 |     89.31 |    0.4471 |
#> |    F| Forward Diff. |    -3505. |     2922. |    -25.49 |
#> |   52|     77091.677 |    -1.069 |     1.807 |     90.85 |
#> |    U|     77091.677 |     3.657 |     90.36 |    0.4692 |
#> |    X|     77091.677 |     3.657 |     90.36 |    0.4692 |
#> |    F| Forward Diff. |    -3967. |     2782. |    -11.86 |
#> |   53|     77085.259 |    -1.066 |     1.815 |     92.67 |
#> |    U|     77085.259 |     3.671 |     90.77 |    0.4783 |
#> |    X|     77085.259 |     3.671 |     90.77 |    0.4783 |
#> |    F| Forward Diff. |    -4138. |     2721. |    -7.067 |
#> |   54|     77085.259 |    -1.066 |     1.815 |     92.67 |
#> |    U|     77085.259 |     3.671 |     90.77 |    0.4783 |
#> |    X|     77085.259 |     3.671 |     90.77 |    0.4783 |
#> calculating covariance matrix
#> done
#> Calculating residuals/tables
#> done
#> Warning: gradient problems with initial estimate and covariance; see $scaleInfo

# }