powered by
VAR regressor matrix, see Section 2.2.3 in Koop and Korobilis (2010)
VAR_regressor_matrix(data, var.order)