== cellmeans: divergences = 0; max R-hat = 1.000; min bulk ESS = 2543
                mean      sd  hdi_2.5%  hdi_97.5%  mcse_mean  mcse_sd  ess_bulk  ess_tail  r_hat
mu_uggFW[0]  125.385  24.958    76.241    175.926      0.320    0.289    6100.0    6907.0    1.0
mu_uggFW[1]  140.259  27.359    84.940    192.197      0.344    0.285    6379.0    7235.0    1.0
mu_uggFW[2]  139.981  25.802    89.201    191.719      0.330    0.302    6167.0    7118.0    1.0
mu_uggFW[3]  122.072  25.349    70.223    170.288      0.331    0.326    5799.0    7095.0    1.0
mu_uggFW[4]   50.862  26.684    -0.492    104.858      0.359    0.293    5540.0    6426.0    1.0
s_ext          0.971   0.695     0.000      2.271      0.013    0.010    2543.0    3317.0    1.0
s_tech         0.008   0.002     0.005      0.012      0.000    0.000    9078.0    7820.0    1.0
b0             0.089   0.005     0.079      0.100      0.000    0.000    4235.0    5506.0    1.0
b1             0.006   0.000     0.005      0.006      0.000    0.000    5752.0    6820.0    1.0

== linear: divergences = 0; max R-hat = 1.000; min bulk ESS = 2703
                mean      sd  hdi_2.5%  hdi_97.5%  mcse_mean  mcse_sd  ess_bulk  ess_tail  r_hat
mu_uggFW[0]  146.772  25.540    94.581    195.584      0.394    0.291    4224.0    6005.0    1.0
mu_uggFW[1]  130.486  21.412    89.727    172.858      0.316    0.236    4604.0    5925.0    1.0
mu_uggFW[2]  114.201  20.337    75.778    155.807      0.283    0.234    5206.0    6585.0    1.0
mu_uggFW[3]   97.915  22.754    52.557    143.049      0.311    0.258    5381.0    6887.0    1.0
mu_uggFW[4]   81.630  27.765    26.945    136.440      0.385    0.315    5193.0    6120.0    1.0
a              7.339   1.277     4.729      9.779      0.020    0.015    4224.0    6005.0    1.0
b             -0.814   0.432    -1.658      0.084      0.007    0.006    4255.0    5124.0    1.0
s_ext          1.542   0.666     0.140      2.805      0.012    0.009    2703.0    2682.0    1.0
s_tech         0.008   0.002     0.005      0.012      0.000    0.000    4765.0    6406.0    1.0
b0             0.089   0.006     0.077      0.100      0.000    0.000    5028.0    5839.0    1.0
b1             0.006   0.000     0.005      0.006      0.000    0.000    5223.0    5595.0    1.0

== hinge: divergences = 0; max R-hat = 1.000; min bulk ESS = 3660
                mean      sd  hdi_2.5%  hdi_97.5%  mcse_mean  mcse_sd  ess_bulk  ess_tail  r_hat
mu_uggFW[0]  134.040  19.210    94.984    170.532      0.257    0.237    5756.0    5825.0    1.0
mu_uggFW[1]  133.959  19.175    94.997    170.532      0.256    0.237    5785.0    5815.0    1.0
mu_uggFW[2]  133.049  19.100    94.984    170.070      0.253    0.234    5914.0    5883.0    1.0
mu_uggFW[3]  120.678  21.655    77.078    161.526      0.272    0.222    6425.0    6995.0    1.0
mu_uggFW[4]   55.595  25.972     5.768    107.946      0.311    0.319    7193.0    7159.0    1.0
a              6.702   0.961     4.749      8.527      0.013    0.012    5756.0    5825.0    1.0
b             -4.201   2.257    -8.986     -0.602      0.031    0.032    5593.0    5449.0    1.0
theta          3.845   0.570     2.706      4.762      0.008    0.010    5115.0    5086.0    1.0
s_ext          0.761   0.516     0.001      1.728      0.008    0.006    3660.0    4797.0    1.0
s_tech         0.008   0.002     0.005      0.011      0.000    0.000    7951.0    7925.0    1.0
b0             0.089   0.006     0.078      0.100      0.000    0.000    5646.0    5959.0    1.0
b1             0.006   0.000     0.005      0.006      0.000    0.000    6394.0    6599.0    1.0

== Contrasts vs T1 (ug/g FW), cell-means model
       median     lo95    hi95  P(<0)  P(< -50)
T2-T1  15.505  -48.429  75.806  0.294     0.024
T3-T1  14.504  -44.814  73.740  0.294     0.018
T4-T1  -3.360  -61.883  54.766  0.552     0.050
T5-T1 -75.299 -131.945 -11.735  0.985     0.824

P(T5 is the lowest treatment mean) = 0.966
P(T5 < every other treatment) = 0.966

Linear slope (ug/g FW per severity step): median -16.3, 95% CrI [-33.7, 1.3], P(slope<0) = 0.968
Hinge post-threshold slope: median -74.5, 95% CrI [-199.1, -21.4], P(<0) = 0.997
Hinge threshold theta: median 3.91, 95% CrI [2.43, 4.66]; P(theta >= 3) = 0.923

== PSIS-LOO comparison (well-level likelihood)
           rank  elpd_loo  p_loo  elpd_diff  weight     se    dse  warning scale
hinge         0    90.385  6.454      0.000     1.0  6.367  0.000    False   log
cellmeans     1    88.587  8.567      1.799     0.0  6.192  0.750    False   log
linear        2    88.142  8.673      2.243     0.0  5.931  1.239    False   log

Between-extract SD: median 17.0 ug/g FW [0.8, 53.8]
Technical (well) SD in concentration units: median 28.3 ug/g FW [18.9, 42.9]

== Assurance (Bayesian expected power) for one-sided T5 < T1 at alpha 0.05
 n_per_treatment  assurance
               2     0.5625
               3     0.7825
               4     0.8500
               5     0.8875
               6     0.8750
               8     0.9125
              10     0.9350