search

nlmixrdevelopment / nlmixr

nlmixr: an R package for population PKPD modeling
https://nlmixrdevelopment.github.io/nlmixr/
GNU General Public License v2.0
115 stars 45 forks source link

Error using plot.dyn.ID #546

Open gcaligari opened 3 years ago

gcaligari commented 3 years ago

Hello,

Thanks for your package.

While using the package I noticed the function plot.dyn.ID does not seem to plot anything and produces the error: Error in for (s in m) {: invalid for() loop sequence. I created the dyn.ID object using the example code from the dynmodel help file. An error for a missing environment was also displayed even though the argument evTable is no longer listed in the help for the function plot.dyn.ID. Hence the argument was included in the example below.

nlmixr version = 2.0.4 R version = 4.0.4

The reproducible code can be found below.

Thanks, Giulia 

ode <- "
      kel = CL/V;
      d/dt(X) = -kel*X;
      C=X/V;
      PRED = C
      "
ode_system <- RxODE::RxODE(model = ode)
model_error_structure <- cp ~ C + add(0.01) + prop(0.01)
inits <- c(CL = 1, V = 10)
control <- nlmixr::dynmodelControl(method = "Nelder-Mead")
ev <- RxODE::eventTable()
ev$add.sampling(c(0, nlmixr::Bolus_1CPT$TIME))
fit <-
  try(nlmixr::dynmodel(
    system = ode_system,
    model = model_error_structure,
    data = nlmixr::Bolus_1CPT,
    inits = inits,
    control = control,
    evTable = ev
  ))
#> |-----+---------------+-----------+-----------+-----------+-----------|
#> |    #| 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 |
#>                                                                                 
#> Warning in sqrt(diag(solve(fit$hessian))): NaNs produced
plot(fit)
#> Error in for (s in m) {: invalid for() loop sequence

Created on 2021-07-19 by the reprex package (v2.0.0)

mattfidler commented 3 years ago

Hi @gcaligari

Thanks for reporting this.

Currently the preferred way to run nlmixr for "dynmodel" is to use 

nlmixr(nlmixrModelFunction, data, est = "dynmodel")

You may also use "focei" and it will run another population only nonlinear regression.