Estimation of parameters in 3-segment (i.e. 2 change-point) regression models with heteroscedastic variances is provided based on both likelihood and hybrid Bayesian approaches, with and without continuity constraints at the change points.
