Hi, I have an error in viewing the bootstrap results. I would like to know what the error indicates and how to rectify it. Below is the code and Error.
Code
BS = infer(Fit,Pumas.Bootstrap(samples=50))
Error
Bootstrap inference results
Successful minimization: true
Likelihood approximation: Pumas.FOCEI
Log-likelihood value: -1540.7111
Number of subjects: 121
Number of parameters: Fixed Optimized
1 7
Observation records: Active Missing
DV: 338 0
Total: 338 0
Error showing value of type Pumas.FittedPumasModelInference{Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:tvka, :tvcl, :tvvc, :tvvp, :tvq, :θ, :Ω, :σ_prop),Tuple{RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},VectorDomain{Array{Int64,1},Array{Int64,1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},RealDomain{Int64,Int64,Float64}}}},Serialization.__deserialized_types__.var"#9#29",Serialization.__deserialized_types__.var"#10#30",Serialization.__deserialized_types__.var"#12#32",Pumas.LinearODE,Serialization.__deserialized_types__.var"#13#33",Serialization.__deserialized_types__.var"#21#41"},Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantCovar{NamedTuple{(:WT, :GROUP, :GRPBMI, :SEX, :BMI, :WTGRP),Tuple{Float64,Int64,Int64,Int64,Float64,Int64}}},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{(:tvka,),Tuple{Float64}},Tuple{},EnsembleThreads}},ParamSet{NamedTuple{(:tvka, :tvcl, :tvvc, :tvvp, :tvq, :θ, :Ω, :σ_prop),Tuple{Pumas.ConstDomain{Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},VectorDomain{Array{Int64,1},Array{Int64,1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},RealDomain{Int64,Int64,Float64}}}}},Pumas.Bootstraps,Float64}:
ERROR: MethodError: no method matching isnan(::Array{Float64,1})
Closest candidates are:
isnan(::Float16) at float.jl:537
isnan(::Missing) at missing.jl:100
isnan(::BigFloat) at mpfr.jl:881
...
Stacktrace:
[1] _quantilesort!(::Array{Array{Float64,1},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!(::Array{Array{Float64,1},1}, ::Float64; sorted::Bool, alpha::Float64, beta::Float64) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\Statistics\src\Statistics.jl:943
[3] quantile(::Array{Array{Float64,1},1}, ::Float64; sorted::Bool, alpha::Float64, beta::Float64) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\Statistics\src\Statistics.jl:1052
[4] quantile(::Array{Array{Float64,1},1}, ::Float64) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\Statistics\src\Statistics.jl:1052
[5] _push_varinfo!(::Array{String,1}, ::Array{Float64,1}, ::Array{Float64,1}, ::Array{Tuple{Float64,Float64},1}, ::String, ::Float64, ::Float64, ::Float64, ::Array{Array{Float64,1},1}) at C:\Users\julia\.julia\packages\Pumas\iIBBr\src\estimation\show.jl:203
[6] _push_varinfo!(::Array{String,1}, ::Array{Float64,1}, ::Array{Float64,1}, ::Array{Tuple{Float64,Float64},1}, ::Symbol, ::Array{Float64,1}, ::Array{Float64,1}, ::Float64, ::Array{Array{Float64,1},1}) at C:\Users\julia\.julia\packages\Pumas\iIBBr\src\estimation\show.jl:183
[7] coeftable(::Pumas.FittedPumasModelInference{Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:tvka, :tvcl, :tvvc, :tvvp, :tvq, :θ, :Ω, :σ_prop),Tuple{RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},VectorDomain{Array{Int64,1},Array{Int64,1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},RealDomain{Int64,Int64,Float64}}}},Serialization.__deserialized_types__.var"#9#29",Serialization.__deserialized_types__.var"#10#30",Serialization.__deserialized_types__.var"#12#32",Pumas.LinearODE,Serialization.__deserialized_types__.var"#13#33",Serialization.__deserialized_types__.var"#21#41"},Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantCovar{NamedTuple{(:WT, :GROUP, :GRPBMI, :SEX, :BMI,
:WTGRP),Tuple{Float64,Int64,Int64,Int64,Float64,Int64}}},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{(:tvka,),Tuple{Float64}},Tuple{},EnsembleThreads}},ParamSet{NamedTuple{(:tvka, :tvcl, :tvvc, :tvvp, :tvq, :θ, :Ω, :σ_prop),Tuple{Pumas.ConstDomain{Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},VectorDomain{Array{Int64,1},Array{Int64,1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},RealDomain{Int64,Int64,Float64}}}}},Pumas.Bootstraps,Float64}) at C:\Users\julia\.julia\packages\Pumas\iIBBr\src\estimation\show.jl:235
[8] show(::IOContext{REPL.Terminals.TTYTerminal}, ::MIME{Symbol("text/plain")}, ::Pumas.FittedPumasModelInference{Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:tvka, :tvcl, :tvvc, :tvvp, :tvq, :θ, :Ω, :σ_prop),Tuple{RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},VectorDomain{Array{Int64,1},Array{Int64,1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},RealDomain{Int64,Int64,Float64}}}},Serialization.__deserialized_types__.var"#9#29",Serialization.__deserialized_types__.var"#10#30",Serialization.__deserialized_types__.var"#12#32",Pumas.LinearODE,Serialization.__deserialized_types__.var"#13#33",Serialization.__deserialized_types__.var"#21#41"},Array{Subject{NamedTuple{(:DV,),Tuple{Array{Union{Missing, Float64},1}}},Pumas.ConstantCovar{NamedTuple{(:WT, :GROUP, :GRPBMI, :SEX, :BMI, :WTGRP),Tuple{Float64,Int64,Int64,Int64,Float64,Int64}}},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{(:tvka,),Tuple{Float64}},Tuple{},EnsembleThreads}},ParamSet{NamedTuple{(:tvka, :tvcl, :tvvc, :tvvp, :tvq, :θ, :Ω, :σ_prop),Tuple{Pumas.ConstDomain{Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},RealDomain{Int64,Int64,Float64},VectorDomain{Array{Int64,1},Array{Int64,1},Array{Float64,1}},PDiagDomain{PDMats.PDiagMat{Float64,Array{Float64,1}}},RealDomain{Int64,Int64,Float64}}}}},Pumas.Bootstraps,Float64}) at C:\Users\julia\.julia\packages\Pumas\iIBBr\src\estimation\show.jl:273
[9] 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
[10] display(::REPL.REPLDisplay, ::Any) at C:\Users\julia\AppData\Local\Julia-1.5.2\share\julia\stdlib\v1.5\REPL\src\REPL.jl:218
[11] display(::Any) at .\multimedia.jl:328
[12] #invokelatest#1 at .\essentials.jl:710 [inlined]
[13] invokelatest at .\essentials.jl:709 [inlined]
[14] 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
[15] 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
[16] (::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
[17] #invokelatest#1 at .\essentials.jl:710 [inlined]
[18] invokelatest at .\essentials.jl:709 [inlined]
[19] 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
[20] 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
[21] (::REPL.var"#38#42"{REPL.LineEditREPL,REPL.REPLBackendRef})() at .\task.jl:356
The number of successful runs are below
Bootstrap inference results
Successful fits: 44 out of 50
No stratification.