expand_TS function

Expand the TS models across the factorial combination of LDA models, formulas, and number of change points