TSLM function

Fit a linear model with time series components