variance_lm function

Variance estimation for linear models with change points