doctest failure in src/examples.md:513-541 ```jldoctest julia> using GLM, RDatasets, StatsBase, DataFrames, Optim julia> trees = DataFrame(dataset("datasets", "trees")); julia> bic_glm(λ) = bic(glm(@formula(Volume ~ Height + Girth), trees, Normal(), PowerLink(λ))); julia> optimal_bic = optimize(bic_glm, -1.0, 1.0); julia> round(optimal_bic.minimizer, digits = 5) # Optimal λ 0.40935 julia> glm(@formula(Volume ~ Height + Girth), trees, Normal(), PowerLink(optimal_bic.minimizer)) # Best model GeneralizedLinearModel Volume ~ 1 + Height + Girth Coefficients: ──────────────────────────────────────────────────────────────────────────── Coef. Std. Error z Pr(>|z|) Lower 95% Upper 95% ──────────────────────────────────────────────────────────────────────────── (Intercept) -1.07586 0.352543 -3.05 0.0023 -1.76684 -0.384892 Height 0.0232172 0.00523331 4.44 <1e-05 0.0129601 0.0334743 Girth 0.242837 0.00922556 26.32 <1e-99 0.224756 0.260919 ──────────────────────────────────────────────────────────────────────────── julia> round(optimal_bic.minimum, digits=5) 156.37638 ``` Subexpression: glm(@formula(Volume ~ Height + Girth), trees, Normal(), PowerLink(optimal_bic.minimizer)) # Best model Evaluated output: GeneralizedLinearModel Volume ~ 1 + Height + Girth Coefficients: ──────────────────────────────────────────────────────────────────────────── Coef. Std. Error z Pr(>|z|) Lower 95% Upper 95% ──────────────────────────────────────────────────────────────────────────── (Intercept) -1.07586 0.352543 -3.05 0.0023 -1.76684 -0.384892 Height 0.0232172 0.00523331 4.44 <1e-05 0.0129601 0.0334743 Girth 0.242837 0.00922555 26.32 <1e-99 0.224756 0.260919 ──────────────────────────────────────────────────────────────────────────── Expected output: GeneralizedLinearModel Volume ~ 1 + Height + Girth Coefficients: ──────────────────────────────────────────────────────────────────────────── Coef. Std. Error z Pr(>|z|) Lower 95% Upper 95% ──────────────────────────────────────────────────────────────────────────── (Intercept) -1.07586 0.352543 -3.05 0.0023 -1.76684 -0.384892 Height 0.0232172 0.00523331 4.44 <1e-05 0.0129601 0.0334743 Girth 0.242837 0.00922556 26.32 <1e-99 0.224756 0.260919 ──────────────────────────────────────────────────────────────────────────── diff = Warning: Diff output requires color. GeneralizedLinearModel Volume ~ 1 + Height + Girth Coefficients: ──────────────────────────────────────────────────────────────────────────── Coef. Std. Error z Pr(>|z|) Lower 95% Upper 95% ──────────────────────────────────────────────────────────────────────────── (Intercept) -1.07586 0.352543 -3.05 0.0023 -1.76684 -0.384892 Height 0.0232172 0.00523331 4.44 <1e-05 0.0129601 0.0334743 Girth 0.242837 0.00922556 0.00922555 26.32 <1e-99 0.224756 0.260919 ────────────────────────────────────────────────────────────────────────────