Issues with NaivePooled() Fittting

Estimation
1

Hi,
I am having some issues for the first time using NaivePooled for fitting data-
I have a model with and without η. When I used FOCEI to fit the model with random effects it all goes well. When I remove random effects to do a NaivePooled fit, I quickly get an error. I am not sure if there a reason this will be the case (devaince did not work either). My fit is attempting to fit two parameters such as my dv value (which is the ratio of two other variables is 1.0 at all times)

thanks,

Ale

0.592398 seconds (914.28 k allocations: 810.475 MiB, 34.59% gc time)
MethodError: no method matching _orth_empirical_bayes!(::Array{Float64,1}, ::PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))}, ::Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}}, ::NamedTuple{(:θ, :σ_prop),Tuple{Array{Float64,1},Float64}}, ::Pumas.NaivePooled)
Closest candidates are:
  _orth_empirical_bayes!(::AbstractArray{T,1} where T, ::PumasModel, ::Subject, ::NamedTuple, !Matched::Union{Pumas.FO, Pumas.FOI, Pumas.LLQuad}, !Matched::Any...; kwargs...) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:208
  _orth_empirical_bayes!(::AbstractArray{T,1} where T, ::PumasModel, ::Subject, ::NamedTuple, !Matched::Union{Pumas.FOCE, Pumas.FOCEI, Pumas.LaplaceI}, !Matched::Any...; kwargs...) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:221

Stacktrace:
 [1] #_orth_empirical_bayes#143(::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}}, ::typeof(Pumas._orth_empirical_bayes), ::PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))}, ::Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}}, ::NamedTuple{(:θ, :σ_prop),Tuple{Array{Float64,1},Float64}}, ::Pumas.NaivePooled) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:197
 [2] _orth_empirical_bayes(::PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))}, ::Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}}, ::NamedTuple{(:θ, :σ_prop),Tuple{Array{Float64,1},Float64}}, ::Pumas.NaivePooled) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:195
 [3] #marginal_nll#151(::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}}, ::typeof(Pumas.marginal_nll), ::PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))}, ::Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}}, ::NamedTuple{(:θ, :σ_prop),Tuple{Array{Float64,1},Float64}}, ::Pumas.NaivePooled) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:339
 [4] marginal_nll(::PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))}, ::Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}}, ::NamedTuple{(:θ, :σ_prop),Tuple{Array{Float64,1},Float64}}, ::Pumas.NaivePooled) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:339
 [5] #marginal_nll#152(::Pumas.ParallelType, ::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}}, ::typeof(Pumas.marginal_nll), ::PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))}, ::Array{Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1}, ::NamedTuple{(:θ, :σ_prop),Tuple{Array{Float64,1},Float64}}, ::Vararg{Any,N} where N) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:351
 [6] marginal_nll(::PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))}, ::Array{Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1}, ::NamedTuple{(:θ, :σ_prop),Tuple{Array{Float64,1},Float64}}, ::Pumas.NaivePooled) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:351
 [7] #deviance#169(::Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}}, ::typeof(deviance), ::PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))}, ::Array{Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1}, ::NamedTuple{(:θ, :σ_prop),Tuple{Array{Float64,1},Float64}}, ::Vararg{Any,N} where N) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:510
 [8] deviance(::PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))}, ::Array{Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1}, ::NamedTuple{(:θ, :σ_prop),Tuple{Array{Float64,1},Float64}}, ::Vararg{Any,N} where N) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:510
 [9] deviance(::Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))},Array{Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},getfield(Pumas, Symbol("##211#212")){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},getfield(Pumas, Symbol("##211#212")){NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.NaivePooled,Array{Array{Any,1},1},Tuple{},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},TransformVariables.TransformTuple{NamedTuple{(:θ, :σ_prop),Tuple{Pumas.ElementArrayTransform{TransformVariables.ScaledShiftedLogistic{Float64},1},TransformVariables.ShiftedExp{true,Float64}}}}}) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\likelihoods.jl:1150
 [10] _print_fit_header(::IOContext{Base.GenericIOBuffer{Array{UInt8,1}}}, ::Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))},Array{Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},getfield(Pumas, Symbol("##211#212")){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},getfield(Pumas, Symbol("##211#212")){NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.NaivePooled,Array{Array{Any,1},1},Tuple{},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},TransformVariables.TransformTuple{NamedTuple{(:θ, :σ_prop),Tuple{Pumas.ElementArrayTransform{TransformVariables.ScaledShiftedLogistic{Float64},1},TransformVariables.ShiftedExp{true,Float64}}}}}) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\show.jl:10
 [11] show(::IOContext{Base.GenericIOBuffer{Array{UInt8,1}}}, ::MIME{Symbol("text/plain")}, ::Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))},Array{Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},getfield(Pumas, Symbol("##211#212")){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},getfield(Pumas, Symbol("##211#212")){NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.NaivePooled,Array{Array{Any,1},1},Tuple{},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},TransformVariables.TransformTuple{NamedTuple{(:θ, :σ_prop),Tuple{Pumas.ElementArrayTransform{TransformVariables.ScaledShiftedLogistic{Float64},1},TransformVariables.ShiftedExp{true,Float64}}}}}) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\Pumas\Twdge\src\estimation\show.jl:23
 [12] limitstringmime(::MIME{Symbol("text/plain")}, ::Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))},Array{Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},getfield(Pumas, Symbol("##211#212")){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},getfield(Pumas, Symbol("##211#212")){NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.NaivePooled,Array{Array{Any,1},1},Tuple{},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},TransformVariables.TransformTuple{NamedTuple{(:θ, :σ_prop),Tuple{Pumas.ElementArrayTransform{TransformVariables.ScaledShiftedLogistic{Float64},1},TransformVariables.ShiftedExp{true,Float64}}}}}) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\IJulia\F1GUo\src\inline.jl:43
 [13] display_mimestring(::MIME{Symbol("text/plain")}, ::Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))},Array{Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},getfield(Pumas, Symbol("##211#212")){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},getfield(Pumas, Symbol("##211#212")){NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.NaivePooled,Array{Array{Any,1},1},Tuple{},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},TransformVariables.TransformTuple{NamedTuple{(:θ, :σ_prop),Tuple{Pumas.ElementArrayTransform{TransformVariables.ScaledShiftedLogistic{Float64},1},TransformVariables.ShiftedExp{true,Float64}}}}}) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\IJulia\F1GUo\src\display.jl:67
 [14] display_dict(::Pumas.FittedPumasModel{PumasModel{ParamSet{NamedTuple{(:θ, :σ_prop),Tuple{VectorDomain{Array{Float64,1},Array{Float64,1},Array{Float64,1}},RealDomain{Float64,TransformVariables.Infinity{true},Float64}}}},getfield(Main, Symbol("##21#28")),getfield(Main, Symbol("##22#29")),getfield(Main, Symbol("##23#30")),ODEProblem{Nothing,Tuple{Nothing,Nothing},false,Nothing,ODEFunction{false,getfield(Main, Symbol("##24#31")),LinearAlgebra.UniformScaling{Bool},Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing,Nothing},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},DiffEqBase.StandardODEProblem},getfield(Main, Symbol("##26#33")),getfield(Main, Symbol("##27#34"))},Array{Subject{NamedTuple{(:Ratio,),Tuple{Array{Union{Missing, Float64},1}}},Nothing,Array{Pumas.Event{Float64,Float64,Float64,Float64,Float64,Float64,Int64},1},Array{Float64,1}},1},Optim.MultivariateOptimizationResults{Optim.BFGS{LineSearches.InitialStatic{Float64},LineSearches.BackTracking{Float64,Int64},getfield(Pumas, Symbol("##211#212")){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},getfield(Pumas, Symbol("##211#212")){NLSolversBase.OnceDifferentiable{Float64,Array{Float64,1},Array{Float64,1}},Array{Float64,1}},Nothing,Optim.Flat}},1},Bool},Pumas.NaivePooled,Array{Array{Any,1},1},Tuple{},Base.Iterators.Pairs{Union{},Union{},Tuple{},NamedTuple{,Tuple{}}},TransformVariables.TransformTuple{NamedTuple{(:θ, :σ_prop),Tuple{Pumas.ElementArrayTransform{TransformVariables.ScaledShiftedLogistic{Float64},1},TransformVariables.ShiftedExp{true,Float64}}}}}) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\IJulia\F1GUo\src\display.jl:96
 [15] #invokelatest#1 at .\essentials.jl:790 [inlined]
 [16] invokelatest at .\essentials.jl:789 [inlined]
 [17] execute_request(::ZMQ.Socket, ::IJulia.Msg) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\IJulia\F1GUo\src\execute_request.jl:112
 [18] #invokelatest#1 at .\essentials.jl:790 [inlined]
 [19] invokelatest at .\essentials.jl:789 [inlined]
 [20] eventloop(::ZMQ.Socket) at C:\Users\awolf-yadlin\.juliapro\JuliaPro_v1.2.0-2\packages\IJulia\F1GUo\src\eventloop.jl:8
 [21] (::getfield(IJulia, Symbol("##15#18")))() at .\task.jl:268
2

It’s a bug in the computation of the marginal likelihood function. It tries to compute the empirical Bayes estimates even though the estimation method is NaivePooled. However, the issue has already been fixed so it will work in the next release of Pumas which we hope will be available pretty soon. I noticed that there is a similar problem with the computation of the standard errors. We’ll try to have that fixed before the release.

3

Thanks for the quick reply!

It’s good to know I didn’t mess up this time :slight_smile:

A