Differential Equations in Pumas

Author

David Widmann

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:

warfarin_pkpd_model = @model begin
    ...

    @dynamics begin
        Depot' = -Ka * Depot
        Central' = Ka * Depot - CL / Vc * Central
        Turnover' =
            rin * (1 + emax * (Central / Vc) / (cp50 + Central / Vc)) - kout * Turnover
    end

    ...
end

The same concise rewriting can be applied in a Pumas @model by defining auxiliary variables (“aliases”) in the @vars block:

warfarin_pkpd_model = @model begin
    ...

    @vars begin
        cp := Central / Vc
        ratein := Ka * Depot
        pd := 1 + emax * cp / (c50 + cp)
    end

    @dynamics begin
        Depot' = -ratein
        Central' = ratein - CL * cp
        Turnover' = rin * pd - kout * Turnover
    end

    ...
end
Tip

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
    cp = Central / Vc
    ratein := Ka * Depot
    pd := 1 + emax * cp / (c50 + cp)
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

warfarin_pkpd_model = @model begin
    @param begin
        # PK parameters
        """
        Clearance (L/h/70kg)
        """
        pop_CL  RealDomain(lower = 0.0, init = 0.134)
        """
        Central Volume L/70kg
        """
        pop_V  RealDomain(lower = 0.0, init = 8.11)
        """
        Absorption time (h)
        """
        pop_tabs  RealDomain(lower = 0.0, init = 0.523)
        """
        Lag time (h)
        """
        pop_lag  RealDomain(lower = 0.0, init = 0.1)
        # PD parameters
        """
        Baseline
        """
        pop_e0  RealDomain(lower = 0.0, init = 100.0)
        """
        Emax
        """
        pop_emax  RealDomain(init = -1.0)
        """
        EC50
        """
        pop_c50  RealDomain(lower = 0.0, init = 1.0)
        """
        Turnover
        """
        pop_tover  RealDomain(lower = 0.0, init = 14.0)
        # Inter-individual variability
        """
          - ΩCL
          - ΩVc
          - ΩTabs
        """
        pk_Ω  PDiagDomain([0.01, 0.01, 0.01])
        """
          - Ωe0
          - Ωemax
          - Ωec50
          - Ωturn
        """
        pd_Ω  PDiagDomain([0.01, 0.01, 0.01, 0.01])
        # Residual variability
        """
        Proportional residual error for drug concentration
        """
        σ_prop  RealDomain(lower = 0.0, init = 0.00752)
        """
        Additive residual error for drug concentration (mg/L)
        """
        σ_add  RealDomain(lower = 0.0, init = 0.0661)
        """
        Additive error for PCA
        """
        σ_fx  RealDomain(lower = 0.0, init = 0.01)
    end

    @random begin
        # mean = 0, covariance = pk_Ω
        pk_η ~ MvNormal(pk_Ω)
        # mean = 0, covariance = pd_Ω
        pd_η ~ MvNormal(pd_Ω)
    end

    @covariates FSZV FSZCL

    @pre begin
        # PK
        CL = FSZCL * pop_CL * exp(pk_η[1])
        Vc = FSZV * pop_V * exp(pk_η[2])
        tabs = pop_tabs * exp(pk_η[3])
        Ka = log(2) / tabs
        # PD
        e0 = pop_e0 * exp(pd_η[1])
        emax = pop_emax * exp(pd_η[2])
        c50 = pop_c50 * exp(pd_η[3])
        tover = pop_tover * exp(pd_η[4])
        kout = log(2) / tover
        rin = e0 * kout
        time = t
    end

    @dosecontrol begin
        lags = (Depot = pop_lag,)
    end

    @init begin
        Turnover = e0
    end

    # aliases for use in @dynamics and @derived
    @vars begin
        cp := Central / Vc
        ratein := Ka * Depot
        pd := 1 + emax * cp / (c50 + cp)
    end

    @dynamics begin
        Depot' = -ratein
        Central' = ratein - CL * cp
        Turnover' = rin * pd - kout * Turnover
    end

    @derived begin
        """
        Warfarin Concentration (mg/L)
        """
        conc ~ @. Normal(cp, sqrt((σ_prop * cp)^2 + σ_add^2))
        """
        PCA
        """
        pca ~ @. Normal(Turnover, σ_fx)
    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.

Important

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:

  1. 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.
  2. 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();
    diffeq_options = (; reltol = 1e-3, abstol = 1e-6),
)
[ 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();
    diffeq_options = (; reltol = 1e-8, abstol = 1e-12),
)
[ 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
-------------------------
Tip

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();
    diffeq_options = (; alg = Rodas5P(), reltol = 1e-8, abstol = 1e-12),
)
[ 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();
    diffeq_options = (; alg = Vern7(), reltol = 1e-8, abstol = 1e-12),
)
[ 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.