Error with Pumas v 2.2

Hi,

I am getting the following error ( the model was running fine previously):

ERROR: MethodError: no method matching _valid_foce_distribution(::Int64)
Closest candidates are:
  _valid_foce_distribution(::AbstractVector{<:Union{LogNormal{<:ForwardDiff.Dual}, Normal{<:ForwardDiff.Dual}, Pumas.FastNormal{<:ForwardDiff.Dual}}}) at /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:2575
  _valid_foce_distribution(::AbstractVector{<:Union{LogNormal{<:AbstractFloat}, Normal{<:AbstractFloat}}}) at /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:2578
  _valid_foce_distribution(::AbstractVector{<:Union{LogNormal{<:AbstractFloat}, Normal{<:AbstractFloat}, Pumas.FastNormal{<:AbstractFloat}}}) at /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:2583
  ...
Stacktrace:
  [1] (::Pumas.var"#533#535"{Val{4}, Val{4}})(ydual_dist::Int64)
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:2414
  [2] map (repeats 5 times)
    @ ./tuple.jl:224 [inlined]
  [3] map(::Function, ::NamedTuple{(:CONC, :conc, :E, :ANC, :NN, :KTR, :EDRUG, :EDRUGC, :FEED), Tuple{Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag, Float64}, Float64, 4}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag, Float64}, Float64, 4}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag, Float64}, Float64, 4}}, Vector{Normal{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag, Float64}, Float64, 4}}}, Int64, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag, Float64}, Float64, 4}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag, Float64}, Float64, 4}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag, Float64}, Float64, 4}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag, Float64}, Float64, 4}}}})
    @ Base ./namedtuple.jl:208
  [4] chunk_mode_jacobian(f::Pumas.var"#537#538"{PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{Float64, Float64, Float64, Float64, Float64, PDMats.PDiagMat{Float64, Vector{Float64}}, PDMats.PDiagMat{Float64, Vector{Float64}}, Float64}}, NamedTuple{(), Tuple{}}}, x::Vector{Float64}, cfg::ForwardDiff.JacobianConfig{ForwardDiff.Tag{Pumas.Tag, Float64}, Float64, 4, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag, Float64}, Float64, 4}}}, #unused#::Val{4})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:2411
  [5] _derived_vηorth_gradient(m::PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, subject::Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, param::NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{Float64, Float64, Float64, Float64, Float64, PDMats.PDiagMat{Float64, Vector{Float64}}, PDMats.PDiagMat{Float64, Vector{Float64}}, Float64}}, vrandeffsorth::Vector{Float64}, diffeq_options::NamedTuple{(), Tuple{}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:2479
  [6] _∂²l∂η²(m::PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, subject::Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, param::NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{Float64, Float64, Float64, Float64, Float64, PDMats.PDiagMat{Float64, Vector{Float64}}, PDMats.PDiagMat{Float64, Vector{Float64}}, Float64}}, vrandeffsorth::Vector{Float64}, approx::Pumas.FOCE, diffeq_options::NamedTuple{(), Tuple{}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:2524
  [7] (::Pumas.var"#467#468"{PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{Float64, Float64, Float64, Float64, Float64, PDMats.PDiagMat{Float64, Vector{Float64}}, PDMats.PDiagMat{Float64, Vector{Float64}}, Float64}}, Pumas.FOCE, NamedTuple{(), Tuple{}}, TransformVariables.TransformTuple{NamedTuple{(:η_1, :η_2), Tuple{Pumas.MvNormalTransform{ZeroMeanDiagNormal{Tuple{Base.OneTo{Int64}}}}, Pumas.MvNormalTransform{ZeroMeanDiagNormal{Tuple{Base.OneTo{Int64}}}}}}}, ParamSet{NamedTuple{(:η_1, :η_2), Tuple{ZeroMeanDiagNormal{Tuple{Base.OneTo{Int64}}}, ZeroMeanDiagNormal{Tuple{Base.OneTo{Int64}}}}}}, Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}})(F::Float64, G::Vector{Float64}, H::Matrix{Float64}, x::Vector{Float64})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:907
  [8] (::NLSolversBase.var"#98#104"{NLSolversBase.InplaceObjective{Nothing, Nothing, Pumas.var"#467#468"{PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{Float64, Float64, Float64, Float64, Float64, PDMats.PDiagMat{Float64, Vector{Float64}}, PDMats.PDiagMat{Float64, Vector{Float64}}, Float64}}, Pumas.FOCE, NamedTuple{(), Tuple{}}, TransformVariables.TransformTuple{NamedTuple{(:η_1, :η_2), Tuple{Pumas.MvNormalTransform{ZeroMeanDiagNormal{Tuple{Base.OneTo{Int64}}}}, Pumas.MvNormalTransform{ZeroMeanDiagNormal{Tuple{Base.OneTo{Int64}}}}}}}, ParamSet{NamedTuple{(:η_1, :η_2), Tuple{ZeroMeanDiagNormal{Tuple{Base.OneTo{Int64}}}, ZeroMeanDiagNormal{Tuple{Base.OneTo{Int64}}}}}}, Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}, Nothing, Nothing}, Float64})(G::Vector{Float64}, H::Matrix{Float64}, x::Vector{Float64})
    @ NLSolversBase /builds/PumasAI/PumasSystemImages-jl/.julia/packages/NLSolversBase/cfJrN/src/objective_types/incomplete.jl:115
  [9] value_gradient_hessian!!(obj::NLSolversBase.TwiceDifferentiable{Float64, Vector{Float64}, Matrix{Float64}, Vector{Float64}}, x::Vector{Float64})
    @ NLSolversBase /builds/PumasAI/PumasSystemImages-jl/.julia/packages/NLSolversBase/cfJrN/src/objective_types/incomplete.jl:138
 [10] initial_state(method::Optim.NewtonTrustRegion{Float64}, options::Optim.Options{Float64, Nothing}, d::NLSolversBase.TwiceDifferentiable{Float64, Vector{Float64}, Matrix{Float64}, Vector{Float64}}, initial_x::Vector{Float64})
    @ Optim /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Optim/rpjtl/src/multivariate/solvers/second_order/newton_trust_region.jl:310
 [11] optimize
    @ /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Optim/rpjtl/src/multivariate/optimize/optimize.jl:36 [inlined]
 [12] _orth_empirical_bayes!(vrandeffsorth::Vector{Float64}, m::PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, _subject::Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, param::NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{Float64, Float64, Float64, Float64, Float64, PDMats.PDiagMat{Float64, Vector{Float64}}, PDMats.PDiagMat{Float64, Vector{Float64}}, Float64}}, approx::Pumas.FOCE, diffeq_options::NamedTuple{(), Tuple{}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:935
 [13] _update_ebes_and_evaluate_marginal_nll!(m::PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, subject::Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, param::NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{Float64, Float64, Float64, Float64, Float64, PDMats.PDiagMat{Float64, Vector{Float64}}, PDMats.PDiagMat{Float64, Vector{Float64}}, Float64}}, vrandeffsorth::Vector{Float64}, vrandeffsorth_tmp::Vector{Float64}, approx::Pumas.FOCE, diffeq_options::NamedTuple{(), Tuple{}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:3008
 [14] _update_ebes_and_evaluate_marginal_nll_threads!(m::PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, population::Vector{Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}, param::NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{Float64, Float64, Float64, Float64, Float64, PDMats.PDiagMat{Float64, Vector{Float64}}, PDMats.PDiagMat{Float64, Vector{Float64}}, Float64}}, vvrandeffs::Vector{Vector{Float64}}, vvrandeffs_tmp::Vector{Vector{Float64}}, approx::Pumas.FOCE, diffeq_options::NamedTuple{(), Tuple{}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:3057
 [15] (::Pumas.var"#581#583"{PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, Vector{Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}, Pumas.FOCE, EnsembleThreads, NamedTuple{(), Tuple{}}, TransformVariables.TransformTuple{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, Pumas.PDiagTransform, Pumas.PDiagTransform, TransformVariables.ShiftedExp{true, Float64}}}}, Vector{Vector{Float64}}, Vector{Vector{Float64}}, Vector{Float64}})(f::Float64, g::Vector{Float64}, _vparam::Vector{Float64})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:3699
 [16] (::NLSolversBase.var"#69#70"{NLSolversBase.InplaceObjective{Nothing, Pumas.var"#581#583"{PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, Vector{Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}, Pumas.FOCE, EnsembleThreads, NamedTuple{(), Tuple{}}, TransformVariables.TransformTuple{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, Pumas.PDiagTransform, Pumas.PDiagTransform, TransformVariables.ShiftedExp{true, Float64}}}}, Vector{Vector{Float64}}, Vector{Vector{Float64}}, Vector{Float64}}, Nothing, Nothing, Nothing}, Float64})(G::Vector{Float64}, x::Vector{Float64})
    @ NLSolversBase /builds/PumasAI/PumasSystemImages-jl/.julia/packages/NLSolversBase/cfJrN/src/objective_types/incomplete.jl:54
 [17] value_gradient!!(obj::NLSolversBase.OnceDifferentiable{Float64, Vector{Float64}, Vector{Float64}}, x::Vector{Float64})
    @ NLSolversBase /builds/PumasAI/PumasSystemImages-jl/.julia/packages/NLSolversBase/cfJrN/src/interface.jl:82
 [18] initial_state(method::Optim.BFGS{LineSearches.InitialStatic{Float64}, LineSearches.BackTracking{Float64, Int64}, Nothing, Float64, Optim.Flat}, options::Optim.Options{Float64, Pumas.var"#551#553"{Pumas.DefaultOptimizeFN{Nothing, NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases), Tuple{Bool, Bool, Bool, Float64, Bool}}}, Pumas.var"#582#584"{Pumas.FOCE, Bool, NLSolversBase.OnceDifferentiable{Float64, Vector{Float64}, Vector{Float64}}, TransformVariables.TransformTuple{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, TransformVariables.ShiftedExp{true, Float64}, Pumas.PDiagTransform, Pumas.PDiagTransform, TransformVariables.ShiftedExp{true, Float64}}}}, Vector{Vector{Float64}}, Vector{Vector{Float64}}}}}, d::NLSolversBase.OnceDifferentiable{Float64, Vector{Float64}, Vector{Float64}}, initial_x::Vector{Float64})
    @ Optim /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Optim/rpjtl/src/multivariate/solvers/first_order/bfgs.jl:94
 [19] optimize
    @ /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Optim/rpjtl/src/multivariate/optimize/optimize.jl:36 [inlined]
 [20] DefaultOptimizeFN
    @ /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:2902 [inlined]
 [21] _fit(m::PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, population::Vector{Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}, param::NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{Float64, Float64, Float64, Float64, Float64, Diagonal{Float64, Vector{Float64}}, Diagonal{Float64, Vector{Float64}}, Float64}}, approx::Pumas.FOCE, ensemblealg::EnsembleThreads, optimize_fn::Pumas.DefaultOptimizeFN{Nothing, NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases), Tuple{Bool, Bool, Bool, Float64, Bool}}}, fixedparamset::ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, fixedparam::NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{Float64, Float64, Float64, Float64, Float64, Diagonal{Float64, Vector{Float64}}, Diagonal{Float64, Vector{Float64}}, Float64}}, checkidentification::Bool, diffeq_options::NamedTuple{(), Tuple{}})


Could you help on what does this error mean ?

Thanks

    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:3785
 [22] __fit
    @ /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:3574 [inlined]
 [23] fit(m::PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, _population::Vector{Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}, param::NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :σ_prop, :Ω_1, :Ω_2), Tuple{Float64, Float64, Float64, Float64, Float64, Float64, Diagonal{Float64, Vector{Float64}}, Diagonal{Float64, Vector{Float64}}}}, approx::Pumas.FOCE; optimize_fn::Pumas.DefaultOptimizeFN{Nothing, NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases), Tuple{Bool, Bool, Bool, Float64, Bool}}}, constantcoef::NamedTuple{(), Tuple{}}, omegas::Tuple{}, ensemblealg::EnsembleThreads, checkidentification::Bool, diffeq_options::NamedTuple{(), Tuple{}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:3542
 [24] fit(m::PumasModel{(tvPOP_CIRC0 = 1, tvPOP_MTT = 1, tvPOP_GAMMA = 1, tvPOP_SLOP = 1, tvPOP_SLOPC = 1, Ω_1 = 3, Ω_2 = 1, σ_prop = 1), 4, (:Central, :PERIPH1, :PERIPH2, :central_carb, :Periph_carb, :Prol, :Transit1, :Transit2, :Transit3, :Circ), ParamSet{NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :Ω_1, :Ω_2, :σ_prop), Tuple{RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, PDiagDomain{PDMats.PDiagMat{Float64, Vector{Float64}}}, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#63#72", Pumas.TimeDispatcher{var"#64#73", var"#65#74"}, Nothing, var"#67#76", ODEProblem{Nothing, Tuple{Nothing, Nothing}, true, Nothing, ODEFunction{true, ModelingToolkit.ODEFunctionClosure{var"#68#77", var"#69#78"}, LinearAlgebra.UniformScaling{Bool}, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Nothing, Vector{Symbol}, Symbol, typeof(SciMLBase.DEFAULT_OBSERVED), Nothing}, Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}}, SciMLBase.StandardODEProblem}, var"#70#79", var"#71#80", ModelingToolkit.ODESystem}, _population::Vector{Subject{NamedTuple{(:ANC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantCovar{NamedTuple{(:BSA, :BILIRUBIN, :SEX, :AGE, :Regimen), Tuple{Float64, Float64, Int64, Int64, String}}}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}, param::NamedTuple{(:tvPOP_CIRC0, :tvPOP_MTT, :tvPOP_GAMMA, :tvPOP_SLOP, :tvPOP_SLOPC, :σ_prop, :Ω_1, :Ω_2), Tuple{Float64, Float64, Float64, Float64, Float64, Float64, Diagonal{Float64, Vector{Float64}}, Diagonal{Float64, Vector{Float64}}}}, approx::Pumas.FOCE)
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/MxXdQ/src/estimation/likelihoods.jl:3521
 [25] top-level scope
    @ ~/data/code/TaxaneFinal/Friburg_updated_clean.jmd:4331

Hi Ahmed,

Most likely, you have a variable in the @vars block that doesn’t depend on ODE variables. Try moving the definition to the @pre block instead. That should resolve the issue.

Best
Andreas

1 Like

Hi Andreas,

Thanks for the reply. Yes that’s true. I have variables that do not depend on ODE variables. The code works now after moving them to @pre block.

Just a note for you that could be helpful from the software development point of view that when I changed the estimation method from FOCEI to laplacian (LaplaceI) the old code worked.

Thanks for letting us know. That is indeed useful to know. It suggests that it might be possible to make this case work. If not, we should improve the error message.

1 Like