-
Notifications
You must be signed in to change notification settings - Fork 2
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Feature request: Identify the data row causing integration issues #219
Comments
Here is another reprex that may help with getting integration issues. I've attached the required file: nlds_trans_sd_pk_chr.csv library(nlmixr2)
#> Loading required package: nlmixr2data
library(nlmixr2lib)
library(tidyverse)
trans_sd_err_chr <- read.csv("../datasets/nlds_trans_sd_pk_chr.csv")
d_mod <-
trans_sd_err_chr %>%
filter(CMPD == "SM") |>
filter(SPEC == "Monkey")
one.cmt.ode <- function() {
ini({
tka <- .5
tcl <- -3.2
tv <- -1
eta.ka ~ 1
eta.cl ~ 2
eta.v ~ 1
add.err <- 0.1
})
model({
ka <- exp(tka + eta.ka)
cl <- exp(tcl + eta.cl)
v <- exp(tv + eta.v)
d/dt(depot) = -ka * depot
d/dt(central) = ka * depot - cl/v * central
cp = central/v
cp ~ add(add.err)
})
}
fit <-
nlmixr2(
one.cmt.ode,
data = d_mod,
est = "focei",
control = foceiControl(print = 10),
table = tableControl(keep = "BW")
)
#> rxode2 2.1.2 using 8 threads (see ?getRxThreads)
#> no cache: create with `rxCreateCache()`
#> 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 | tka | tcl | tv | add.err |
#> |.....................| o1 | o2 | o3 |...........|
#> |-----+---------------+-----------+-----------+-----------+-----------|
#> | 10| 1128701.9 | -0.2605 | -1.132 | 0.7888 | 26.03 |
#> |.....................| 1.329 | 2.453 | 1.969 |...........|
#> | U| 1128701.9 | -0.5224 | -3.332 | -0.2588 | 1.373 |
#> |.....................| 1.329 | 2.659 | 1.969 |...........|
#> | X| 1128701.9 | 0.5931 | 0.03572 | 0.7720 | 1.373 |
#> |.....................| 1.329 | 2.659 | 1.969 |...........|
#> | M| Mixed Diff. | 229.6 | 169.6 | 123.8 |-4.120e+04 |
#> |.....................| -9159. | -5103. | -8285. |...........|
#> |-----+---------------+-----------+-----------+-----------+-----------|
#> | 20| 158567.34 | 15.27 | -15.43 | 30.59 | 140.3 |
#> |.....................| 3.162 | 0.6900 | 0.6687 |...........|
#> | U| 158567.34 | 15.01 | -17.63 | 29.54 | 7.088 |
#> |.....................| 3.162 | 0.5623 | 0.6687 |...........|
#> | X| 158567.34 | 3.309e+06 | 2.203e-08 | 6.741e+12 | 7.088 |
#> |.....................| 3.162 | 0.5623 | 0.6687 |...........|
#> | M| Mixed Diff. | 4518. | 2032. | 6698. | 268.6 |
#> |.....................| 6.443e+04 | 3.883e+04 | 1.825e+05 |...........|
#> |-----+---------------+-----------+-----------+-----------+-----------|
#> | 30| 21880.156 | -3.871 | -31.17 | -18.06 | 234.5 |
#> |.....................| 0.6687 | 0.6900 | 0.6687 |...........|
#> | U| 21880.156 | -4.133 | -33.37 | -19.10 | 11.80 |
#> |.....................| 0.6687 | 0.5623 | 0.6687 |...........|
#> | X| 21880.156 | 0.01604 | 3.231e-15 | 5.048e-09 | 11.80 |
#> |.....................| 0.6687 | 0.5623 | 0.6687 |...........|
#> Theta reset (zero/bad gradient values); Switch to bobyqa
#> | 40| 2.1196987e+08 | 0.5619 | -1.000 | 0.04762 | 0.5714 |
#> |.....................| 1.000 | 0.9973 | 1.000 |...........|
#> | U| 2.1196987e+08 | 0.3000 | -3.200 | -1.000 | 0.1000 |
#> |.....................| 1.000 | 0.9278 | 1.000 |...........|
#> | X| 2.1196987e+08 | 1.350 | 0.04076 | 0.3679 | 0.1000 |
#> |.....................| 1.000 | 0.9278 | 1.000 |...........|
#> | 50| 83685704. | 0.7545 | -1.001 | 0.3902 | 1.754 |
#> |.....................| 1.082 | 1.424 | 0.6687 |...........|
#> | U| 83685704 | 0.4926 | -3.201 | -0.6574 | 0.1592 |
#> |.....................| 1.082 | 1.436 | 0.6687 |...........|
#> | X| 83685704 | 1.637 | 0.04074 | 0.5182 | 0.1592 |
#> |.....................| 1.082 | 1.436 | 0.6687 |...........|
#> | 60| 76954437. | 0.8477 | -0.9145 | -0.02850 | 1.891 |
#> |.....................| 1.172 | 1.292 | 0.6687 |...........|
#> | U| 76954437 | 0.5858 | -3.115 | -1.076 | 0.1660 |
#> |.....................| 1.172 | 1.278 | 0.6687 |...........|
#> | X| 76954437 | 1.796 | 0.04440 | 0.3409 | 0.1660 |
#> |.....................| 1.172 | 1.278 | 0.6687 |...........|
#> | 70| 53501118. | 0.4267 | -0.8903 | -0.09208 | 2.552 |
#> |.....................| 1.140 | 1.639 | 0.9031 |...........|
#> | U| 53501118 | 0.1648 | -3.090 | -1.140 | 0.1990 |
#> |.....................| 1.140 | 1.691 | 0.9031 |...........|
#> | X| 53501118 | 1.179 | 0.04549 | 0.3199 | 0.1990 |
#> |.....................| 1.140 | 1.691 | 0.9031 |...........|
#> | 80| 43152133. | 0.7668 | -0.8650 | -0.1221 | 3.004 |
#> |.....................| 1.013 | 1.793 | 1.349 |...........|
#> | U| 43152133 | 0.5049 | -3.065 | -1.170 | 0.2216 |
#> |.....................| 1.013 | 1.874 | 1.349 |...........|
#> | X| 43152133 | 1.657 | 0.04665 | 0.3105 | 0.2216 |
#> |.....................| 1.013 | 1.874 | 1.349 |...........|
#> | 90| 36202130. | 0.4390 | -0.5840 | -0.05015 | 3.411 |
#> |.....................| 1.138 | 1.482 | 1.633 |...........|
#> | U| 36202130 | 0.1771 | -2.784 | -1.098 | 0.2420 |
#> |.....................| 1.138 | 1.505 | 1.633 |...........|
#> | X| 36202130 | 1.194 | 0.06179 | 0.3336 | 0.2420 |
#> |.....................| 1.138 | 1.505 | 1.633 |...........|
#> | 100| 23143646. | 1.050 | -0.6733 | -0.2346 | 4.625 |
#> |.....................| 0.8718 | 1.765 | 2.564 |...........|
#> | U| 23143646 | 0.7877 | -2.873 | -1.282 | 0.3027 |
#> |.....................| 0.8718 | 1.840 | 2.564 |...........|
#> | X| 23143646 | 2.198 | 0.05651 | 0.2774 | 0.3027 |
#> |.....................| 0.8718 | 1.840 | 2.564 |...........|
#> | 110| 13036310. | 0.8172 | 1.623 | -1.273 | 6.641 |
#> |.....................| 1.171 | 1.609 | 2.662 |...........|
#> | U| 13036310 | 0.5553 | -0.5774 | -2.321 | 0.4035 |
#> |.....................| 1.171 | 1.655 | 2.662 |...........|
#> | X| 13036310 | 1.742 | 0.5614 | 0.09820 | 0.4035 |
#> |.....................| 1.171 | 1.655 | 2.662 |...........|
#> | 120| 12163780. | 1.245 | 1.438 | -0.8737 | 6.924 |
#> |.....................| 0.9442 | 1.663 | 2.521 |...........|
#> | U| 12163780 | 0.9833 | -0.7617 | -1.921 | 0.4176 |
#> |.....................| 0.9442 | 1.719 | 2.521 |...........|
#> | X| 12163780 | 2.673 | 0.4669 | 0.1464 | 0.4176 |
#> |.....................| 0.9442 | 1.719 | 2.521 |...........|
#> | 130| 10589363. | 1.178 | 1.635 | -0.8318 | 7.525 |
#> |.....................| 0.7773 | 1.620 | 2.649 |...........|
#> | U| 10589363 | 0.9156 | -0.5652 | -1.879 | 0.4477 |
#> |.....................| 0.7773 | 1.668 | 2.649 |...........|
#> | X| 10589363 | 2.498 | 0.5682 | 0.1527 | 0.4477 |
#> |.....................| 0.7773 | 1.668 | 2.649 |...........|
#> | 140| 9341667.7 | 1.063 | 1.669 | -0.9767 | 8.107 |
#> |.....................| 1.056 | 1.488 | 2.845 |...........|
#> | U| 9341667.7 | 0.8006 | -0.5314 | -2.024 | 0.4768 |
#> |.....................| 1.056 | 1.512 | 2.845 |...........|
#> | X| 9341667.7 | 2.227 | 0.5878 | 0.1321 | 0.4768 |
#> |.....................| 1.056 | 1.512 | 2.845 |...........|
#> | 150| 7707929.4 | 0.7560 | 1.207 | -1.167 | 9.077 |
#> |.....................| 1.763 | 1.025 | 3.081 |...........|
#> | U| 7707929.4 | 0.4941 | -0.9930 | -2.215 | 0.5253 |
#> |.....................| 1.763 | 0.9603 | 3.081 |...........|
#> | X| 7707929.4 | 1.639 | 0.3705 | 0.1092 | 0.5253 |
#> |.....................| 1.763 | 0.9603 | 3.081 |...........|
#> | 160| 6389900.4 | 0.6964 | 1.373 | -1.122 | 10.11 |
#> |.....................| 1.816 | 0.7973 | 2.866 |...........|
#> | U| 6389900.4 | 0.4345 | -0.8266 | -2.169 | 0.5768 |
#> |.....................| 1.816 | 0.6900 | 2.866 |...........|
#> | X| 6389900.4 | 1.544 | 0.4375 | 0.1142 | 0.5768 |
#> |.....................| 1.816 | 0.6900 | 2.866 |...........|
#> | 170| 2494006.0 | 2.325 | 1.347 | -0.3678 | 17.03 |
#> |.....................| 0.6687 | 0.7973 | 2.051 |...........|
#> | U| 2494006 | 2.063 | -0.8534 | -1.415 | 0.9227 |
#> |.....................| 0.6687 | 0.6900 | 2.051 |...........|
#> | X| 2494006 | 7.867 | 0.4260 | 0.2428 | 0.9227 |
#> |.....................| 0.6687 | 0.6900 | 2.051 |...........|
#> | 180| 1764636.4 | -0.2012 | 2.867 | -1.487 | 20.51 |
#> |.....................| 2.229 | 1.001 | 1.632 |...........|
#> | U| 1764636.4 | -0.4631 | 0.6672 | -2.535 | 1.097 |
#> |.....................| 2.229 | 0.9321 | 1.632 |...........|
#> | X| 1764636.4 | 0.6293 | 1.949 | 0.07928 | 1.097 |
#> |.....................| 2.229 | 0.9321 | 1.632 |...........|
#> | 190| 1248876.5 | -1.032 | 3.998 | -1.431 | 24.64 |
#> |.....................| 3.162 | 2.280 | 0.6687 |...........|
#> | U| 1248876.5 | -1.294 | 1.798 | -2.479 | 1.303 |
#> |.....................| 3.162 | 2.453 | 0.6687 |...........|
#> | X| 1248876.5 | 0.2741 | 6.037 | 0.08384 | 1.303 |
#> |.....................| 3.162 | 2.453 | 0.6687 |...........|
#> | 200| 982766.06 | 0.09229 | 8.013 | -5.963 | 27.99 |
#> |.....................| 0.7925 | 1.918 | 0.7041 |...........|
#> | U| 982766.06 | -0.1696 | 5.813 | -7.010 | 1.471 |
#> |.....................| 0.7925 | 2.022 | 0.7041 |...........|
#> | X| 982766.06 | 0.8440 | 334.5 | 0.0009027 | 1.471 |
#> |.....................| 0.7925 | 2.022 | 0.7041 |...........|
#> | 210| 449041.29 | -5.444 | 21.51 | -9.569 | 42.67 |
#> |.....................| 1.951 | 2.280 | 0.6687 |...........|
#> | U| 449041.29 | -5.705 | 19.31 | -10.62 | 2.205 |
#> |.....................| 1.951 | 2.453 | 0.6687 |...........|
#> | X| 449041.29 | 0.003328 | 2.443e+08 | 2.451e-05 | 2.205 |
#> |.....................| 1.951 | 2.453 | 0.6687 |...........|
#> | 220| 267287.85 | -23.53 | 47.71 | -2.544 | 56.03 |
#> |.....................| 0.6687 | 2.280 | 2.017 |...........|
#> | U| 267287.85 | -23.79 | 45.51 | -3.592 | 2.873 |
#> |.....................| 0.6687 | 2.453 | 2.017 |...........|
#> | X| 267287.85 | 4.657e-11 | 5.841e+19 | 0.02754 | 2.873 |
#> |.....................| 0.6687 | 2.453 | 2.017 |...........|
#> | 230| 218167.37 | -22.17 | 45.99 | 1.487 | 61.49 |
#> |.....................| 2.184 | 2.280 | 1.762 |...........|
#> | U| 218167.37 | -22.43 | 43.79 | 0.4395 | 3.146 |
#> |.....................| 2.184 | 2.453 | 1.762 |...........|
#> | X| 218167.37 | 1.811e-10 | 1.039e+19 | 1.552 | 3.146 |
#> |.....................| 2.184 | 2.453 | 1.762 |...........|
#> | 240| 178557.98 | -23.54 | 47.38 | 1.395 | 68.32 |
#> |.....................| 1.090 | 1.938 | 1.877 |...........|
#> | U| 178557.98 | -23.80 | 45.18 | 0.3478 | 3.487 |
#> |.....................| 1.090 | 2.047 | 1.877 |...........|
#> | X| 178557.98 | 4.606e-11 | 4.198e+19 | 1.416 | 3.487 |
#> |.....................| 1.090 | 2.047 | 1.877 |...........|
#> | 250| 141491.18 | -31.37 | 46.08 | 1.718 | 78.04 |
#> |.....................| 1.451 | 1.309 | 2.237 |...........|
#> | U| 141491.18 | -31.63 | 43.88 | 0.6706 | 3.973 |
#> |.....................| 1.451 | 1.298 | 2.237 |...........|
#> | X| 141491.18 | 1.836e-14 | 1.135e+19 | 1.955 | 3.973 |
#> |.....................| 1.451 | 1.298 | 2.237 |...........|
#> | 260| 69926.666 | -58.79 | 33.26 | 3.934 | 111.1 |
#> |.....................| 3.162 | 0.7973 | 0.6687 |...........|
#> | U| 69926.666 | -59.05 | 31.06 | 2.886 | 5.626 |
#> |.....................| 3.162 | 0.6900 | 0.6687 |...........|
#> | X| 69926.666 | 2.257e-26 | 3.093e+13 | 17.92 | 5.626 |
#> |.....................| 3.162 | 0.6900 | 0.6687 |...........|
#> | 270| 62316.580 | -55.71 | 26.53 | -0.3657 | 118.0 |
#> |.....................| 2.720 | 0.7973 | 0.9972 |...........|
#> | U| 62316.58 | -55.97 | 24.33 | -1.413 | 5.972 |
#> |.....................| 2.720 | 0.6900 | 0.9972 |...........|
#> | X| 62316.58 | 4.904e-25 | 3.697e+10 | 0.2433 | 5.972 |
#> |.....................| 2.720 | 0.6900 | 0.9972 |...........|
#> | 280| 61870.951 | -55.90 | 26.64 | 0.01046 | 118.4 |
#> |.....................| 2.773 | 0.9519 | 0.7602 |...........|
#> | U| 61870.951 | -56.16 | 24.44 | -1.037 | 5.989 |
#> |.....................| 2.773 | 0.8738 | 0.7602 |...........|
#> | X| 61870.951 | 4.063e-25 | 4.103e+10 | 0.3545 | 5.989 |
#> |.....................| 2.773 | 0.8738 | 0.7602 |...........|
#> | 290| 61138.208 | -55.88 | 26.24 | 0.4106 | 119.1 |
#> |.....................| 2.400 | 0.8713 | 0.6687 |...........|
#> | U| 61138.208 | -56.14 | 24.04 | -0.6370 | 6.027 |
#> |.....................| 2.400 | 0.7780 | 0.6687 |...........|
#> | X| 61138.208 | 4.147e-25 | 2.764e+10 | 0.5289 | 6.027 |
#> |.....................| 2.400 | 0.7780 | 0.6687 |...........|
#> | 300| 59705.452 | -56.33 | 25.97 | 1.015 | 120.6 |
#> |.....................| 2.352 | 0.8115 | 0.6687 |...........|
#> | U| 59705.452 | -56.59 | 23.77 | -0.03259 | 6.103 |
#> |.....................| 2.352 | 0.7069 | 0.6687 |...........|
#> | X| 59705.452 | 2.650e-25 | 2.096e+10 | 0.9679 | 6.103 |
#> |.....................| 2.352 | 0.7069 | 0.6687 |...........|
#> | 310| 57009.078 | -56.61 | 23.16 | 0.1715 | 123.6 |
#> |.....................| 1.168 | 1.451 | 0.6687 |...........|
#> | U| 57009.078 | -56.88 | 20.96 | -0.8761 | 6.253 |
#> |.....................| 1.168 | 1.467 | 0.6687 |...........|
#> | X| 57009.078 | 1.993e-25 | 1.262e+09 | 0.4164 | 6.253 |
#> |.....................| 1.168 | 1.467 | 0.6687 |...........|
#> | 320| 53014.222 | -58.99 | 25.25 | 0.8209 | 128.6 |
#> |.....................| 3.162 | 0.7973 | 0.6687 |...........|
#> | U| 53014.222 | -59.25 | 23.05 | -0.2267 | 6.501 |
#> |.....................| 3.162 | 0.6900 | 0.6687 |...........|
#> | X| 53014.222 | 1.850e-26 | 1.021e+10 | 0.7972 | 6.501 |
#> |.....................| 3.162 | 0.6900 | 0.6687 |...........|
#> | 330| 28703.464 | -74.30 | 27.02 | 7.037 | 185.6 |
#> |.....................| 3.162 | 0.7973 | 0.6687 |...........|
#> | U| 28703.464 | -74.56 | 24.82 | 5.989 | 9.353 |
#> |.....................| 3.162 | 0.6900 | 0.6687 |...........|
#> | X| 28703.464 | 4.145e-33 | 6.039e+10 | 399.0 | 9.353 |
#> |.....................| 3.162 | 0.6900 | 0.6687 |...........|
#> | 340| 33672.783 | -73.58 | 24.33 | 7.012 | 186.2 |
#> |.....................| 1.951 | 0.7973 | 0.9682 |...........|
#> | U| 33672.783 | -73.84 | 22.13 | 5.964 | 9.381 |
#> |.....................| 1.951 | 0.6900 | 0.9682 |...........|
#> | X| 33672.783 | 8.511e-33 | 4.101e+09 | 389.2 | 9.381 |
#> |.....................| 1.951 | 0.6900 | 0.9682 |...........|
#> | 350| 28325.976 | -73.33 | 24.39 | 6.709 | 186.0 |
#> |.....................| 1.926 | 0.8045 | 0.6796 |...........|
#> | U| 28325.976 | -73.59 | 22.19 | 5.662 | 9.373 |
#> |.....................| 1.926 | 0.6985 | 0.6796 |...........|
#> | X| 28325.976 | 1.098e-32 | 4.315e+09 | 287.6 | 9.373 |
#> |.....................| 1.926 | 0.6985 | 0.6796 |...........|
#> | 360| 28305.241 | -73.30 | 24.37 | 6.710 | 186.1 |
#> |.....................| 1.899 | 0.7973 | 0.6782 |...........|
#> | U| 28305.241 | -73.57 | 22.17 | 5.663 | 9.374 |
#> |.....................| 1.899 | 0.6900 | 0.6782 |...........|
#> | X| 28305.241 | 1.124e-32 | 4.240e+09 | 287.9 | 9.374 |
#> |.....................| 1.899 | 0.6900 | 0.6782 |...........|
#> | 370| 28414.845 | -73.31 | 24.37 | 6.711 | 186.0 |
#> |.....................| 1.915 | 0.7973 | 0.6785 |...........|
#> | U| 28414.845 | -73.57 | 22.17 | 5.664 | 9.374 |
#> |.....................| 1.915 | 0.6900 | 0.6785 |...........|
#> | X| 28414.845 | 1.116e-32 | 4.263e+09 | 288.3 | 9.374 |
#> |.....................| 1.915 | 0.6900 | 0.6785 |...........|
#> | 380| 28554.872 | -73.31 | 24.37 | 6.714 | 186.0 |
#> |.....................| 1.915 | 0.7973 | 0.6786 |...........|
#> | U| 28554.872 | -73.57 | 22.17 | 5.666 | 9.374 |
#> |.....................| 1.915 | 0.6900 | 0.6786 |...........|
#> | X| 28554.872 | 1.117e-32 | 4.260e+09 | 288.9 | 9.374 |
#> |.....................| 1.915 | 0.6900 | 0.6786 |...........|
#> | 390| 28271.957 | -73.31 | 24.37 | 6.713 | 186.0 |
#> |.....................| 1.914 | 0.7973 | 0.6791 |...........|
#> | U| 28271.957 | -73.57 | 22.17 | 5.666 | 9.374 |
#> |.....................| 1.914 | 0.6900 | 0.6791 |...........|
#> | X| 28271.957 | 1.117e-32 | 4.257e+09 | 288.8 | 9.374 |
#> |.....................| 1.914 | 0.6900 | 0.6791 |...........|
#> | 400| 28273.180 | -73.31 | 24.37 | 6.713 | 186.0 |
#> |.....................| 1.914 | 0.7974 | 0.6790 |...........|
#> | U| 28273.18 | -73.57 | 22.17 | 5.666 | 9.374 |
#> |.....................| 1.914 | 0.6901 | 0.6790 |...........|
#> | X| 28273.18 | 1.117e-32 | 4.258e+09 | 288.8 | 9.374 |
#> |.....................| 1.914 | 0.6901 | 0.6790 |...........|
#> calculating covariance matrix
#> done
#> → Calculating residuals/tables
#> ✔ done
#> → compress origData in nlmixr2 object, save 72992
#> → compress parHistData in nlmixr2 object, save 32792
fit
#> ── nlmixr² FOCEi (outer: bobyqa) ──
#>
#> OBJF AIC BIC Log-likelihood Condition#(Cov) Condition#(Cor)
#> FOCEi 28271.73 29388.45 29419.23 -14687.23 165141.4 1
#>
#> ── Time (sec fit$time): ──
#>
#> setup optimize covariance table other
#> elapsed 0.019 0.198 0.198 0.05 6.265
#>
#> ── Population Parameters (fit$parFixed or fit$parFixedDf): ──
#>
#> Est. SE %RSE Back-transformed(95%CI) BSV(CV%) Shrink(SD)%
#> tka -73.6 175 238 1.12e-32 (1.46e-181, 8.55e+116) 27.8 1.74%
#> tcl 22.2 175 789 4.26e+09 (5.57e-140, 3.25e+158) 901. 24.3%
#> tv 5.67 0.43 7.6 289 (124, 671) 1050. -47.5%
#> add.err 9.37 9.37
#>
#> Covariance Type (fit$covMethod): s+
#> No correlations in between subject variability (BSV) matrix
#> Full BSV covariance (fit$omega) or correlation (fit$omegaR; diagonals=SDs)
#> Distribution stats (mean/skewness/kurtosis/p-value) available in fit$shrink
#> Information about run found (fit$runInfo):
#> • tolerances (atol/rtol) were increased (after 4 bad solves) for some difficult ODE solving during the optimization. can control with foceiControl(stickyRecalcN=) consider increasing sigdig/atol/rtol changing initial estimates or changing the structural model
#> • gradient problems with initial estimate and covariance; see $scaleInfo
#> • using S matrix to calculate covariance, can check sandwich or R matrix with $covRS and $covR
#> • S matrix had problems solving for some subject and parameters
#> • bad solve during optimization
#> • last objective function was not at minimum, possible problems in optimization
#> • ETAs were reset to zero during optimization; (Can control by foceiControl(resetEtaP=.))
#> • initial ETAs were nudged; (can control by foceiControl(etaNudge=., etaNudge2=))
#> • Hessian reset during optimization; (can control by foceiControl(resetHessianAndEta=.))
#> • thetas were reset during optimization because of a zero gradient
#> Censoring (fit$censInformation): No censoring
#> Minimization message (fit$message):
#> Normal exit from bobyqa
#>
#> ── Fit Data (object fit is a modified tibble): ──
#> # A tibble: 600 × 23
#> ID TIME DV PRED RES WRES IPRED IRES IWRES CPRED CRES
#> <fct> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
#> 1 37 0.083 125. -3.28e-19 125. 13.3 3.48e-12 125. 13.3 -4.61e-13 125.
#> 2 37 0.5 31.5 3.67e-40 31.5 3.36 2.48e-41 31.5 3.36 2.29e-41 31.5
#> 3 37 1 56.6 2.36e-41 56.6 6.04 2.48e-41 56.6 6.04 2.36e-41 56.6
#> # ℹ 597 more rows
#> # ℹ 12 more variables: CWRES <dbl>, eta.ka <dbl>, eta.cl <dbl>, eta.v <dbl>,
#> # depot <dbl>, central <dbl>, ka <dbl>, cl <dbl>, v <dbl>, tad <dbl>,
#> # dosenum <dbl>, BW <dbl> Created on 2024-04-23 with reprex v2.1.0 |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
When getting an error like the following (and others like it), it would help if the error message identified the specific subject, time point, and row which caused the issue.
Ideally, the subject identification would be from the original subject identifier and not from a translated identifier (in other words, I know that identifiers are numeric-only in rxode2, but if they come in as character or factor variables (or maybe doubles), then, having them show as the source version would be preferred.
Here are some data and a model with the issue to give a use case.
It should work like the following (not quite a reprex as I'm traveling now):
anon_data.RDS.txt
anon_model.RDS.txt
The text was updated successfully, but these errors were encountered: