REVIEW 5 major objections 6 minor 16 references
Horizon Activation Mapping for Neural Networks in Time Series Forecasting
T0 review · 5 major / 6 minor · reviewed 2026-08-03 · deepseek-v4-flash
Pith's one-line read A grad-CAM-style tool maps which forecast timesteps a network learns to emphasize, across model families.
desk verdict A novel masking-based visualization for time-series forecasting, but the normalization and baseline assumptions undermine the central claims. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The core object is the pair of HAM curves: for a horizon of length H, causal subseries [0→h] and anti-causal subseries [h→H] are masked in the loss, and the average gradient norm over the whole training set is computed for every h. A line of proportionality connecting (0,0) to (H, G) serves as a baseline for uniform activity, where G is the maximum norm average across subseries. Derived quantities include the gradient equivariant point (where causal and anti-causal norms are equal), signed area curves measuring deviation from proportionality, interpolated area plots for comparing different horizon sizes, and difference plots d(t) = [g(0→t) − g(t→H)]/max|·|.
What would settle it
Train a model with known per-timestep importance, such as a linear model with preset weights or a transformer with determinable attention. Compute HAM curves for that model and check whether the gradient equivariant point and the deviation from the proportionality line correlate with the known importance profile. If the curves are identical across random and trained models, or ignore the known importance structure, the interpretive claims of HAM would be falsified.
Extended reading notes
Core claim
HAM computes, for each subseries of the forecast horizon, the average gradient norm of the loss when that subseries is masked before backpropagation. Causal mode masks timesteps after a cut point; anti-causal mode masks timesteps before it. Plotting these norm averages against the cut point yields two curves whose shape, intersection (the 'gradient equivariant point'), and deviation from a line of proportionality indicate the model's relative emphasis on shorter versus longer horizon subseries. The paper's central claim is that these curves are not just mathematical artifacts but reflect meaningful learning behavior: they change reproducibly with batch size, dropout, early stopping, data spl
Load-bearing premise
The gradient norm of the masked loss, computed by scaling per-timestep contributions by the full horizon length, faithfully measures how much the model 'activates' for each subseries—rather than merely reflecting the increasing size of the cumulative sum.
Editorial extensions
If this is right
- HAM provides a single visualization that works across fundamentally different model families (MLP, attention, state-space, diffusion), making cross-architecture comparison of gradient behavior possible without modifying the models.
- Batch-size comparisons in HAM suggest that converged models exhibit a polynomial saturation in gradient norm averages as batch size increases, which could be used to estimate the effect of batch size from a few runs.
- Dropout is shown to amplify gradient magnitudes, especially for longer horizon subseries, aligning with improved multivariate performance—offering a diagnostic for regularizer effects.
- HAM can inform validation set selection by comparing relative variability across timesteps between candidate validation splits.
- Horizon-size HAM patterns reveal family-specific inductive biases, such as NHITS's linear-combination assumption and SpaceTime's exponential weighting of longer horizons, which could guide architecture selection.
- The study of early stopping shows that post-stopping epochs reduce gradient magnitudes monotonically, suggesting HAM could be used to detect when a model has left a local minimum.
Reading between the lines
- Because the causal and anti-causal curves are cumulative sums over subseries, they are mathematically forced to be monotonically increasing and decreasing, respectively; the 'line of proportionality' comparison therefore needs a null model (e.g., random weights or a purely linear model) to calibrate whether observed deviations exceed what any model would produce.
- A natural validation experiment is to compare HAM's gradient-norm curves against known per-timestep importance signals—for instance, attention weights in a transformer or learned weights in a linear model—to test whether HAM tracks ground-truth emphasis or merely reflects gradient scaling.
- The paper's suggestion of layer-wise HAM is underdeveloped but could be extended to study learning dynamics within a single network, potentially revealing which layers specialize on short versus long horizon subseries.
- For diffusion models, the consistently larger-than-proportional norms for small subseries hint that noise-injected training creates a characteristic gradient signature; testing whether this signature predicts denoising performance could make HAM a useful probe for generative forecasting models.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper introduces Horizon Activation Mapping (HAM), a grad-CAM-inspired visualization technique for time series forecasting models. HAM computes causal and anti-causal masked-loss gradient norms for horizon subseries and plots them with a 'line of proportionality' and derived area/difference plots. The authors apply HAM to several model families (NHITS, CycleNet, N-Linear, FEDformer, Pyraformer, SpaceTime, Multi-Resolution DDPM) trained on ETTm2, and use visual patterns to draw conclusions about model selection, validation-set choice, early stopping, batch size, architectural variations, and cross-family comparisons. The code is released.
Significance. If HAM were a valid layer-agnostic interpretability tool, it could be a useful complement to error metrics for comparing diverse forecasting architectures. The paper's strength is its breadth: it applies a single visualization recipe across MLP, attention, SSM, and diffusion models and releases code, which would facilitate reproducibility. However, the central quantity in Eq. (1) is a scaled cumulative-sum norm, not a per-subseries average of activation magnitudes, and the 'proportionality' baseline is anchored by the observed maximum. These issues make the reported monotonic trends and model-family differences potentially mathematical artifacts rather than model-behavior signals. The paper provides no ground-truth validation linking HAM curves to actual model emphasis. The significance of the claimed utility is therefore not established.
major comments (5)
- [Eq. (1), §2] The masked loss in Eq. (1) divides the cumulative loss by the full horizon length H, not by the length of the masked subseries. Consequently, the causal-mode value at subseries 0→h is ||Σ_{t≤h} ∇θ l_t|| / H, i.e., a scaled norm of a partial sum. The 'gradient norm average' is thus not an average over the subseries; it is a cumulative quantity. Under this definition, a linear proportionality curve is only expected if per-timestep gradient vectors are parallel and equal in norm. In realistic networks, gradient vectors from different timesteps can cancel or be orthogonal, causing the cumulative norm to grow sublinearly even under perfectly uniform per-timestep activity. The monotone increases/decreases of the causal and anti-causal curves, their 'equivariant points,' and their proximity to the proportionality line could therefore reflect the geometry of vector summation rather than the mode
- [Abstract, §1, §3–§4] The central claim—that HAM enables 'granular model selection, validation set choices and comparisons across different neural network model families'—is not validated against any ground truth. There are no ablation studies, no synthetic tasks with known horizon importance, no comparison with alternative attribution methods, and no quantitative measure of whether HAM curves correlate with model behavior. All conclusions are drawn from visual inspection of single training runs without error bars or statistical tests (e.g., Figs. 2–9). For instance, the batch-size claim in §3.2 is based on three configurations and no variance estimate. Without a validation experiment, the paper's empirical conclusions are unfalsifiable and the utility claim is unsupported.
- [§2, 'line of proportionality', Eq. (2)] The line of proportionality is defined using G = max over all subseries of the gradient norm average. This makes the baseline self-referential: every curve intersects or touches the line at its own maximum by construction. The 'areas' computed in Eq. (2) are deviations from an anchor that depends on the observed data, so 'closer to the line' statements (e.g., §3.3, §4.1) are not calibrated across models or runs. A meaningful uniform-activity baseline should be independent of the observed maximum, for example a line based on the expected cumulative norm under a stated assumption about per-step gradients. As written, the proportionality line cannot support claims that a model's activity is 'more uniform' or 'more concentrated' than another's.
- [Abstract vs. §3.2] The abstract states that batch-size differences 'seem to indicate potential for existence of an exponential approximation across them per epoch relative to each other,' while §3.2 and Figure 4 describe 'a polynomially saturating trend' and 'polynomial approximation across batch sizes.' These are contradictory. Moreover, the claim is based on only three batch sizes (500, 2000, 4000) with no repeats, and the interpolated plots use a single model each. The text should commit to one functional form and provide evidence (e.g., residual analysis or multiple seeds) before making this claim.
- [§3.6 (CycleNet cycle queue)] In the architectural-change study, the authors plot 'CycleNet’s cycle queue’s norms ... instead of all the layers’ to prevent too much averaging.' This changes the quantity being visualized from the full-model gradient norm defined in Eq. (1) to a submodule's gradient norm. The paper does not specify how the partial gradient is computed (e.g., whether the loss is masked with respect to all parameters or only the queue parameters), nor does it discuss whether HAM's interpretation remains valid for a submodule. This ambiguity makes the corresponding conclusions about cycle length effects difficult to interpret.
minor comments (6)
- [§2, Eq. (1)] Notation is inconsistent: the loss L is written as a scalar average over h, but the masked loss then appears as a vector or sum. The precise meaning of the '·' product and the index ranges should be clarified. Also, the text says masks are defined over h∈[1,H], but for anti-causal mode the empty/zero cases are not handled explicitly.
- [§3.2, Fig. 4] Figure 4 caption says 'polynomially saturating trend' while the text calls it 'polynomial curve'; this inconsistency should be resolved.
- [§3.3, Fig. 5] The four early-stopping runs are presented as a time series (epochs 20,30,40,50) but no stopping criterion is described. Please specify how 'after early stopping' is defined and what the epochs correspond to.
- [§3.5, Fig. 7] The text says 'High UseFul Load' but the ETTm2 dataset's variable is usually 'HUFL'; please correct the typo.
- [§4.1, Table 3] Area values are mentioned in the text ('anti-causal areas converge to 1') but the table reports only MSE/MAE and plots. It would help to report the numerical area values or point to a supplementary table.
- [General] There are several typos and grammatical issues (e.g., 'alongwith', 'Differences plots over 4 models' in Fig. 5, 'Interpolated Norms w.r.t. Lines of Proportionality' in Fig. 8). A careful proofread is needed.
Circularity Check
HAM's construction is independent, but the 'line of proportionality' benchmark is defined from the observed maximum G, so uniform-distribution claims are self-referential; central model-selection claim is not circular but is under-validated.
-
self definitional
[Section 2, definition of line of proportionality after Eq. (1); used in Sections 3.2–3.5 and 4.1]
"A line of proportionality connecting y=0 and y=G is used in both causal and anti-causal modes corresponding to a uniform distribution of gradient norm averages by the subseries, where G= max ∀i(µθ(∇θL(fθ(x(g)), y(g), i))) represents the maximum gradient norm average across all the subseries."
G is the maximum of the very causal/anti-causal gradient-norm-average curves that the line is used to interpret. The line through (0,0) and (H,G) is a chord determined by the observed curve's own maximum, not an independent uniform-activation baseline. Statements like 'saturating closer to lines of proportionality' (Sec. 3.2) describe the cumulative partial-sum norm's shape relative to its own scale. Also, a linear cumulative norm equals uniform per-timestep gradients only if per-timestep gradient vectors are collinear; Eq. (1) plots the norm of a sum, so that correspondence is an extra assumption, not implied by the definition.
full rationale
No significant circularity in the central construction. HAM is a new visualization computed from trained-model gradients via Eq. (1); the causal/anti-causal curves are not fitted to any target, and no parameter is renamed as a prediction. There are no load-bearing self-citations: the cited NHITS, SpaceTime and MRL works are external, and the paper does not invoke a self-authored uniqueness theorem. The main interpretive claims (proximity to the line of proportionality, gradient equivariant points, area curves) are descriptive statistics of the same cumulative-norm curves; the paper treats these shapes as evidence about model behavior without ground-truth validation, which is a validation/correctness weakness rather than a circular derivation. The one mildly circular component is the line of proportionality: because G is the observed maximum of the curves being interpreted, 'uniform distribution' is not an independent baseline. This affects the strength of some qualitative readings but does not force the central model-selection claim, so it is scored as a minor self-referential element (2/10).
Assumptions & free parameters
assumptions (4)
- domain assumption The gradient norm of the masked loss reflects the model's emphasis on each subseries.
- domain assumption ETTm2 dataset is representative and sufficient for cross-model comparisons.
- domain assumption The line of proportionality (straight line from (0,0) to (H,G)) represents a uniform distribution of gradient norm averages.
- domain assumption NHITS's neural approximation theorem and SpaceTime's exponential autoregressive properties are applicable to explain HAM trends.
Cite this review
Pith. "Pith review of Horizon Activation Mapping for Neural Networks in Time Series Forecasting." pith.science (2026). https://pith.science/paper/REZR4ZY5
@misc{pith2026260102094,
author = {Pith},
title = {Pith review of: Horizon Activation Mapping for Neural Networks in Time Series Forecasting},
year = {2026},
howpublished = {\url{https://pith.science/paper/REZR4ZY5}},
note = {Machine review of arXiv:2601.02094}
}
read the original abstract
Neural networks for time series forecasting have relied on error metrics and architecture-specific interpretability approaches for model selection that don't apply across models of different families. To interpret forecasting models agnostic to the types of layers across state-of-the-art model families, we introduce Horizon Activation Mapping (HAM), a visual interpretability technique inspired by grad-CAM that uses gradient norm averages to study the horizon's subseries where grad-CAM studies attention maps over image data. We introduce causal and anti-causal modes to calculate gradient update norm averages across subseries at every timestep and lines of proportionality signifying uniform distributions of the norm averages. Optimization landscape studies with respect to changes in batch sizes, early stopping, train-val-test splits, architectural choices, univariate forecasting and dropouts are studied with respect to performances and subseries in HAM. Interestingly, batch size based differences in activities seem to indicate potential for existence of an exponential approximation across them per epoch relative to each other. Multivariate forecasting models including MLP-based CycleNet, N-Linear, N-HITS, self attention-based FEDformer, Pyraformer, SSM-based SpaceTime and diffusion-based Multi-Resolution DDPM over different horizon sizes trained over the ETTm2 dataset are used for HAM plots in this study. NHITS' neural approximation theorem and SpaceTime's exponential autoregressive activities have been attributed to trends in HAM plots over their training, validation and test sets. In general, HAM can be used for granular model selection, validation set choices and comparisons across different neural network model families.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[1]
In: CCCG (2000)
Boland, R.P., Urrutia, J.: Polygon area problems. In: CCCG (2000)
2000
-
[2]
Hans and V.A
Challu, C., Olivares, K.G., Oreshkin, B.N., Ramirez, F.G., Canseco, M.M., Dubrawski, A.: Nhits: Neural hierarchical interpolation for time series forecast- 12 K. Hans and V.A. Kandappan ing. In: Proceedings of the AAAI conference on artificial intelligence. vol. 37, pp. 6989–6997 (2023)
2023
-
[3]
In: 1st ICML Workshop on Foundation Models for Structured Data (2025), https: //openreview.net/forum?id=eVGS73kXEu
Cohen, B., Khwaja, E., Doubli, Y., Lemaachi, S., Lettieri, C., Masson, C., Mic- cinilli, H., Ramé, E., Ren, Q., Rostamizadeh, A., du Terrail, J.O., Toon, A.M., Wang, K., Xie, S., Xu, Z., Zhukova, V., Asker, D., Talwalkar, A., Abou-Amal, O.: Toto: An open time series foundation model optimized for observability. In: 1st ICML Workshop on Foundation Models f...
2025
-
[4]
Contreras,F.:Cuttingpolygonsandaproblemonilluminationofstages.University of Ottawa (Canada) (1998)
1998
-
[5]
In: Forty-first International Conference on Machine Learning (2024), https://openreview.net/forum?id=jn2iTJas6h
Das, A., Kong, W., Sen, R., Zhou, Y.: A decoder-only foundation model for time- series forecasting. In: Forty-first International Conference on Machine Learning (2024), https://openreview.net/forum?id=jn2iTJas6h
2024
-
[6]
Kusupati, A., Bhatt, G., Rege, A., Wallingford, M., Sinha, A., Ramanujan, V., Howard-Snyder, W., Chen, K., Kakade, S., Jain, P., et al.: Matryoshka representa- tionlearning.AdvancesinNeuralInformationProcessingSystems35,30233–30249 (2022)
2022
-
[7]
Advances in Neural Infor- mation Processing Systems37, 106315–106345 (2024)
Lin, S., Lin, W., Hu, X., Wu, W., Mo, R., Zhong, H.: Cyclenet: Enhancing time series forecasting through modeling periodic patterns. Advances in Neural Infor- mation Processing Systems37, 106315–106345 (2024)
2024
-
[8]
In: International conference on learning representations (2021)
Liu, S., Yu, H., Liao, C., Li, J., Lin, W., Liu, A.X., Dustdar, S.: Pyraformer: Low-complexity pyramidal attention for long-range time series modeling and fore- casting. In: International conference on learning representations (2021)
2021
Show all 16 references
-
[9]
arXiv preprint arXiv:2412.12953 (2024)
Reuss, M., Pari, J., Agrawal, P., Lioutikov, R.: Efficient diffusion transformer policies with mixture of expert denoisers for multitask learning. arXiv preprint arXiv:2412.12953 (2024)
2024 arXiv
-
[10]
In: Proceedings of the IEEE international conference on computer vision
Selvaraju, R.R., Cogswell, M., Das, A., Vedantam, R., Parikh, D., Batra, D.: Grad- cam: Visual explanations from deep networks via gradient-based localization. In: Proceedings of the IEEE international conference on computer vision. pp. 618–626 (2017)
2017
-
[11]
In: The Twelfth International Conference on Learning Representations (2024)
Shen, L., Chen, W., Kwok, J.: Multi-resolution diffusion models for time series forecasting. In: The Twelfth International Conference on Learning Representations (2024)
2024
-
[12]
Zeng, A., Chen, M., Zhang, L., Xu, Q.: Are transformers effective for time se- ries forecasting? In: Proceedings of the AAAI conference on artificial intelligence. vol. 37, pp. 11121–11128 (2023)
2023
-
[13]
arXiv preprint arXiv:2303.09489 (2023)
Zhang, M., Saab, K.K., Poli, M., Dao, T., Goel, K., Ré, C.: Effectively model- ing time series with simple discrete state spaces. arXiv preprint arXiv:2303.09489 (2023)
2023 arXiv
-
[14]
net/forum?id=kQDzqIkXLO
Zhang, Y., Zeng, Q., Zhang, Y., Xu, Z., Zheng, M., Gao, C., Jiang, M., Zheng, Z.: Are time series foundation models ready for zero-shot forecasting? In: 1st ICML Workshop on Foundation Models for Structured Data (2025), https://openreview. net/forum?id=kQDzqIkXLO
2025
-
[15]
In: Proceed- ings of the AAAI conference on artificial intelligence
Zhou, H., Zhang, S., Peng, J., Zhang, S., Li, J., Xiong, H., Zhang, W.: Informer: Beyond efficient transformer for long sequence time-series forecasting. In: Proceed- ings of the AAAI conference on artificial intelligence. vol. 35, pp. 11106–11115 (2021)
2021
-
[16]
In: International conference on machine learning
Zhou, T., Ma, Z., Wen, Q., Wang, X., Sun, L., Jin, R.: Fedformer: Frequency en- hanced decomposed transformer for long-term series forecasting. In: International conference on machine learning. pp. 27268–27286. PMLR (2022)
2022
Reviewed August 3, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.