“…Two statically loaded independent I(1) common factors (F j , j = 1, 2), the first with a drift δ = 0.01, shape long-term growth over time (with T = 50, 100, 150) and pervasive dependence across N spatial units, hereafter regions (with N = 50, 100, 150). Short-term and local fluctuations of the error at time t for unit i, ε t,i , are governed by a STARMA(1,1,1) process (that is, an ARMA(1,1) augmented with a spatial lag term at time (t−1); see, e.g, Dai and Billard, 1998). Thus, the value of the variable Y at time t for region i is generated as follows:…”