Hello,
I met this error when fitting the dataset using a two-compartment model. The error occurred when I fix the intercompartment clearance, Q, meaning I do not estimate the variability for Q. However, the same model can run through if I do estimate Ω for Q .
The error message is shown below. Could anyone help me out to interpret the error message and how to fix it? Thank you very much.
BoundsError: attempt to access 0-element Array{Float64,1} at index [1]
getindex at array.jl:728 [inlined]
#solve_tr_subproblem!#36(::Float64, ::Int64, ::typeof(Optim.solve_tr_subproblem!), ::Array{Float64,1}, ::Array{Float64,2}, ::Float64, ::Array{Float64,1}) at newton_trust_region.jl:105
solve_tr_subproblem! at newton_trust_region.jl:94 [inlined]
update_state!(::NLSolversBase.TwiceDifferentiable{Float64,Array{Float64,1},Array{Float64,2},Array{Float64,1}}, ::Optim.NewtonTrustRegionState{Array{Float64,1},Float64,Array{Float64,1}}, ::Optim.NewtonTrustRegion{Float64}) at newton_trust_region.jl:314
optimize(::NLSolversBase.TwiceDifferentiable{Float64,Array{Float64,1},Array{Float64,2},Array{Float64,1}}, ::Array{Float64,1}, ::Optim.NewtonTrustRegion{Float64}, ::Optim.Options{Float64,Nothing}, ::Optim.NewtonTrustRegionState{Array{Float64,1},Float64,Array{Float64,1}}) at optimize.jl:57
optimize at optimize.jl:33 [inlined]
#_orth_empirical_bayes!#145(::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}}, ::typeof(Pumas._orth_empirical_bayes!), ::Array{Float64,1}, ::PumasModel{ParamSet{NamedTuple{(:tvcl, :tvQ, :tvv, :tvvt, :Ω, :σ²_add, :σ²_prop),Tuple{RealDomain{Int64,TransformVariables.Infinity{true},Float64},RealDomain{Int64,TransformVariables.Infinity{true},Float64},RealDomain{Int64,TransformVariables.Infinity{true},Float64},RealDomain{Int64,TransformVariables.Infinity{true},Float64},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},RealDomain{Float64,TransformVariables.Infinity{true},Float64},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##49#54")),getfield(Main, Symbol("##50#55")),getfield(Main, Symbol("##51#56")),Central1Periph1,getfield(Main, Symbol("##52#57")),getfield(Main, Symbol("##53#58"))}, ::Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},NamedTuple{(:WT, :CRCLM_final, :ALB),Tuple{Float64,Float64,Float64}},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}}, ::NamedTuple{(:tvcl, :tvQ, :tvv, :tvvt, :Ω, :σ²_add, :σ²_prop),Tuple{Float64,Float64,Float64,Float64,PDMats.PDiagMat{Float64,Array{Float64,1}},Float64,Float64}}, ::Pumas.FOCEI) at likelihoods.jl:248
_orth_empirical_bayes! at likelihoods.jl:221 [inlined]
#219 at likelihoods.jl:1079 [inlined]
mapfoldl_impl(::getfield(Pumas, Symbol("##219#225")){NamedTuple{(:tvcl, :tvQ, :tvv, :tvvt, :Ω, :σ²_add, :σ²_prop),Tuple{Float64,Float64,Float64,Float64,PDMats.PDiagMat{Float64,Array{Float64,1}},Float64,Float64}},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},PumasModel{ParamSet{NamedTuple{(:tvcl, :tvQ, :tvv, :tvvt, :Ω, :σ²_add, :σ²_prop),Tuple{RealDomain{Int64,TransformVariables.Infinity{true},Float64},RealDomain{Int64,TransformVariables.Infinity{true},Float64},RealDomain{Int64,TransformVariables.Infinity{true},Float64},RealDomain{Int64,TransformVariables.Infinity{true},Float64},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},RealDomain{Float64,TransformVariables.Infinity{true},Float64},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##49#54")),getfield(Main, Symbol("##50#55")),getfield(Main, Symbol("##51#56")),Central1Periph1,getfield(Main, Symbol("##52#57")),getfield(Main, Symbol("##53#58"))},Pumas.FOCEI,Tuple{}}, ::typeof(Base.add_sum), ::NamedTuple{(:init,),Tuple{Float64}}, ::Base.Iterators.Zip{Tuple{Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},NamedTuple{(:WT, :CRCLM_final, :ALB),Tuple{Float64,Float64,Float64}},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1},Array{Array{Float64,1},1},Array{Array{Float64,1},1}}}, ::Tuple{Int64,Int64,Int64}) at reduce.jl:49
mapfoldl_impl(::Function, ::Function, ::NamedTuple{,Tuple{}}, ::Base.Iterators.Zip{Tuple{Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},NamedTuple{(:WT, :CRCLM_final, :ALB),Tuple{Float64,Float64,Float64}},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1},Array{Array{Float64,1},1},Array{Array{Float64,1},1}}}) at reduce.jl:61
#mapfoldl#188 at reduce.jl:72 [inlined]
mapfoldl at reduce.jl:72 [inlined]
#mapreduce#192 at reduce.jl:208 [inlined]
mapreduce at reduce.jl:208 [inlined]
sum(::Function, ::Base.Iterators.Zip{Tuple{Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},NamedTuple{(:WT, :CRCLM_final, :ALB),Tuple{Float64,Float64,Float64}},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1},Array{Array{Float64,1},1},Array{Array{Float64,1},1}}}) at reduce.jl:403
(::getfield(Pumas, Symbol("##218#224")){Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},PumasModel{ParamSet{NamedTuple{(:tvcl, :tvQ, :tvv, :tvvt, :Ω, :σ²_add, :σ²_prop),Tuple{RealDomain{Int64,TransformVariables.Infinity{true},Float64},RealDomain{Int64,TransformVariables.Infinity{true},Float64},RealDomain{Int64,T...