Some useful tools for the cross-sectional forecast reconciliation of a linearly constrained (e.g., hierarchical/grouped) multiple time series.
hts_tools(C, h = 1, Ut, nb, sparse = TRUE)
A list of five elements:
C
(n n_b) cross-sectional (contemporaneous) aggregation matrix.
S
(n n_b) cross-sectional (contemporaneous) summing matrix,
S = [arrayc C
I_n_barray].
Ut
(n_a n) zero constraints cross-sectional (contemporaneous) kernel matrix. If the hierarchy admits a structural representation U' = [I \ -C]
n
Number of variables n_a + n_b.
na
Number of upper level variables.
nb
Number of bottom level variables.
(n_a n_b) cross-sectional (contemporaneous) matrix mapping the bottom level series into the higher level ones.
Forecast horizon (default is 1
).
Zero constraints cross-sectional (contemporaneous) kernel matrix
(U'y = 0) spanning the null space valid
for the reconciled forecasts. It can be used instead of parameter
C
, but nb
is needed if
U' [I \ -C]. If the hierarchy
admits a structural representation, U' has dimension
(n_a n).
Number of bottom time series; if C
is present, nb
and Ut
are not used.
Option to return sparse matrices (default is TRUE
).
Other utilities:
Cmatrix()
,
FoReco2ts()
,
agg_ts()
,
arrange_hres()
,
commat()
,
ctf_tools()
,
lcmat()
,
oct_bounds()
,
residuals_matrix()
,
score_index()
,
shrink_estim()
,
thf_tools()
# One level hierarchy (na = 1, nb = 2)
obj <- hts_tools(C = matrix(c(1, 1), 1), sparse = FALSE)
Run the code above in your browser using DataLab