= @model begin
warfarin_pkpd_model ...
@dynamics begin
' = -Ka * Depot
Depot' = Ka * Depot - CL / Vc * Central
Central' =
Turnover* (1 + emax * (Central / Vc) / (cp50 + Central / Vc)) - kout * Turnover
rin end
...
end
Differential Equations in Pumas
1 Introduction
Pumas automatically chooses a differential equation solver that is suitable for the simulation or estimation of the dynamical system of the NLME (Nonlinear Mixed Effects) model at hand. This default solver is the preferred choice and optimized for most users and use cases. Nevertheless, in some cases the performance-accuracy trade-off can be improved by adjusting the tolerances, or possibly even the algorithm, of the differential equation solver.
In this tutorial, the Warfarin PK/PD model is used to demonstrate how to configure the differential equation solver.
2 Learning Goals
- Observe the utility of the
@vars
block of a Pumas model with respect to storing dynamic variables associated with differential equations - Understand the main differences between common differential equation solvers for nonlinear dynamical systems
- Learn how to adjust the algorithm and the tolerances of the differential equation solver
3 Warfarin PK/PD Model
We return to the Warfarin PK/PD model. Its dynamical system consists of three states, \(\operatorname{Depot}\), \(\operatorname{Central}\), and \(\operatorname{Turnover}\), whose dynamics are governed by the ordinary differential equations:
\[ \begin{aligned} \operatorname{Depot}'(t) &= - \operatorname{Ka} \operatorname{Depot}(t),\\ \operatorname{Central}'(t) &= \operatorname{Ka} \operatorname{Depot}(t) - \frac{\operatorname{CL}}{\operatorname{Vc}} \operatorname{Central}(t),\\ \operatorname{Turnover}'(t) &= \operatorname{rin} (1 + \operatorname{emax} \frac{\operatorname{Central}(t) / \operatorname{Vc}}{\operatorname{c50} + \operatorname{Central}(t)/\operatorname{Vc}}) - \operatorname{kout} \operatorname{Turnover}(t) \end{aligned} \]
with PK parameters \(\operatorname{Ka}\) (absorption rate), \(\operatorname{CL}\) (clearance), and \(\operatorname{Vc}\) (volume of distribution) and PD parameters \(\operatorname{rin}\), \(\operatorname{emax}\), \(\operatorname{c50}\), and \(\operatorname{kout}\).
The dynamical system can be written more concisely by introducing auxiliary variables for repeated expressions:
\[ \begin{aligned} \operatorname{Depot}'(t) &= -\operatorname{ratein}(t),\\ \operatorname{Central}'(t) &= \operatorname{ratein}(t) - \operatorname{CL} \operatorname{cp}(t),\\ \operatorname{Turnover}'(t) &= \operatorname{rin} \operatorname{pd}(t) - \operatorname{kout} \operatorname{Turnover}(t) \end{aligned} \]
with influx rate \(\operatorname{ratein}(t) := \operatorname{Ka} \operatorname{Depot}(t)\), concentration \(\operatorname{cp}(t) := \operatorname{Central}(t) / \operatorname{Vc}\), and \(\operatorname{pd}(t) := 1 + \operatorname{emax} \frac{\operatorname{cp}(t)}{\operatorname{c50} + \operatorname{cp}(t)}\).
4 Auxiliary Variables in @vars
In Pumas, dynamical systems are defined in the @dynamics
block inside of the @model
definition. For instance, the dynamical system of the Warfarin PK/PD model can be implemented as follows:
The same concise rewriting can be applied in a Pumas @model
by defining auxiliary variables (“aliases”) in the @vars
block:
= @model begin
warfarin_pkpd_model ...
@vars begin
:= Central / Vc
cp := Ka * Depot
ratein := 1 + emax * cp / (c50 + cp)
pd end
@dynamics begin
' = -ratein
Depot' = ratein - CL * cp
Central' = rin * pd - kout * Turnover
Turnoverend
...
end
The walrus operator (:=
) ensures that the aliases do not show up in the simulation output of the model. However, if you would like to access an alias in the simulation output, you should define the alias with =
. For instance, if you want to obtain concentration cp
as part of the simulation output, you can change the @vars
block to
@vars begin
= Central / Vc
cp := Ka * Depot
ratein := 1 + emax * cp / (c50 + cp)
pd end
5 Differential Equation Solvers
The differential equation in the Warfarin model is non-linear, as detected by Pumas (“Dynamical system type: Nonlinear ODE”):
using Pumas
= @model begin
warfarin_pkpd_model @param begin
# PK parameters
"""
Clearance (L/h/70kg)
"""
∈ RealDomain(lower = 0.0, init = 0.134)
pop_CL """
Central Volume L/70kg
"""
∈ RealDomain(lower = 0.0, init = 8.11)
pop_V """
Absorption time (h)
"""
∈ RealDomain(lower = 0.0, init = 0.523)
pop_tabs """
Lag time (h)
"""
∈ RealDomain(lower = 0.0, init = 0.1)
pop_lag # PD parameters
"""
Baseline
"""
∈ RealDomain(lower = 0.0, init = 100.0)
pop_e0 """
Emax
"""
∈ RealDomain(init = -1.0)
pop_emax """
EC50
"""
∈ RealDomain(lower = 0.0, init = 1.0)
pop_c50 """
Turnover
"""
∈ RealDomain(lower = 0.0, init = 14.0)
pop_tover # Inter-individual variability
"""
- ΩCL
- ΩVc
- ΩTabs
"""
∈ PDiagDomain([0.01, 0.01, 0.01])
pk_Ω """
- Ωe0
- Ωemax
- Ωec50
- Ωturn
"""
∈ PDiagDomain([0.01, 0.01, 0.01, 0.01])
pd_Ω # Residual variability
"""
Proportional residual error for drug concentration
"""
∈ RealDomain(lower = 0.0, init = 0.00752)
σ_prop """
Additive residual error for drug concentration (mg/L)
"""
∈ RealDomain(lower = 0.0, init = 0.0661)
σ_add """
Additive error for PCA
"""
∈ RealDomain(lower = 0.0, init = 0.01)
σ_fx end
@random begin
# mean = 0, covariance = pk_Ω
~ MvNormal(pk_Ω)
pk_η # mean = 0, covariance = pd_Ω
~ MvNormal(pd_Ω)
pd_η end
@covariates FSZV FSZCL
@pre begin
# PK
= FSZCL * pop_CL * exp(pk_η[1])
CL = FSZV * pop_V * exp(pk_η[2])
Vc = pop_tabs * exp(pk_η[3])
tabs = log(2) / tabs
Ka # PD
= pop_e0 * exp(pd_η[1])
e0 = pop_emax * exp(pd_η[2])
emax = pop_c50 * exp(pd_η[3])
c50 = pop_tover * exp(pd_η[4])
tover = log(2) / tover
kout = e0 * kout
rin = t
time end
@dosecontrol begin
= (Depot = pop_lag,)
lags end
@init begin
= e0
Turnover end
# aliases for use in @dynamics and @derived
@vars begin
:= Central / Vc
cp := Ka * Depot
ratein := 1 + emax * cp / (c50 + cp)
pd end
@dynamics begin
' = -ratein
Depot' = ratein - CL * cp
Central' = rin * pd - kout * Turnover
Turnoverend
@derived begin
"""
Warfarin Concentration (mg/L)
"""
~ @. Normal(cp, sqrt((σ_prop * cp)^2 + σ_add^2))
conc """
PCA
"""
~ @. Normal(Turnover, σ_fx)
pca end
end
PumasModel
Parameters: pop_CL, pop_V, pop_tabs, pop_lag, pop_e0, pop_emax, pop_c50, pop_tover, pk_Ω, pd_Ω, σ_prop, σ_add, σ_fx
Random effects: pk_η, pd_η
Covariates: FSZV, FSZCL
Dynamical system variables: Depot, Central, Turnover
Dynamical system type: Nonlinear ODE
Derived: conc, pca
Observed: conc, pca
Pumas approximates the solution of the differential equation with a numerical differential equation solver. Generally, one distinguishes between solvers for stiff and non-stiff differential equations.
5.1 Stiff vs. Non-Stiff Systems
A key distinction among numerical solvers is whether they are designed for stiff or non-stiff differential equations:
Non-Stiff Differential Equations: These systems exhibit relatively moderate changes in their variables. Standard non-stiff solvers can efficiently approximate solutions of these systems.
Stiff Differential Equations: These systems contain rapidly changing components alongside more slowly varying dynamics. Non-stiff solvers typically perform poorly on stiff systems, as they may require exceedingly small step sizes to maintain numerical stability. Specialized stiff solvers are therefore employed to handle the sharp gradients and large timescale differences without compromising accuracy.
6 Pumas’s Automatic Solver Selection
By default, Pumas adopts a hybrid approach with automatic stiffness detection to switch between stiff and non-stiff solvers as needed.
6.1 Model Simulation
- Default Solvers: Rosenbrock23 (stiff) and Tsit5 (non-stiff)
- Tolerances: Relative tolerance \(1 \times 10^{-3}\) and absolute tolerance \(1 \times 10^{-6}\)
- Rationale: These higher (less stringent) tolerances permit faster simulations while maintaining sufficient accuracy for exploratory and predictive modeling.
6.2 Model Fitting
- Default Solvers: Rodas5P (stiff) and Vern7 (non-stiff)
- Tolerances: Relative tolerance \(10^{-8}\) and absolute tolerance \(10^{-12}\)
- Rationale: These lower (more stringent) tolerances ensure high precision during parameter estimation, which is critical for matching the model’s predictions to observed data.
The default solvers and tolerances are recommended for most users in most instances. If desired, however, it is possible to adjust these settings with the diffeq_options
keyword argument.
Computation time decreases as tolerances are increased. However, higher tolerances come at the cost of a less strict error control, and hence generally a less accurate solution.
6.3 Adjusting the Tolerances
The absolute and relative tolerance of the solver can be specified with abstol
and reltol
.
7 Absolute and Relative Tolerances
When employing a numerical solver, it is necessary to specify how accurately the solution should be computed. This precision is controlled by two key parameters:
Absolute Tolerance \((\text{abstol})\)
- Interpreted as the maximum allowable error when the solution values are near zero.
- Ensures that numerical approximations stay within a reasonable bound, preventing physically impossible outcomes (e.g., negative concentrations) or excessive drift at small scales.
- For instance, an absolute tolerance of \(10^{-6}\) means the solver attempts to keep the absolute error below \(10^{-6}\) whenever the solution magnitude is close to zero.
Relative Tolerance \((\text{reltol})\)
- Enforces the number of correct digits throughout the simulation, effectively controlling error relative to the current scale of the solution.
- For example, a relative tolerance of \(10^{-3}\) implies the solver aims for three correct decimal places (i.e., the solution is accurate to within 0.1% of its current magnitude).
- As the solution grows or shrinks, the solver adjusts its time-step size and internal computations to maintain this relative accuracy.
Sometimes decreasing tolerances can help to reduce numerical problems, e.g. to keep solutions non-negative that are mathematically guaranteed to be non-negative. Additionally, the choice of tolerances can be motivated by the application of the numerical solution: For plotting a less accurate solution, and hence larger tolerances, might be tolerable, whereas typically for model fitting a more accurate solution, and hence smaller tolerances, are beneficial.
This can be demonstrated when fitting the Warfarin model with an example dataset: Optimization fails with large tolerances of 1e-3 (relative) and 1e-6 (absolute)
fit(
warfarin_pkpd_model,
pop,init_params(warfarin_pkpd_model),
FOCE();
= (; reltol = 1e-3, abstol = 1e-6),
diffeq_options )
[ Info: Checking the initial parameter values.
[ Info: The initial negative log likelihood and its gradient are finite. Check passed.
Iter Function value Gradient norm
0 3.004020e+06 5.745853e+06
* time: 0.036778926849365234
1 4.816669e+05 8.316173e+05
* time: 2.5890910625457764
2 3.598373e+05 6.060150e+05
* time: 4.041985988616943
3 1.682124e+05 2.709072e+05
* time: 5.025707960128784
4 9.168914e+04 1.446832e+05
* time: 6.04728102684021
5 4.791138e+04 6.671463e+04
* time: 6.980186939239502
6 2.909104e+04 3.426663e+04
* time: 7.950915098190308
7 1.843737e+04 1.678953e+04
* time: 8.913994073867798
8 1.282447e+04 1.021463e+04
* time: 9.825294971466064
9 9.723361e+03 9.291650e+03
* time: 10.753484010696411
10 7.635613e+03 8.325311e+03
* time: 11.570909023284912
11 6.243129e+03 7.337397e+03
* time: 12.43933892250061
12 5.279957e+03 6.304542e+03
* time: 13.235424041748047
13 4.422123e+03 5.001962e+03
* time: 14.043703079223633
14 3.744748e+03 3.550400e+03
* time: 14.780565023422241
15 3.409674e+03 2.461547e+03
* time: 15.548466920852661
16 3.312049e+03 1.856430e+03
* time: 16.269973039627075
17 3.300559e+03 1.648439e+03
* time: 17.003747940063477
18 3.299489e+03 1.599460e+03
* time: 17.72153091430664
19 3.298436e+03 1.567520e+03
* time: 18.457744121551514
20 3.294921e+03 1.501363e+03
* time: 19.170675039291382
21 3.286641e+03 1.407664e+03
* time: 19.913458108901978
22 3.264637e+03 1.259077e+03
* time: 20.630484104156494
23 3.209594e+03 1.043534e+03
* time: 21.32778000831604
24 3.072117e+03 7.454677e+02
* time: 22.049246072769165
25 2.739794e+03 3.818063e+02
* time: 22.678839921951294
26 1.949883e+03 2.246644e+02
* time: 23.424750089645386
27 1.722434e+03 2.224562e+02
* time: 25.341300010681152
28 1.479355e+03 2.042452e+02
* time: 27.792783975601196
29 1.310583e+03 1.531868e+02
* time: 30.72933006286621
30 1.224675e+03 1.756826e+02
* time: 31.43550205230713
31 1.204841e+03 2.738440e+02
* time: 31.999047994613647
32 1.196432e+03 2.403217e+02
* time: 32.6150860786438
33 1.189131e+03 2.301186e+02
* time: 33.17582893371582
34 1.177881e+03 2.083030e+02
* time: 33.80488991737366
35 1.176989e+03 2.201993e+02
* time: 34.38284206390381
36 1.176899e+03 2.230826e+02
* time: 35.01331090927124
37 1.176882e+03 2.233531e+02
* time: 35.57376194000244
38 1.176795e+03 2.238567e+02
* time: 36.187248945236206
39 1.176609e+03 2.236914e+02
* time: 36.76012706756592
40 1.176090e+03 2.210941e+02
* time: 37.37300491333008
41 1.174841e+03 2.112741e+02
* time: 37.949657917022705
42 1.171973e+03 1.825298e+02
* time: 38.557745933532715
43 1.166737e+03 1.210559e+02
* time: 39.13270711898804
44 1.160698e+03 4.450697e+01
* time: 39.746416091918945
45 1.157926e+03 6.411275e+01
* time: 40.31949710845947
46 1.157537e+03 6.121762e+01
* time: 40.92166495323181
47 1.157520e+03 5.867223e+01
* time: 41.48825192451477
48 1.157519e+03 5.814523e+01
* time: 42.06284809112549
49 1.157515e+03 5.686007e+01
* time: 42.61605906486511
50 1.157505e+03 5.497418e+01
* time: 43.20294690132141
51 1.157477e+03 5.129724e+01
* time: 43.77531695365906
52 1.157408e+03 4.463171e+01
* time: 44.34547996520996
53 1.157239e+03 3.279117e+01
* time: 44.89983606338501
54 1.156869e+03 3.681526e+01
* time: 45.47150707244873
55 1.156230e+03 3.198754e+01
* time: 46.036349058151245
56 1.155563e+03 4.576438e+01
* time: 46.61383891105652
57 1.155264e+03 4.832159e+01
* time: 47.16887903213501
58 1.155212e+03 4.052175e+01
* time: 47.7317750453949
59 1.155207e+03 3.797125e+01
* time: 48.277997970581055
60 1.155206e+03 3.778452e+01
* time: 48.83329105377197
61 1.155201e+03 3.723003e+01
* time: 49.383882999420166
62 1.155190e+03 3.642375e+01
* time: 49.94686007499695
63 1.155159e+03 3.487535e+01
* time: 50.49422812461853
64 1.155079e+03 3.209060e+01
* time: 51.052669048309326
65 1.154875e+03 2.670227e+01
* time: 51.645885944366455
66 1.154375e+03 2.366165e+01
* time: 52.18662095069885
67 1.153278e+03 5.499525e+01
* time: 52.78028607368469
68 1.151457e+03 7.903833e+01
* time: 53.30179810523987
69 1.149737e+03 6.771227e+01
* time: 53.874305963516235
70 1.148983e+03 5.387708e+01
* time: 54.38840293884277
71 1.148872e+03 4.683950e+01
* time: 54.95218896865845
72 1.148868e+03 4.512036e+01
* time: 55.458534955978394
73 1.148866e+03 4.489706e+01
* time: 56.01118612289429
74 1.148857e+03 4.426314e+01
* time: 56.51303195953369
75 1.148838e+03 4.321161e+01
* time: 57.05819892883301
76 1.148785e+03 4.095232e+01
* time: 57.55684995651245
77 1.148652e+03 3.780572e+01
* time: 58.104042053222656
78 1.148325e+03 3.474239e+01
* time: 58.60068106651306
79 1.147599e+03 3.678957e+01
* time: 59.152535915374756
80 1.146330e+03 5.134277e+01
* time: 59.66378998756409
81 1.144979e+03 6.444063e+01
* time: 60.2073540687561
82 1.144325e+03 7.574751e+01
* time: 60.712568044662476
83 1.144210e+03 6.968296e+01
* time: 61.25182795524597
84 1.144204e+03 6.648234e+01
* time: 61.75245213508606
85 1.144202e+03 6.560999e+01
* time: 62.28234505653381
86 1.144194e+03 6.305667e+01
* time: 62.772743940353394
87 1.144176e+03 5.956959e+01
* time: 63.307044982910156
88 1.144127e+03 5.307828e+01
* time: 63.811546087265015
89 1.144005e+03 4.218776e+01
* time: 64.3437271118164
90 1.143708e+03 4.120261e+01
* time: 64.86069798469543
91 1.143074e+03 3.675516e+01
* time: 65.39756393432617
92 1.142052e+03 3.855439e+01
* time: 65.91438698768616
93 1.141101e+03 5.193197e+01
* time: 66.44793796539307
94 1.140701e+03 4.606858e+01
* time: 66.9569981098175
95 1.140646e+03 4.857736e+01
* time: 67.49195313453674
96 1.140643e+03 4.817596e+01
* time: 68.00139594078064
97 1.140642e+03 4.793197e+01
* time: 68.52508807182312
98 1.140637e+03 4.723245e+01
* time: 69.03045296669006
99 1.140626e+03 4.619695e+01
* time: 69.55674290657043
100 1.140597e+03 4.418761e+01
* time: 70.0640001296997
101 1.140523e+03 4.050986e+01
* time: 70.58925604820251
102 1.140334e+03 3.335334e+01
* time: 71.10701894760132
103 1.139892e+03 2.982120e+01
* time: 71.63840007781982
104 1.138994e+03 4.360237e+01
* time: 72.15890002250671
105 1.137687e+03 5.246892e+01
* time: 72.6933479309082
106 1.136679e+03 4.825586e+01
* time: 73.21530795097351
107 1.136391e+03 4.467499e+01
* time: 73.74992108345032
108 1.136364e+03 3.961591e+01
* time: 74.26277899742126
109 1.136362e+03 3.951402e+01
* time: 74.78474307060242
110 1.136360e+03 3.943358e+01
* time: 75.32679295539856
111 1.136355e+03 3.923965e+01
* time: 75.80602192878723
112 1.136344e+03 3.894674e+01
* time: 76.3524649143219
113 1.136314e+03 3.841518e+01
* time: 76.83977103233337
114 1.136235e+03 3.748824e+01
* time: 77.39240407943726
115 1.136029e+03 3.577210e+01
* time: 77.884202003479
116 1.135499e+03 3.252332e+01
* time: 78.42985892295837
117 1.134182e+03 2.630844e+01
* time: 78.94658613204956
118 1.131238e+03 4.869163e+01
* time: 79.58505606651306
119 1.126345e+03 6.537856e+01
* time: 80.15365791320801
120 1.121879e+03 5.123373e+01
* time: 80.78502798080444
121 1.119645e+03 2.457316e+01
* time: 81.33928608894348
122 1.119031e+03 2.124721e+01
* time: 81.93651103973389
123 1.118958e+03 1.857508e+01
* time: 82.49152398109436
124 1.118953e+03 1.896595e+01
* time: 83.08413195610046
125 1.118952e+03 1.849021e+01
* time: 83.62071204185486
126 1.118952e+03 1.852941e+01
* time: 84.21458005905151
127 1.118951e+03 1.856812e+01
* time: 84.78116011619568
128 1.118951e+03 1.862403e+01
* time: 85.3782889842987
129 1.118949e+03 1.871484e+01
* time: 85.92496800422668
130 1.118946e+03 1.886332e+01
* time: 86.5058970451355
131 1.118939e+03 1.910217e+01
* time: 87.04735898971558
132 1.118922e+03 1.948050e+01
* time: 87.61697292327881
133 1.118879e+03 2.006903e+01
* time: 88.15827798843384
134 1.118769e+03 2.097106e+01
* time: 88.73619198799133
135 1.118486e+03 2.231801e+01
* time: 89.28660297393799
136 1.117755e+03 3.049039e+01
* time: 89.877769947052
137 1.115883e+03 5.345487e+01
* time: 90.44304895401001
138 1.111258e+03 8.489274e+01
* time: 91.02758312225342
139 1.101149e+03 7.643749e+01
* time: 91.59094595909119
140 1.090612e+03 5.802659e+01
* time: 92.2185070514679
141 1.087971e+03 5.776554e+01
* time: 92.82703804969788
142 1.087205e+03 1.148234e+02
* time: 93.4706621170044
143 1.085558e+03 4.724872e+01
* time: 94.12859010696411
144 1.085271e+03 2.018848e+01
* time: 94.77615809440613
145 1.085124e+03 2.004201e+01
* time: 95.44430208206177
146 1.085089e+03 2.025616e+01
* time: 96.04406809806824
147 1.085081e+03 2.049682e+01
* time: 96.6955509185791
148 1.085073e+03 2.086075e+01
* time: 97.31198000907898
149 1.085070e+03 2.095706e+01
* time: 97.95375108718872
150 1.085063e+03 2.094978e+01
* time: 98.54541397094727
151 1.085062e+03 2.093855e+01
* time: 99.22301006317139
152 1.085062e+03 2.093439e+01
* time: 99.9742591381073
153 1.085062e+03 2.093335e+01
* time: 100.71732807159424
154 1.085062e+03 2.093252e+01
* time: 101.56392908096313
155 1.085062e+03 2.093227e+01
* time: 102.41550993919373
156 1.085062e+03 2.093208e+01
* time: 103.22258806228638
157 1.085062e+03 2.093187e+01
* time: 104.08704996109009
158 1.085062e+03 2.093181e+01
* time: 104.99790096282959
159 1.085062e+03 2.093175e+01
* time: 105.91559791564941
160 1.085062e+03 2.093172e+01
* time: 106.81795907020569
161 1.085062e+03 2.046653e+01
* time: 107.48654913902283
162 1.085061e+03 2.045920e+01
* time: 108.22480893135071
163 1.085056e+03 2.021786e+01
* time: 108.81096601486206
164 1.085055e+03 2.018381e+01
* time: 109.50923705101013
165 1.085043e+03 1.997498e+01
* time: 110.09617209434509
166 1.085029e+03 1.961695e+01
* time: 110.73259806632996
167 1.084997e+03 1.929634e+01
* time: 111.3284900188446
168 1.084881e+03 1.855033e+01
* time: 111.95556402206421
169 1.084633e+03 1.799616e+01
* time: 112.55737090110779
170 1.083831e+03 1.856893e+01
* time: 113.18822407722473
171 1.081240e+03 2.069096e+01
* time: 113.8402681350708
172 1.072465e+03 3.262436e+01
* time: 114.65537810325623
173 1.072446e+03 6.686587e+01
* time: 115.4879469871521
174 1.068493e+03 7.417042e+01
* time: 116.35348296165466
175 1.067469e+03 7.109936e+01
* time: 117.29630708694458
176 1.065854e+03 2.774005e+01
* time: 118.06751894950867
177 1.065215e+03 1.328875e+01
* time: 118.86322808265686
178 1.065128e+03 1.390400e+01
* time: 119.61015701293945
179 1.065120e+03 1.383172e+01
* time: 120.376699924469
180 1.065120e+03 1.384995e+01
* time: 121.12983798980713
FittedPumasModel
Successful minimization: false
Likelihood approximation: FOCE
Likelihood Optimizer: BFGS
Dynamical system type: Nonlinear ODE
Solver(s):(OrdinaryDiffEq.Vern7,OrdinaryDiffEq.Rodas5P)
Log-likelihood value: -1065.1201
Number of subjects: 31
Number of parameters: Fixed Optimized
0 18
Observation records: Active Missing
conc: 239 47
pca: 225 61
Total: 464 108
-------------------------
Estimate
-------------------------
pop_CL 0.13142
pop_V 7.9732
pop_tabs 0.57389
pop_lag 0.83565
pop_e0 96.508
pop_emax -1.0982
pop_c50 1.7031
pop_tover 14.804
pk_Ω₁,₁ 0.078665
pk_Ω₂,₂ 0.031456
pk_Ω₃,₃ 0.083202
pd_Ω₁,₁ 0.00316
pd_Ω₂,₂ 0.0023756
pd_Ω₃,₃ 0.025823
pd_Ω₄,₄ 0.013808
σ_prop 0.013696
σ_add 0.84173
σ_fx 3.3937
-------------------------
but succeeds with lower tolerances of 1e-8 (relative) and 1e-12 (absolute):
fit(
warfarin_pkpd_model,
pop,init_params(warfarin_pkpd_model),
FOCE();
= (; reltol = 1e-8, abstol = 1e-12),
diffeq_options )
[ Info: Checking the initial parameter values.
[ Info: The initial negative log likelihood and its gradient are finite. Check passed.
Iter Function value Gradient norm
0 2.998927e+06 5.740466e+06
* time: 0.00011706352233886719
1 4.804005e+05 8.280312e+05
* time: 1.5192348957061768
2 3.595691e+05 6.035843e+05
* time: 2.83803391456604
3 1.691071e+05 2.702705e+05
* time: 4.2914719581604
4 9.193152e+04 1.451040e+05
* time: 5.5921080112457275
5 4.768539e+04 6.630202e+04
* time: 6.934375047683716
6 2.904205e+04 3.416919e+04
* time: 8.278393983840942
7 1.840153e+04 1.671742e+04
* time: 9.582650899887085
8 1.280989e+04 1.022268e+04
* time: 10.996672868728638
9 9.713287e+03 9.298485e+03
* time: 12.30501103401184
10 7.628127e+03 8.330656e+03
* time: 13.64023208618164
11 6.236278e+03 7.340367e+03
* time: 14.880173921585083
12 5.271958e+03 6.303773e+03
* time: 16.160606861114502
13 4.411047e+03 4.994130e+03
* time: 17.385520935058594
14 3.732936e+03 3.538184e+03
* time: 18.640158891677856
15 3.399273e+03 2.450045e+03
* time: 19.945761919021606
16 3.302987e+03 1.848476e+03
* time: 21.230284929275513
17 3.291788e+03 1.642864e+03
* time: 22.455607891082764
18 3.290743e+03 1.594765e+03
* time: 23.669301986694336
19 3.289687e+03 1.562641e+03
* time: 24.95719003677368
20 3.286186e+03 1.496632e+03
* time: 26.265692949295044
21 3.277917e+03 1.402865e+03
* time: 27.66506791114807
22 3.255959e+03 1.254297e+03
* time: 28.98529291152954
23 3.200980e+03 1.038621e+03
* time: 30.30991506576538
24 3.063521e+03 7.403605e+02
* time: 31.462276935577393
25 2.730316e+03 3.768701e+02
* time: 32.64803194999695
26 1.934186e+03 2.246121e+02
* time: 33.711530923843384
27 1.704399e+03 2.219826e+02
* time: 38.09420990943909
28 1.471183e+03 2.025826e+02
* time: 45.78375005722046
29 1.316789e+03 1.554714e+02
* time: 52.24638295173645
30 1.226363e+03 1.775666e+02
* time: 53.75022101402283
31 1.206626e+03 2.780980e+02
* time: 55.33333396911621
32 1.197830e+03 2.430550e+02
* time: 59.57191705703735
33 1.190342e+03 2.327905e+02
* time: 62.08125901222229
34 1.178783e+03 2.105809e+02
* time: 64.50521302223206
35 1.177853e+03 2.228585e+02
* time: 66.45746088027954
36 1.177763e+03 2.256116e+02
* time: 68.1512839794159
37 1.177745e+03 2.259359e+02
* time: 69.66983103752136
38 1.177670e+03 2.264633e+02
* time: 71.21446800231934
39 1.177500e+03 2.264574e+02
* time: 72.701730966568
40 1.177033e+03 2.244136e+02
* time: 74.37466597557068
41 1.175894e+03 2.160226e+02
* time: 75.89445996284485
42 1.173240e+03 1.905912e+02
* time: 77.49400186538696
43 1.168189e+03 1.332783e+02
* time: 78.92913603782654
44 1.161835e+03 5.289952e+01
* time: 80.445631980896
45 1.158453e+03 6.595387e+01
* time: 81.81252598762512
46 1.157897e+03 6.411419e+01
* time: 83.19129300117493
47 1.157870e+03 6.098509e+01
* time: 84.42896509170532
48 1.157868e+03 6.032325e+01
* time: 85.76575899124146
49 1.157864e+03 5.920812e+01
* time: 86.92808103561401
50 1.157856e+03 5.740225e+01
* time: 88.0699999332428
51 1.157832e+03 5.406795e+01
* time: 89.14641809463501
52 1.157774e+03 4.797879e+01
* time: 90.25668597221375
53 1.157628e+03 3.658732e+01
* time: 91.34145092964172
54 1.157299e+03 3.707874e+01
* time: 92.45667600631714
55 1.156690e+03 3.452947e+01
* time: 93.57003903388977
56 1.155959e+03 4.269004e+01
* time: 94.67274188995361
57 1.155557e+03 5.000322e+01
* time: 95.79495692253113
58 1.155472e+03 4.219349e+01
* time: 96.80988001823425
59 1.155463e+03 3.895049e+01
* time: 97.91653108596802
60 1.155462e+03 3.872099e+01
* time: 99.04453897476196
61 1.155458e+03 3.823491e+01
* time: 100.3542640209198
62 1.155450e+03 3.747370e+01
* time: 101.41949796676636
63 1.155426e+03 3.608648e+01
* time: 102.50969004631042
64 1.155365e+03 3.361118e+01
* time: 103.50373101234436
65 1.155207e+03 2.894678e+01
* time: 104.5814778804779
66 1.154814e+03 2.504284e+01
* time: 105.58042001724243
67 1.153910e+03 4.541126e+01
* time: 106.65499997138977
68 1.152230e+03 7.426044e+01
* time: 107.64266586303711
69 1.150286e+03 7.473044e+01
* time: 108.68848896026611
70 1.149188e+03 5.605967e+01
* time: 109.65592694282532
71 1.148944e+03 4.801123e+01
* time: 110.65061688423157
72 1.148931e+03 4.659478e+01
* time: 111.6021020412445
73 1.148929e+03 4.635180e+01
* time: 112.60002589225769
74 1.148922e+03 4.568579e+01
* time: 113.56269907951355
75 1.148907e+03 4.466128e+01
* time: 114.59013986587524
76 1.148864e+03 4.249742e+01
* time: 115.56250286102295
77 1.148758e+03 3.824840e+01
* time: 116.58168697357178
78 1.148490e+03 3.552459e+01
* time: 117.53972887992859
79 1.147881e+03 3.038181e+01
* time: 118.55775094032288
80 1.146736e+03 4.467281e+01
* time: 119.50921106338501
81 1.145347e+03 5.929158e+01
* time: 120.49515104293823
82 1.144535e+03 7.572784e+01
* time: 121.46721291542053
83 1.144354e+03 7.034986e+01
* time: 122.4580090045929
84 1.144341e+03 6.591584e+01
* time: 123.40814089775085
85 1.144340e+03 6.496353e+01
* time: 124.41140294075012
86 1.144333e+03 6.268267e+01
* time: 125.37297487258911
87 1.144320e+03 5.950121e+01
* time: 126.34827303886414
88 1.144281e+03 5.368613e+01
* time: 127.30502200126648
89 1.144186e+03 4.393956e+01
* time: 128.32886385917664
90 1.143949e+03 4.102324e+01
* time: 129.29173707962036
91 1.143423e+03 3.752005e+01
* time: 130.29346108436584
92 1.142497e+03 3.278884e+01
* time: 131.2672619819641
93 1.141489e+03 5.217814e+01
* time: 132.26393795013428
94 1.140969e+03 4.310032e+01
* time: 133.22649693489075
95 1.140872e+03 4.759545e+01
* time: 134.21879506111145
96 1.140867e+03 4.727034e+01
* time: 135.17387199401855
97 1.140866e+03 4.704365e+01
* time: 136.14086389541626
98 1.140862e+03 4.637923e+01
* time: 137.09379386901855
99 1.140853e+03 4.542582e+01
* time: 138.0620620250702
100 1.140829e+03 4.357378e+01
* time: 139.016352891922
101 1.140767e+03 4.023060e+01
* time: 139.99046397209167
102 1.140610e+03 3.376323e+01
* time: 140.95349097251892
103 1.140237e+03 2.857028e+01
* time: 141.9597909450531
104 1.139451e+03 3.895699e+01
* time: 142.94640803337097
105 1.138223e+03 5.064036e+01
* time: 143.93478202819824
106 1.137140e+03 4.763314e+01
* time: 144.8974449634552
107 1.136763e+03 4.604968e+01
* time: 145.88860893249512
108 1.136720e+03 3.993479e+01
* time: 146.86068892478943
109 1.136717e+03 3.890565e+01
* time: 147.87167406082153
110 1.136716e+03 3.883645e+01
* time: 148.91256308555603
111 1.136712e+03 3.865517e+01
* time: 149.86452198028564
112 1.136702e+03 3.838833e+01
* time: 150.878112077713
113 1.136674e+03 3.789981e+01
* time: 151.81833386421204
114 1.136603e+03 3.705558e+01
* time: 152.83747386932373
115 1.136418e+03 3.549940e+01
* time: 153.7962989807129
116 1.135940e+03 3.257136e+01
* time: 154.90574288368225
117 1.134743e+03 2.697248e+01
* time: 155.8784420490265
118 1.132014e+03 4.863416e+01
* time: 156.92252802848816
119 1.127203e+03 6.901202e+01
* time: 157.87802386283875
120 1.122363e+03 5.794311e+01
* time: 158.88915586471558
121 1.119759e+03 2.906766e+01
* time: 159.84126496315002
122 1.118972e+03 2.137843e+01
* time: 160.84685492515564
123 1.118870e+03 1.909639e+01
* time: 161.80524706840515
124 1.118863e+03 1.887024e+01
* time: 162.8046269416809
125 1.118862e+03 1.852691e+01
* time: 163.75453400611877
126 1.118861e+03 1.864783e+01
* time: 164.77575588226318
127 1.118861e+03 1.862870e+01
* time: 165.75831508636475
128 1.118860e+03 1.866000e+01
* time: 166.7681589126587
129 1.118859e+03 1.869625e+01
* time: 167.68520188331604
130 1.118856e+03 1.878717e+01
* time: 168.68680906295776
131 1.118850e+03 1.894283e+01
* time: 169.65949606895447
132 1.118834e+03 1.922180e+01
* time: 170.68194007873535
133 1.118794e+03 1.968060e+01
* time: 171.69221997261047
134 1.118693e+03 2.042353e+01
* time: 172.80190587043762
135 1.118431e+03 2.156998e+01
* time: 173.87150692939758
136 1.117751e+03 2.596600e+01
* time: 174.96889090538025
137 1.116007e+03 4.716393e+01
* time: 176.03343200683594
138 1.111677e+03 7.744992e+01
* time: 177.14787602424622
139 1.102024e+03 7.496268e+01
* time: 178.2622299194336
140 1.090649e+03 6.246756e+01
* time: 179.45081901550293
141 1.087843e+03 8.941834e+01
* time: 180.59681606292725
142 1.087400e+03 1.188539e+02
* time: 181.76632595062256
143 1.085563e+03 6.082023e+01
* time: 182.83661198616028
144 1.085185e+03 1.947062e+01
* time: 184.1010708808899
145 1.085039e+03 1.958695e+01
* time: 185.38184595108032
146 1.084989e+03 1.989918e+01
* time: 186.84682393074036
147 1.084983e+03 2.009395e+01
* time: 191.64635586738586
148 1.084973e+03 2.065564e+01
* time: 200.11242699623108
149 1.084971e+03 2.080010e+01
* time: 202.01024389266968
150 1.084971e+03 2.084232e+01
* time: 203.72682404518127
151 1.084970e+03 2.084506e+01
* time: 205.44128894805908
152 1.084970e+03 2.084321e+01
* time: 207.03729009628296
153 1.084968e+03 2.084737e+01
* time: 208.7056610584259
154 1.084963e+03 2.088045e+01
* time: 210.17149090766907
155 1.084954e+03 2.098376e+01
* time: 211.607764005661
156 1.084935e+03 2.121619e+01
* time: 212.90796208381653
157 1.084904e+03 2.160158e+01
* time: 214.2358798980713
158 1.084864e+03 2.208328e+01
* time: 215.61956691741943
159 1.084807e+03 2.260004e+01
* time: 216.9590950012207
160 1.084700e+03 2.319870e+01
* time: 218.29315996170044
161 1.084457e+03 2.382124e+01
* time: 219.66228008270264
162 1.083882e+03 2.406095e+01
* time: 221.03799200057983
163 1.082589e+03 2.849219e+01
* time: 222.4134500026703
164 1.081660e+03 3.379635e+01
* time: 223.90836906433105
165 1.078576e+03 4.983497e+01
* time: 225.3506679534912
166 1.075004e+03 6.870527e+01
* time: 226.7394859790802
167 1.069377e+03 6.309176e+01
* time: 227.99089789390564
168 1.067978e+03 3.811944e+01
* time: 229.22918891906738
169 1.065604e+03 1.825143e+01
* time: 230.40376496315002
170 1.065251e+03 1.739631e+01
* time: 231.6250729560852
171 1.065088e+03 1.384781e+01
* time: 232.80305290222168
172 1.065086e+03 1.382386e+01
* time: 233.98204803466797
173 1.065086e+03 1.383866e+01
* time: 235.16159296035767
174 1.065086e+03 1.384608e+01
* time: 236.44933986663818
175 1.065086e+03 1.385678e+01
* time: 237.65214586257935
176 1.065084e+03 1.387273e+01
* time: 238.8369550704956
177 1.065081e+03 1.389858e+01
* time: 240.02831506729126
178 1.065073e+03 1.393897e+01
* time: 241.38331389427185
179 1.065052e+03 1.400288e+01
* time: 242.56688404083252
180 1.064997e+03 1.842052e+01
* time: 243.7803089618683
181 1.064853e+03 3.031555e+01
* time: 245.1414818763733
182 1.064477e+03 4.894991e+01
* time: 246.64757800102234
183 1.063505e+03 7.692533e+01
* time: 248.13937401771545
184 1.061069e+03 1.122192e+02
* time: 249.58867192268372
185 1.056257e+03 1.328367e+02
* time: 251.27166295051575
186 1.052013e+03 1.408303e+02
* time: 253.18860602378845
187 1.047713e+03 1.506293e+02
* time: 263.16789388656616
188 1.039359e+03 9.039176e+01
* time: 266.12340903282166
189 1.037334e+03 8.945366e+01
* time: 268.6719880104065
190 1.035766e+03 2.509488e+01
* time: 271.0166530609131
191 1.035510e+03 5.323250e+00
* time: 273.2789189815521
192 1.035401e+03 3.326584e+00
* time: 275.2801830768585
193 1.035202e+03 1.011867e+01
* time: 277.16745591163635
194 1.035130e+03 9.089057e+00
* time: 278.9494700431824
195 1.035055e+03 3.941605e+00
* time: 280.68632197380066
196 1.035048e+03 3.160942e+00
* time: 282.3372440338135
197 1.035046e+03 3.164670e+00
* time: 283.9616279602051
198 1.035046e+03 3.160486e+00
* time: 285.61218094825745
199 1.035046e+03 3.159428e+00
* time: 287.09520196914673
200 1.035046e+03 3.159138e+00
* time: 288.50277495384216
201 1.035046e+03 3.158158e+00
* time: 289.97307300567627
202 1.035046e+03 3.156875e+00
* time: 291.45791602134705
203 1.035046e+03 3.154403e+00
* time: 292.9741539955139
204 1.035045e+03 3.150100e+00
* time: 294.4060890674591
205 1.035043e+03 3.141854e+00
* time: 295.824914932251
206 1.035038e+03 3.125468e+00
* time: 297.27209186553955
207 1.035025e+03 3.887763e+00
* time: 298.5941939353943
208 1.034991e+03 6.349435e+00
* time: 299.88971495628357
209 1.034902e+03 1.025429e+01
* time: 301.0621929168701
210 1.034678e+03 1.614807e+01
* time: 302.25922894477844
211 1.034143e+03 2.389258e+01
* time: 303.4212498664856
212 1.033018e+03 3.065325e+01
* time: 304.75169587135315
213 1.031157e+03 2.988569e+01
* time: 305.9808940887451
214 1.029913e+03 2.000449e+01
* time: 307.4042270183563
215 1.029264e+03 1.806624e+01
* time: 308.8183898925781
216 1.028739e+03 1.876674e+01
* time: 310.2684009075165
217 1.028529e+03 1.788313e+01
* time: 311.79619097709656
218 1.028407e+03 1.869928e+01
* time: 313.29398488998413
219 1.028318e+03 2.103440e+01
* time: 314.76568388938904
220 1.028259e+03 2.278116e+01
* time: 316.195494890213
221 1.028193e+03 2.168720e+01
* time: 317.61125802993774
222 1.028177e+03 2.174690e+01
* time: 319.0970690250397
223 1.028175e+03 2.195554e+01
* time: 320.6408689022064
224 1.028175e+03 2.193442e+01
* time: 322.37642002105713
225 1.028174e+03 2.202037e+01
* time: 324.6540448665619
226 1.028174e+03 2.201903e+01
* time: 326.3532609939575
227 1.028173e+03 2.200004e+01
* time: 328.0789530277252
228 1.028171e+03 2.197102e+01
* time: 329.66717505455017
229 1.028164e+03 2.191369e+01
* time: 331.3471488952637
230 1.028148e+03 2.181626e+01
* time: 332.8448438644409
231 1.028106e+03 2.163664e+01
* time: 334.3984749317169
232 1.027999e+03 3.038649e+01
* time: 335.91544699668884
233 1.027721e+03 5.035972e+01
* time: 337.41603088378906
234 1.027021e+03 8.114702e+01
* time: 338.7693510055542
235 1.025404e+03 1.233505e+02
* time: 340.1757719516754
236 1.022343e+03 1.646448e+02
* time: 341.80963706970215
237 1.017556e+03 1.679173e+02
* time: 344.42520093917847
238 1.013261e+03 5.600035e+01
* time: 346.33879709243774
239 1.012953e+03 2.691504e+01
* time: 348.19176506996155
240 1.012758e+03 1.329779e+01
* time: 349.7655589580536
241 1.012704e+03 2.847118e+00
* time: 351.3634760379791
242 1.012699e+03 1.333414e+00
* time: 352.96707797050476
243 1.012697e+03 9.979337e-01
* time: 354.5908889770508
244 1.012697e+03 9.969545e-01
* time: 355.9790198802948
245 1.012697e+03 9.971341e-01
* time: 357.3837299346924
246 1.012697e+03 9.971972e-01
* time: 358.797394990921
247 1.012697e+03 9.972237e-01
* time: 360.2021849155426
248 1.012697e+03 9.972927e-01
* time: 361.5180449485779
249 1.012697e+03 9.973804e-01
* time: 362.78346586227417
250 1.012697e+03 9.975286e-01
* time: 364.24108386039734
251 1.012697e+03 9.977461e-01
* time: 365.9457528591156
252 1.012697e+03 9.980635e-01
* time: 367.9481599330902
253 1.012697e+03 9.984743e-01
* time: 369.49633288383484
254 1.012696e+03 9.988765e-01
* time: 370.99748396873474
255 1.012694e+03 9.988363e-01
* time: 372.4875409603119
256 1.012690e+03 1.474525e+00
* time: 373.90662002563477
257 1.012678e+03 2.326520e+00
* time: 375.34226298332214
258 1.012649e+03 3.510529e+00
* time: 376.7298250198364
259 1.012590e+03 4.737310e+00
* time: 378.16983103752136
260 1.012497e+03 4.897808e+00
* time: 379.5339090824127
261 1.012416e+03 2.875911e+00
* time: 380.86388993263245
262 1.012389e+03 9.301211e-01
* time: 382.17883491516113
263 1.012386e+03 8.804644e-01
* time: 383.4761519432068
264 1.012386e+03 8.547350e-01
* time: 384.7260539531708
265 1.012386e+03 8.516621e-01
* time: 385.8806688785553
266 1.012386e+03 8.511780e-01
* time: 386.99943590164185
267 1.012386e+03 8.489977e-01
* time: 388.10281205177307
268 1.012386e+03 8.463566e-01
* time: 389.2003860473633
269 1.012386e+03 8.415231e-01
* time: 390.3088130950928
270 1.012386e+03 8.340037e-01
* time: 391.40973806381226
271 1.012386e+03 8.215489e-01
* time: 392.5239210128784
272 1.012386e+03 8.013257e-01
* time: 393.6287319660187
273 1.012385e+03 7.680260e-01
* time: 394.74793100357056
274 1.012383e+03 7.407512e-01
* time: 395.8580870628357
275 1.012377e+03 7.469071e-01
* time: 396.9811329841614
276 1.012364e+03 1.103988e+00
* time: 398.10390400886536
277 1.012330e+03 1.710981e+00
* time: 399.2443630695343
278 1.012257e+03 2.412723e+00
* time: 400.4257559776306
279 1.012125e+03 2.705165e+00
* time: 401.75014901161194
280 1.011948e+03 1.883597e+00
* time: 403.2552399635315
281 1.011809e+03 8.215973e-01
* time: 404.56367588043213
282 1.011775e+03 9.067203e-01
* time: 405.92310404777527
283 1.011774e+03 9.565073e-01
* time: 407.2017879486084
284 1.011774e+03 9.628263e-01
* time: 408.5160689353943
285 1.011774e+03 9.654767e-01
* time: 409.72981905937195
286 1.011774e+03 9.654769e-01
* time: 411.0573809146881
287 1.011774e+03 9.656366e-01
* time: 412.14817690849304
288 1.011774e+03 9.656622e-01
* time: 413.2797749042511
289 1.011774e+03 9.656624e-01
* time: 414.52914690971375
290 1.011774e+03 9.656627e-01
* time: 415.9566168785095
291 1.011774e+03 9.660395e-01
* time: 417.411838054657
292 1.011774e+03 9.661257e-01
* time: 420.7361218929291
293 1.011774e+03 9.663115e-01
* time: 422.9264678955078
294 1.011774e+03 9.662851e-01
* time: 424.89053201675415
295 1.011773e+03 9.654773e-01
* time: 426.63135409355164
296 1.011772e+03 9.623124e-01
* time: 428.21404695510864
297 1.011769e+03 9.521889e-01
* time: 429.60455107688904
298 1.011763e+03 9.239616e-01
* time: 431.01576685905457
299 1.011746e+03 9.242954e-01
* time: 432.3588728904724
300 1.011711e+03 1.257052e+00
* time: 433.70727801322937
301 1.011651e+03 1.361278e+00
* time: 434.99998593330383
302 1.011591e+03 9.076916e-01
* time: 436.30766892433167
303 1.011568e+03 3.084752e-01
* time: 437.6138319969177
304 1.011565e+03 2.946404e-01
* time: 438.9141249656677
305 1.011565e+03 2.821273e-01
* time: 440.64502000808716
306 1.011565e+03 2.766260e-01
* time: 447.7602529525757
307 1.011565e+03 2.778010e-01
* time: 450.07408809661865
308 1.011565e+03 2.830281e-01
* time: 452.05622005462646
309 1.011565e+03 2.842619e-01
* time: 453.7743148803711
310 1.011565e+03 2.861093e-01
* time: 455.3576180934906
311 1.011565e+03 2.861138e-01
* time: 456.8638188838959
312 1.011565e+03 2.867641e-01
* time: 458.39902687072754
313 1.011565e+03 2.867756e-01
* time: 460.0616879463196
314 1.011565e+03 2.884589e-01
* time: 461.60845589637756
315 1.011565e+03 2.897335e-01
* time: 463.0596649646759
316 1.011565e+03 2.927744e-01
* time: 464.4553608894348
317 1.011565e+03 2.967424e-01
* time: 465.7813789844513
318 1.011565e+03 3.027347e-01
* time: 467.1630620956421
319 1.011564e+03 3.099827e-01
* time: 468.45212602615356
320 1.011563e+03 3.155826e-01
* time: 469.80787897109985
321 1.011560e+03 3.082844e-01
* time: 471.32473707199097
322 1.011555e+03 2.625121e-01
* time: 472.89328503608704
323 1.011548e+03 2.136271e-01
* time: 474.3493719100952
324 1.011542e+03 1.165033e-01
* time: 475.70230507850647
325 1.011538e+03 3.914224e-02
* time: 477.2863130569458
326 1.011537e+03 3.083888e-02
* time: 478.6780700683594
327 1.011536e+03 2.766364e-02
* time: 480.2767000198364
328 1.011536e+03 1.873262e-02
* time: 481.8227870464325
329 1.011536e+03 5.453951e-03
* time: 483.27135491371155
330 1.011536e+03 4.882951e-03
* time: 485.43862295150757
331 1.011535e+03 4.331879e-03
* time: 487.00595903396606
332 1.011535e+03 2.711268e-03
* time: 488.7913589477539
333 1.011535e+03 4.474256e-04
* time: 490.4663028717041
FittedPumasModel
Successful minimization: true
Likelihood approximation: FOCE
Likelihood Optimizer: BFGS
Dynamical system type: Nonlinear ODE
Solver(s):(OrdinaryDiffEq.Vern7,OrdinaryDiffEq.Rodas5P)
Log-likelihood value: -1011.5354
Number of subjects: 31
Number of parameters: Fixed Optimized
0 18
Observation records: Active Missing
conc: 239 47
pca: 225 61
Total: 464 108
-------------------------
Estimate
-------------------------
pop_CL 0.13137
pop_V 8.0255
pop_tabs 0.42635
pop_lag 0.92336
pop_e0 96.688
pop_emax -1.0628
pop_c50 1.5016
pop_tover 14.055
pk_Ω₁,₁ 0.050744
pk_Ω₂,₂ 0.021666
pk_Ω₃,₃ 0.97884
pd_Ω₁,₁ 0.0028828
pd_Ω₂,₂ 7.3726e-9
pd_Ω₃,₃ 0.15184
pd_Ω₄,₄ 0.011913
σ_prop 0.088292
σ_add 0.36362
σ_fx 3.5336
-------------------------
It is not recommended to decrease tolerances below 1e-14.
7.1 Changing the Algorithm
Usually, it should not be necessary to adjust the differential equation solver. If you change the solver, you should follow the guidelines in the SciML documentation that explains which solvers are the most efficient at the desired tolerance level.
For instance, if it is known that a differential equation is stiff, a stiff solver such as Rosenbrock23
at high tolerances or Rodas5P
at low tolerances could be a possible alternative to the default auto-switching solver:
# Fitting with stiff solver Rodas5P at low tolerances (relative: 1e-8, absolute: 1e-12)
fit(
warfarin_pkpd_model,
pop,init_params(warfarin_pkpd_model),
FOCE();
= (; alg = Rodas5P(), reltol = 1e-8, abstol = 1e-12),
diffeq_options )
[ Info: Checking the initial parameter values.
[ Info: The initial negative log likelihood and its gradient are finite. Check passed.
Iter Function value Gradient norm
0 2.998927e+06 5.740466e+06
* time: 0.0004229545593261719
1 4.804005e+05 8.280311e+05
* time: 20.43834400177002
2 3.595691e+05 6.035842e+05
* time: 29.471380949020386
3 1.691071e+05 2.702705e+05
* time: 37.48148512840271
4 9.193152e+04 1.451040e+05
* time: 55.84230709075928
5 4.768538e+04 6.630202e+04
* time: 66.06435704231262
6 2.904205e+04 3.416919e+04
* time: 73.87840008735657
7 1.840152e+04 1.671742e+04
* time: 81.54669404029846
8 1.280989e+04 1.022268e+04
* time: 95.6118631362915
9 9.713285e+03 9.298485e+03
* time: 104.83492994308472
10 7.628127e+03 8.330656e+03
* time: 116.6948390007019
11 6.236278e+03 7.340367e+03
* time: 125.74649715423584
12 5.271958e+03 6.303773e+03
* time: 133.56066608428955
13 4.411047e+03 4.994130e+03
* time: 143.04981112480164
14 3.732935e+03 3.538183e+03
* time: 155.91171598434448
15 3.399272e+03 2.450045e+03
* time: 165.12851905822754
16 3.302987e+03 1.848475e+03
* time: 172.9835181236267
17 3.291788e+03 1.642864e+03
* time: 180.44070601463318
18 3.290743e+03 1.594765e+03
* time: 190.4156060218811
19 3.289686e+03 1.562640e+03
* time: 201.36314702033997
20 3.286186e+03 1.496632e+03
* time: 209.5697090625763
21 3.277917e+03 1.402865e+03
* time: 222.28419995307922
22 3.255958e+03 1.254296e+03
* time: 231.74119114875793
23 3.200979e+03 1.038621e+03
* time: 239.39262199401855
24 3.063520e+03 7.403602e+02
* time: 251.87802696228027
25 2.730315e+03 3.768701e+02
* time: 262.83934712409973
26 1.934186e+03 2.246121e+02
* time: 271.126501083374
27 1.704399e+03 2.219826e+02
* time: 297.37584495544434
28 1.471184e+03 2.025825e+02
* time: 352.5652551651001
29 1.316789e+03 1.554715e+02
* time: 390.85428404808044
30 1.226363e+03 1.775721e+02
* time: 399.8102309703827
31 1.206626e+03 2.781030e+02
* time: 408.9569561481476
32 1.197830e+03 2.430603e+02
* time: 416.63732504844666
33 1.190342e+03 2.327966e+02
* time: 423.44104409217834
34 1.178784e+03 2.105878e+02
* time: 435.07282996177673
35 1.177853e+03 2.228653e+02
* time: 443.9634311199188
36 1.177764e+03 2.256182e+02
* time: 451.4774830341339
37 1.177746e+03 2.259425e+02
* time: 459.1620099544525
38 1.177670e+03 2.264699e+02
* time: 476.0657501220703
39 1.177500e+03 2.264641e+02
* time: 485.15774512290955
40 1.177034e+03 2.244207e+02
* time: 493.3200659751892
41 1.175895e+03 2.160307e+02
* time: 501.2916090488434
42 1.173241e+03 1.906017e+02
* time: 514.9100179672241
43 1.168190e+03 1.332917e+02
* time: 523.6566340923309
44 1.161836e+03 5.291022e+01
* time: 531.5969970226288
45 1.158453e+03 6.595471e+01
* time: 540.7448670864105
46 1.157897e+03 6.411640e+01
* time: 553.9287431240082
47 1.157870e+03 6.098679e+01
* time: 561.7589581012726
48 1.157868e+03 6.032481e+01
* time: 568.8653681278229
49 1.157864e+03 5.920991e+01
* time: 575.6833369731903
50 1.157856e+03 5.740417e+01
* time: 582.4387221336365
51 1.157832e+03 5.407031e+01
* time: 591.1384210586548
52 1.157774e+03 4.798188e+01
* time: 598.6969130039215
53 1.157628e+03 3.659183e+01
* time: 606.7312109470367
54 1.157299e+03 3.707862e+01
* time: 619.2901420593262
55 1.156690e+03 3.453183e+01
* time: 628.3015520572662
56 1.155959e+03 4.268473e+01
* time: 636.4888861179352
57 1.155557e+03 5.000309e+01
* time: 644.342481136322
58 1.155472e+03 4.219411e+01
* time: 664.7153830528259
59 1.155463e+03 3.894982e+01
* time: 672.343031167984
60 1.155462e+03 3.872028e+01
* time: 686.8097219467163
61 1.155458e+03 3.823435e+01
* time: 698.6388101577759
62 1.155450e+03 3.747328e+01
* time: 708.0100240707397
63 1.155426e+03 3.608641e+01
* time: 718.1292691230774
64 1.155365e+03 3.361168e+01
* time: 728.1509051322937
65 1.155207e+03 2.894844e+01
* time: 736.8533210754395
66 1.154814e+03 2.504332e+01
* time: 753.8598899841309
67 1.153910e+03 4.539953e+01
* time: 785.2044761180878
68 1.152231e+03 7.425240e+01
* time: 849.2698309421539
69 1.150287e+03 7.473669e+01
* time: 898.6862740516663
70 1.149189e+03 5.606074e+01
* time: 945.3819961547852
71 1.148944e+03 4.801310e+01
* time: 990.8594880104065
72 1.148932e+03 4.659470e+01
* time: 1039.267070055008
73 1.148930e+03 4.635167e+01
* time: 1071.9612309932709
74 1.148923e+03 4.568576e+01
* time: 1108.8116471767426
75 1.148907e+03 4.466140e+01
* time: 1154.0720889568329
76 1.148865e+03 4.249800e+01
* time: 1180.8980050086975
77 1.148758e+03 3.824645e+01
* time: 1208.1219000816345
78 1.148491e+03 3.552369e+01
* time: 1236.714872121811
79 1.147881e+03 3.037355e+01
* time: 1264.948559999466
80 1.146737e+03 4.466673e+01
* time: 1284.4998350143433
81 1.145347e+03 5.927932e+01
* time: 1305.3711149692535
82 1.144536e+03 7.572710e+01
* time: 1330.274708032608
83 1.144354e+03 7.035223e+01
* time: 1374.2716090679169
84 1.144341e+03 6.591571e+01
* time: 1390.9949271678925
85 1.144340e+03 6.496308e+01
* time: 1402.0454921722412
86 1.144333e+03 6.268335e+01
* time: 1411.9175870418549
87 1.144320e+03 5.950296e+01
* time: 1427.5938081741333
88 1.144281e+03 5.369030e+01
* time: 1436.283518075943
89 1.144186e+03 4.394749e+01
* time: 1445.050843000412
90 1.143949e+03 4.102392e+01
* time: 1453.3583381175995
91 1.143423e+03 3.752367e+01
* time: 1462.9932429790497
92 1.142498e+03 3.276758e+01
* time: 1472.3590321540833
93 1.141490e+03 5.217551e+01
* time: 1480.7130789756775
94 1.140969e+03 4.309381e+01
* time: 1488.746691942215
95 1.140872e+03 4.759588e+01
* time: 1497.1185879707336
96 1.140867e+03 4.727133e+01
* time: 1516.8032920360565
97 1.140866e+03 4.704459e+01
* time: 1542.2029631137848
98 1.140862e+03 4.638036e+01
* time: 1581.6669120788574
99 1.140853e+03 4.542717e+01
* time: 1604.6622800827026
100 1.140829e+03 4.357571e+01
* time: 1617.0542480945587
101 1.140767e+03 4.023368e+01
* time: 1629.1456451416016
102 1.140611e+03 3.376885e+01
* time: 1641.0703761577606
103 1.140237e+03 2.856928e+01
* time: 1649.546592950821
104 1.139452e+03 3.894107e+01
* time: 1658.2586181163788
105 1.138224e+03 5.063688e+01
* time: 1668.774108171463
106 1.137140e+03 4.762519e+01
* time: 1678.1966819763184
107 1.136763e+03 4.605051e+01
* time: 1687.0871460437775
108 1.136720e+03 3.993288e+01
* time: 1697.618327140808
109 1.136717e+03 3.890566e+01
* time: 1706.8516991138458
110 1.136716e+03 3.883647e+01
* time: 1715.1650550365448
111 1.136712e+03 3.865523e+01
* time: 1724.0706040859222
112 1.136702e+03 3.838845e+01
* time: 1733.4759459495544
113 1.136674e+03 3.790005e+01
* time: 1746.6718690395355
114 1.136604e+03 3.705604e+01
* time: 1758.4053790569305
115 1.136418e+03 3.550032e+01
* time: 1768.1316890716553
116 1.135940e+03 3.257323e+01
* time: 1776.9608490467072
117 1.134744e+03 2.697619e+01
* time: 1785.4685339927673
118 1.132016e+03 4.861352e+01
* time: 1793.8143169879913
119 1.127206e+03 6.900228e+01
* time: 1801.833261013031
120 1.122365e+03 5.795238e+01
* time: 1810.062390089035
121 1.119760e+03 2.908163e+01
* time: 1818.3642930984497
122 1.118972e+03 2.137897e+01
* time: 1826.815092086792
123 1.118870e+03 1.909827e+01
* time: 1835.1649680137634
124 1.118863e+03 1.886936e+01
* time: 1844.1320791244507
125 1.118862e+03 1.852746e+01
* time: 1856.5950450897217
126 1.118861e+03 1.864804e+01
* time: 1866.5084581375122
127 1.118861e+03 1.862831e+01
* time: 1875.7426221370697
128 1.118860e+03 1.865947e+01
* time: 1884.8379831314087
129 1.118859e+03 1.869546e+01
* time: 1893.4879431724548
130 1.118856e+03 1.878621e+01
* time: 1901.8185760974884
131 1.118850e+03 1.894171e+01
* time: 1910.514722108841
132 1.118834e+03 1.922062e+01
* time: 1920.2844569683075
133 1.118794e+03 1.967952e+01
* time: 1929.663906097412
134 1.118693e+03 2.042289e+01
* time: 1938.819440126419
135 1.118431e+03 2.157038e+01
* time: 1950.001897096634
136 1.117752e+03 2.595717e+01
* time: 1964.6537470817566
137 1.116008e+03 4.715215e+01
* time: 1975.397437095642
138 1.111680e+03 7.743613e+01
* time: 1985.6807510852814
139 1.102031e+03 7.499119e+01
* time: 1994.522628068924
140 1.090653e+03 6.250170e+01
* time: 2003.9429359436035
141 1.087845e+03 8.960496e+01
* time: 2013.5378630161285
142 1.087404e+03 1.190607e+02
* time: 2023.2112259864807
143 1.085564e+03 6.087787e+01
* time: 2033.3335530757904
144 1.085185e+03 1.947487e+01
* time: 2043.5302469730377
145 1.085039e+03 1.959012e+01
* time: 2053.359142065048
146 1.084989e+03 1.990081e+01
* time: 2064.4666249752045
147 1.084983e+03 2.009415e+01
* time: 2075.745949983597
148 1.084973e+03 2.065532e+01
* time: 2085.8533611297607
149 1.084971e+03 2.079993e+01
* time: 2095.0596799850464
150 1.084971e+03 2.084259e+01
* time: 2103.9194509983063
151 1.084970e+03 2.084552e+01
* time: 2112.93279004097
152 1.084970e+03 2.084403e+01
* time: 2122.2051179409027
153 1.084968e+03 2.084868e+01
* time: 2132.152698993683
154 1.084963e+03 2.088252e+01
* time: 2149.3556940555573
155 1.084954e+03 2.098670e+01
* time: 2162.3499059677124
156 1.084935e+03 2.121929e+01
* time: 2173.0724289417267
157 1.084905e+03 2.160241e+01
* time: 2182.630905151367
158 1.084865e+03 2.207894e+01
* time: 2191.1150369644165
159 1.084808e+03 2.259104e+01
* time: 2200.0709969997406
160 1.084701e+03 2.318755e+01
* time: 2209.0549890995026
161 1.084458e+03 2.380805e+01
* time: 2218.8124191761017
162 1.083881e+03 2.403001e+01
* time: 2228.604233980179
163 1.082583e+03 2.865879e+01
* time: 2239.9191510677338
164 1.081739e+03 3.341843e+01
* time: 2253.2874200344086
165 1.078499e+03 5.054188e+01
* time: 2265.629585981369
166 1.075168e+03 6.966535e+01
* time: 2277.107300043106
167 1.069624e+03 6.533378e+01
* time: 2287.227961063385
168 1.068098e+03 3.856272e+01
* time: 2300.397463083267
169 1.065571e+03 1.611062e+01
* time: 2314.5375130176544
170 1.065230e+03 1.527933e+01
* time: 2324.6764891147614
171 1.065088e+03 1.385503e+01
* time: 2335.1522631645203
172 1.065086e+03 1.382742e+01
* time: 2345.117233991623
173 1.065086e+03 1.383663e+01
* time: 2354.6991741657257
174 1.065086e+03 1.384273e+01
* time: 2363.707862138748
175 1.065086e+03 1.385109e+01
* time: 2373.4430961608887
176 1.065084e+03 1.386455e+01
* time: 2382.745407104492
177 1.065082e+03 1.388591e+01
* time: 2392.2373571395874
178 1.065074e+03 1.391915e+01
* time: 2403.7202920913696
179 1.065053e+03 1.397116e+01
* time: 2417.62606716156
180 1.065000e+03 1.742031e+01
* time: 2429.4155411720276
181 1.064862e+03 2.855619e+01
* time: 2440.538100004196
182 1.064503e+03 4.585022e+01
* time: 2451.8257830142975
183 1.063581e+03 7.120893e+01
* time: 2463.167657136917
184 1.061297e+03 1.015565e+02
* time: 2473.319419145584
185 1.056901e+03 1.172079e+02
* time: 2482.9640481472015
186 1.052919e+03 1.304962e+02
* time: 2491.9052131175995
187 1.048870e+03 1.445092e+02
* time: 2501.2800829410553
188 1.037225e+03 8.004893e+01
* time: 2511.557326078415
189 1.035989e+03 6.689301e+01
* time: 2522.5835371017456
190 1.035708e+03 3.351073e+01
* time: 2533.8712990283966
191 1.035220e+03 1.023182e+01
* time: 2548.2996430397034
192 1.035148e+03 4.588384e+00
* time: 2560.0548300743103
193 1.035107e+03 3.198226e+00
* time: 2570.0095920562744
194 1.035086e+03 3.196151e+00
* time: 2580.793991088867
195 1.035067e+03 3.160387e+00
* time: 2590.7326350212097
196 1.035054e+03 3.147012e+00
* time: 2600.224065065384
197 1.035047e+03 3.144947e+00
* time: 2609.702609062195
198 1.035046e+03 3.157031e+00
* time: 2624.5390470027924
199 1.035046e+03 3.159475e+00
* time: 2642.155135154724
200 1.035046e+03 3.160059e+00
* time: 2657.7700049877167
201 1.035046e+03 3.159779e+00
* time: 2669.1894719600677
202 1.035046e+03 3.158994e+00
* time: 2683.2103910446167
203 1.035046e+03 3.157781e+00
* time: 2692.936823129654
204 1.035046e+03 3.155534e+00
* time: 2702.3234000205994
205 1.035045e+03 3.151677e+00
* time: 2712.068423986435
206 1.035044e+03 3.144628e+00
* time: 2722.602807998657
207 1.035041e+03 3.131349e+00
* time: 2732.3776140213013
208 1.035033e+03 3.104858e+00
* time: 2741.1890740394592
209 1.035012e+03 3.049096e+00
* time: 2749.8483719825745
210 1.034957e+03 3.798101e+00
* time: 2760.982553958893
211 1.034816e+03 6.029607e+00
* time: 2772.9031381607056
212 1.034472e+03 9.154140e+00
* time: 2784.46430516243
213 1.033706e+03 1.258550e+01
* time: 2796.1049270629883
214 1.032296e+03 1.427962e+01
* time: 2807.5622651576996
215 1.030400e+03 1.353250e+01
* time: 2831.049805164337
216 1.029562e+03 1.451097e+01
* time: 2851.1893751621246
217 1.028793e+03 1.534582e+01
* time: 2866.5833320617676
218 1.028722e+03 2.309928e+01
* time: 2880.0622301101685
219 1.028409e+03 1.971948e+01
* time: 2891.0739810466766
220 1.028326e+03 1.936501e+01
* time: 2900.0354540348053
221 1.028234e+03 2.040540e+01
* time: 2920.857339143753
222 1.028219e+03 2.201826e+01
* time: 2935.3177139759064
223 1.028180e+03 2.177994e+01
* time: 2947.1054689884186
224 1.028176e+03 2.186278e+01
* time: 2968.3688700199127
225 1.028175e+03 2.201733e+01
* time: 2980.8037869930267
226 1.028174e+03 2.197587e+01
* time: 2990.280543088913
227 1.028174e+03 2.195893e+01
* time: 2999.281928062439
228 1.028172e+03 2.186809e+01
* time: 3015.903007030487
229 1.028169e+03 2.178132e+01
* time: 3031.9204399585724
230 1.028160e+03 2.162024e+01
* time: 3054.154732942581
231 1.028139e+03 2.139041e+01
* time: 3066.4955101013184
232 1.028082e+03 2.100359e+01
* time: 3077.6246061325073
233 1.027936e+03 2.808708e+01
* time: 3087.8173410892487
234 1.027561e+03 4.544763e+01
* time: 3110.7782230377197
235 1.026628e+03 7.187426e+01
* time: 3118.881999015808
236 1.024510e+03 1.065337e+02
* time: 3126.5078880786896
237 1.020621e+03 1.329824e+02
* time: 3134.26610994339
238 1.015709e+03 1.045562e+02
* time: 3141.806561946869
239 1.014365e+03 2.918128e+01
* time: 3150.7529051303864
240 1.013219e+03 1.803992e+01
* time: 3160.667901992798
241 1.012886e+03 1.450533e+01
* time: 3169.6189131736755
242 1.012748e+03 9.039639e+00
* time: 3185.847612142563
243 1.012704e+03 2.255849e+00
* time: 3194.1536099910736
244 1.012698e+03 1.107778e+00
* time: 3202.0567860603333
245 1.012697e+03 9.968755e-01
* time: 3209.287098169327
246 1.012697e+03 9.970599e-01
* time: 3216.6433091163635
247 1.012697e+03 9.971842e-01
* time: 3225.622558116913
248 1.012697e+03 9.972169e-01
* time: 3233.1766381263733
249 1.012697e+03 9.972169e-01
* time: 3241.4692890644073
250 1.012697e+03 9.973243e-01
* time: 3255.4805409908295
251 1.012697e+03 9.973078e-01
* time: 3263.834365129471
252 1.012697e+03 9.972048e-01
* time: 3271.2622430324554
253 1.012697e+03 9.970536e-01
* time: 3278.3843700885773
254 1.012697e+03 9.967402e-01
* time: 3285.1638259887695
255 1.012697e+03 9.961565e-01
* time: 3292.202826023102
256 1.012696e+03 9.949478e-01
* time: 3301.139825105667
257 1.012694e+03 1.312115e+00
* time: 3314.0078580379486
258 1.012690e+03 2.118030e+00
* time: 3340.3247499465942
259 1.012679e+03 3.352490e+00
* time: 3362.176973104477
260 1.012651e+03 5.076526e+00
* time: 3396.4928760528564
261 1.012594e+03 6.891495e+00
* time: 3412.8716809749603
262 1.012501e+03 7.239074e+00
* time: 3443.629009962082
263 1.012418e+03 4.479611e+00
* time: 3459.2118260860443
264 1.012390e+03 1.267149e+00
* time: 3495.185975074768
265 1.012386e+03 8.840304e-01
* time: 3528.58144903183
266 1.012386e+03 8.549097e-01
* time: 3565.8253951072693
267 1.012386e+03 8.516869e-01
* time: 3588.4442241191864
268 1.012386e+03 8.516869e-01
* time: 3617.1038789749146
269 1.012386e+03 8.516869e-01
* time: 3657.302843093872
270 1.012386e+03 8.516817e-01
* time: 3688.437945127487
271 1.012386e+03 8.516817e-01
* time: 3708.840376138687
272 1.012386e+03 8.516817e-01
* time: 3740.1402111053467
273 1.012386e+03 8.516817e-01
* time: 3756.975336074829
274 1.012386e+03 8.516817e-01
* time: 3787.3287370204926
275 1.012386e+03 8.516817e-01
* time: 3801.0230989456177
FittedPumasModel
Successful minimization: true
Likelihood approximation: FOCE
Likelihood Optimizer: BFGS
Dynamical system type: Nonlinear ODE
Solver(s): OrdinaryDiffEq.Rodas5P
Log-likelihood value: -1012.3861
Number of subjects: 31
Number of parameters: Fixed Optimized
0 18
Observation records: Active Missing
conc: 239 47
pca: 225 61
Total: 464 108
--------------------------
Estimate
--------------------------
pop_CL 0.13136
pop_V 8.0217
pop_tabs 0.41997
pop_lag 0.92474
pop_e0 96.692
pop_emax -1.0629
pop_c50 1.5003
pop_tover 14.06
pk_Ω₁,₁ 0.05162
pk_Ω₂,₂ 0.019945
pk_Ω₃,₃ 1.0369
pd_Ω₁,₁ 0.0025468
pd_Ω₂,₂ 0.00037903
pd_Ω₃,₃ 0.15787
pd_Ω₄,₄ 0.012052
σ_prop 0.088654
σ_add 0.36167
σ_fx 3.5399
--------------------------
On the other hand, if it is known that a differential equation is non-stiff (this might be difficult to guarantee for all admissible parameter values), a non-stiff solver such as Tsit5
at high tolerances or Vern7
at low tolerances could be an alternative to the default solver:
# Fitting with the non-stiff solver Vern7 at low tolerances (relative: 1e-8, absolute: 1e-12)
fit(
warfarin_pkpd_model,
pop,init_params(warfarin_pkpd_model),
FOCE();
= (; alg = Vern7(), reltol = 1e-8, abstol = 1e-12),
diffeq_options )
[ Info: Checking the initial parameter values.
[ Info: The initial negative log likelihood and its gradient are finite. Check passed.
Iter Function value Gradient norm
0 2.998927e+06 5.740466e+06
* time: 0.00010085105895996094
1 4.804005e+05 8.280312e+05
* time: 1.5635318756103516
2 3.595691e+05 6.035843e+05
* time: 2.973449945449829
3 1.691071e+05 2.702705e+05
* time: 4.2098469734191895
4 9.193152e+04 1.451040e+05
* time: 5.552165985107422
5 4.768539e+04 6.630202e+04
* time: 6.8227269649505615
6 2.904205e+04 3.416919e+04
* time: 8.140636920928955
7 1.840153e+04 1.671742e+04
* time: 9.441663026809692
8 1.280989e+04 1.022268e+04
* time: 10.77936601638794
9 9.713287e+03 9.298485e+03
* time: 12.059754848480225
10 7.628127e+03 8.330656e+03
* time: 13.41556692123413
11 6.236278e+03 7.340367e+03
* time: 14.714256048202515
12 5.271958e+03 6.303773e+03
* time: 16.078125953674316
13 4.411047e+03 4.994130e+03
* time: 17.400704860687256
14 3.732936e+03 3.538184e+03
* time: 18.949026823043823
15 3.399273e+03 2.450045e+03
* time: 20.62395191192627
16 3.302987e+03 1.848476e+03
* time: 25.080665826797485
17 3.291788e+03 1.642864e+03
* time: 28.05945897102356
18 3.290743e+03 1.594765e+03
* time: 30.176079988479614
19 3.289687e+03 1.562641e+03
* time: 32.12293887138367
20 3.286186e+03 1.496632e+03
* time: 33.85771083831787
21 3.277917e+03 1.402865e+03
* time: 35.5943169593811
22 3.255959e+03 1.254297e+03
* time: 37.18429493904114
23 3.200980e+03 1.038621e+03
* time: 38.81414198875427
24 3.063521e+03 7.403605e+02
* time: 40.30588984489441
25 2.730316e+03 3.768701e+02
* time: 41.80408501625061
26 1.934186e+03 2.246121e+02
* time: 43.00061798095703
27 1.704399e+03 2.219826e+02
* time: 47.89810395240784
28 1.471183e+03 2.025826e+02
* time: 56.518481969833374
29 1.316789e+03 1.554714e+02
* time: 66.16471004486084
30 1.226363e+03 1.775666e+02
* time: 68.11983394622803
31 1.206626e+03 2.780980e+02
* time: 69.6958019733429
32 1.197830e+03 2.430550e+02
* time: 71.36149883270264
33 1.190342e+03 2.327905e+02
* time: 72.73438596725464
34 1.178783e+03 2.105809e+02
* time: 74.20251488685608
35 1.177853e+03 2.228585e+02
* time: 76.28807783126831
36 1.177763e+03 2.256116e+02
* time: 78.41123604774475
37 1.177745e+03 2.259359e+02
* time: 80.35065484046936
38 1.177670e+03 2.264633e+02
* time: 82.22239804267883
39 1.177500e+03 2.264574e+02
* time: 84.12630200386047
40 1.177033e+03 2.244136e+02
* time: 85.90185689926147
41 1.175894e+03 2.160226e+02
* time: 87.9045569896698
42 1.173240e+03 1.905912e+02
* time: 89.87055397033691
43 1.168189e+03 1.332783e+02
* time: 91.62784004211426
44 1.161835e+03 5.289952e+01
* time: 95.45163083076477
45 1.158453e+03 6.595387e+01
* time: 101.9491708278656
46 1.157897e+03 6.411419e+01
* time: 103.81520986557007
47 1.157870e+03 6.098509e+01
* time: 105.44981002807617
48 1.157868e+03 6.032325e+01
* time: 107.01936388015747
49 1.157864e+03 5.920812e+01
* time: 108.46431493759155
50 1.157856e+03 5.740225e+01
* time: 109.97215604782104
51 1.157832e+03 5.406795e+01
* time: 111.35535287857056
52 1.157774e+03 4.797879e+01
* time: 113.28410601615906
53 1.157628e+03 3.658732e+01
* time: 116.53428292274475
54 1.157299e+03 3.707874e+01
* time: 119.05746793746948
55 1.156690e+03 3.452947e+01
* time: 121.14047002792358
56 1.155959e+03 4.269004e+01
* time: 123.06517696380615
57 1.155557e+03 5.000322e+01
* time: 124.99146604537964
58 1.155472e+03 4.219349e+01
* time: 126.7311680316925
59 1.155463e+03 3.895049e+01
* time: 128.3339970111847
60 1.155462e+03 3.872099e+01
* time: 129.4409909248352
61 1.155458e+03 3.823491e+01
* time: 130.61320304870605
62 1.155450e+03 3.747370e+01
* time: 131.68385696411133
63 1.155426e+03 3.608648e+01
* time: 132.77817296981812
64 1.155365e+03 3.361118e+01
* time: 133.86697483062744
65 1.155207e+03 2.894678e+01
* time: 135.00072503089905
66 1.154814e+03 2.504284e+01
* time: 136.139466047287
67 1.153910e+03 4.541126e+01
* time: 137.23268485069275
68 1.152230e+03 7.426044e+01
* time: 138.36527681350708
69 1.150286e+03 7.473044e+01
* time: 139.43518090248108
70 1.149188e+03 5.605967e+01
* time: 140.59311890602112
71 1.148944e+03 4.801123e+01
* time: 141.67724895477295
72 1.148931e+03 4.659478e+01
* time: 142.79204893112183
73 1.148929e+03 4.635180e+01
* time: 143.93104100227356
74 1.148922e+03 4.568579e+01
* time: 145.03991103172302
75 1.148907e+03 4.466128e+01
* time: 146.3771059513092
76 1.148864e+03 4.249742e+01
* time: 148.11872696876526
77 1.148758e+03 3.824840e+01
* time: 149.93944883346558
78 1.148490e+03 3.552459e+01
* time: 151.6986379623413
79 1.147881e+03 3.038181e+01
* time: 153.57162594795227
80 1.146736e+03 4.467281e+01
* time: 155.25954794883728
81 1.145347e+03 5.929158e+01
* time: 156.83591294288635
82 1.144535e+03 7.572784e+01
* time: 157.95326590538025
83 1.144354e+03 7.034986e+01
* time: 159.19690585136414
84 1.144341e+03 6.591584e+01
* time: 160.30735087394714
85 1.144340e+03 6.496353e+01
* time: 161.812805891037
86 1.144333e+03 6.268267e+01
* time: 164.6780068874359
87 1.144320e+03 5.950121e+01
* time: 167.5503158569336
88 1.144281e+03 5.368613e+01
* time: 169.54465889930725
89 1.144186e+03 4.393956e+01
* time: 171.03845691680908
90 1.143949e+03 4.102324e+01
* time: 172.5093879699707
91 1.143423e+03 3.752005e+01
* time: 173.9874210357666
92 1.142497e+03 3.278884e+01
* time: 175.46158981323242
93 1.141489e+03 5.217814e+01
* time: 176.93660688400269
94 1.140969e+03 4.310032e+01
* time: 178.29149794578552
95 1.140872e+03 4.759545e+01
* time: 179.76309490203857
96 1.140867e+03 4.727034e+01
* time: 181.105073928833
97 1.140866e+03 4.704365e+01
* time: 182.47545981407166
98 1.140862e+03 4.637923e+01
* time: 183.747820854187
99 1.140853e+03 4.542582e+01
* time: 185.09347581863403
100 1.140829e+03 4.357378e+01
* time: 186.2819948196411
101 1.140767e+03 4.023060e+01
* time: 187.4135389328003
102 1.140610e+03 3.376323e+01
* time: 188.5290129184723
103 1.140237e+03 2.857028e+01
* time: 189.61854004859924
104 1.139451e+03 3.895699e+01
* time: 190.73187899589539
105 1.138223e+03 5.064036e+01
* time: 191.79952001571655
106 1.137140e+03 4.763314e+01
* time: 192.93470001220703
107 1.136763e+03 4.604968e+01
* time: 193.96098589897156
108 1.136720e+03 3.993479e+01
* time: 195.05673694610596
109 1.136717e+03 3.890565e+01
* time: 196.10263895988464
110 1.136716e+03 3.883645e+01
* time: 197.18500590324402
111 1.136712e+03 3.865517e+01
* time: 198.26260900497437
112 1.136702e+03 3.838833e+01
* time: 199.3391568660736
113 1.136674e+03 3.789981e+01
* time: 200.49608302116394
114 1.136603e+03 3.705558e+01
* time: 201.55509996414185
115 1.136418e+03 3.549940e+01
* time: 202.68112897872925
116 1.135940e+03 3.257136e+01
* time: 203.74488401412964
117 1.134743e+03 2.697248e+01
* time: 204.99595403671265
118 1.132014e+03 4.863416e+01
* time: 206.21633887290955
119 1.127203e+03 6.901202e+01
* time: 207.85340690612793
120 1.122363e+03 5.794311e+01
* time: 211.80425691604614
121 1.119759e+03 2.906766e+01
* time: 213.65995287895203
122 1.118972e+03 2.137843e+01
* time: 215.50017499923706
123 1.118870e+03 1.909639e+01
* time: 217.05376482009888
124 1.118863e+03 1.887024e+01
* time: 218.55808997154236
125 1.118862e+03 1.852691e+01
* time: 219.98240685462952
126 1.118861e+03 1.864783e+01
* time: 221.39966082572937
127 1.118861e+03 1.862870e+01
* time: 222.65939688682556
128 1.118860e+03 1.866000e+01
* time: 223.93079900741577
129 1.118859e+03 1.869625e+01
* time: 225.24359893798828
130 1.118856e+03 1.878717e+01
* time: 226.48283791542053
131 1.118850e+03 1.894283e+01
* time: 227.70876693725586
132 1.118834e+03 1.922180e+01
* time: 228.8762059211731
133 1.118794e+03 1.968060e+01
* time: 230.08033084869385
134 1.118693e+03 2.042353e+01
* time: 231.23646998405457
135 1.118431e+03 2.156998e+01
* time: 232.43370485305786
136 1.117751e+03 2.596600e+01
* time: 233.56890487670898
137 1.116007e+03 4.716393e+01
* time: 234.72252488136292
138 1.111677e+03 7.744992e+01
* time: 235.90466785430908
139 1.102024e+03 7.496268e+01
* time: 237.0937340259552
140 1.090649e+03 6.246756e+01
* time: 238.3414590358734
141 1.087843e+03 8.941834e+01
* time: 239.5851218700409
142 1.087400e+03 1.188539e+02
* time: 240.97220993041992
143 1.085563e+03 6.082023e+01
* time: 242.4430480003357
144 1.085185e+03 1.947062e+01
* time: 244.05774903297424
145 1.085039e+03 1.958695e+01
* time: 245.42364597320557
146 1.084989e+03 1.989918e+01
* time: 247.16838383674622
147 1.084983e+03 2.009395e+01
* time: 248.5596878528595
148 1.084973e+03 2.065564e+01
* time: 250.0434489250183
149 1.084971e+03 2.080010e+01
* time: 251.4054868221283
150 1.084971e+03 2.084232e+01
* time: 252.79676985740662
151 1.084970e+03 2.084506e+01
* time: 254.0913689136505
152 1.084970e+03 2.084321e+01
* time: 255.40042400360107
153 1.084968e+03 2.084737e+01
* time: 256.6979010105133
154 1.084963e+03 2.088045e+01
* time: 257.9676058292389
155 1.084954e+03 2.098376e+01
* time: 259.2608778476715
156 1.084935e+03 2.121619e+01
* time: 260.49714493751526
157 1.084904e+03 2.160158e+01
* time: 261.7877678871155
158 1.084864e+03 2.208328e+01
* time: 263.0279948711395
159 1.084807e+03 2.260004e+01
* time: 264.3532350063324
160 1.084700e+03 2.319870e+01
* time: 265.5825810432434
161 1.084457e+03 2.382124e+01
* time: 266.9293019771576
162 1.083882e+03 2.406095e+01
* time: 268.20076084136963
163 1.082589e+03 2.849219e+01
* time: 269.86068892478943
164 1.081660e+03 3.379635e+01
* time: 271.4090678691864
165 1.078576e+03 4.983497e+01
* time: 273.0225110054016
166 1.075004e+03 6.870527e+01
* time: 274.5193569660187
167 1.069377e+03 6.309176e+01
* time: 275.9228479862213
168 1.067978e+03 3.811944e+01
* time: 277.26304388046265
169 1.065604e+03 1.825143e+01
* time: 278.61276483535767
170 1.065251e+03 1.739631e+01
* time: 279.9821689128876
171 1.065088e+03 1.384781e+01
* time: 281.32171082496643
172 1.065086e+03 1.382386e+01
* time: 282.9433870315552
173 1.065086e+03 1.383866e+01
* time: 284.6018970012665
174 1.065086e+03 1.384608e+01
* time: 286.16933703422546
175 1.065086e+03 1.385678e+01
* time: 287.58597588539124
176 1.065084e+03 1.387273e+01
* time: 289.0493450164795
177 1.065081e+03 1.389858e+01
* time: 290.45479702949524
178 1.065073e+03 1.393897e+01
* time: 291.8324179649353
179 1.065052e+03 1.400288e+01
* time: 293.222718000412
180 1.064997e+03 1.842052e+01
* time: 294.5791208744049
181 1.064853e+03 3.031555e+01
* time: 295.9338729381561
182 1.064477e+03 4.894991e+01
* time: 297.2636408805847
183 1.063505e+03 7.692533e+01
* time: 298.64283299446106
184 1.061069e+03 1.122192e+02
* time: 300.0239908695221
185 1.056257e+03 1.328367e+02
* time: 301.56972885131836
186 1.052013e+03 1.408303e+02
* time: 303.2595419883728
187 1.047713e+03 1.506293e+02
* time: 304.965430021286
188 1.039359e+03 9.039176e+01
* time: 306.5875189304352
189 1.037334e+03 8.945366e+01
* time: 308.29118394851685
190 1.035766e+03 2.509488e+01
* time: 310.01715302467346
191 1.035510e+03 5.323250e+00
* time: 311.8450210094452
192 1.035401e+03 3.326584e+00
* time: 313.63545203208923
193 1.035202e+03 1.011867e+01
* time: 315.8587238788605
194 1.035130e+03 9.089057e+00
* time: 318.66000390052795
195 1.035055e+03 3.941605e+00
* time: 320.7836048603058
196 1.035048e+03 3.160942e+00
* time: 322.8379089832306
197 1.035046e+03 3.164670e+00
* time: 324.80938482284546
198 1.035046e+03 3.160486e+00
* time: 326.60200095176697
199 1.035046e+03 3.159428e+00
* time: 328.3718559741974
200 1.035046e+03 3.159138e+00
* time: 330.21579003334045
201 1.035046e+03 3.158158e+00
* time: 331.9061498641968
202 1.035046e+03 3.156875e+00
* time: 333.6177718639374
203 1.035046e+03 3.154403e+00
* time: 335.3565728664398
204 1.035045e+03 3.150100e+00
* time: 336.99276304244995
205 1.035043e+03 3.141854e+00
* time: 338.61231899261475
206 1.035038e+03 3.125468e+00
* time: 340.17678904533386
207 1.035025e+03 3.887763e+00
* time: 341.84859800338745
208 1.034991e+03 6.349435e+00
* time: 343.4519729614258
209 1.034902e+03 1.025429e+01
* time: 345.3474099636078
210 1.034678e+03 1.614807e+01
* time: 348.23253297805786
211 1.034143e+03 2.389258e+01
* time: 351.4320168495178
212 1.033018e+03 3.065325e+01
* time: 353.58906292915344
213 1.031157e+03 2.988569e+01
* time: 355.58328795433044
214 1.029913e+03 2.000449e+01
* time: 357.50524282455444
215 1.029264e+03 1.806624e+01
* time: 359.45595598220825
216 1.028739e+03 1.876674e+01
* time: 361.33223485946655
217 1.028529e+03 1.788313e+01
* time: 363.2868709564209
218 1.028407e+03 1.869928e+01
* time: 365.07710003852844
219 1.028318e+03 2.103440e+01
* time: 366.791934967041
220 1.028259e+03 2.278116e+01
* time: 368.49730682373047
221 1.028193e+03 2.168720e+01
* time: 370.1747679710388
222 1.028177e+03 2.174690e+01
* time: 371.8954529762268
223 1.028175e+03 2.195554e+01
* time: 373.45725989341736
224 1.028175e+03 2.193442e+01
* time: 375.025593996048
225 1.028174e+03 2.202037e+01
* time: 376.589900970459
226 1.028174e+03 2.201903e+01
* time: 378.12916588783264
227 1.028173e+03 2.200004e+01
* time: 379.7699429988861
228 1.028171e+03 2.197102e+01
* time: 381.72052001953125
229 1.028164e+03 2.191369e+01
* time: 384.04214000701904
230 1.028148e+03 2.181626e+01
* time: 385.9039158821106
231 1.028106e+03 2.163664e+01
* time: 387.77518582344055
232 1.027999e+03 3.038649e+01
* time: 389.5375349521637
233 1.027721e+03 5.035972e+01
* time: 391.33608388900757
234 1.027021e+03 8.114702e+01
* time: 393.0492420196533
235 1.025404e+03 1.233505e+02
* time: 394.6975269317627
236 1.022343e+03 1.646448e+02
* time: 396.37359285354614
237 1.017556e+03 1.679173e+02
* time: 397.99962186813354
238 1.013261e+03 5.600035e+01
* time: 399.59919691085815
239 1.012953e+03 2.691504e+01
* time: 401.2294578552246
240 1.012758e+03 1.329779e+01
* time: 403.0844578742981
241 1.012704e+03 2.847118e+00
* time: 406.4183418750763
242 1.012699e+03 1.333414e+00
* time: 416.36782598495483
243 1.012697e+03 9.979337e-01
* time: 418.6901168823242
244 1.012697e+03 9.969545e-01
* time: 420.8819649219513
245 1.012697e+03 9.971341e-01
* time: 422.6343219280243
246 1.012697e+03 9.971972e-01
* time: 424.5509388446808
247 1.012697e+03 9.972237e-01
* time: 426.53851103782654
248 1.012697e+03 9.972927e-01
* time: 428.41429686546326
249 1.012697e+03 9.973804e-01
* time: 430.20265889167786
250 1.012697e+03 9.975286e-01
* time: 431.8866000175476
251 1.012697e+03 9.977461e-01
* time: 433.5195748806
252 1.012697e+03 9.980635e-01
* time: 435.15022897720337
253 1.012697e+03 9.984743e-01
* time: 436.83702397346497
254 1.012696e+03 9.988765e-01
* time: 438.42469596862793
255 1.012694e+03 9.988363e-01
* time: 440.04763984680176
256 1.012690e+03 1.474525e+00
* time: 441.7132430076599
257 1.012678e+03 2.326520e+00
* time: 443.29159593582153
258 1.012649e+03 3.510529e+00
* time: 444.8700189590454
259 1.012590e+03 4.737310e+00
* time: 446.3869159221649
260 1.012497e+03 4.897808e+00
* time: 447.9421339035034
261 1.012416e+03 2.875911e+00
* time: 449.47058391571045
262 1.012389e+03 9.301211e-01
* time: 451.22642993927
263 1.012386e+03 8.804644e-01
* time: 452.809378862381
264 1.012386e+03 8.547350e-01
* time: 454.4557628631592
265 1.012386e+03 8.516621e-01
* time: 456.1050579547882
266 1.012386e+03 8.511780e-01
* time: 457.7878530025482
267 1.012386e+03 8.489977e-01
* time: 459.3069999217987
268 1.012386e+03 8.463566e-01
* time: 460.8405690193176
269 1.012386e+03 8.415231e-01
* time: 462.47891902923584
270 1.012386e+03 8.340037e-01
* time: 464.2725908756256
271 1.012386e+03 8.215489e-01
* time: 467.33807492256165
272 1.012386e+03 8.013257e-01
* time: 469.4641718864441
273 1.012385e+03 7.680260e-01
* time: 471.569787979126
274 1.012383e+03 7.407512e-01
* time: 473.7246849536896
275 1.012377e+03 7.469071e-01
* time: 475.8183100223541
276 1.012364e+03 1.103988e+00
* time: 477.8015630245209
277 1.012330e+03 1.710981e+00
* time: 479.49283385276794
278 1.012257e+03 2.412723e+00
* time: 481.2047188282013
279 1.012125e+03 2.705165e+00
* time: 482.78704595565796
280 1.011948e+03 1.883597e+00
* time: 484.4337248802185
281 1.011809e+03 8.215973e-01
* time: 486.1502878665924
282 1.011775e+03 9.067203e-01
* time: 487.87961602211
283 1.011774e+03 9.565073e-01
* time: 489.38594698905945
284 1.011774e+03 9.628263e-01
* time: 490.93089604377747
285 1.011774e+03 9.654767e-01
* time: 492.50087881088257
286 1.011774e+03 9.654769e-01
* time: 494.37740182876587
287 1.011774e+03 9.656366e-01
* time: 498.4528298377991
288 1.011774e+03 9.656622e-01
* time: 502.6207528114319
289 1.011774e+03 9.656624e-01
* time: 505.4759418964386
290 1.011774e+03 9.656627e-01
* time: 507.77419686317444
291 1.011774e+03 9.660395e-01
* time: 509.9345908164978
292 1.011774e+03 9.661257e-01
* time: 511.86144399642944
293 1.011774e+03 9.663115e-01
* time: 513.6751699447632
294 1.011774e+03 9.662851e-01
* time: 515.5480098724365
295 1.011773e+03 9.654773e-01
* time: 517.3696918487549
296 1.011772e+03 9.623124e-01
* time: 519.0805349349976
297 1.011769e+03 9.521889e-01
* time: 520.6302828788757
298 1.011763e+03 9.239616e-01
* time: 522.2563509941101
299 1.011746e+03 9.242954e-01
* time: 523.7630770206451
300 1.011711e+03 1.257052e+00
* time: 525.2770118713379
301 1.011651e+03 1.361278e+00
* time: 526.8750038146973
302 1.011591e+03 9.076916e-01
* time: 528.4926760196686
303 1.011568e+03 3.084752e-01
* time: 529.9842410087585
304 1.011565e+03 2.946404e-01
* time: 531.4863159656525
305 1.011565e+03 2.821273e-01
* time: 533.0531280040741
306 1.011565e+03 2.766260e-01
* time: 534.5547268390656
307 1.011565e+03 2.778010e-01
* time: 536.0586438179016
308 1.011565e+03 2.830281e-01
* time: 537.6060948371887
309 1.011565e+03 2.842619e-01
* time: 539.1675560474396
310 1.011565e+03 2.861093e-01
* time: 540.7683999538422
311 1.011565e+03 2.861138e-01
* time: 542.404373884201
312 1.011565e+03 2.867641e-01
* time: 544.0829439163208
313 1.011565e+03 2.867756e-01
* time: 545.9445879459381
314 1.011565e+03 2.884589e-01
* time: 547.6148238182068
315 1.011565e+03 2.897335e-01
* time: 549.1516959667206
316 1.011565e+03 2.927744e-01
* time: 550.7983009815216
317 1.011565e+03 2.967424e-01
* time: 552.3457388877869
318 1.011565e+03 3.027347e-01
* time: 553.8705599308014
319 1.011564e+03 3.099827e-01
* time: 555.4642078876495
320 1.011563e+03 3.155826e-01
* time: 557.1347358226776
321 1.011560e+03 3.082844e-01
* time: 558.9035758972168
322 1.011555e+03 2.625121e-01
* time: 560.5107469558716
323 1.011548e+03 2.136271e-01
* time: 562.1413838863373
324 1.011542e+03 1.165033e-01
* time: 563.7222368717194
325 1.011538e+03 3.914224e-02
* time: 565.3413758277893
326 1.011537e+03 3.083888e-02
* time: 566.9747710227966
327 1.011536e+03 2.766364e-02
* time: 568.6530740261078
328 1.011536e+03 1.873262e-02
* time: 570.711532831192
329 1.011536e+03 5.453951e-03
* time: 572.6278629302979
330 1.011536e+03 4.882951e-03
* time: 575.0212609767914
331 1.011535e+03 4.331879e-03
* time: 576.7694199085236
332 1.011535e+03 2.711268e-03
* time: 578.4879369735718
333 1.011535e+03 4.474256e-04
* time: 580.1902258396149
FittedPumasModel
Successful minimization: true
Likelihood approximation: FOCE
Likelihood Optimizer: BFGS
Dynamical system type: Nonlinear ODE
Solver(s): OrdinaryDiffEq.Vern7
Log-likelihood value: -1011.5354
Number of subjects: 31
Number of parameters: Fixed Optimized
0 18
Observation records: Active Missing
conc: 239 47
pca: 225 61
Total: 464 108
-------------------------
Estimate
-------------------------
pop_CL 0.13137
pop_V 8.0255
pop_tabs 0.42635
pop_lag 0.92336
pop_e0 96.688
pop_emax -1.0628
pop_c50 1.5016
pop_tover 14.055
pk_Ω₁,₁ 0.050744
pk_Ω₂,₂ 0.021666
pk_Ω₃,₃ 0.97884
pd_Ω₁,₁ 0.0028828
pd_Ω₂,₂ 7.3726e-9
pd_Ω₃,₃ 0.15184
pd_Ω₄,₄ 0.011913
σ_prop 0.088292
σ_add 0.36362
σ_fx 3.5336
-------------------------
8 Concluding Remarks
In this tutorial, you have seen how to adjust the tolerances and the algorithm of the differential solver. Usually, the default differential equation solver in Pumas is an efficient choice. To reduce numerical issues, sometimes it can be helpful to decrease the default tolerances.