Two LH stages, no environment, non-linear Bayesian fit, heterogeneous residual variance. R code with knitr.