calcregions function

Fit segmented regression models for all combinations of breakpoints