FourierGNN: Rethinking Multivariate Time Series Forecasting from a Pure Graph Perspective

Multivariate time series (MTS) forecasting has shown great importance in numerous industries. Current state-of-the-art graph neural network (GNN)-based forecasting methods usually require both graph networks (e.g., GCN) and temporal networks (e.g., LSTM) to capture inter-series (spatial) dynamics and intra-series (temporal) dependencies, respectively. However, the uncertain compatibility of the two networks puts an extra burden on handcrafted model designs. Moreover, the separate spatial and temporal modeling naturally violates the unified spatiotemporal inter-dependencies in real world, which largely hinders the forecasting performance. To overcome these problems, we explore an interesting direction of directly applying graph networks and rethink MTS forecasting from a pure graph perspective. We first define a novel data structure, hypervariate graph, which regards each series value (regardless of variates or timestamps) as a graph node, and represents sliding windows as space-time fully-connected graphs. This perspective considers spatiotemporal dynamics unitedly and reformulates classic MTS forecasting into the predictions on hypervariate graphs. Then, we propose a novel architecture Fourier Graph Neural Network (FourierGNN) by stacking our proposed Fourier Graph Operator (FGO) to perform matrix multiplications in Fourier space. FourierGNN accommodates adequate expressiveness and achieves much lower complexity, which can effectively and efficiently accomplish the forecasting. Besides, our theoretical analysis reveals FGO's equivalence to graph convolutions in the time domain, which further verifies the validity of FourierGNN. Extensive experiments on seven datasets have demonstrated our superior performance with higher efficiency and fewer parameters compared with state-of-the-art methods.

Paper

References (39)

Scroll for more · 27 remaining

Similar papers

Peer review

Reviewer yaJh7/10 · confidence 4/52023-07-01

Summary

The paper presents FourierGNN for multivariate time series forecasting from a pure graph perspective, performing matrix multiplications in Fourier space, which has not been investigated so far. They design a new hypervariate graph structure to consider spatiotemporal dynamics unitedly and reformulate the graph operations on the hypervariate graph in Fourier space. Extensive experiments have demonstrated superior performance with higher efficiency and fewer parameters over state-of-the-art methods.

Strengths

1. The paper is overall well-written. I think the contribution and novelty of the proposed work are clearly presented. 2. The idea of directly applying graph networks for multivariate time series forecasting seems novel and exciting. The hypervariate graph structure can encode spatiotemporal dependencies unitedly, and the visualizations in the experiments have also verified its advantages. 3. The authors argue that performing the multiplication between the input and the proposed FGO in Fourier space is equivalent to a graph convolution in the time domain, while the multiplications in Fourier space have lower complexity. The argument is proved by the authors and seems very meaningful. I think this argument provides a new path for conducting graph operations in Fourier space. 4. The effectiveness is validated through extensive experiments across seven real-world datasets. Sufficient analysis including efficiency analysis and ablation study is also provided.

Weaknesses

1. In the efficiency analysis section, the authors assess the overall efficiency of the forecasting process. To specifically evaluate the efficiency of FourierGNN, it would be more appropriate to focus on comparing the number of parameters involved in the graph operations with those of the baselines. 2. The MAPE (Mean Absolute Percentage Error) results on some datasets (e.g., Traffic, COVID-19) do not achieve state-of-the-art performance.

Questions

FourierGNN demonstrates impressive performance across multiple datasets, with notable success observed on the COVID dataset. Could you explain the reasons that contribute to the superior performance of FourierGNN on the COVID dataset compared to other datasets?

Rating

7: Accept: Technically solid paper, with high impact on at least one sub-area, or moderate-to-high impact on more than one areas, with good-to-excellent evaluation, resources, reproducibility, and no unaddressed ethical considerations.

Confidence

4: You are confident in your assessment, but not absolutely certain. It is unlikely, but not impossible, that you did not understand some parts of the submission or that you are unfamiliar with some pieces of related work.

Soundness

3 good

Presentation

3 good

Contribution

3 good

Limitations

N/A

Reviewer ipTP7/10 · confidence 4/52023-07-05

Summary

This paper has studied a popular and important problem, i.e., modeling the intricate spatial and temporal dependencies among multivariate time series for accurate forecasting. To overcome the main limitation that existing works always separately model spatial and temporal, this work designs a new hypervariate graph structure to encode the spatiotemporal dynamics unitedly, and proposes a novel FourierGNN to learn the spatiotemporal dependencies. Extensive experiments on seven real-world datasets demonstrate FourierGNN achieves good performances in both accuracy and efficiency compared with state-of-the-art methods.

Strengths

This paper presents a novel formulation from a pure graph perspective to model spatiotemporal dependencies for MTS forecasting. The work is interesting and original, which is different from previous GNN-based methods that typically model spatial and temporal dependencies with distinct graph and temporal networks. The novel pure graph modeling is straightforward yet quite meaningful for MTS forecasting, which brings up good inspiring insights. This paper proposes a graph neural network, namely FourierGNN, to cooperate with the pure graph formulation and provides a theoretical guarantee of effectiveness. This paper seems to have solid, extensive and diverse experiments, which are conducted on seven real-world different datasets. The extensive experimental results have clearly demonstrated the obvious performance improvement of FourierGNN in terms of accuracy and efficiency over state-of-the-art MTS forecasting models.

Weaknesses

I have two concerns. My first concern is about the hypervariate graph structure. Although the authors stated that all variables at all timestamps can encode high-resolution relationships, it also may introduce some redundant information or unwanted correlations, such as some variables that are far apart in time. How can this be properly handled? Besides, I find that the core operations of FourierGNN are conducted in the Fourier space, so I am curious whether transferring the hypervariate graph structure into the Fourier space help to learn dependencies. In addition, another concern is whether FourierGNN can also be applied to other domains apart from the multivariate forecasting. Maybe it can be extended to other graph tasks. Minor typos: In the caption of Figure 2, Given the hypervariate graph $\mathcal{G}=(X^\mathcal{G}_t, A^\mathcal{G}_t)$ should be $\mathcal{G}_t=(X^\mathcal{G}_t, A^\mathcal{G}_t)$

Questions

1.Besides MTS forecasting, can FourierGNN be applied to other domains as well? 2.How does the proposed model handle the redundant information in the hypervariate graph structure?

Rating

7: Accept: Technically solid paper, with high impact on at least one sub-area, or moderate-to-high impact on more than one areas, with good-to-excellent evaluation, resources, reproducibility, and no unaddressed ethical considerations.

Confidence

4: You are confident in your assessment, but not absolutely certain. It is unlikely, but not impossible, that you did not understand some parts of the submission or that you are unfamiliar with some pieces of related work.

Soundness

4 excellent

Presentation

3 good

Contribution

4 excellent

Limitations

n/a

Reviewer ndBv3/10 · confidence 4/52023-07-06

Summary

The paper addresses the time series forecasting problem. The authors propose a model that represents each scalar observation as a node in a (fully connected) graph, encodes the nodes and finally regresses the future observations on all previous encoded observations (via a fully connected layer). As encoding they propose a discrete fourier transform followed by several node-wise fully connected layers and finally an inverse fourier transform. In experiments on 6 datasets they show improvements over several baselines.

Strengths

s1. interesting problem: time series forecasting s2. interesting approach: parametrizing transformations in fourier space s3. results showing improvements over many baselines

Weaknesses

w1. a main aspect of the proposed model is not clear. w2. missed recent related work. w3. experimental protocol is unclear and deviates from the literature. w4. for a key component of the proposed model, there is no ablation study demonstrating its impact.

Questions

w1. a main aspect of the proposed model is not clear. - the authors describe the fourier graph layer as applying the discrete fourier transform on the graph features X\in\R^{nd\times 1} of a complete graph, explicitly not the graph fourier transform. Which dimension is the time dimension for this fourier transform? The node dimension of X contains nd many entries, one for each time point and channel. While one can compute a fourier transform also over this vector, what is it semantics? How does it capture time varying information? And would this not crucially depend on how on is ordering these nd many entries? w2. missed recent related work. - esp. PatchTST [Nie et al. 2023] and dlinear [Zeng et al. 2023] are well known model outperforming your baselines by a good margin. w3. experimental protocol is unclear and deviates from the literature. - what is the forecasting horizon \tau you are using? - the experiments report only a single number per dataset and error measure, while in the literature usually results for a portfolio of forecasting horizons are reported (e.g., 96, 192, 336, 720, e.g., the Fedformer). Why do you deviate from this standard? This way it also is not possible to compare your results with the results published in those papers. - usually results are also reported for further datasets such as ETTm2, Exchange, Weather and ILI. How does your model compare to those published results? w4. for a key component of the proposed model, there is no ablation study demonstrating its impact. - while there are many papers that simply operate on the fourier spectrum of a time series, one key component of the model of the authors seems to be that after these operations they move back into the time domain by the inverse fourier transform. What is the impact of this back transformation? references: - Nie, Yuqi, Nam H. Nguyen, Phanwadee Sinthong, and Jayant Kalagnanam. “A Time Series Is Worth 64 Words: Long-Term Forecasting with Transformers.” arXiv, March 5, 2023. https://doi.org/10.48550/arXiv.2211.14730. - Zeng, Ailing, Muxi Chen, Lei Zhang, and Qiang Xu. “Are Transformers Effective for Time Series Forecasting?” In AAAI, 2023.

Rating

3: Reject: For instance, a paper with technical flaws, weak evaluation, inadequate reproducibility and incompletely addressed ethical considerations.

Confidence

4: You are confident in your assessment, but not absolutely certain. It is unlikely, but not impossible, that you did not understand some parts of the submission or that you are unfamiliar with some pieces of related work.

Soundness

3 good

Presentation

3 good

Contribution

2 fair

Limitations

n/a

Reviewer G8af6/10 · confidence 5/52023-07-08

Summary

In this paper, the authors study a problem in GNN-based multivariate time series (MTS) forecasting, i.e. modeling spatial correlations and temporal dependencies in the same time. In particular, the authors do not follow previous works on regarding the input as T graphs and capturing temporal dependencies between graphs by temporal networks, but propose to formulate T graphs as a hypervariate graph(a pure fully-connected graph with NT nodes). Specifically, to deal with the huge node amount NT in a hypervariate graph, the authors propose the FourierGNN which utilize matrix multiplications in the Fourier space of graphs to decrease the quadratic computational complexity to quasi linear one. Empirical results on seven datasets show the effectiveness of this method.

Strengths

1. The authors provide a novel view of graph-based MTS forecasting called hypervariate graph. 2. The authors propose Fourier Graph Operator and FourierGNN for efficient convolution on hypervariate graph. 3. In the authors’ experiments, FourierGNN obtains superior performance.

Weaknesses

1. The detail of the experiments is not clear. The authors run previous models on different datasets but do not show whether some hyperparameters, such as feature dim, change with the dataset. 2. The authors claim that FourierGNN is efficient for hypervariate graph, but do not provide other networks’ performance such as GCN on hypervariate graph. 3. The models compared with Efficiency Analysis are too old. 4. The authors do not explain the reason for the Green kernel. The Green kernel is a strong assumpition. In line 161, k_{ij} = A_{ij} W and k_{ij} = k_{i-j} results A_{j+k}_{j} equals for any j, which strictly restricts the topology of the hypervariate graph. Besides, there is no visualization of A in each layer to check whether the green kernel assumption is satisfied during training.

Questions

1. What’s the performance of FourierGNN on longer prediction lengths such as {24, 36, 48, 60} in FEDformer? 2. What’s the meaning of “the same sparsity pattern” in line 196? 3. How to determine the value of A_i in equation (7), A_0 is identity and what about A_1, A_2?

Rating

6: Weak Accept: Technically solid, moderate-to-high impact paper, with no major concerns with respect to evaluation, resources, reproducibility, ethical considerations.

Confidence

5: You are absolutely certain about your assessment. You are very familiar with the related work and checked the math/other details carefully.

Soundness

3 good

Presentation

3 good

Contribution

3 good

Limitations

Please check the weakness

Reviewer ndBv2023-08-18

answer to rebuttal

Dear authors, thanks for your extensive answers and additional experiments. My concerns w1 and w4 you resolved. About the other two: w2. PatchTST and dlinear - The numbers you report for e.g. PatchTST deviate from the published numbers, e.g., for dataset weather, horizon 96, MAE - you report 0.034 (table 2 in your rebuttal pdf), - the PatchTST paper reports 0.198. Can you explain these differences? It would be more convincing to reproduce the experimental settings of the baseline papers (and their numbers). w3. experimental protocol unclear and deviates from the literature. 1. for the main experiment in tab. 1: what is \tau? 1 ? 2. what does "follow the experimental settings in short-term forecasting baselines, like LSTNet and StemGNN" mean exactly? - The LSTNet paper reports two different error measures, RSE and corr, so one cannot compare numbers directly. - The StemGNN paper reports RMSE for datasets Solar and Electricity, but they report different numbers than you do (their table 2): - Solar: they report 0.07, you report 0.222 - Electricity: they report 0.06, you report 0.101 Are you using a vastly different split? Or how are these differences explained? 3. Do any of the numbers in your table 1 coincide with some published numbers in the baseline papers? And if so, could you mark them, say with a star?

Authorsrebuttal2023-08-19

Thanks for your feedback (1/2)

Dear Reviewer ndBv, We greatly appreciate your feedback, and thanks for your careful readings. We would like to clarify our experimental settings and the two points you mentioned. **About experimental settings** In the literature on short-term forecasting, previous models have employed diversified experimental settings in their experiments. - They use **different normalization methods**. For example, LSTNet normalizes raw data row by row using the maximum absolute value, short for max-abs normalization; GraphWaveNet uses Z-score normalization; StemGNN uses Z-score normalization for some datasets and min-max normalization for other datasets. - They use **different data splitting ratios**. For example, LSTNet uses 6:2:2; DCRNN uses 7:2:1; MTGNN and StemGNN use 6:2:2 for some datasets and 7:2:1 for other datasets. - They use **different prediction lengths**. For example, GraphWaveNet uses \{3,6,12\}; STGCN uses \{3,6,9\}; LSTNet uses \{3,6,12,24\}; AGCRN uses \{12\}; MTGNN uses \{1, 12\}; StemGNN uses different prediction lengths for different datasets, such as \{3,12,28\}. It is important to note that due to the significant variations in experimental settings among different baselines, we did not directly replicate the results reported in the baseline papers for our paper. In our work, 1. We first unify the experimental settings to guarantee a fair and more convincing comparison. Specifically, we - adopt the min-max normalization and the splitting ratio of 7:2:1 for all datasets, - fix the input length 12 and the output length 12 in Table 1, and - set the input length 12 and the output length \{3,6,9,12\} for the multi-step forecasting. An exception is made for COVID-19 where we adopt the ratio of 6:2:2 because the number of samples in COVID-19 is too small (355). Please refer to **Lines 240-244** for more details. 2. We then re-run **all baselines under the above settings on all datasets for both the experiments presented in our paper and those conducted during the rebuttal phase**. Despite the substantial workload and large time consumption, we firmly believe that employing a unified experimental framework and reproducing the results ensures a fairer comparison, ultimately contributing to the advancement of short-term forecasting. We believe our dedicated efforts in conducting these experiments coincide with your rigorous attitude and expectations toward our experimental settings/results. Furthermore, due to the different experimental settings, **the results of the baselines reported in our paper may differ from those reported in their original papers**. This discrepancy is particularly prominent for baselines that utilize different normalization methods, such as LSTnet using max-abs normalization as default, StemGNN using Z-score or min-max normalization, and our settings using min-max normalization. This disparity is expected, as data normalized using different methods can have varying scales, while most multivariate time series (MTS) forecasting methods evaluate results on normalized data.

Authorsrebuttal2023-08-19

Thanks for your feedback (2/2)

Thanks for your valuable time and careful review. **W2** **The numbers you report for e.g. PatchTST deviate from the published numbers. Can you explain these differences?** The differences are attributed to the fact that we adopt a different normalization method (**min-max normalization**) in the experiment to that (**standard normalization**) in the PatchTST paper. - As we stated in **Lines 241-242**, all datasets are normalized using the **min-max normalization**. Keeping consistent with the setting, we also performed the **min-max normalization** on the four datasets, i.e., Exchange, Weather, ETTh1, and ETTm1. Accordingly, we rerun PatchTST, DLinear, and our FourierGNN on the normalized data with the same input length of 96. The corresponding results are reported for comparison in **Table 2 in the attached PDF**. - In addition, we aim to investigate the long-term forecasting performance of our FourierGNN compared with two SOTA long-term forecasting performances. Under the same experimental settings, the results in Table 2 in the rebuttal PDF verify that our FourierGNN performs worse than the baselines under long-term forecasting settings, which confirms our expectation that our FourierGNN is designed for short-term forecasting and unsuitable for long-term forecasting. - Honestly, comparing FourierGNN with PatchTST and DLinear under the settings of the baseline papers would be more appropriate. We reperformed the experiments and rerun our FourierGNN under the same normalization settings of the PatchTST paper. We report the corresponding results in the following table, where we can achieve the same conclusion that FourierGNN does not perform as well as PatchTST and DLinear under large horizons, i.e., the long-term forecasting settings. ETTh1 | Model| Metric|96 | 192 | 336 | 720 | |:--------|:----|:----|:----|:----|:----| |PatchTST|MAE| 0.370| 0.413| 0.422| 0.447| | |MSE|0.400 | 0.429| 0.440| 0.468| |DLinear|MAE| 0.375| 0.405| 0.439| 0.472| | |MSE|0.399 | 0.416| 0.443| 0.490| |FourierGNN|MAE| 0.433| 0.486| 0.507| 0.564| | |MSE| 0.451| 0.553| 0.594| 0.663| **W3** **1. for the main experiment in tab. 1: what is \tau? 1 ?** The $\tau$ in Table 1 is 12, and we will mark this in the caption of Table 1. **2. what does "follow the experimental settings in short-term forecasting baselines, like LSTNet and StemGNN" mean exactly?** The sentence "follow the experimental settings in short-term forecasting baselines, like LSTNet and StemGNN" in the rebuttal means that we adopt $\tau\in\{3,6,9,12\}$ in our experiments, following the two papers. **2.1 The LSTNet paper reports two different error measures, RSE and corr, so one cannot compare numbers directly.** Since the metrics, the datasets, and the experimental settings of LSTNet are different from our paper and some other baseline papers, we re-run the source code of LSTNet according to our experimental settings and datasets. To ensure a comprehensive assessment, we employed three evaluation metrics in MTS forecasting, namely MAE, RMSE, and MAPE. These evaluation metrics are also commonly utilized in MTS forecasting research, as exemplified by the baselines StemGNN and AGCRN. **2.2 The StemGNN paper reports RMSE for datasets Solar and Electricity, but they report different numbers than you do (their table 2)** As aforementioned, the different numbers are derived from the different experimental settings between the StemGNN paper and ours. In the StemGNN paper, for Solar and Electricity datasets, the input length is 24 and the prediction length is 3; while in our settings, both the input length and prediction length are 12. Meanwhile, StemGNN uses Z-score normalization, whereas our model employs min-max normalization. **3. Do any of the numbers in your table 1 coincide with some published numbers in the baseline papers? And if so, could you mark them, say with a star?** As explained in our responses to the experimental settings, in our work, we adopt the unified experiments settings for all methods for a fair comparison. Since the experimental settings are different in each baseline paper, we re-run all baseline source codes and report the result numbers carefully. Correspondingly, the results of the baselines reported in our paper may differ from those reported in their original papers. We believe we could advance short-term time series forecasting by using such unified settings for fair comparisons. Thanks again for your feedback. Hope we have addressed all your concerns. Authors

Program Chairsdecision2023-09-21

Decision

Accept (poster)

© 2026 NYSGPT2525 LLC