ERROR: ArgumentError: quantiles are undefined in presence of NaNs

sai_matcha · January 22, 2021, 12:31pm

I tried Bootstrap with below code

BS1 = infer(R130_FINALMODEL_result_p,Pumas.Bootstrap(samples = 1000))

it gives me below error

julia> BS1
Bootstrap inference results

Successful minimization:                     false

Likelihood approximation:              Pumas.FOCEI
Log-likelihood value:                   -419.04306
Number of subjects:                            134
Number of parameters:         Fixed      Optimized
                                  1             11
Observation records:         Active        Missing
    DV:                         356             49
    Total:                      356             49

Error showing value of type Pumas.FittedPumasModelInference{Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}},Serialization.__deserialized_types__.var"#533#563",Serialization.__deserialized_types__.var"#534#564",Serialization.__deserialized_types__.var"#536#566",Pumas.LinearODE,Serialization.__deserialized_types__.var"#537#567",Serialization.__deserialized_types__.var"#550#580"},Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantInterpolationStructArray{Array{Float64,1},StructArrays.StructArray{NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Float64,Float64,Float64,Int64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64}},1,NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Int64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1}}},Int64},Symbol},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Nothing,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat},Float64,Array{Float64,1},Float64,Float64,Array{Optim.OptimizationState{Float64,Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.FOCEI,Array{Array{Float64,1},1},NamedTuple{(:optimize_fn, :constantcoef, :omegas, :ensemblealg),Tuple{Pumas.DefaultOptimizeFN{Nothing,NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases),Tuple{Bool,Bool,Bool,Float64,Bool}}},NamedTuple{(:tvbio,),Tuple{Float64}},Tuple{},EnsembleThreads}},ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},Pumas.ConstDomain{Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}}},Pumas.Bootstraps,Float64}:
ERROR: ArgumentError: quantiles are undefined in presence of NaNs
Stacktrace:
 [1] _quantilesort!(::Array{Float64,1}, ::Bool, ::Float64, ::Float64) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\Statistics\src\Statistics.jl:960
 [2] #quantile!#52 at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\Statistics\src\Statistics.jl:943 [inlined]
 [3] #quantile#55 at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\Statistics\src\Statistics.jl:1052 [inlined]
 [4] quantile at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\Statistics\src\Statistics.jl:1052 [inlined]
 [5] _push_varinfo!(::Array{String,1}, ::Array{Float64,1}, ::Array{Float64,1}, ::Array{Tuple{Float64,Float64},1}, ::Symbol, ::Float64, ::Float64, ::Float64, ::Array{Float64,1}) at C:\Users\julia\.julia\packages\Pumas\iIBBr\src\estimation\show.jl:203
 [6] coeftable(::Pumas.FittedPumasModelInference{Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}},Serialization.__deserialized_types__.var"#533#563",Serialization.__deserialized_types__.var"#534#564",Serialization.__deserialized_types__.var"#536#566",Pumas.LinearODE,Serialization.__deserialized_types__.var"#537#567",Serialization.__deserialized_types__.var"#550#580"},Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantInterpolationStructArray{Array{Float64,1},StructArrays.StructArray{NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Float64,Float64,Float64,Int64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64}},1,NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Int64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1}}},Int64},Symbol},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Nothing,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat},Float64,Array{Float64,1},Float64,Float64,Array{Optim.OptimizationState{Float64,Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.FOCEI,Array{Array{Float64,1},1},NamedTuple{(:optimize_fn, :constantcoef, :omegas, :ensemblealg),Tuple{Pumas.DefaultOptimizeFN{Nothing,NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases),Tuple{Bool,Bool,Bool,Float64,Bool}}},NamedTuple{(:tvbio,),Tuple{Float64}},Tuple{},EnsembleThreads}},ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},Pumas.ConstDomain{Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}}},Pumas.Bootstraps,Float64}) at C:\Users\julia\.julia\packages\Pumas\iIBBr\src\estimation\show.jl:235
 [7] show(::IOContext{REPL.Terminals.TTYTerminal}, ::MIME{Symbol("text/plain")}, ::Pumas.FittedPumasModelInference{Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}},Serialization.__deserialized_types__.var"#533#563",Serialization.__deserialized_types__.var"#534#564",Serialization.__deserialized_types__.var"#536#566",Pumas.LinearODE,Serialization.__deserialized_types__.var"#537#567",Serialization.__deserialized_types__.var"#550#580"},Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantInterpolationStructArray{Array{Float64,1},StructArrays.StructArray{NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Float64,Float64,Float64,Int64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64}},1,NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, 
:SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Int64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1}}},Int64},Symbol},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Nothing,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat},Float64,Array{Float64,1},Float64,Float64,Array{Optim.OptimizationState{Float64,Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.FOCEI,Array{Array{Float64,1},1},NamedTuple{(:optimize_fn, :constantcoef, :omegas, :ensemblealg),Tuple{Pumas.DefaultOptimizeFN{Nothing,NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases),Tuple{Bool,Bool,Bool,Float64,Bool}}},NamedTuple{(:tvbio,),Tuple{Float64}},Tuple{},EnsembleThreads}},ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},Pumas.ConstDomain{Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}}},Pumas.Bootstraps,Float64}) at C:\Users\julia\.julia\packages\Pumas\iIBBr\src\estimation\show.jl:273
 [8] display(::REPL.REPLDisplay, ::MIME{Symbol("text/plain")}, ::Any) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\REPL\src\REPL.jl:214
 [9] display(::REPL.REPLDisplay, ::Any) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\REPL\src\REPL.jl:218
 [10] display(::Any) at .\multimedia.jl:328
 [11] #invokelatest#1 at .\essentials.jl:710 [inlined]
 [12] invokelatest at .\essentials.jl:709 [inlined]
 [13] print_response(::IO, ::Any, ::Bool, ::Bool, ::Any) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\REPL\src\REPL.jl:238
 [14] print_response(::REPL.AbstractREPL, ::Any, ::Bool, ::Bool) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\REPL\src\REPL.jl:223
 [15] (::REPL.var"#do_respond#54"{Bool,Bool,Atom.var"#246#247",REPL.LineEditREPL,REPL.LineEdit.Prompt})(::Any, ::Any, ::Any) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\REPL\src\REPL.jl:822     
 [16] #invokelatest#1 at .\essentials.jl:710 [inlined]
 [17] invokelatest at .\essentials.jl:709 [inlined]
 [18] run_interface(::REPL.Terminals.TextTerminal, ::REPL.LineEdit.ModalInterface, ::REPL.LineEdit.MIState) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\REPL\src\LineEdit.jl:2355
 [19] run_frontend(::REPL.LineEditREPL, ::REPL.REPLBackendRef) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\REPL\src\REPL.jl:1144
 [20] (::REPL.var"#38#42"{REPL.LineEditREPL,REPL.REPLBackendRef})() at .\task.jl:356

andreasnoack · January 22, 2021, 12:48pm

So it seems that some of the bootstrapped parameter values are NaNs. Do the parameter estimates look reasonable?

sai_matcha · January 22, 2021, 1:17pm

Yeah, they are reasonable compared to literature

julia> R130_FINALMODEL_result_p
FittedPumasModel

Successful minimization:                     false

Likelihood approximation:              Pumas.FOCEI
Log-likelihood value:                   -419.04306
Number of subjects:                            134
Number of parameters:         Fixed      Optimized
                                  1             11
Observation records:         Active        Missing
    DV:                         356             49
    Total:                      356             49

--------------------
          Estimate
--------------------
tvka       0.60938
tvk23      4.3642
tvvc     155.68
tvvp       0.49669
tvQ        3.6573
tvCL       3.7717
tvbio      0.63
θ₁         0.63565
θ₂         1.3489
θ₃        -0.59024
θ₄         5.6256
Ω1₁,₁      0.10003
Ω2₁,₁      0.23318
Ω2₂,₁      0.11298
Ω2₂,₂      0.40715
σ1         0.26796
σ2         0.089607
--------------------

andreasnoack · January 22, 2021, 4:21pm

Could you please try to take a look inside BS1.samples to see if any of the results there look completely off?

sai_matcha · January 23, 2021, 4:16am

R151_Bootstrap_1000 is serialization file of BS1 . here is my result.


julia> R151_Bootstrap_1000.fpm
FittedPumasModel

Successful minimization:                     false

Likelihood approximation:              Pumas.FOCEI
Log-likelihood value:                   -419.04306
Number of subjects:                            134
Number of parameters:         Fixed      Optimized
                                  1             11
Observation records:         Active        Missing
    DV:                         356             49
    Total:                      356             49

--------------------
          Estimate
--------------------
tvka       0.60938
tvk23      4.3642
tvvc     155.68
tvvp       0.49669
tvQ        3.6573
tvCL       3.7717
tvbio      0.63
θ₁         0.63565
θ₂         1.3489
θ₃        -0.59024
θ₄         5.6256
Ω1₁,₁      0.10003
Ω2₁,₁      0.23318
Ω2₂,₁      0.11298
Ω2₂,₂      0.40715
σ1         0.26796
σ2         0.089607
--------------------


julia> 

julia> R151_Bootstrap_1000.vcov
Bootstrap inference results

Successful fits: 932 out of 1000
No stratification.

julia> R151_Bootstrap_1000.level
0.95
julia> R151_Bootstrap_1000.samples
ERROR: type FittedPumasModelInference has no field samples
Stacktrace:
 [1] getproperty(::Pumas.FittedPumasModelInference{Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}},Serialization.__deserialized_types__.var"#533#563",Serialization.__deserialized_types__.var"#534#564",Serialization.__deserialized_types__.var"#536#566",Pumas.LinearODE,Serialization.__deserialized_types__.var"#537#567",Serialization.__deserialized_types__.var"#550#580"},Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantInterpolationStructArray{Array{Float64,1},StructArrays.StructArray{NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Float64,Float64,Float64,Int64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64}},1,NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Int64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1}}},Int64},Symbol},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Nothing,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat},Float64,Array{Float64,1},Float64,Float64,Array{Optim.OptimizationState{Float64,Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.FOCEI,Array{Array{Float64,1},1},NamedTuple{(:optimize_fn, :constantcoef, :omegas, :ensemblealg),Tuple{Pumas.DefaultOptimizeFN{Nothing,NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases),Tuple{Bool,Bool,Bool,Float64,Bool}}},NamedTuple{(:tvbio,),Tuple{Float64}},Tuple{},EnsembleThreads}},ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},Pumas.ConstDomain{Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}}},Pumas.Bootstraps,Float64}, ::Symbol) at .\Base.jl:33

I got only three types in BS1 result.

julia> R151_Bootstrap_1000.
fpm   level  vcov

andreasnoack · January 23, 2021, 8:10am

Sorry. It’s R151_Bootstrap_1000.vcov.samples

lør. 23. jan. 2021 kl. 05.27 skrev saikumar matcha via Pumas <pumas@discoursemail.com>:

sai_matcha · January 23, 2021, 10:31am

julia> R151_Bootstrap_1000.vcov.samples
ERROR: type Bootstraps has no field samples
Stacktrace:
 [1] getproperty(::Pumas.Bootstraps, ::Symbol) at .\Base.jl:33

andreasnoack · January 23, 2021, 10:51am

Should have been

R151_Bootstrap_1000.vcov.fits

sai_matcha · January 24, 2021, 9:43am

julia> R151_Bootstrap_1000.vcov.fits
1000-element Array{Any,1}:
 Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}},Serialization.__deserialized_types__.var"#533#563",Serialization.__deserialized_types__.var"#534#564",Serialization.__deserialized_types__.var"#536#566",Pumas.LinearODE,Serialization.__deserialized_types__.var"#537#567",Serialization.__deserialized_types__.var"#550#580"},Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantInterpolationStructArray{Array{Float64,1},StructArrays.StructArray{NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Float64,Float64,Float64,Int64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64}},1,NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Int64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1}}},Int64},Symbol},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Nothing,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat},Float64,Array{Float64,1},Float64,Float64,Array{Optim.OptimizationState{Float64,Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.FOCEI,Array{Array{Float64,1},1},NamedTuple{(:optimize_fn, :constantcoef, :omegas, :ensemblealg),Tuple{Pumas.DefaultOptimizeFN{Nothing,NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases),Tuple{Bool,Bool,Bool,Float64,Bool}}},NamedTuple{(:tvbio,),Tuple{Float64}},Tuple{},EnsembleThreads}},ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},Pumas.ConstDomain{Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}}}(PumasModel
  Parameters: tvka, tvk23, tvvc, tvvp, tvQ, tvCL, tvbio, θ, Ω1, Ω2, σ1, σ2
  Random effects: η1, η2
  Covariates: WT, AGEYRS, PMA, ALBD, SCRD, FFM
  Dynamical variables: Depot1, Prodrug, Central, Peripheral
  Derived: DV
  Observed: DV
, Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantInterpolationStructArray{Array{Float64,1},StructArrays.StructArray{NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Float64,Float64,Float64,Int64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64}},1,NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Int64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1}}},Int64},Symbol},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Nothing,Int64},1},Array{Float64,1}}[Subject
  ID: 141
  Events: 18
  Observables: DV: (n=3)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 132
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 115
  Events: 4
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 109
  Events: 6
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 9
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 27
  Events: 23
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 32
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 19
  Events: 60
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 81
  Events: 12
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 81
  Events: 12
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
  …  Subject
  ID: 145
  Events: 4
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 76
  Events: 22
  Observables: DV: (n=2)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 32
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 178
  Events: 32
  Observables: DV: (n=2)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 10
  Events: 284
  Observables: DV: (n=26)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 138
  Events: 146
  Observables: DV: (n=22)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 79
  Events: 10
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 44
  Events: 244
  Observables: DV: (n=30)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 58
  Events: 45
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 100
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
],  * Status: success

 * Candidate solution
    Final objective value:     4.883529e+02

 * Found with
    Algorithm:     BFGS

 * Convergence measures
    |x - x'|               = 5.02e+01 ≰ 0.0e+00
    |x - x'|/|x'|          = 1.87e-01 ≰ 0.0e+00
    |f(x) - f(x')|         = 6.30e-06 ≰ 0.0e+00
    |f(x) - f(x')|/|f(x')| = 1.29e-08 ≰ 0.0e+00
    |g(x)|                 = 6.56e-04 ≤ 1.0e-03

 * Work counters
    Seconds run:   60362  (vs limit Inf)
    Iterations:    59
    f(x) calls:    66
    ∇f(x) calls:   60
, Pumas.FOCEI(), [[0.0, -0.015650381923024877, 0.009838617598326157], [0.0, -0.432686580301894, -0.05005274548782389], [0.0, -0.047407067225383945, -0.034312958787973656], [0.0, 0.4233781317856444, 0.06906012192622052], [0.0, -0.7136180797168394, -0.3402191649861466], [-9.294381626030422e-5, -0.22069108414155655, -0.44103183725036926], [0.0, -0.09752964384166136, -0.04292447542042066], [0.0, 0.607359915545319, 1.2384802420890586], [0.0, -0.8470263217295175, -0.5855326583798955], [0.0, -0.8470263217295175, -0.5855326583798955]  …  [0.0, -0.347574558594557, -0.020980438701487265], [0.0, 0.23125835415338833, 0.177253459833323], [0.0, -0.09752964384166136, -0.04292447542042066], [0.0, -0.4020886538362563, -0.6212875156894225], [-0.0002615962488434546, 0.7659829569943182, 0.36356913732343676], [0.0, 0.09874257054205653, -0.01693923829294349], [0.0, -0.0017434243099462008, 0.1670536644742547], [0.0, -1.5178487764111648, 0.09124207401417189], [2.9446366668327095e-5, 0.2691225195308215, 0.5065726752024755], [0.0, 0.4041959953264322, 0.033807242334137594]], (optimize_fn = Pumas.DefaultOptimizeFN{Nothing,NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases),Tuple{Bool,Bool,Bool,Float64,Bool}}}(nothing, (show_trace = false, store_trace = true, extended_trace = false, g_tol = 0.001, allow_f_increases = true)), constantcoef = (tvbio = 0.63,), omegas = , ensemblealg = EnsembleThreads()), ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},Pumas.ConstDomain{Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}}((tvka = RealDomain{Float64,Int64,Float64}(0.1, 5, 2.55), tvk23 = RealDomain{Float64,Int64,Float64}(0.1, 20, 10.05), tvvc = RealDomain{Float64,Int64,Float64}(0.1, 250, 125.05), tvvp = RealDomain{Float64,Int64,Float64}(0.1, 250, 125.05), tvQ = RealDomain{Float64,Int64,Float64}(0.1, 250, 125.05), tvCL = RealDomain{Float64,Int64,Float64}(0.1, 15, 7.55), tvbio = Pumas.ConstDomain{Float64}(0.63), θ = VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}}([TransformVariables.Infinity{false}(), TransformVariables.Infinity{false}(), TransformVariables.Infinity{false}(), TransformVariables.Infinity{false}()], [TransformVariables.Infinity{true}(), TransformVariables.Infinity{true}(), TransformVariables.Infinity{true}(), TransformVariables.Infinity{true}()], [0.0, 0.0, 0.0, 0.0]), Ω1 = PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}}(PDMats.PDiagMat{Float64,Array{Float64,1}}(1, [1.0], [1.0])), Ω2 = PSDDomain{Array{Float64,2}}([1.0 0.0; 0.0 1.0]), σ1 = RealDomain{Float64,Int64,Float64}(0.01, 1, 0.505), σ2 = RealDomain{Float64,Int64,Float64}(0.01, 1, 0.505))))
 Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}},Serialization.__deserialized_types__.var"#533#563",Serialization.__deserialized_types__.var"#534#564",Serialization.__deserialized_types__.var"#536#566",Pumas.LinearODE,Serialization.__deserialized_types__.var"#537#567",Serialization.__deserialized_types__.var"#550#580"},Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantInterpolationStructArray{Array{Float64,1},StructArrays.StructArray{NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Float64,Float64,Float64,Int64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64}},1,NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Int64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1}}},Int64},Symbol},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Nothing,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat},Float64,Array{Float64,1},Float64,Float64,Array{Optim.OptimizationState{Float64,Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.FOCEI,Array{Array{Float64,1},1},NamedTuple{(:optimize_fn, :constantcoef, :omegas, :ensemblealg),Tuple{Pumas.DefaultOptimizeFN{Nothing,NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases),Tuple{Bool,Bool,Bool,Float64,Bool}}},NamedTuple{(:tvbio,),Tuple{Float64}},Tuple{},EnsembleThreads}},ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},Pumas.ConstDomain{Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}}}(PumasModel
  Parameters: tvka, tvk23, tvvc, tvvp, tvQ, tvCL, tvbio, θ, Ω1, Ω2, σ1, σ2
  Random effects: η1, η2
  Covariates: WT, AGEYRS, PMA, ALBD, SCRD, FFM
  Dynamical variables: Depot1, Prodrug, Central, Peripheral
  Derived: DV
  Observed: DV
, Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantInterpolationStructArray{Array{Float64,1},StructArrays.StructArray{NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Float64,Float64,Float64,Int64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64}},1,NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Int64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1}}},Int64},Symbol},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Nothing,Int64},1},Array{Float64,1}}[Subject
  ID: 141
  Events: 18
  Observables: DV: (n=3)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 140
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 24
  Events: 10
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 39
  Events: 18
  Observables: DV: (n=4)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 38
  Events: 3
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 137
  Events: 4
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 118
  Events: 60
  Observables: DV: (n=3)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 159
  Events: 4
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 85
  Events: 6
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 87
  Events: 26
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
  …  Subject
  ID: 32
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 119
  Events: 64
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 1
  Events: 27
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 11
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 76
  Events: 22
  Observables: DV: (n=2)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 75
  Events: 50
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 19
  Events: 60
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 150
  Events: 10
  Observables: DV: (n=2)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 167
  Events: 467
  Observables: DV: (n=9)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 142
  Events: 14
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
],  * Status: success

 * Candidate solution
    Final objective value:     3.662277e+02

 * Found with
    Algorithm:     BFGS

 * Convergence measures
    |x - x'|               = 1.07e-14 ≰ 0.0e+00
    |x - x'|/|x'|          = 3.61e-16 ≰ 0.0e+00
    |f(x) - f(x')|         = 0.00e+00 ≤ 0.0e+00
    |f(x) - f(x')|/|f(x')| = 0.00e+00 ≤ 0.0e+00
    |g(x)|                 = 2.82e-03 ≰ 1.0e-03

 * Work counters
    Seconds run:   306  (vs limit Inf)
    Iterations:    94
    f(x) calls:    274
    ∇f(x) calls:   95
, Pumas.FOCEI(), [[0.0, 0.020300959011351156, -0.7867710744179646], [0.0, 0.21067056231250836, 0.0201150099330132], [0.0, 1.1436539243475392, 0.5496052794755483], [0.0, 0.5539517209779848, 0.33161694656473967], [0.0, -0.2732292356004539, -0.023917439828655875], [0.0, -0.0744593522493071, -0.008920016230495976], [0.0, -0.9217794148000251, -0.9830046144540607], [0.0, 
0.14290193786418, 0.019153104557669432], [0.0, 0.3611048592378587, 0.02415500742680584], [0.0, -0.10015483619381738, -0.35430345151100967]  …  [0.0, -0.3046892612197491, -0.2287284304293414], [0.0, 0.2911609608623102, 0.9455987780005921], [-0.0035385777409184817, -0.3328968213765026, -1.1266935087556684], [0.0, 0.26808051112551556, 0.3900902117199193], [0.0, 0.2915722860844069, 0.2183919635078862], [0.0, 0.1159508473401361, 0.552486453937084], [0.0, 0.2636355235654636, 0.7010504799729986], [0.0, 0.4034262577418525, 0.29284839130536383], [0.0, 0.16802428031527789, 0.40831499211674965], [0.0, -1.0861130225156763, -1.3444853437819235]], (optimize_fn = Pumas.DefaultOptimizeFN{Nothing,NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases),Tuple{Bool,Bool,Bool,Float64,Bool}}}(nothing, (show_trace = false, store_trace = true, extended_trace = false, g_tol = 0.001, allow_f_increases = true)), constantcoef = (tvbio = 0.63,), omegas = , ensemblealg = EnsembleThreads()), ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},Pumas.ConstDomain{Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}}((tvka = RealDomain{Float64,Int64,Float64}(0.1, 5, 2.55), tvk23 = RealDomain{Float64,Int64,Float64}(0.1, 20, 10.05), tvvc = RealDomain{Float64,Int64,Float64}(0.1, 250, 125.05), tvvp = RealDomain{Float64,Int64,Float64}(0.1, 250, 125.05), tvQ = RealDomain{Float64,Int64,Float64}(0.1, 250, 125.05), tvCL = RealDomain{Float64,Int64,Float64}(0.1, 15, 7.55), tvbio = Pumas.ConstDomain{Float64}(0.63), θ = VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}}([TransformVariables.Infinity{false}(), TransformVariables.Infinity{false}(), TransformVariables.Infinity{false}(), TransformVariables.Infinity{false}()], [TransformVariables.Infinity{true}(), TransformVariables.Infinity{true}(), TransformVariables.Infinity{true}(), TransformVariables.Infinity{true}()], [0.0, 0.0, 0.0, 0.0]), Ω1 = PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}}(PDMats.PDiagMat{Float64,Array{Float64,1}}(1, [1.0], [1.0])), Ω2 = PSDDomain{Array{Float64,2}}([1.0 0.0; 0.0 1.0]), σ1 = RealDomain{Float64,Int64,Float64}(0.01, 1, 0.505), σ2 = RealDomain{Float64,Int64,Float64}(0.01, 1, 0.505))))
 Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}},Serialization.__deserialized_types__.var"#533#563",Serialization.__deserialized_types__.var"#534#564",Serialization.__deserialized_types__.var"#536#566",Pumas.LinearODE,Serialization.__deserialized_types__.var"#537#567",Serialization.__deserialized_types__.var"#550#580"},Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantInterpolationStructArray{Array{Float64,1},StructArrays.StructArray{NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Float64,Float64,Float64,Int64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64}},1,NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Int64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1}}},Int64},Symbol},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Nothing,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat},Float64,Array{Float64,1},Float64,Float64,Array{Optim.OptimizationState{Float64,Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},Pumas.var"#339#340"{NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.FOCEI,Array{Array{Float64,1},1},NamedTuple{(:optimize_fn, :constantcoef, :omegas, :ensemblealg),Tuple{Pumas.DefaultOptimizeFN{Nothing,NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases),Tuple{Bool,Bool,Bool,Float64,Bool}}},NamedTuple{(:tvbio,),Tuple{Float64}},Tuple{},EnsembleThreads}},ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},Pumas.ConstDomain{Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}}}(PumasModel
  Parameters: tvka, tvk23, tvvc, tvvp, tvQ, tvCL, tvbio, θ, Ω1, Ω2, σ1, σ2
  Random effects: η1, η2
  Covariates: WT, AGEYRS, PMA, ALBD, SCRD, FFM
  Dynamical variables: Depot1, Prodrug, Central, Peripheral
  Derived: DV
  Observed: DV
, Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantInterpolationStructArray{Array{Float64,1},StructArrays.StructArray{NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Float64,Float64,Float64,Int64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64,Float64}},1,NamedTuple{(:WT, :INITWT, :DFWT, :M, :SCR, :INITBMI, :DFBMI, :INITFFM, :DFFFM, :SCRD, :GA, :PMA, :AGEYRS, :BMI, :FFM, :ALBD),Tuple{Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Int64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1},Array{Float64,1}}},Int64},Symbol},Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Nothing,Int64},1},Array{Float64,1}}[Subject
  ID: 64
  Events: 24
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 10
  Events: 284
  Observables: DV: (n=26)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 154
  Events: 182
  Observables: DV: (n=6)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 42
  Events: 134
  Observables: DV: (n=11)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 147
  Events: 24
  Observables: DV: (n=5)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 50
  Events: 112
  Observables: DV: (n=3)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 46
  Events: 6
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 125
  Events: 32
  Observables: DV: (n=2)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 163
  Events: 34
  Observables: DV: (n=6)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 127
  Events: 14
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
  …  Subject
  ID: 12
  Events: 78
  Observables: DV: (n=4)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 42
  Events: 134
  Observables: DV: (n=11)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 30
  Events: 54
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 123
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 20
  Events: 116
  Observables: DV: (n=3)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 139
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 139
  Events: 2
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
, Subject
  ID: 79
  Events: 10
  Observables: DV: (n=1)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD

sai_matcha · January 24, 2021, 9:45am

It runs with out error. But full result was more than 90000 characters , so posted part of result here.

sai_matcha · January 24, 2021, 9:47am

last part of result is here

, Subject
  ID: 110
  Events: 18
  Observables: DV: (n=2)
  Covariates: WT, INITWT, DFWT, M, SCR, INITBMI, DFBMI, INITFFM, DFFFM, SCRD, GA, PMA, AGEYRS, BMI, FFM, ALBD
],  * Status: success

 * Candidate solution
    Final objective value:     3.564276e+02

 * Found with
    Algorithm:     BFGS

 * Convergence measures
    |x - x'|               = 1.45e+00 ≰ 0.0e+00
    |x - x'|/|x'|          = 4.98e-02 ≰ 0.0e+00
    |f(x) - f(x')|         = 2.22e-05 ≰ 0.0e+00
    |f(x) - f(x')|/|f(x')| = 6.24e-08 ≰ 0.0e+00
    |g(x)|                 = 9.89e-04 ≤ 1.0e-03

 * Work counters
    Seconds run:   135  (vs limit Inf)
    Iterations:    39
    f(x) calls:    45
    ∇f(x) calls:   40
, Pumas.FOCEI(), [[0.0, -1.3343096726444401, -0.3879317448798752], [0.0, -0.13856329522081576, -0.4501662966865892], [0.0, -0.0015110934736817946, 0.013616456183023363], [0.0, -0.31943276531723325, -0.004736795927651099], [0.0, 0.7084906805844491, 1.2834857695405288], [0.0, 0.014655756892671462, -0.431601843966014], [0.0, 0.770117805590958, 0.8021511071938359], [0.0, -0.42281084834308086, -0.6494343029389286], [0.0, -0.44344626690204453, -0.020515835850542147], [0.0, 0.022791106049660153, 0.019492673311563905]  …  [0.0, -0.381971125125141, -0.2289072632317178], [0.0, -0.14204753130313827, -0.16554534318431682], [0.0, -0.05943901675511194, -0.08282606269223304], [0.0, 0.11470140250140486, 0.11771524755109197], [0.0, -0.22692754086872974, -0.014537212946267664], [0.0, 0.5975236195853255, 1.1098743582213246], [0.0, -1.3343096726444401, -0.3879317448798752], [0.0, 0.18055845502900378, 0.054906084482718], [0.0, 
-0.6062520817017301, -0.019837116155745235], [0.0, 0.4283838744646378, 0.57655493597566]], (optimize_fn = Pumas.DefaultOptimizeFN{Nothing,NamedTuple{(:show_trace, :store_trace, :extended_trace, :g_tol, :allow_f_increases),Tuple{Bool,Bool,Bool,Float64,Bool}}}(nothing, (show_trace = false, store_trace = true, extended_trace = false, g_tol = 0.001, allow_f_increases = 
true)), constantcoef = (tvbio = 0.63,), omegas = , ensemblealg = EnsembleThreads()), ParamSet{NamedTuple{(:tvka, :tvk23, :tvvc, :tvvp, :tvQ, :tvCL, :tvbio, :θ, :Ω1, :Ω2, :σ1, :σ2),Tuple{RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64},Pumas.ConstDomain{Float64},VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},PSDDomain{Array{Float64,2}},RealDomain{Float64,Int64,Float64},RealDomain{Float64,Int64,Float64}}}}((tvka = RealDomain{Float64,Int64,Float64}(0.1, 
5, 2.55), tvk23 = RealDomain{Float64,Int64,Float64}(0.1, 20, 10.05), tvvc = RealDomain{Float64,Int64,Float64}(0.1, 250, 125.05), tvvp = RealDomain{Float64,Int64,Float64}(0.1, 250, 125.05), tvQ = RealDomain{Float64,Int64,Float64}(0.1, 250, 125.05), tvCL = RealDomain{Float64,Int64,Float64}(0.1, 15, 7.55), tvbio = Pumas.ConstDomain{Float64}(0.63), θ = VectorDomain{Array{TransformVariables.Infinity{false},1},Array{TransformVariables.Infinity{true},1},Array{Float64,1}}([TransformVariables.Infinity{false}(), TransformVariables.Infinity{false}(), TransformVariables.Infinity{false}(), TransformVariables.Infinity{false}()], [TransformVariables.Infinity{true}(), TransformVariables.Infinity{true}(), TransformVariables.Infinity{true}(), TransformVariables.Infinity{true}()], [0.0, 0.0, 0.0, 0.0]), Ω1 = PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}}(PDMats.PDiagMat{Float64,Array{Float64,1}}(1, [1.0], [1.0])), Ω2 
= PSDDomain{Array{Float64,2}}([1.0 0.0; 0.0 1.0]), σ1 = RealDomain{Float64,Int64,Float64}(0.01, 1, 0.505), σ2 = RealDomain{Float64,Int64,Float64}(0.01, 1, 0.505))))
 InexactError(:trunc, Int64, Inf)

andreasnoack · January 24, 2021, 2:17pm

I suppose one or some of the bootstrap fits fail for some reason but it’s difficult to know exactly what is going on. You could try to browse through the 1000 bootstrap fits and see if any of them look odd.

Topic		Replies	Views
Error in Displaying Bootstrap results bootstrap	0	369	March 5, 2021
Infer() error after fitting data Estimation	12	904	March 5, 2020
Error while using infer(res) Estimation	8	982	March 12, 2020
Error with Bootstrap bootstrap	3	191	November 17, 2022
Fields of SIR infer Pumas issues	1	242	November 8, 2021

ERROR: ArgumentError: quantiles are undefined in presence of NaNs

Related Topics