## Enhancing Multistep Prediction of Multivariate Market Indices Using Weighted Optical Reservoir Computing

Fang Wang , Ting Bu , Yuping Huang a b a,b

a Department of Physics, Stevens Institute of Technology, Hoboken, 07030, New Jersey, USA b Quantum Computing Inc., Hoboken, 07030, New Jersey, USA

## Abstract

We propose and experimentally demonstrate an innovative stock index prediction method using a weighted optical reservoir computing system. We construct fundamental market data combined with macroeconomic data and technical indicators to capture the broader behavior of the stock market. Our approach shows significant higher performance than state-of-the-art methods such as linear regression, decision trees, and neural network architectures including long short-term memory. It captures well the market's high volatility and nonlinear behaviors despite limited data, demonstrating great potential for real-time, parallel, multi-dimensional data processing and predictions.

Keywords: Optical reservoir computing, Parallel data processing, Stock index prediction, LSTM, Regression, Multivariate Forecasting

## 1. Introduction

Over the past decades, machine learning (ML) and deep learning have emerged as a transformative tool in the field of market prediction, particularly for predicting market indices [1]. The primary challenge to generate accurate prediction is aggregating diverse information for the benchmarks and constructing a reliable model. Techniques including regression and classification, as well as more advanced methods like deep learning and ensemble models, offer sophisticated methods for analyzing data to forecast stock and

Email address: yhuang5@stevens.edu (Yuping Huang)

market trends, [2, 3, 4, 5, 6, 7, 8, 9, 10, 11]. Especially, long short-term memory (LSTM) networks, a specialized form of recurrent neural networks (RNNs), have gained significant attention in the realm of stock prediction due to their ability to model temporal dependencies and sequence data effectively [12, 13, 14, 15, 16, 17, 18]. Together, ML models and LSTM can incorporate diverse data sources, including historical prices, technical indicators [19], and macroeconomic indicators to generate comprehensive and nuanced forecasts.

However, existing ML research on market predictions has suffered several limitations. One significant drawback is that most methods require a large amount of dense historical data for training. This requirement can be resource-intensive and may not always be feasible. Additionally, these models often struggle to accurately capture market peaks and volatility, leading to less reliable forecasts during periods of rapid market changes. As such, these limitations hinder the effectiveness of stock prediction models in providing timely and precise insights, particularly under volatile market conditions.

Recently, reservoir computing (RC) has gained prominence as a streamlined yet effective neural network architecture. It has been demonstrated to accurately predict time-series data with a limited number of training samples [20, 21, 22]. This architecture holds significant potential as an alternative for traditional, complex neural network structures in processing financial datasets, particularly for real-time data processing when storage capacity for large datasets is limited and/or time-to-prediction is sensitive. An RC system comprises a reservoir that maps inputs into a high-dimensional space and a readout that analyzes patterns from these high-dimensional states [23]. In this system, the reservoir architecture remains fixed, and only the readout is trained using simple methods such as linear regression and classification [24]. Therefore, reservoir computing offers fast learning and low training costs compared with other RNNs. Additionally, the static nature of the reservoir facilitates hardware implementation using a variety of optical systems, substrates, and devices. As such, optical reservoir computing has garnered increasing attention across diverse research fields[23, 25, 26, 27, 22]. These RC systems can be categorized into three main types: free-space optical RC, fiber-based temporal RC, and integrated photonic RC. In comparison, the free-space optical RC offers significant advantages in large parallelism and system optimization flexibility.

In this paper, we introduce and demonstrate a free-space optical RC system for predicting daily stock index prices. Utilizing a spatial light modulator (SLM), the temporal data from benchmarks is transformed into high-

Figure 1: Experiment diagram for stock benchmark prediction to compare spatial reservoir computing and competitive machine learning models.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00652v1-figure-5c02c953-2390-4601-821b-82a990471595.jpg)

dimensional optical signals in both spatial and temporal domains for highquality benchmark predictions. Owing to an innovative design of the mapping method and the large number of pixels on the modulator, the system is flexible to implement the training samples with different feature sizes simultaneously. Therefore, the information can be processed in parallel in the optical domain without compromising the time for signals' interconnection or nonlinear effect as in the temporal RC [21]. Furthermore, this large-scale parallelism allows multiple-step prediction, by translating temporal correlation in the time-series data to spatial coherence in the free-space optics, and conveniently encoding and reading out in the spatial domain.

Weselect seven popular stock indices and gather 500 historical data points for each. Additionally, we incorporate seven macroeconomic and technical analysis features as input variables for the model training. We further optimize the optical RC system by incorporating feature correlation as weights. Our comparative analyses with popular ML pricing models indicate that our system outperforms most ML models for both one-step and multi-step predictions. Both ML and LSTM models serve as significant reference points in experimental financial research and act as essential benchmarks when compared with our optical RC experiment. Unlike previous research requiring extensive historical data, our system efficiently achieves parallel real-time predictions using no more than 500 data points with multi-dimensional features. Notably, the inclusion of feature correlation weights significantly enhances prediction performance by capturing signals and peaks instantly. The feature-weighted optimized RC demonstrates substantial potential for confidential, low-latency edge processing and practical pricing strategies.

## 2. Model

## 2.1. Experiment diagram

Figure 1 illustrates the experimental procedures described in this paper. Initially, we download historical trading datasets from stock benchmarks and perform basic data and feature engineering. Subsequently, we apply optical RC and various popular ML techniques to predict future stock index prices, respectively. We compare different machine learning models, utilizing crossvalidation and hyperparameter tuning to select the ML model with the lowest error for each benchmark. Finally, we evaluate our models and compare the predictions from optical RC with the best-performing ML model for each benchmark.

## 3. Datasets preparation

This study examines seven stock indices frequently analyzed in investment strategies and portfolio management, as referenced in [11]. The primary datasets consist of weekday historical data from May 2023 to September 2023, sourced from Yahoo Finance. These datasets include crucial pricing details, specifically the daily closing prices. Comprehensive ticker symbols and index information are provided in Table I. The closing prices of all the stock indices are inherently noisy, non-linear, and chaotic.

Table 1: Global market indices used in this paper.

| Ticker   | Index                                        | Region   |
|----------|----------------------------------------------|----------|
| VTHR     | Vanguard Russell 3000 ETF                    | US       |
| ^N225    | Nikkei 225                                   | Japan    |
| UKX.L    | The Financial Times Stock Exchange 100 Index | UK       |
| ^NYA     | The NYSE Composite                           | US       |
| ^IXIC    | NASDAQ Composite                             | US       |
| ^HSI     | HANG SENG INDEX                              | Chine    |
| MME=F    | MSCI Emerging Markets Index Fut              |          |

Due to the complexity of the benchmarks, a prediction model based solely on the closing price cannot gather all the necessary information for accurate predictions. Consequently, feature selection is a crucial part of the prediction enhancement. However, it is challenging to extract features from financial data. On the one hand, limiting the features can hinder the predictive model's performance. On the other hand, including all available features

from the financial market can lead to a dimensionally complex and difficultto-interpret model. Additionally, collinearity among multiple variables can negatively impact the model's performance. Some studies select the features rely solely on technical indicators [28, 29, 30, 31], while others incorporate historical data and macroeconomic indicators [32, 33, 34, 35, 36].

Table 2: Features used for stock prediction in this paper.

| Fundamental data closing price closing price(t-1)                | Source Yahoo Finance Yahoo Finance   | Abbreviation Close Close(t-1)   |
|------------------------------------------------------------------|--------------------------------------|---------------------------------|
| Macroeconomic data Cboe volatility index Effective Federal Funds | Source Yahoo Finance                 | Abbreviation VIX EFFR           |
| Consumer sentiment index                                         | FRED FRED                            | UMSCENT                         |
| Rate                                                             |                                      |                                 |
| US dollar index                                                  | Yahoo Finance                        | DXYNYB                          |
| Technical analysis data                                          |                                      | Abbreviation                    |
|                                                                  | Source                               |                                 |
| Moving average convergence divergence Average true range         | Calculated Calculated                | MACD ATR                        |
| Relative strength index                                          | Calculated                           | RSI                             |

Figure 2: (a) An overview of one stock index sequence ˆ NYA with training and testing parts. (b) A feature correlation matrix heatmap for ˆ NYA

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00652v1-figure-ff57366f-2fd6-4de9-8eb0-e866a4732d4b.jpg)

In the above considerations, we select seven main features either from macroeconomic indicators or technical indicators as the prediction inputs, as listed in Table 2. These features are chosen based on their high correlation

with the closing price [37]. As macroeconomic data can significantly impact the stock market, we consider the Cboe Volatility Index (VIX), the Interest Rate (EFFR), the Consumer Sentiment Index (UMCSENT), and the US Dollar Index (DXYNYB). The technical indicators include Moving Average Convergence Divergence (MACD), Average True Range (ATR), and Relative Strength Index (RSI). The details of the features are listed in Table 2. The correlation matrix among features is illustrated in Figure 2(b). These values are depicted based on the intensity of the color, indicating the strength of the relationship between the given variables. In the following RC prediction, we add the bottom row of the correlation weight (CW) matrix as extra input weights W cor . For each index, we used 500 days as the training data and 100 days as the testing data. Each day, the training input is the previous day's closing price and the seven selected features.

## 3.1. Optical RC

Our spatial optical RC follows the basic working principle of conventional RC, which can be described by

$$x ( t ) = f [ \alpha x ( t - 1 ) W _ { r e s } + \beta u ( t ) W _ { i n } ],$$

where x t ( ) represents the reservoir state at time t . f is the nonlinear function. α and β are the feedback gain and input gain, which will be optimized for the best performance. W res and W in are the reservoir nodes' interconnection and input weight, respectively. Due to this updating rule of reservoir states, the current state x t ( ) inherently contains information from the preceding state x t ( -1) . This endows reservoir computers with the crucial property of short-term memory, a feature for effective time-series data prediction. In the training stage, the input is injected sequentially into the reservoir, and the resulting reservoir states x for all time instances within the range [1, T] are recorded. Subsequently, the system's output weights ( W out ) can be determined by minimizing the error between the predicted results ( ˆ = y W out x ) and their corresponding true time series ( y ). The prediction performance for spatial RC and ML models are evaluated using the normalized root mean square error (NRMSE) between the true time series y and the prediction results ˆ y on the testing dataset

$$\text{NRME} ( y, \hat { y } ) = \frac { \sqrt { \sum _ { i = 0 } ^ { N - 1 } ( y _ { i } - \hat { y } _ { i } ) ^ { 2 } } } { \sqrt { \sum _ { i = 0 } ^ { N - 1 } ( y _ { i } - \bar { y } ) ^ { 2 } } },$$

Figure 3: The configuration of the spatial reservoir computer. SLM: spatial light modulator, BS: beam splitter, HWP: half-wave plate, QWP: quarter-wave plate, PBS: polarizing beam splitter, PM: power meter

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00652v1-figure-5ce52965-67d5-4f8e-831f-75a65c69ab4b.jpg)

where ¯ y is the mean value of the true time series.

In our optical setup, a continuous-wave 775 nm laser diode serves as the light source, as illustrated in Fig. 3. The beam waist is adjusted to 1 mm after collimation. To match the effective area on the spatial light modulator (SLM), a 10 × beam expander is used to enlarge the beam size. A power meter is implemented to capture the transmitted light passing through the polarizing beam splitter (PBS) for real-time monitoring of system stability during training and testing. Before the beam splitter (BS), we employ a half-wave plate to rotate the light's polarization, thereby maximizing the modulation effect on the SLM. Following modulation, the light is focused, and a camera is positioned in the time-fourier domain [38] to align the sensor area with the beam size for subsequent information processing.

The inset in Fig.3 illustrates the data flow diagram in our spatial RC. One key here is the combination of input values and the feedback state for spatial light modulation. Here, we use a similar mapping method as the "linear-combination mapping" in our previous paper [20]. The difference is with how the overall mapping area is split into several input portals. Instead of equally dividing the SLM into n sections for n input values, we repeat and stack the inputs until all the modulation areas on the SLM are filled. In the case where n = 4 , as shown in Fig.4, the 4 input values are repeated and

Figure 4: The mapping method to combine multi-feature inputs and the camera state for spatial light modulation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00652v1-figure-78a20c5b-3be0-4a89-a3f0-86f50eb32c92.jpg)

stacked sequentially, each with a height and length of 20 SLM pixels. The mapped input is multiplied by an input scale ( β ) and a random input weight, while the captured camera state is rescaled with a feedback gain ( α ). The combination is then linearly rescaled again to match with the 0 -2 π phase modulation on SLM. The nonlinearity of the system is achieved through the camera intensity reading and the intensity saturation rate. Similar to Eq.(1), our spatial RC can be simplified as

$$x ( t ) = f \{ [ \alpha x ( t - 1 ) + \beta W _ { i n } \odot W _ { c o r } \odot u ( t ) ] W _ { r e s } \}.$$

Here W cor is a designed weight for the features of each stock index benchmark. f {} and W res are the accumulated nonlinear effect and reservoir weight, respectively, through the whole system.

Thanks to the large number of pixels on SLM (400 × 400), it can easily function as a multi-input portal, allowing a large amount of information to be encoded into the system at once. In this way, the temporal correlation among the time-series data is translated into spatial coherence of the optical field, whose dynamics implement a complex-number kernel function. Therefore, besides evaluating the 1-step prediction, where the result contains only the next step value, we also assess the system's performance for multi-step ahead prediction using more inputs. In this case, u t ( ) is the vector containing all the feature information for different time steps. Specifically, we present examples of multi-step prediction with four-step and ten-step ahead predictions, labeled as 4-step and 10-step, respectively. Taking 4-step prediction as an example, the input vector u t ( ) has 32 values, corresponding to one closing

price and 7 features per time step. Accordingly, the prediction horizon is 4 steps containing 4 days of closing prices beyond the most recent input. This encoding method exploits the spatial correlation of the free space optics that occurs naturally in our setup.

## 3.2. Multiple ML and LSTM models

This paper uses LSTM and various regression methods to predict each benchmark [37, 39]. The LSTM architecture includes an internal memory that acts like a local variable, storing information as it processes a sequence of inputs. It uniquely handles outputs by recycling the production from one time-step as the input for the subsequent time-step[40]. This recycling process enables the LSTM to make informed predictions by considering both the current input and the output from the immediate preceding timestep. Besides LSTM, we also compared multiple regression models for stock index prediction. The regression methods used contain Lasso, ElasticNet, RandomForestRegressor, LinearRegression, Ridge, SGDRegressor, KNeighborsRegressor, DecisionTreeRegressor, and BayesianRidge[41, 42, 43].

For multi-step prediction, the data preparation for the time-series sequence with features follows the same approach as in optical RC in the previous subsection. The regression methods remain consistent with the one-step prediction, while the model dimensions in the LSTM are adjusted to align with the prediction horizon.

## 4. Results and analysis

Figure 5 (a) and (b) are two examples of 1-step prediction results for stock indices ˆ NYA and ˆ N225, respectively, showing the ground truth, the prediction results from RC with and without CW. Also shown are the prediction results of the best ML model, which is the Lasso regression model for both ˆ NYA and ˆ N225. The NRMSEs of RC prediction for ˆ NYA with and without CW are 0.104 and 0.108, respectively, whereas the best ML regression result is 0.558, which is over five times worse than the RC prediction. To compare the overall prediction performance with the results for ˆ N225 in Fig. 5(b), the prediction results for ˆ NYA show a small difference between predictions with and without CW. The difference for ˆ N225 between the two is much larger, as illustrated in Fig.5(b), where the NRMSE is 0.688 without

Figure 5: True values versus 1-step prediction results from optical RC with and without feature correlation weights, contrast to the best-performing ML results, shown for ˆ NYA (a) and ˆ N225 (b).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00652v1-figure-c1ec5b34-4a90-4be8-acc4-a7f93a0a4425.jpg)

CW, and only 0.236 with CW. This indicates that the added feature correlation weights can significantly improve the prediction performance for this benchmark.

In contrast, the predictive performance of ML models is considerably inferior to that of the RC, with the most accurate ML model achieving an NRMSE of 1.309 for ˆ N225. This is clear in Fig. 5, where both the best ML curves for ˆ NYA and ˆ N225 in Fig. 5 (a) and (b) lag behind the price movement trend and deviate significantly from the true values. Conversely, our optical RC predictions capture the stock volatility and nonlinear trends more accurately. These predictions align better with the real-time series stock prices than the best ML regression predictions.

To quantify their performance difference, Figure 6 lists the error reduction (ER) between the best ML model and optical RC predictions for seven stock indices. The ER is calculated using the formula

$$\text{ER} = 1 - \frac { \text{NRMSE} _ { R C } } { \text{NRMSE} _ { M L } },$$

where NRMSE RC and NRMSE ML are the NRMSE of the RC and the bestperforming ML methods, respectively. The above ER measures the performance improvement by RC compared to the best ML model. A near 100% ER indicates negligible relative error achieved by RC. As shown, the RC model significantly outperforms the ML model, with all ER values being positive and more than half of the results exceeding 70% for the different stock indices. Interestingly, for the stock indices ˆ N225, ˆ HSI, MME = F, and VTHR,

Figure 6: The error reduction between results from the best ML and RC models for different stock indices.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00652v1-figure-92bcc3be-b64e-469a-b0a6-3b86c6fad260.jpg)

the RC model shows superior performance in 1-step predictions compared to multi-step predictions. For other indices, the RC model achieves comparable or higher ER values in multi-step predictions than in 1-step predictions.

To gain insights into prediction performances between RC and ML models, we plot the prediction results of the two models for ˆ NYA in Fig. 7(a) and (b), respectively. Fig.7(a) shows the prediction outcomes from our RC model. The different-step prediction lines fit closely to the ground truth, with only a slight offset around a few peaks for 10-step prediction. However, the results of the ML model results in Figure 7(b) are more dispersed to compare the results for RC prediction. Additionally, the red rectangular area in Figure 7(b) highlights a 'delayed' effect in the multi-step prediction. Larger prediction dimensions result in greater drift from the ground truth. Consequently, as the forecasting horizon extends, prediction accuracy diminishes, leading to greater deviations from actual outcomes.

## 5. Conclusion

In conclusion, we have established an optical reservoir computing system in free space that effectively captures broader market behavior and demonstrates superior performance in stock index prediction. Compared to the

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00652v1-figure-65daf992-026c-446d-b471-803be6bd4aab.jpg)

2023-10

Date

Figure 7: 1-step and multi-step predictions for stock index ˆ NYA by (a) RC and (b) ML models

existing machine learning methods and long short-term memory neural networks, our system offloads intense computation to an optical system with low energy consumption, showing significant advantages in handling high volatility and nonlinear market behaviors using limited data. This method promises a cost-effective approach to real-time, parallel, multi-dimensional data processing and predictions, with broad applications in financial forecasting. The next steps include moving more data processing into the optical domain, for high-speed, high-quality, and versatile applications in financial analysis.

## References

- [1] A. W. Lo, What is an index?, Available at SSRN 2672755 (2015).
- [2] S. Shen, H. Jiang, T. Zhang, Stock market forecasting using machine learning algorithms, Department of Electrical Engineering, Stanford University, Stanford, CA (2012) 1-5.
- [3] C. K.-S. Leung, R. K. MacKinnon, Y. Wang, A machine learning approach for stock price prediction, in: Proceedings of the 18th international database engineering &amp; applications symposium, 2014, pp. 274277.
- [4] P. Soni, Y. Tewari, D. Krishnan, Machine learning approaches in stock price prediction: A systematic review, in: Journal of Physics: Conference Series, Vol. 2161, IOP Publishing, 2022, p. 012065.

| [5]   | B. M. Henrique, V. A. Sobreiro, H. Kimura, Literature review: Ma- chine learning techniques applied to financial market prediction, Expert Systems with Applications 124 (2019) 226-251.                                                        |
|-------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [6]   | W. Jiang, Applications of deep learning in stock market prediction: re- cent progress, Expert Systems with Applications 184 (2021) 115537.                                                                                                      |
| [7]   | M. Vijh, D. Chandola, V. A. Tikkiwal, A. Kumar, Stock closing price prediction using machine learning techniques, Procedia computer science 167 (2020) 599-606.                                                                                 |
| [8]   | R. Singh, S. Srivastava, Stock prediction using deep learning, Multime- dia Tools and Applications 76 (2017) 18569-18584.                                                                                                                       |
| [9]   | M. Obthong, N. Tantisantiwong, W. Jeamwatthanachai, G. Wills, A survey on machine learning for stock price prediction: Algorithms and techniques (2020).                                                                                        |
| [10]  | M. Nikou, G. Mansourfar, J. Bagherzadeh, Stock price prediction using deep learning algorithm and its comparison with machine learning al- gorithms, Intelligent Systems in Accounting, Finance and Management 26 (4) (2019) 164-174.           |
| [11]  | M. M. Kumbure, C. Lohrmann, P. Luukka, J. Porras, Machine learning techniques and data for stock market forecasting: A literature review, Expert Systems with Applications 197 (2022) 116659.                                                   |
| [12]  | D. M. Nelson, A. C. Pereira, R. A. De Oliveira, Stock market's price movement prediction with lstm neural networks, in: 2017 International joint conference on neural networks (IJCNN), Ieee, 2017, pp. 1419-1426.                              |
| [13]  | A. Moghar, M. Hamiche, Stock market prediction using lstm recurrent neural network, Procedia computer science 170 (2020) 1168-1173.                                                                                                             |
| [14]  | Y. Wang, Y. Liu, M. Wang, R. Liu, Lstm model optimization on stock price forecasting, in: 2018 17th international symposium on dis- tributed computing and applications for business engineering and science (dcabes), IEEE, 2018, pp. 173-177. |
| [15]  | S. Mehtab, J. Sen, A. Dutta, Stock price prediction using machine learning and lstm-based deep learning models, in: Machine Learning                                                                                                            |

|      | and Metaheuristics Algorithms, and Applications: Second Symposium, SoMMA 2020, Chennai, India, October 14-17, 2020, Revised Selected Papers 2, Springer, 2021, pp. 88-106.                                                                                                  |
|------|-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [16] | P. Ghosh, A. Neufeld, J. K. Sahoo, Forecasting directional movements of stock prices for intraday trading using lstm and random forests, Finance Research Letters 46 (2022) 102280.                                                                                         |
| [17] | M. A. I. Sunny, M. M. S. Maswood, A. G. Alharbi, Deep learning-based stock price prediction using lstm and bi-directional lstm model, in: 2020 2nd novel intelligent and leading emerging sciences conference (NILES), IEEE, 2020, pp. 87-92.                               |
| [18] | S. Selvin, R. Vinayakumar, E. Gopalakrishnan, V. K. Menon, K. Soman, Stock price prediction using lstm, rnn and cnn-sliding window model, in: 2017 international conference on advances in computing, communi- cations and informatics (icacci), IEEE, 2017, pp. 1643-1647. |
| [19] | R. Dash, P. K. Dash, Ahybrid stock trading framework integrating tech- nical analysis with machine learning techniques, The Journal of Finance and Data Science 2 (1) (2016) 42-57.                                                                                         |
| [20] | T. Bu, H. Zhang, S. Kumar, M. Jin, P. Kumar, Y. Huang, Efficient optical reservoir computing for parallel data processing, Optics Letters 47 (15) (2022) 3784-3787.                                                                                                         |
| [21] | P. Kumar, M. Jin, T. Bu, S. Kumar, Y.-P. Huang, Efficient reservoir computing using field programmable gate array and electro-optic mod- ulation, OSA Continuum 4 (3) (2021) 1086-1098.                                                                                     |
| [22] | M. Rafayelyan, J. Dong, Y. Tan, F. Krzakala, S. Gigan, Large-scale op- tical reservoir computing for spatiotemporal chaotic systems prediction, Physical Review X 10 (4) (2020) 041037.                                                                                     |
| [23] | B. Schrauwen, D. Verstraeten, J. Van Campenhout, An overview of reservoir computing: theory, applications and implementations, in: Pro- ceedings of the 15th european symposium on artificial neural networks. p. 471-482 2007, 2007, pp. 471-482.                          |

| [24]   | G. Tanaka, T. Yamane, J. B. Héroux, R. Nakane, N. Kanazawa, S. Takeda, H. Numata, D. Nakano, A. Hirose, Recent advances in physi- cal reservoir computing: A review, Neural Networks 115 (2019) 100-123.                |
|--------|-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [25]   | F. Duport, B. Schneider, A. Smerieri, M. Haelterman, S. Massar, All- optical reservoir computing, Optics express 20 (20) (2012) 22783-22795.                                                                            |
| [26]   | G. Van der Sande, D. Brunner, M. C. Soriano, Advances in photonic reservoir computing, Nanophotonics 6 (3) (2017) 561-576.                                                                                              |
| [27]   | J. Dong, M. Rafayelyan, F. Krzakala, S. Gigan, Optical reservoir com- puting using multiple light scattering for chaotic systems prediction, IEEE Journal of Selected Topics in Quantum Electronics 26 (1) (2019) 1-12. |
| [28]   | C.-H. Park, S. H. Irwin, What do we know about the profitability of technical analysis?, Journal of Economic surveys 21 (4) (2007) 786-826.                                                                             |
| [29]   | R. D. Edwards, J. Magee, W. C. Bassetti, Technical analysis of stock trends, CRC press, 2018.                                                                                                                           |
| [30]   | A. W. Lo, H. Mamaysky, J. Wang, Foundations of technical analysis: Computational algorithms, statistical inference, and empirical imple- mentation, The journal of finance 55 (4) (2000) 1705-1765.                     |
| [31]   | J. J. Murphy, Technical analysis of the financial markets: A comprehen- sive guide to trading methods and applications, Penguin, 1999.                                                                                  |
| [32]   | I. K. Nti, A. F. Adekoya, B. A. Weyori, A systematic review of fun- damental and technical analysis of stock market predictions, Artificial Intelligence Review 53 (4) (2020) 3007-3057.                                |
| [33]   | J. Dow, G. Gorton, Stock market efficiency and economic efficiency: Is there a connection?, The Journal of Finance 52 (3) (1997) 1087-1129.                                                                             |
| [34]   | N.-F. Chen, R. Roll, S. A. Ross, Economic forces and the stock market, Journal of business (1986) 383-403.                                                                                                              |
| [35]   | R. Morck, A. Shleifer, R. W. Vishny, M. Shapiro, J. M. Poterba, The stock market and investment: is the market a sideshow?, Brookings papers on economic Activity 1990 (2) (1990) 157-215.                              |

| [36]   | R. F. Engle, E. Ghysels, B. Sohn, Stock market volatility and macroeco- nomic fundamentals, Review of Economics and Statistics 95 (3) (2013) 776-797.                                                    |
|--------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [37]   | H. N. Bhandari, B. Rimal, N. R. Pokhrel, R. Rimal, K. R. Dahal, R. K. Khatri, Predicting stock market index using lstm, Machine Learning with Applications 9 (2022) 100320.                              |
| [38]   | S. Kumar, T. Bu, H. Zhang, I. Huang, Y. Huang, Robust and efficient single-pixel image classification with nonlinear optics, Optics Letters 46 (8) (2021) 1848-1851.                                     |
| [39]   | M. Sangiorgio, F. Dercole, Robustness of lstm neural networks for multi- step forecasting of chaotic time series, Chaos, Solitons & Fractals 139 (2020) 110045.                                          |
| [40]   | S. Hochreiter, J. Schmidhuber, Long short-term memory, Neural com- putation 9 (8) (1997) 1735-1780.                                                                                                      |
| [41]   | J. Li, W. Chen, Forecasting macroeconomic time series: Lasso-based approaches and their forecast combinations with dynamic factor models, International Journal of Forecasting 30 (4) (2014) 996-1015.   |
| [42]   | S. Basak, S. Kar, S. Saha, L. Khaidem, S. R. Dey, Predicting the di- rection of stock market prices using tree-based classifiers, The North American Journal of Economics and Finance 47 (2019) 552-567. |
| [43]   | B. Qian, K. Rasheed, Stock market prediction with multiple classifiers, Applied Intelligence 26 (2007) 25-33.                                                                                            |
Fang Wang Ting Bu Yuping Huang
Enhancing Multistep Prediction of Multivariate Market Indices Using Weighted Optical Reservoir Computing
We propose and experimentally demonstrate an innovative stock index
prediction method using a weighted optical reservoir computing system. We
construct fundamental market data combined with macroeconomic data and
technical indicators to capture the broader behavior of the stock market. Our
approach shows significant higher performance than state-of-the-art methods
such as linear regression, decision trees, and neural network architectures
including long short-term memory. It captures well the market's high volatility
and nonlinear behaviors despite limited data, demonstrating great potential for
real-time, parallel, multi-dimensional data processing and predictions.