Fitting Error with Pumas 2.1

Estimation
1

Hello, I’m getting an error with fitting with FOCE following the transition to Pumas 2.1. Does anyone know what may cause this error? I pasted part of the error below.

ERROR: MethodError: no method matching unique_times_sorted(::Vector{Pumas.Event{T1, T2, Float64, T4, Float64, Float64, Int64} where {T1, T2, T4}}, ::Float64)
Closest candidates are:
  unique_times_sorted(::AbstractVector{E}, ::Any) where {T1, T2, E<:(Pumas.Event{T1, T2})} at /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/Ozwn0/src/simulate_methods/utils.jl:27
Stacktrace:
  [1] _build_analytical_problem(prob::Pumas.LinearODE, u0::LabelledArrays.LArray{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, 1, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, (:Depot, :Transit1, :Transit2, :Transit3, :Transit4, :Transit5, :rispC, :rispP)}, tspan::Tuple{Float64, Float64}, events::Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, covariates::Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, col::Pumas.Pre{Nothing, var"#71#99"{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}}}, NamedTuple{(:ηka1, :ηka2, :ηkrel, :ηV, :ηktr, :ηkrrp, :ηkrpr, :κka1, :κka2, :κkrel, :κV, :κktr, :κkrrp, :κkrpr), Tuple{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}}}, Subject{NamedTuple{(:dv_rispC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}}, dcp::Pumas.Pre{Nothing, var"#73#101"{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}}}, NamedTuple{(:ηka1, :ηka2, :ηkrel, :ηV, :ηktr, :ηkrrp, :ηkrpr, :κka1, :κka2, :κkrel, :κV, :κktr, :κkrrp, :κkrpr), Tuple{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}}}, Subject{NamedTuple{(:dv_rispC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/Ozwn0/src/simulate_methods/analytical.jl:174
  [2] #_problem#1332
    @ /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/Ozwn0/src/simulate_methods/analytical.jl:41 [inlined]
  [3] _problem(model::PumasModel{(tvtr_ratio = 1, tvka1 = 1, tvka2 = 1, tvkrel = 1, tvV = 1, tvktr = 1, tvkrrp = 1, tvkrpr = 1, ωka1_BSV = 1, ωka2_BSV = 1, ωkrel_BSV = 1, ωV_BSV = 1, ωktr_BSV = 1, ωkrrp_BSV = 1, ωkrpr_BSV = 1, ω²ka1_BOV = 1, ω²ka2_BOV = 1, ω²krel_BOV = 1, ω²V_BOV = 1, ω²ktr_BOV = 1, ω²krrp_BOV = 1, ω²krpr_BOV = 1, σ_add = 1, σ_prop = 1), 49, ParamSet{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#69#97", var"#70#98", var"#72#100", var"#74#102", Pumas.LinearODE, var"#75#103", var"#86#114", ModelingToolkit.ODESystem}, subject::Subject{NamedTuple{(:dv_rispC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, col::Pumas.Pre{Nothing, var"#71#99"{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}}}, NamedTuple{(:ηka1, :ηka2, :ηkrel, :ηV, :ηktr, :ηkrrp, :ηkrpr, :κka1, :κka2, :κkrel, :κV, :κktr, :κkrrp, :κkrpr), Tuple{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}}}, Subject{NamedTuple{(:dv_rispC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}}, dcp::Pumas.Pre{Nothing, var"#73#101"{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}}}, NamedTuple{(:ηka1, :ηka2, :ηkrel, :ηV, :ηktr, :ηkrrp, :ηkrpr, :κka1, :κka2, :κkrel, :κV, :κktr, :κkrrp, :κkrpr), Tuple{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag2, Float64}, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}, 49}}}}, Subject{NamedTuple{(:dv_rispC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}}}; diffeq_options::NamedTuple{(:abstol, :reltol, :alg, :saveat), Tuple{Float64, Float64, CompositeAlgorithm{Tuple{Vern7, Rodas5{0, true, DefaultLinSolve, Val{:forward}, true}}, AutoSwitch{Vern7, Rodas5{0, true, DefaultLinSolve, Val{:forward}, true}, Rational{Int64}, Int64}}, Vector{Float64}}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/Ozwn0/src/models/model_api.jl:828
  [4] #_derived#266
    @ /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/Ozwn0/src/models/model_api.jl:871 [inlined]
  [5] _derived_vηorth_gradient(m::PumasModel{(tvtr_ratio = 1, tvka1 = 1, tvka2 = 1, tvkrel = 1, tvV = 1, tvktr = 1, tvkrrp = 1, tvkrpr = 1, ωka1_BSV = 1, ωka2_BSV = 1, ωkrel_BSV = 1, ωV_BSV = 1, ωktr_BSV = 1, ωkrrp_BSV = 1, ωkrpr_BSV = 1, ω²ka1_BOV = 1, ω²ka2_BOV = 1, ω²krel_BOV = 1, ω²V_BOV = 1, ω²ktr_BOV = 1, ω²krrp_BOV = 1, ω²krpr_BOV = 1, σ_add = 1, σ_prop = 1), 49, ParamSet{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#69#97", var"#70#98", var"#72#100", var"#74#102", Pumas.LinearODE, var"#75#103", var"#86#114", ModelingToolkit.ODESystem}, subject::Subject{NamedTuple{(:dv_rispC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, param::NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}}}, vrandeffsorth::Vector{Float64}, diffeq_options::NamedTuple{(), Tuple{}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/Ozwn0/src/estimation/likelihoods.jl:1830
  [6] _∂²l∂η²(m::PumasModel{(tvtr_ratio = 1, tvka1 = 1, tvka2 = 1, tvkrel = 1, tvV = 1, tvktr = 1, tvkrrp = 1, tvkrpr = 1, ωka1_BSV = 1, ωka2_BSV = 1, ωkrel_BSV = 1, ωV_BSV = 1, ωktr_BSV = 1, ωkrrp_BSV = 1, ωkrpr_BSV = 1, ω²ka1_BOV = 1, ω²ka2_BOV = 1, ω²krel_BOV = 1, ω²V_BOV = 1, ω²ktr_BOV = 1, ω²krrp_BOV = 1, ω²krpr_BOV = 1, σ_add = 1, σ_prop = 1), 49, ParamSet{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#69#97", var"#70#98", var"#72#100", var"#74#102", Pumas.LinearODE, var"#75#103", var"#86#114", ModelingToolkit.ODESystem}, subject::Subject{NamedTuple{(:dv_rispC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, param::NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}}}, vrandeffsorth::Vector{Float64}, approx::Pumas.FOCE, diffeq_options::NamedTuple{(), Tuple{}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/Ozwn0/src/estimation/likelihoods.jl:1870
  [7] _marginal_nll(model::PumasModel{(tvtr_ratio = 1, tvka1 = 1, tvka2 = 1, tvkrel = 1, tvV = 1, tvktr = 1, tvkrrp = 1, tvkrpr = 1, ωka1_BSV = 1, ωka2_BSV = 1, ωkrel_BSV = 1, ωV_BSV = 1, ωktr_BSV = 1, ωkrrp_BSV = 1, ωkrpr_BSV = 1, ω²ka1_BOV = 1, ω²ka2_BOV = 1, ω²krel_BOV = 1, ω²V_BOV = 1, ω²ktr_BOV = 1, ω²krrp_BOV = 1, ω²krpr_BOV = 1, σ_add = 1, σ_prop = 1), 49, ParamSet{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#69#97", var"#70#98", var"#72#100", var"#74#102", Pumas.LinearODE, var"#75#103", var"#86#114", ModelingToolkit.ODESystem}, subject::Subject{NamedTuple{(:dv_rispC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, param::NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}}}, vrandeffsorth::Vector{Float64}, approx::Pumas.FOCE, diffeq_options::NamedTuple{(), Tuple{}})
    @ Pumas /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/Ozwn0/src/estimation/likelihoods.jl:1112
  [8] #465
    @ /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/Ozwn0/src/estimation/likelihoods.jl:1480 [inlined]
  [9] chunk_mode_gradient(f::Pumas.var"#465#466"{PumasModel{(tvtr_ratio = 1, tvka1 = 1, tvka2 = 1, tvkrel = 1, tvV = 1, tvktr = 1, tvkrrp = 1, tvkrpr = 1, ωka1_BSV = 1, ωka2_BSV = 1, ωkrel_BSV = 1, ωV_BSV = 1, ωktr_BSV = 1, ωkrrp_BSV = 1, ωkrpr_BSV = 1, ω²ka1_BOV = 1, ω²ka2_BOV = 1, ω²krel_BOV = 1, ω²V_BOV = 1, ω²ktr_BOV = 1, ω²krrp_BOV = 1, ω²krpr_BOV = 1, σ_add = 1, σ_prop = 1), 49, ParamSet{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#69#97", var"#70#98", var"#72#100", var"#74#102", Pumas.LinearODE, var"#75#103", var"#86#114", ModelingToolkit.ODESystem}, Subject{NamedTuple{(:dv_rispC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, Vector{Float64}, Pumas.FOCE, TransformVariables.TransformTuple{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, TransformVariables.ShiftedExp{true, Float64}}}}, NamedTuple{(), Tuple{}}}, x::Vector{Float64}, cfg::ForwardDiff.GradientConfig{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1, Vector{ForwardDiff.Dual{ForwardDiff.Tag{Pumas.Tag1, Float64}, Float64, 1}}})
    @ ForwardDiff /builds/PumasAI/PumasSystemImages-jl/.julia/packages/ForwardDiff/PBzup/src/gradient.jl:150
 [10] gradient
    @ /builds/PumasAI/PumasSystemImages-jl/.julia/packages/ForwardDiff/PBzup/src/gradient.jl:21 [inlined]
 [11] gradient
    @ /builds/PumasAI/PumasSystemImages-jl/.julia/packages/ForwardDiff/PBzup/src/gradient.jl:17 [inlined]
 [12] _∂ℓᵐ∂θ
    @ /builds/PumasAI/PumasSystemImages-jl/.julia/packages/Pumas/Ozwn0/src/estimation/likelihoods.jl:1488 [inlined]
 [13] _marginal_nll_gradient!(g::Vector{Float64}, model::PumasModel{(tvtr_ratio = 1, tvka1 = 1, tvka2 = 1, tvkrel = 1, tvV = 1, tvktr = 1, tvkrrp = 1, tvkrpr = 1, ωka1_BSV = 1, ωka2_BSV = 1, ωkrel_BSV = 1, ωV_BSV = 1, ωktr_BSV = 1, ωkrrp_BSV = 1, ωkrpr_BSV = 1, ω²ka1_BOV = 1, ω²ka2_BOV = 1, ω²krel_BOV = 1, ω²V_BOV = 1, ω²ktr_BOV = 1, ω²krrp_BOV = 1, ω²krpr_BOV = 1, σ_add = 1, σ_prop = 1), 49, ParamSet{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, RealDomain{Float64, TransformVariables.Infinity{true}, Float64}}}}, var"#69#97", var"#70#98", var"#72#100", var"#74#102", Pumas.LinearODE, var"#75#103", var"#86#114", ModelingToolkit.ODESystem}, subject::Subject{NamedTuple{(:dv_rispC,), Tuple{Vector{Union{Missing, Float64}}}}, Pumas.ConstantInterpolationStructArray{Vector{Float64}, StructArrays.StructVector{NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{Int64, Int64, Int64, Int64, String, Int64, String, Int64, Float64}}, NamedTuple{(:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio), Tuple{FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, Vector{Int64}, Vector{Int64}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{String, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Int64, 1, Tuple{Base.OneTo{Int64}}}, FillArrays.Fill{Float64, 1, Tuple{Base.OneTo{Int64}}}}}, Int64}, Symbol}, Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}, Vector{Float64}}, vparam::Vector{Float64}, vrandeffsorth::Vector{Float64}, approx::Pumas.FOCE, trf::TransformVariables.TransformTuple{NamedTuple{(:tvtr_ratio, :tvka1, :tvka2, :tvkrel, :tvV, :tvktr, :tvkrrp, :tvkrpr, :ωka1_BSV, :ωka2_BSV, :ωkrel_BSV, :ωV_BSV, :ωktr_BSV, :ωkrrp_BSV, :ωkrpr_BSV, :ω²ka1_BOV, :ω²ka2_BOV, :ω²krel_BOV, :ω²V_BOV, :ω²ktr_BOV, :ω²krrp_BOV, :ω²krpr_BOV, :σ_add, :σ_prop), NTuple{24, TransformVariables.ShiftedExp{true, Float64}}}}, diffeq_options::NamedTuple{(), Tuple{}})
2

This can be considered a bug but I’m surprised that you’ve ended up hitting it. I’d be interested in knowing how your event vector has been created and how it looks like.

3

What do you mean by event vector for fitting? I’m only familiar with that for simulations.

Below is a snippet of the data being read for 1 subject, if that helps.

id	time	cp_rispC	dv_rispC	amt	evid	cmt	rate	ab_NSubj	iii	occ	period	sequence	seq_n	formulation	isT	ab_rep	route	input_tr_ratio
1	0			1.20E+08	1	1	0	24	672	1	0	R	1	R	0	1	ev	0.8
1	669.6	2.080088146	1.954009289		0		0	24	672	1	0	R	1	R	0	1	ev	0.8
1	672			1.20E+08	1	1	0	24	672	2	0	R	1	R	0	1	ev	0.8
1	673	3.235381981	4.484136298		0		0	24	672	2	0	R	1	R	0	1	ev	0.8
1	674	4.816475708	3.472602778		0		0	24	672	2	0	R	1	R	0	1	ev	0.8
1	675	5.85096083	2.751162352		0		0	24	672	2	0	R	1	R	0	1	ev	0.8
1	676	6.519690401	3.17529513		0		0	24	672	2	0	R	1	R	0	1	ev	0.8
1	678	7.204335389	7.293956245		0		0	24	672	2	0	R	1	R	0	1	ev	0.8
1	680	7.434155329	7.52093271		0		0	24	672	2	0	R	1	R	0	1	ev	0.8
1	684	7.401034967	11.17459771		0		0	24	672	2	0	R	1	R	0	1	ev	0.8
1	792	11.37795264	9.156762942		0		0	24	672	2	0	R	1	R	0	1	ev	0.8
1	984	14.17490782	17.52771548		0		0	24	672	2	0	R	1	R	0	1	ev	0.8
1	1341.6	0.857686897	0.471644758		0		0	24	672	2	0	R	1	R	0	1	ev	0.8
1	1344			1.20E+08	1	1	0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1344.01	0.673314662	0.756680178		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1345	3.700926326	4.429653848		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1346	5.681894551	6.404575023		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1347	6.961907863	6.53283577		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1348	7.78802078	7.145519556		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1350	8.659063985	7.612304123		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1352	9.004468356	9.409678086		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1356	9.133882314	11.06442957		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1368	8.736108666	6.364671985		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1416	12.85764127	17.41568415		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1464	24.14330165	21.64670178		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1512	28.06632868	29.37534411		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1560	24.34758762	18.86051545		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1608	18.25638238	13.12160606		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1656	12.78742282	15.35400217		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1704	8.695304511	10.37775034		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1752	5.842714026	8.106173498		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	1824	3.194635647	4.479504925		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	2013.6	0.647220616	0.432121245		0		0	24	672	3	0	R	1	R	0	1	ev	0.8
1	2016			1.20E+08	1	1	0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2016.01	1.883535154	2.235690785		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2017	4.084963131	4.813674779		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2018	5.554552965	4.808905669		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2019	6.530653388	4.959932862		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2020	7.184790666	7.46639685		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2022	7.931571591	5.730254662		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2024	8.290954644	9.322676354		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2028	8.566600018	2.399716749		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2088	13.65877684	4.162259947		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2136	26.26853545	40.4375542		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2184	31.73679423	18.93940135		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2232	29.12368476	20.71986468		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2328	17.92153821	14.94544852		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2424	9.918042142	12.1114475		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2496	6.301593243	8.496131221		0		0	24	672	4	1	R	1	R	0	1	ev	0.8
1	2685.6	1.902834598	2.239710654		0		0	24	672	4	1	R	1	R	0	1	ev	0.8

And here is what I ran for read_pumas:

i_BE_sim_pop = read_pumas(data,
                            id                      =   :id,
                            time                    =   :time,
                            amt                     =   :amt,
                            cmt                     =   :cmt,
                            ii                      =   :iii,                        
                            observations            =   [:dv_rispC],
                            route                   =   :route,
                            covariates              =   [:ab_rep, :ab_NSubj, :occ, :period, :sequence, :seq_n, :formulation, :isT, :input_tr_ratio],
                            covariates_direction    =   :left,
                            evid                    =   :evid)

And the fitting:

i_reest_fit =   @time fit(runin_pl_abbrev_BE_ML_model,
                            i_BE_sim_pop,
                            reest_params,
                            Pumas.FOCE(),)
                            ensemblealg = EnsembleThreads())
4

I also encountered an error with I believe it was read_pumas, in which it said I wasn’t allowed to have ‘ii’ as the name of a covariate. Is this a bug as well, or is this intended?

5

I’d be interested in how i_BE_sim_pop[1].events look like. Does it say where {T1,T2,T4} when printed. If not could you try to find a subject that looks like that?

6

This is intended behavior. The ii name is a reserved name so you’d have to create a copy of the column under a different name if you’d like to use it as a covariate.

7

I’m getting the following:

6-element Vector{Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}}:
 Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}(1.2e8, 0.0, 1, 0.0, Inf, 0.0, 0.0, NCA.EV, 1, 0, 1)
 Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}(1.2e8, 672.0, 1, 0.0, Inf, 0.0, 672.0, NCA.EV, 1, 0, 1)
 Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}(1.2e8, 1344.0, 1, 0.0, Inf, 0.0, 1344.0, NCA.EV, 1, 0, 1)
 Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}(1.2e8, 2016.0, 1, 0.0, Inf, 0.0, 2016.0, NCA.EV, 1, 0, 1)
 Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}(1.2e8, 2688.0, 1, 0.0, Inf, 0.0, 2688.0, NCA.EV, 1, 0, 1)
 Pumas.Event{Float64, Float64, Float64, Float64, Float64, Float64, Int64}(1.2e8, 3360.0, 1, 0.0, Inf, 0.0, 3360.0, NCA.EV, 1, 0, 1)
8

Okay. It can’t be that subject that causes the issue. Maybe try

filter(t -> eltype(t.events) != eltype(i_BE_sim_pop[1].events), i_BE_sim_pop)

and see what shows up.

9 
Population
  Subjects: 0
10

Do you have a dose control block in your model?

11

Yes, I am using it for bioavailability.

12

Would you be able to show us the the dose control block?

man. 16. maj 2022 kl. 15.40 skrev Donald Lee via Pumas <notifications@pumas.discoursemail.com>:

13 
@dosecontrol begin
        bioav       =   (Depot = ((isT == 1) & (occ > 3)) ? tvtr_ratio : 1,)
    end
14

Thanks. We’ll try to debug this issue. Meanwhile, there is a chance that you can avoid the error by changing the dosecontrol block to

@dosecontrol begin
    bioav       =   (Depot = ((isT == 1) & (occ > 3)) ? tvtr_ratio : one(tvtr_ratio),)
end