Skip to main content

Simple Structure of Sectoral Sales

Simple Structure of Sectoral Sales

Assuming a structure for the time series of sectoral sales involves defining it as a sum of time series that partially compose it. For instance, fitting sectoral sales with a generic spts_{pt} results in:

Sp=β^pst+rptS_p = \hat \beta_p s_t + r_{pt}

This leads to an expression for aggregate variance, which is the sum of elements of the cross covariance matrix:

var(X)=ijcov(Si,Sj)=ijcov(β^ist+rit,β^jst+rjt)var(X) = \sum_{ij} cov(S_i, S_j) = \sum_{ij} cov(\hat \beta_i s_t + r_{it}, \hat \beta_j s_t + r_{jt})

Here, β^p\hat \beta_p are elements of β=(stTst)1stTY\bm \beta = (s_t^T s_t)^{-1} s_t^T Y, with YRT×PY \in \mathbb{R}^{T \times P} containing the observed parts' time series and sts_t as a column vector used for fitting. It can be z-standardized so that cov(st,st)=1cov(s_t, s_t) = 1 (see Figure 4).

Note: A specific example is fitting (OLS) using a factor of linear time evolution st=ts_t = t, then: Spt=β^pt+rptS_{pt} = \hat \beta_{p} t + r_{pt}, where βp\beta_{p} is as described and rptr_{pt} are the residuals of the fit.

The elements of the cross covariance matrix are expressed as:

cov(β^ist+rit,β^jst+rjt)=cov(β^ist,β^jst)+cov(β^ist,rjt)+cov(rit,β^jst)+cov(rit,rjt)\begin{split} cov(\hat \beta_i s_t + r_{it}, \hat \beta_j s_t + r_{jt}) = & cov(\hat \beta_i s_t, \hat \beta_j s_t) + cov(\hat \beta_i s_t, r_{jt}) \\ & + cov(r_{it}, \hat \beta_j s_t) + cov(r_{it}, r_{jt}) \end{split}

Each of the P×PP \times P elements is composed of four components. Generally, a structure of sectoral sales made from NN terms results in N×NN \times N terms forming each element of the cross covariance matrix.

The β^\hat \beta coefficients are arranged into a column vector β^\bm{\hat \beta} of length PP, and the standard deviation of residuals σ\sigma in the column vector σ^\bm{\hat \sigma} of length PP. The matrix product β^β^T\bm{\hat \beta} \bm{\hat \beta}^T creates a P×PP \times P matrix (an outer product), and the cross covariance matrix among sectors is expressed as:

C=β^β^TC(st,st)+β^σTC(st,ϵt)+σβ^TC(ϵt,st)+σσTC(ϵt,ϵt)\begin{split} C = & \bm{\hat \beta} \bm{\hat \beta}^T \cdot C(s_t, s_t) + \bm{\hat \beta} \bm{\sigma}^T \cdot C(s_t, \epsilon_{t}) \\ & + \bm{\sigma} \bm{\hat \beta}^T \cdot C(\epsilon_{t}, s_t) + \bm{\sigma} \bm{\sigma}^T \cdot C(\epsilon_{t}, \epsilon_{t}) \end{split}

where C(a,b)C({a}, {b}) denotes the P×PP \times P cross covariance matrix computed from vectors aa and bb. They are multiplied with the outer product matrices element-wise.

Structure of Cross Covariance Matrix = Outer Product × Net Covariance

Figure 4: Structure of cross covariance matrix (N×NN \times N blocks with elements made of outer product times net covariance). This illustration is from a simplified toy example where sectoral time series are decomposed into dependence as mean time series (fitted by OLS) and residuals.

This structure of the cross covariance matrix as an element-wise product of an outer product matrix and a net cross covariance is illustrated in Figure 4.

In this context, aggregate variance comprises a comovement contribution and an idiosyncrasies contribution, apart from cross terms that often do not significantly contribute.

Note: If the time series of the parts are independent, the expectation is E[var[X]]=kσk2E[var[X]] = \sum_k \sigma^2_k. The cross covariance terms are, on average, zero, although it is likely that none of them is truly null. They will add up to the uncertainty with which we can estimate σ2\sigma^2 as kσk2\sum_k \sigma^2_k. We maintain cross covariance terms even when their total expectation is zero.