add comment

JuliaStats · May 20, 2022 · 70d2537 · 70d2537
1 parent 2a1ccc4
commit 70d2537
Showing 1 changed file with 5 additions and 0 deletions.
diff --git a/src/linearmixedmodel.jl b/src/linearmixedmodel.jl
@@ -444,6 +444,11 @@ function fit!(
         val = try
             objective(updateL!(setθ!(m, x)))
         catch ex
+            # This can happen when the optimizer drifts into an area where
+            # there isn't enough shrinkage. Why finitial? Generally, it will
+            # be the (near) worst case scenario value, so the optimizer won't
+            # view it as an optimum. Using Inf messes up the quadratic
+            # approximation in BOBYQA.
             ex isa PosDefException || rethrow()
             iter == 0 && rethrow()
             m.optsum.finitial