For the withing-customer errors (εijk)k, we explored different covariance structures: compound symmetry, first-order autoregressive, first-order autoregressive moving average, Toeplitz, as well as heterogeneous autoregressive and heterogeneous Toeplitz, of which the latter two allowed for unequal variance (View the MathML source) across weeks.