special_xreg function

Helper special for producing a model matrix of exogenous regressors