lagfit function

Fits a piecewise glm model with lags