How to estimate uncertainty in a LinearMixedModel?

Specific Domains Statistics
4

That produces an interval for a single Cut & Clarity parameter, but I was looking for an interval for every value of Cut & Clarity. I want to make a coefficient plot, to show the credible ranges for each value of Cut and Clarity .

pbdf = DataFrame(pb.allpars)
combine(groupby(pbdf, [:type, :group, :names]), :value => shortestcovint => :interval)

#= 3×4 DataFrame
 Row │ type    group          names        interval           
     │ String  String?        String?      Tuple…             
─────┼────────────────────────────────────────────────────────
   1 │ β       missing        (Intercept)  (3600.75, 4104.61)
   2 │ σ       Cut & Clarity  (Intercept)  (552.975, 906.98)
   3 │ σ       residual       missing      (3892.35, 3939.2) =#