Available online at www.HighTechJournal.org HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 743 ISSN: 2723-9535 SARIMAX–GARCH Model to Forecast Composite Index with Inflation Rate and Exchange Rate Factors M. Fariz Fadillah Mardianto 1* , Elly Pusporani 1 , Diana Ulya 1, I Kadek Pasek Kusuma Adi Putra 1 , Rico Ramadhan 1 1 Faculty of Science and Technology, Airlangga University, Surabaya 60115, Indonesia. Received 15 April 2024; Revised 10 August 2024; Accepted 16 August 2024; Published 01 September 2024 Abstract Investors should consider the Indonesia Composite Index (ICI) as a key indicator before making investment decisions, as it reflects the performance of industries and the broader economic growth. In Indonesia, the ICI exhibits fluctuating movements, making accurate forecasting essential for understanding the country's economic conditions, which are closely tied to capital flows, growth, and tax revenues. This study aims to forecast the ICI using the SARIMAX-GARCH model, incorporating macroeconomic factors such as the inflation rate and exchange rate. The findings reveal that both variables significantly impact the ICI, with the model achieving a Mean Absolute Percentage Error (MAPE) of 0.952% for training data and 5.233% for test data. The model's performance is supported by an RΒ² value of 0.9782 and a Mean Squared Error (MSE) of 0.0003. This research not only improves the accuracy of ICI forecasts but also supports Indonesia's 8th Sustainable Development Goal (SDG) for decent work and economic growth. Keywords: Indonesia Composite Index; SARIMAX-GARCH; Inflation Rate;, Exchange Rate; Sustainable Development Goals. 1. Introduction 1.1. Indonesia Composite Index The Indonesia Composite Index (ICI) is a statistical measurement used to determine changes in the share prices of all companies listed in the capital market at a certain time compared to the base year [1]. The development of the ICI not only reflects the performance of a country's companies or industries but can also be considered a broader fundamental indicator of national economic health [2]. The movement of the ICI in Indonesia tends to fluctuate. In 2022, the ICI closed down 0.14% to 6,850.61 in the last trade. However, on a year-to-date (YTD) basis, the index increased by 4.09% throughout 2022, though this growth was not as high as in the previous period [3]. The volatile movement of the ICI is closely related to macroeconomic conditions. Macroeconomic factors are elements outside the company that can affect performance, either directly or indirectly. Inflation and exchange rates are macroeconomic variables that are known to affect the capital market [4]. 1.2. Inflation Rate Inflation is generally defined as a continuous increase in prices that applies broadly or impacts other goods and services [5]. As a key macroeconomic indicator, inflation significantly influences various aspects of the economy, * Corresponding author: m.fariz.fadillah.m@fst.unair.ac.id http://dx.doi.org/10.28991/HIJ-2024-05-03-014 ➒ This is an open access article under the CC-BY license (https://creativecommons.org/licenses/by/4.0/). Β© Authors retain all copyrights. https://creativecommons.org/licenses/by/4.0/ https://orcid.org/0000-0002-4541-4552 https://orcid.org/0000-0001-7948-3719 https://orcid.org/0009-0008-7856-7302 https://orcid.org/0000-0002-2565-7944 HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 744 including the stock market. Understanding the relationship between inflation and the ICI on the Indonesia Stock Exchange (IDX) is crucial for assessing its implications for the wider economy [6]. Studies such as Amalia (2016) show that inflation has a significant positive effect on the composite stock price index [7]. However, this finding contrasts with research conducted by Yudha (2016), which suggests that inflation has a significant negative effect on the composite stock price index [8]. Typically, a stable or low inflation rate reflects a healthy economy, boosting investor confidence. Conversely, high inflation can raise concerns about the government's ability to maintain economic stability, potentially reducing the ICI [9]. 1.3. Exchange Rate The exchange rate, specifically the comparison of the value of the rupiah with foreign currencies, is a crucial factor in determining the cost of goods and services across borders [10]. Research by Tampubolon (2021) indicates that the exchange rate has a significant negative effect on the ICI [11]. This study explains that if the rupiah depreciates, it can increase operational costs, thereby reducing company profits and affecting the ICI negatively. In contrast, Hasanudin (2021) demonstrates that a stronger rupiah can lead to better ICI performance [12]. The exchange rate used is the middle rate of the rupiah against the US dollar. A weakening dollar against the rupiah, coupled with less promising alternative investments, might lead investors to prefer holding dollars, further influencing the ICI [13]. 1.4. Vision Statement The urgency of this research lies in the need to improve the accuracy of ICI forecasts, an essential indicator of Indonesia's economic health. The SARIMAX-GARCH method, which extends the basic ARIMA model by including exogenous variables such as inflation and exchange rates while accounting for seasonal effects, offers a significant advancement. The SARIMAX model provides the ability to predict time-series observations and incorporate the influence of exogenous variables on the response variable. However, when applied to economic data, the SARIMAX model can encounter issues with non-constant variance. The GARCH model addresses this by managing variance non- uniformity in time series data, resulting in a robust prediction model [14]. Given the global economic uncertainties, such as those caused by the COVID-19 pandemic, this research is particularly urgent as it aims to provide more accurate forecasts of the ICI. This enhanced predictive capability supports better investment decisions and policymaking, directly contributing to the achievement of Indonesia's 8th Sustainable Development Goal (SDG), which targets decent work and economic growth. 2. Research Methodology 2.1. Autoregressive Integrated Moving Average (ARIMA) The Autoregressive Integrated Moving Average (ARIMA) method, also known as the Box-Jenkins method, was introduced by George Box and Gwilym Jenkins in 1970 [15]. ARIMA represents an Autoregressive Moving Average (ARMA) model that lacks stationarity. While the ARIMA method excels in short-term forecasting, its accuracy diminishes for long-term projections, often resulting in forecasted values remaining relatively constant over extended periods. The ARIMA model comprises three processes: autoregressive, integration, and moving average, represented by the order autoregre. An essential assumption of ARIMA is the requirement of stationarity in variation. To address non- stationarity in the data, a differencing process is employed to make the data stationary, with the number of differencing steps denoted as (𝑑). The fundamental principle of time series posits that the current observation (𝑧𝑑). is influenced by one or more preceding observations (𝑧𝑑). The general expression for the autoregressive integrated moving average model is denoted as ARIMA(𝑝, 𝑑, π‘ž) [16]. Ο†p(B)(1 βˆ’ B)dZt = ΞΈq(B)Ξ΅t (1) where, 𝐡 is denoted as backshift operation and πœ‘π‘(𝐡) is denoted as autoregressive (AR) backshift that follows second equation below. πœ‘π‘(𝐡) = (1 βˆ’ πœ‘1𝐡 βˆ’ β‹― βˆ’ πœ‘π‘π΅π‘) (2) while, πœƒπ‘ž(𝐡) is denoted as moving average (MA) backshift operator as follows. πœƒπ‘ž(𝐡) = (1 βˆ’ πœƒ1𝐡 βˆ’ β‹― βˆ’ πœƒπ‘žπ΅π‘ž) (3) with (1 βˆ’ 𝐡)𝑑𝑍𝑑 as time series that is stationary at the 𝑑-th differencing. This process is denoted with ARIMA (𝑝, 𝑑, π‘ž). ARIMA is usually estimated with maximum likelihood (MLE) method or conditional sum squares (CSS) method to find the estimation of all parameters. However, the selection of the number of parameters 𝑝 and π‘ž in the HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 745 ARIMA model is based on the patterns observed in the autocorrelation function (ACF) and the partial autocorrelation function (PACF) as referred to Table 1 [17]. Table 1. Determine 𝒑 and 𝒒 in ARIMA Model Model ACF PACF ARIMA (𝒑, 𝒅, 𝟎) Dies down Drop off after lag π‘ž ARIMA (𝟎, 𝒅, 𝒒) Drop off after lag π‘ž Dies down ARIMA (𝒑, 𝒅, 𝒒) Dies down (until lag π‘ž is still different from zero) Dies down (until lag 𝑝 is still different from zero) 2.2. Seasonal Autoregressive Integrated Moving Average (SARIMA) A time series may occasionally exhibit seasonal phenomena that recur at specific intervals. The shortest time interval for these recurring phenomena is known as the seasonal period [18]. Seasonal periods commonly used include 1 month, 3 months, 4 months, 6 months, and 12 months or annually. In general, a seasonal ARIMA model is expressed as ARIMA(𝑝 𝑑 π‘ž)(𝑃 𝐷 𝑄)𝑆, where 𝑑 is the nonseasonal differencing order, 𝐷 is the seasonal differencing order, and 𝑆 is the seasonal period. The general form of the seasonal ARIMA has been stated in Equation 4 as follows. βˆ…π‘(𝐡) Φ𝑃(𝐡𝑆)(1 βˆ’ 𝐡)𝑑(1 βˆ’ 𝐡𝑆)𝐷𝑍�̇� = πœƒπ‘ž(𝐡) Ξ˜π‘„(𝐡𝑆) π‘Žπ‘‘ (4) where, βˆ…π‘(𝐡) ∢ (1 βˆ’ βˆ…1𝐡 βˆ’ βˆ…2𝐡2 βˆ’ β‹― βˆ’ βˆ…p𝐡𝑝) polynomial of non seasonal AR(𝑝) πœƒπ‘ž(𝐡) ∢ (1 βˆ’ πœƒ1𝐡 βˆ’ πœƒ2𝐡2 βˆ’ β‹― βˆ’ πœƒqπ΅π‘ž) polynomial of non seasonal MA(π‘ž) Φ𝑃(𝐡𝑆) ∢ (1 βˆ’ Ξ¦1𝐡𝑆 βˆ’ Ξ¦2𝐡2𝑆 … βˆ’ Ξ¦P𝐡𝑃𝑆) polynomial of seasonal AR(𝑃) Ξ˜π‘„(𝐡𝑆) ∢ (1 βˆ’ Θ1𝐡𝑆 βˆ’ Θ2𝐡2𝑆 … βˆ’ ΘQ𝐡𝑄𝑆) polynomial of seasonal MA(𝑄) 2.3. Autoregressive Integrated Moving Average with Exogenous Variable (SARIMAX) The ARIMAX model, an extension of the ARIMA model, incorporates additional or exogenous variables that are considered to have a significant impact on the data, thereby enhancing the accuracy of forecasting [19]. In the field of forecasting, there are other variables that are postulated to influence the model. The existence of these influential variables can lead to significant fluctuations in observation values, which recur over distinct time periods. As a result, a specialized model is necessary for forecasting under these conditions. The first step in ARIMAX modeling involves testing the stationarity of the exogenous variables [20]. According to Simms et al. (2022) [21], the general form of the ARIMAX (𝒑, 𝒅, 𝒒) model can be expressed with the following equation. 𝑍𝑑 = 𝛽1𝑋1,𝑑 + 𝛽2𝑋2,𝑑 + β‹― + 𝛽𝑝𝑋𝑝,𝑑 + πœƒπ‘ž(𝐡) πœ‘π‘(𝐡)(1 βˆ’ 𝐡)𝑑 πœ€π‘‘ (5) The exogenous variables 𝑋𝑖,𝑑 (𝑖 = 1,2, … , 𝑑) in Equation 4 are presented under the condition of stationarity with each variable have their unknown parameter 𝛽. Extending this further, the ARIMAX model equation, when considering a stochastic trend and seasonality, is referred to as the SARIMAX model. The SARIMAX model represents an enhancement of the ARIMAX model, as it incorporates considerations for seasonal factors. The typical expression for the SARIMAX model is generally articulated as follows [22]. 𝑍𝑑 = 𝛽1𝑋1,𝑑 + 𝛽2𝑋2,𝑑 + β‹― + 𝛽𝑝𝑋𝑝,𝑑 + πœƒπ‘ž(𝐡)Ξ˜π‘„(𝐡𝑆) πœ‘π‘(𝐡)Φ𝑃(𝐡𝑆)(1 βˆ’ 𝐡)𝑑(1 βˆ’ 𝐡𝑆)𝐷 πœ€π‘‘ (6) To determine the appropriateness of the obtained model, a diagnostic check is necessary. It is possible that the results of time series modeling may yield several models with all significant parameters. In such cases, the residuals should meet the white noise assumption and be normally distributed [23]. 2.4. Generalized Autoregressive Conditional Heterocedasticity (GARCH) The term β€˜volatility’ is commonly used to describe the β€˜volatile’ behavior observed in financial markets. Volatility has become significant in both financial theory and practice, playing a pivotal role in areas such as risk management and portfolio selection. In statistical analyses, volatility is typically quantified using variance or standard deviation. In 1982, Engle successfully introduced a volatility model for financial time series data, known as the Autoregressive Conditional Heteroscedasticity (ARCH) model. Furthermore, in 1986, Bollerslev developed a more adaptable volatility model, termed the Generalized Autoregressive Conditional Heteroscedasticity (GARCH) [24]. When the ARCH effect test indicates significance for a time series, it becomes feasible to estimate the ARCH model and concurrently derive an estimate of volatility, denoted as πœŽπ‘‘, based on historical information. In practical applications, HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 746 the lag count, p, is often substantial, leading to a considerable number of parameters being estimated in the model. In 2002, Bollerslev, Zivot, and Wang introduced a more concise model, replacing the AR model with the subsequent formulation [25]. πœŽπ‘‘ 2 = 𝛼0 + βˆ‘ π›Όπ‘–πœŽπ‘‘βˆ’π‘– 2 + 𝑝 𝑖=1 βˆ‘ π›½π‘—πœŽπ‘—βˆ’π‘‘ 2 π‘ž 𝑗=1 (7) Ensuring that all coefficients 𝛼𝑖 > 0 (π‘€β„Žπ‘’π‘Ÿπ‘’ 𝑖 = 0,1, … , 𝑝) and 𝛽𝑗 > 0 (π‘€β„Žπ‘’π‘Ÿπ‘’ 𝑗 = 1, … , 𝑝) are positive is necessary to maintain the positivity of the conditional variance πœŽπ‘‘ 2. The equation presented above, coupled with the stationary time series equation (π‘Ÿπ‘‘), is recognized as the Generalized Autoregressive Conditional Heteroscedasticity or GARCH(𝑝, π‘ž) model. In the given equation, when 𝛽𝑗 > 0 (π‘€β„Žπ‘’π‘Ÿπ‘’ 𝑗 = 1, … , 𝑝), the GARCH model transforms into an Autoregressive Conditional Heteroscedasticity or ARCH model. 2.5. Goodness of Fit Model In modeling to generate predictions, goodness-of-fit measures are used to assess how accurately a model can estimate true values. The following are several techniques that can be employed to measure the goodness-of-fit of a model. 1) Mean Absolute Percentage Error (MAPE) The Mean Absolute Percentage Error (MAPE) is a statistical measure employed in the evaluation of the accuracy of predictions or forecasts within a model. MAPE quantifies the average error in percentage terms between actual and predicted values [26], providing insight into the extent to which forecast errors compare to the actual data values. The statistical formula for MAPE is presented in Equation 8: MAPE = βˆ‘ | 𝑍𝑑 βˆ’ Ẑ𝑑 Zt | Γ— 100% 𝑛 t=1 (8) 𝑍𝑑 indicates the actual value and οΏ½Μ‚οΏ½ 𝑑 indicates the predicted value. The smaller the percentage error value in MAPE, the higher the accuracy level of the forecasting results. Table 2 provides an interpretation of MAPE values [27]. Table 2. Interpretation of MAPE Values MAPE Value Ranges Interpretation < 10% The ranges indicate that the model has a high degree of forecasting accuracy. 10% ≀ 𝑀𝐴𝑃𝐸 ≀ 20% The ranges indicate that the model's forecasting accuracy is good. 20% < 𝑀𝐴𝑃𝐸 ≀ 50% The ranges indicate that the model has sufficient predictive power. > 50% The ranges indicate that the model's predictive power is at its lowest. 2) Mean Squared Error (MSE) The Root Mean Square Error (RMSE) serves as an error metric utilized to evaluate forecasting methods and measure the precision of a model’s predictive outcomes. The Mean Squared Error (MSE), which is the average of squared errors, quantifies the extent of discrepancies between the model-predicted values and the actual values [28]. Consequently, a lower RMSE value signifies a higher level of accuracy in the model’s predictions of actual values. The computation formula for RMSE is delineated in Equation 9 (𝑛 indicates sum of data). MSE = βˆ‘ (𝑍𝑑 βˆ’ �̂�𝑑) 2𝑛 t=1 𝑛 (9) 3) R Squared (R2) – Coefficients of Determination According to Ghozali [29], a low coefficient of determination indicates a limited capacity of independent variables to explain the dependent variable. Conversely, a value approaching 1, distinct from 0, suggests that the independent variables have the potential to provide all the necessary information for predicting the dependent variable. The formula for calculating this coefficient is outlined below. R2 = 1 βˆ’ βˆ‘ (𝑍𝑑 βˆ’ �̂�𝑑) 2𝑛 𝑑=1 βˆ‘ (𝑍𝑑 βˆ’ οΏ½Μ…οΏ½)2𝑛 𝑑=1 (10) HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 747 The larger the value of the coefficient of determination, the better the estimation results for the model. However, a coefficient value very close to 1 can lead to overfitting in the estimation results. Therefore, additional evaluations such as MSE and MAPE, as explained earlier, are needed. 2.6. Data Sources and Research Variables This research employs a quantitative approach, focusing on the analysis of time series data using SARIMAX- GARCH methods. The data for the Indonesia Composite Index (ICI) was sourced from the Yahoo Finance website, while the data for Inflation and Exchange Rate was obtained from the Bank Indonesia website. The study utilizes monthly data spanning from January 2015 to March 2023. The research data is bifurcated into two segments: training data and testing data. The training data, which is used to construct the model, comprises data from January 2015 to September 2022. Conversely, the testing data, used to gauge the accuracy of the model, includes data from October 2022 to March 2023. In this context, the 𝑑 index on each variable signifies time, and the significance level of Ξ± is set at 0.1. The variables employed in this research are delineated in Table 3. As an additional note, the analysis will be conducted in RStudio for model analysis and Minitab for visualization. Table 3. Research Variables Variable Description 𝒁𝒕 Indonesia Composite Index π‘ΏπŸπ’• Inflation Rate π‘ΏπŸπ’• Exchange Rate (Dollars to Rupiah) 2.7. Data Analysis Stages The procedure or stages of the analysis method in this study are systematically presented in the flow chart of Figure 1 as follows. Figure 1. Data Analysis Stages Chart 3. Results and Discussion 3.1. Time Series Plots and Descriptive Statistics Descriptive statistics serve as the initial observation to examine the characteristics of each variable, namely the Indonesia Composite Index (ICI), inflation rate, and the dollar-to-rupiah exchange rate. Figure 2 presents the time series plots of the ICI, inflation rate, and the Indonesian dollar-to-rupiah exchange rate on a monthly basis from January 2015 to September 2022. Determination of research variables Checking data stationarity with the ADF Test Determination of the best ARIMA model Determination of the best SARIMAX model Checking the heteroscedasticity of the model GARCH modeling for model residuals Heteroscedasticity check Model assumption checking (normality, multicollinearity, and white noise) ICI prediction for the period September 2022-March 2023 Checking the goodness of prediction HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 748 (a) (b) (c) Figure 2. (a) Time Series of ICI, (b) Time Series of Inflation, (c) Time Series of Exchange Rate The observation plot reveals fluctuations in the ICI, inflation rate, and dollar-to-rupiah exchange rate. The plot of the dollar-to-rupiah exchange rate exhibits an upward trend, beginning in January 2015 at Rp 12,688.00 (also the minimum value) and ending at Rp 15,323.00. Similarly, the ICI plot demonstrates an upward trend amidst the fluctuating curve. The minimum ICI was 4223.91 in September 2015, and the maximum value of ICI was 7228.91, which occurred in April 2022. However, the year-on-year inflation rate in Indonesia displayed a downward trend, starting at 6.96% and ending at 5.95% in September 2022. Table 4 provides a summary of the data in the form of descriptive statistics for each variable. Table 4. Summary Results of Descriptive Statistics Variables Average Variance Standard Deviation Median Range Minimum Maximum ICI 5777.7 530225.4 728.2 5936.4 3005 4223.9 7228.9 Inflation Rate 3402 2.336 1.529 3.23 5.94 1.32 7.26 Exchange Rate 14077 421791 649 14154 3761 12688 16449 Descriptive statistics yield seven indicator values derived from measures of data centering and deviation. The standard deviation of the ICI is 728.2, indicating that, on average, the ICI deviates from its parameters by 728.2. The inflation rate has a unit value compared to the other two variables, the ICI and the dollar-to-rupiah exchange rate, which are in the thousands. This suggests an imbalance in the variance of all variables. Additionally, the selling exchange rate has a high standard deviation of 649. Based on the descriptive statistical analysis, a logarithmic transformation is performed to minimize the range of variation of the three variables [30]. As a result, variables with a large variance are drastically reduced. The results of the descriptive statistics can be seen in Table 5. Table 5. Summary Results of Descriptive Statistics (Logarithm Transformation) Variables Average Variance Standard Deviation Median Range Minimum Maximum ICI 8.65 0.0166 0.1287 8.35 0.5373 8.35 8.89 Inflation Rate 1.13 0.2001 0.4473 0.28 1.7047 0.28 1.98 Exchange Rate 9.55 0.0021 0.0459 9.45 0.2596 9.45 9.71 HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 749 3.2. Data Stationarity Stationarity is a prerequisite in classical time series modeling, such as ARIMA, ARIMAX, or SARIMAX. Stationarity in data implies stability over time, with a constant mean and covariance. Graphically, data with fluctuations typically does not yield stationary data, as it often lacks a constant mean or covariance. The Augmented Dickey-Fuller (ADF) Test is the statistical test used to assess the stationarity of time series data [31]. The null hypothesis for the test posits that the data is not stationary, while the alternative hypothesis asserts stationarity. Stationarity is achieved if the ADF test yields a probability value (p-value) less than the specified significance level (𝜢 = 0.1). Table 6 presents the ADF test results for all the variables. Regrettably, the initial trial did not satisfy the stationarity condition. To address the issue of data non-stationarity, a transformation in the form of differencing the data is required. As an experiment, differencing was performed with a lag of 1 on each data point. Subsequently, the transformed data was retested with the ADF test. The results of this test are displayed in Table 6. Table 6. ADF Test Results Transformation Variables P-value Decision Not transformed with first differencing Indonesia Composite Index 0.4573 Data is not stationary Inflation Rate 0.9683 Data is not stationary Exchange Rate 0.1157 Data is not stationary Transformed with first differencing Indonesia Composite Index 0.000 Data is stationary Inflation Rate 0.041 Data is stationary Exchange Rate 0.000 Data is stationary The ADF test results indicate that the ICI, inflation rate, and dollar-to-rupiah exchange rate are stationary following logarithmic transformation and first differencing. Consequently, modeling with ARIMAX can proceed, as the data has fulfilled the stationarity assumption. 3.3. Autocorrelation Function (ACF) and Partial Autocorrelation Function (PACF) ACF and PACF are two primary indicators that describe how the data is correlated with time. ACF is used to measure the correlation between the previous observation and the current observation in the time series data (for MA order determination), while PACF measures the direct correlation between two observations in the time series after accounting for the influence between them (for AR order determination) [32]. Figure 3 illustrates the ACF and PACF plots of the ICI data that has been transformed to stationarity by differencing. (a) (b) Figure 3. (a) ACF Plot for Differencing 1, (b) PACF Plot for Differencing 1 The results reveal minimal lags in both the ACF and PACF. However, lag-9 on the ACF plot slightly crosses the red line of significance. This suggests that the ICI data is likely correlated with its own data every 9 periods. Additionally, the significant lag that is quite distant is indicative of the possibility that the ICI exhibits a seasonal trend every 9 periods or every 9 months. Figure 4 presents the ACF and PACF plot of the ICI data that has been differenced for 9 periods. HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 750 (a) (b) Figure 4. (a) ACF Plot for Differencing 9, (b) PACF Plot for Differencing 9 Based on the results of the ACF and PACF plots with a seasonal trend every 9 periods, it is observed that the ACF plot increasingly shows its significance at lag-9 and several lags emerge in the PACF. Therefore, based on these plots, the potential SARIMAX models to be used are SARIMAX(0,1,0)(1,1,0)9, SARIMAX(0,1,0)(0,1,1)9, SARIMAX(0,1,0) (1,1,1)9, SARIMAX(0,1,0)(2,1,0)9, and SARIMAX(0,1,0)(2,1,1)9. 3.4. SARIMAX Model Estimation The ACF and PACF plots have suggested several potential SARIMAX models for estimation. In SARIMAX modeling, it is crucial that exogenous variables, which have achieved stationarity, do not exhibit multicollinearity [33]. Accordingly, Table 7 presents the results of the multicollinearity test, conducted using the Variance Inflation Factor (VIF) method, for each exogenous variable [33]. Table 7. Multicollinearity Test for Exogenous Variables Variables VIF Conclusion Difference 1 Inflation 1.00 There is no multicollinearity between the inflation rate and the exchange rate. Difference 1 Exchange Rate 1.00 The absence of multicollinearity suggests that the exogenous variables satisfy the assumptions and should be incorporated into the SARIMA model, thereby transforming it into a SARIMAX model. Table 8 provides details about the five potential SARIMAX models discussed earlier in the β€œACF and PACF Plots” section. In addition to meeting the assumptions of the SARIMAX model, such as residual normality and white noise residual, the Shapiro-Wilk test will be used for normality testing, and the Ljung-Box Test will be employed for white noise testing [34]. Table 8. Summary of ICI Model Estimation with SARIMAX SARIMAX Model Order P-value Residual Normality Test P-value Residual White Noise Test Parameter Significance (0,1,0)(1,1,0)9 0.1286 (Normal) 0.1104 (White Noise) All Parameters Significant (0,1,0)(0,1,1)9 0.0051 (Not Normal) 0.1602 (White Noise) All Parameters Significant (0,1,0)(1,1,1)9 0.0013 (Not Normal) 0.1691 (White Noise) All Parameters Significant (0,1,0)(2,1,0)9 0.0128 (Not Normal) 0.0071 (Not White Noise) All Parameters Significant (0,1,0)(2,1,1)9 0.1502 (Normal) 0.0027 (Not White Noise) SAR-18 Not Significant Based on Table 8, the SARIMAX model that fulfills all assumptions is SARIMAX(0,1,0)(1,1,0)9. Table 9 provides a more complete estimate of the SARIMAX(0,1,0)(1,1,0)9 model with its estimated parameters. HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 751 Table 9. Summary of ICI Model Estimation with SARIMAX Parameter Estimated Coefficient Standard Deviation of Coefficient P-value Seasonal AR (οΏ½Μ‚οΏ½πŸ—) -0.556 0.0893 0.000 Difference 1 Inflation Rate (οΏ½Μ‚οΏ½πŸ) -0.0089 0.0305 0.000 Difference 1 Exchange Rate (οΏ½Μ‚οΏ½πŸ) -0.4053 0.1083 0.000 Thus, the estimated SARIMAX(0,1,0)(1,1,0)9 model will be used in the SARIMAX-GARCH estimation. The mathematical equation of the logarithmically transformed ICI model estimation with SARIMAX(0,1,0)(1,1,0)9 is as follows: 𝑍�̂� βˆ— = βˆ’0.0089𝑋1𝑑 βˆ— βˆ’ 0.4053𝑋2𝑑 βˆ— + πœ€οΏ½Μ‚οΏ½ [(1 + 0.556 𝐡9)(1 βˆ’ 𝐡9)(1 βˆ’ 𝐡)] (11) With, 𝑍�̂� βˆ— = 𝑙𝑛(�̂�𝑑), while 𝑋1𝑑 βˆ— is the first differentiation of the inflation rate and 𝑋2𝑑 βˆ— is the first differentiation of the exchange rate. Another interpretation is that all of the exogenous variables are exhibiting a significant negative effect on the ICI over time. The negative coefficient in the first differencing of the inflation rate suggests that a larger shock in uncertain inflation rates leads to a decrease in the performance of the ICI. Regrettably, this also applies to the exchange rate differences. The coefficient of the first difference exchange rate, which is larger in a negative sense compared to the inflation rate, indicates a greater decrease in the ICI due to the inconsistency of changes in the exchange rate. Therefore, both inflation rates and exchange rates have a significant impact and play important roles in the performance of the Indonesia Composite Index. Based on the equation obtained, Figure 5 shows the comparison graph between the actual ICI’s value with the ICI’s estimated value with the SARIMAX(0,1,0)(1,1,0)9 model. The red line shows the estimation of ICI with SARIMAX model and the blue line represents the actual ICI’s values. Figure 5. Indonesia Composite Index with SARIMAX Estimation The plot of the estimation with the SARIMAX(0,1,0)(1,1,0)9 model is approaching the pattern with very precise to the original data. To look how good is the estimation, Table 10 provides criteria for the goodness of the SARIMAX model in the form of the R2 value, MSE, and the value of MAPE in sample. Table 10. Summary of Goodnes of Fit from SARIMAX Criteria Value Description π‘πŸ 0.8785 Scores in the good category MSE 0.0021 Small MSE value MAPE in sample 3.265% MAPE is in the excellent category HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 752 3.5. Heteroscedasticity Detection for SARIMAX Residuals In the SARIMAX model, it is understood that the error follows a normal distribution with a mean of 0 and a homogeneous variance 𝝈𝟐. However, economic data generally exhibits high volatility, leading to a deviation from the assumption of homogeneous variance [35]. To detect heteroscedasticity, the squared error of the SARIMAX(0,1,0) (1,1,0)9 estimation is employed as a representative estimator of the residuals 𝝈𝟐 of the residuals [36]. So, in Figure 6 it will shows the ACF and PACF of the squared residual from SARIMAX(0,1,0)(1,1,0)9. (a) (b) Figure 6. (a) ACF Plot for SARIMAX Squared Residuals, (b) PACF Plot for SARIMAX Squared Residuals It can be observed that there is a lag out, indicating that the squared error has a correlation with its own data and previous residuals. This results in the variance of the error not being homogeneous. Therefore, the integration of the GARCH model into the SARIMAX model will address the issue of heteroscedasticity. 3.6. GARCH Estimation Autoregressive Conditionally Heteroscedasticity (ARCH) is a model designed for the variance residuals of time series models that exhibit cases of heteroscedasticity. The ARCH model emphasizes that the variance of the error depends on the squared residual (variance estimator) of the preceding period. However, ARCH does not detect heteroscedasticity based solely on variability. The Generalized-ARCH (GARCH) model accommodates such variability. Thus, in addition to calculating the squared residuals from the previous period, GARCH enables the estimation of past variability. In the preceding discussion, it was identified that the SARIMAX(0,1,0)(1,1,0)9 estimation exhibits heteroscedastic residuals. Consequently, ARCH/GARCH variance modeling is necessitated. Based on the ACF and PACF plots of squared residuals, the number of lags that emerge on the ACF determines the ARCH(𝑝), and the number of lags that emerge on the PACF determines the GARCH(𝑝, π‘ž) model [37]. Therefore, the potential variance model estimates, based on Figure 6, are ARCH(1) (equivalent to GARCH(1,0)) or GARCH(1,1). Table 11 presents the estimation results of both models using the residual data from SARIMAX(0,1,0)(1,1,0)9. Table 11. Summary of GARCH Model Estimation Model Parameter Significance GARCH(1,0) All Parameters Significant GARCH(1,1) Insignificant Variability Parameter In Table 11, GARCH(1,0) has model simplicity and has parameters that are all significant to the model. Therefore, to overcome the case of heteroscedasticity, GARCH(1,0) will become addition in SARIMAX(0,1,0)(1,1,0)9. Equation 12 shows the equation of residual variance as follows. πœŽοΏ½Μ‚οΏ½ 2 = 0.0011983 + 0.4378342πœ€οΏ½Μ‚οΏ½βˆ’1 2 (12) with πœŽοΏ½Μ‚οΏ½ 2 is the estimated variance to solve the heteroscedasticity and πœ€οΏ½Μ‚οΏ½βˆ’1 2 is the estimated square residual 1 lag previous period from SARIMAX(0,1,0)(1,1,0)9. Therefore, Equation 13 shows how to have the estimated residuals based on GARCH(1,0) model [38]. πœ€οΏ½Μ‚οΏ½ = �̂�𝑑�̂�𝑑 (13) With, �̂�𝑑 is the standardization of the estimated SARIMAX(0,1,0)(1,1,0)9 error and is assumed to be normally distributed with mean 0 and variance 1 and �̂�𝑑 is the square root of estimated variance that was shown in Equation 12. HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 753 3.7. SARIMAX-GARCH Estimation The combination of Equations 11 and 13 will be formed into SARIMAX(0,1,0)(1,1,0)9-GARCH(1,0) which can be mathematically written as follows: 𝑍�̂� βˆ— = βˆ’0.0089𝑋1𝑑 βˆ— βˆ’ 0.4053𝑋2𝑑 βˆ— + �̂�𝑑�̂�𝑑 [(1 + 0.556 𝐡9)(1 βˆ’ 𝐡9)(1 βˆ’ 𝐡)] (14) By changing the �̂�𝒕 with Equation 12 and returning 𝒁𝒕 to its original form (inverse of logarithmic transformation), it will have the final SARIMAX(0,1,0)(1,1,0)9-GARCH(1,0) estimation as follows: �̂�𝑑 = exp [βˆ’0.0089𝑋1𝑑 βˆ— βˆ’ 0.4053𝑋2𝑑 βˆ— + οΏ½Μ‚οΏ½π‘‘βˆš0.0011983 + 0.4378342πœ€οΏ½Μ‚οΏ½βˆ’1 2 [(1 + 0.556 𝐡9)(1 βˆ’ 𝐡9)(1 βˆ’ 𝐡)] ] (15) The results of this equation are remarkable, rendering the estimation of the Indonesia Composite Index highly similar to the original value, as depicted in Figure 7. Figure 7. Indonesia Composite Index with SARIMAX-GARCH Estimation Based on the plot with SARIMAX(0,1,0)(1,1,0)9 - GARCH(1,0) in blue, it appears to closely match the black graph (actual ICI value) and provides a better fit than the red line (SARIMAX(0,1,0)(1,1,0)9 model without the GARCH model). This is because the residual variance does not fluctuate, resulting in a smaller and more precise residual value for estimation. Table 12 presents the goodness of fit from the SARIMAX-GARCH model. Table 12. Summary of Goodness of Fit from SARIMAX-GARCH Criteria Value Description π‘πŸ 0.9782 Scores in the excellent category MSE 0.0003 MSE Value Very Small MAPE in sample 0.952% MAPE is in the excellent category The results given by SARIMAX(0,1,0)(1,1,0)9-GARCH(1,0) compared to only SARIMAX(0,1,0)(1,1,0)9 in Table 10 have better criteria, such as 𝑅2 higher, smaller MSE, and very small MAPE. The purpose of incorporating the GARCH model into SARIMAX is to eliminate heteroscedastic effects in the errors or residuals. Figure 8 displays the heteroscedasticity test results of the SARIMAX(0,1,0)(1,1,0)9-GARCH(1,0) model, with the ACF and PACF plot once again. The results from Figure 8 reveals that there are no lags outside of the ACF or PACF plots. This indicates that the squared errors (variance estimators) are independent and have no correlation with their variance or variability. Thus, the heteroscedasticity of errors in the SARIMAX(0,1,0)(1,1,0)9 model has been addressed by incorporating the GARCH(1,0) model. HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 754 (a) (b) Figure 8. (a) ACF for SARIMAX-GARCH Squared Residuals, (b) PACF for SARIMAX-GARCH Squared Residuals 3.8. Forecasting Indonesia Composite Index Forecasting has the purpose of knowing the prediction of an event over the next few periods. In Table 13, we present the forecast data of Indonesia Composite Index for the next 6 month periods (October 2022 - March 2023) and compare it with actual values by using SARIMAX(0,1,0)(1,1,0)9 - GARCH(1,0) estimation. Table 13. Summary of ICI Prediction with SARIMAX(0,1,0)(1,1,0)9-GARCH(1,0) Period Actual ICI Predicted Lower Limit Average Prediction Predicted Upper Limit October 2022 7098.89 6270.19 7444.60 8838.98 November 2022 7081.31 6121.24 7412.32 8975.72 December 2022 6850.62 5929.73 7289.77 8961.75 January 2023 6839.34 5833.37 7263.27 9043.67 February 2023 6843.24 5733.58 7219.61 9090.79 March 2023 6708.93 5471.73 6960.18 8853.52 The results of the ICI prediction tend to exhibit a decreasing trend, which is consistent with the actual ICI data. Although the ICI prediction slightly deviates on a monthly average, it remains within the permissible prediction interval. The Mean Absolute Percentage Error (MAPE) for the prediction is 5.233%, attesting to its exceptional predictive accuracy. Therefore, it can be concluded that the Indonesia Composite Index (ICI) can be effectively predicted using time series SARIMAX-GARCH models, with inflation rate and exchange rate as its exogenous variables. Figure 9, derived from Table 13, illustrates the time series plot and demonstrates the model’s proficiency in predicting highly uncertain and fluctuating data, such as the Indonesia Composite Index. Figure 9. Indonesia Composite Index Prediction with SARIMAX-GARCH Estimation HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 755 3.9. A Brief of Discussion This research focuses on predicting the Indonesia Composite Index (ICI), a vital indicator for investors and a broader measure of economic growth in Indonesia. The ICI mirrors the performance of the country’s industries and is closely linked to capital flows, economic growth, and state tax revenues. Accurate forecasting of the ICI is crucial for understanding the country’s economic conditions, particularly during periods of economic instability such as the COVID-19 pandemic and recent global economic disruptions. The significance of the Composite Stock Price Index has been explored by Rosyadi et al. (2020), who examined the relationship between the Jakarta Composite Index (JCI) and Indonesia’s economic growth. Their findings indicated a positive relationship and significant feedback in the long run [39]. Other analyses, including recent studies by Haryanto et al. (2023) and Sutrisno (2024), underscore the importance of predicting the JCI given its reflection of the Exchange Market Index in Indonesia [40, 41]. While numerous studies have focused on ICI prediction, many overlook critical macroeconomic factors. For instance, Devianto et al. (2020) used artificial neural networks and nonparametric MARS models to predict the ICI, achieving commendable results with a Mean Absolute Percentage Error (MAPE) below 10% [42]. However, these models did not account for the volatile economic conditions that can drastically influence the ICI, such as significant policy changes or global economic shocks. This research introduces the SARIMAX-GARCH model, which incorporates macroeconomic factors as exogenous variables, specifically the inflation rate and exchange rate. Based on Equations 14 and 15, the negative coefficients for both first order of difference in inflation rate and exchange rate were indicating that increases in either the inflation rate or exchange rate lead to a decrease in the Composite Index. Specifically, the infla tion rate has a coefficient of -0.0089, suggesting that changes (differences) in the inflation rate have a smaller impact on the Composite Index, while the exchange rate, with a coefficient of -0.4053, exerts a more significant negative influence. This emphasizes that fluctuations in the exchange rate are more strongly correlated with movements in the Composite Index than inflation. Recent studies, such as Wijaya et al. (2023) and Pratama (2024), have also highlighted the importance of these variables in predicting financial indices [43, 44]. The SARIMAX-GARCH model’s ability to account for these factors provides a more accurate forecasting model for the ICI, enabling more informed decision-making in both the public and private sectors. This can potentially foster more stable economic growth and job creation, particularly in times of economic uncertainty. Given the economic context of the period studied, characterized by fluctuating global markets and domestic policy shifts, the SARIMAX-GARCH model's robustness in capturing these dynamics is critical. The model’s accuracy and its ability to incorporate significant macroeconomic factors make it a valuable tool for investors and policymakers alike. Policymakers can use this model to anticipate market movements, design responsive economic policies, and enhance overall economic stability. For future research, the VAR-GARCH model could be explored as an alternative approach. The VAR-GARCH model offers the advantage of capturing the dynamic interrelationships between multiple time series and accounting for volatility clustering across different variables, making it particularly useful for studying the interconnectedness of macroeconomic factors and financial markets. This approach could provide deeper insights into the systemic risks and spillover effects in emerging markets, where economic conditions are highly variable. 4. Conclusion Based on the analysis, it can be concluded that both the inflation rate and the exchange rate against the US dollar significantly influence the Indonesia Composite Index (ICI). The negative coefficients found in the first differencing of these macroeconomic factors indicate that an increase in uncertainty surrounding inflation and exchange rates tends to reduce the ICI's performance. Specifically, a larger spike in inflation or a depreciation of the rupiah against the dollar is associated with a decrease in the ICI. The SARIMAX-GARCH model used in this study has shown a high level of accuracy, as evidenced by an in-sample Mean Absolute Percentage Error (MAPE) of 0.952%, a relatively low Mean Squared Error (MSE) of 0.0003, and a high RΒ² value of 0.9782. For out-of-sample data, the MAPE was 5.233%, demonstrating the model’s robustness in capturing the dynamics of the ICI. This model not only improves predictive accuracy but also offers a deeper understanding of the intricate relationships between these critical economic indicators. Given these findings, it is imperative that investors in Indonesia closely monitor inflation rates and exchange rates, as these factors have proven to be crucial determinants of market performance. For policymakers, this study highlights the need to refine monetary policies that address macroeconomic variables, particularly inflation and exchange rates, to ensure economic stability. The significant correlations identified between the ICI, inflation, and exchange rates suggest that macroeconomic management is essential for sustaining market confidence and growth. Moreover, this research contributes to the broader goal of achieving the 8th Sustainable Development Goal (SDG), which focuses on promoting decent work and economic growth. By providing a more accurate and reliable forecasting model for the ICI, this study supports informed decision-making across both public and private sectors. This, in turn, could foster a more stable economic environment in Indonesia, facilitating job creation and sustainable growth in the long term. HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 756 5. Declarations 5.1. Author Contributions Conceptualization, M.M., E.P., and I.P.; methodology, R.R.; software, I.P. and D.U.; validation, M.M. and E.P.; formal analysis, I.P., D.U., and R.R.; investigation, M.M. and E.P.; resources, D.U.; data curation, R.R.; writingβ€” original draft preparation, I.P. and E.P.; writingβ€”review and editing, M.M., D.U., and R.R.; visualization, I.P. and D.U.; supervision, M.M. and E.P.; project administration, R.R.; funding acquisition, M.M. and E.P. All authors have read and agreed to the published version of the manuscript. 5.2. Data Availability Statement The data presented in this study are available in: β€’ Indonesia Composite Index Data: https://finance.yahoo.com; β€’ Inflation Rate Data: https://www.bi.go.id/id/statistik/indikator/data-inflasi.aspx; β€’ Exchange Rate Data: https://www.bi.go.id/id/statistik/informasi-kurs/transaksi-bi/default.aspx. 5.3. Funding The authors received no financial support for the research, authorship, and/or publication of this article. 5.4. Acknowledgements Authors express deepest gratitude to Bank Indonesia and Yahoo for their convenience in providing public data, as well as to Universitas Airlangga with all its facilities that have supported of writing this article. 5.5. Institutional Review Board Statement Not applicable. 5.6. Informed Consent Statement Not applicable. 5.7. Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. 6. References [1] Nugraha, N. M., Novan, D., & Nugraha, S. (2020). The Influence of Macroeconomic Factors on the Volatility of Composite Price Stock Index: A Study on the Indonesia Stock Exchange. International Journal of Psychosocial Rehabilitation, 24(01), 2507–2513. [2] Sumaryoto, Nurfarkhana, A., & Anita, T. (2021). The Impact of Money Supply and the Inflation Rate on Indonesia Composite Index: Case Study in Indonesia Stock Exchange 2008-2017. Business and Accounting Research International Journal, 5(2), 196– 213. [3] Republika. (2024). JCI throughout 2022 Strengthens 4.09 Percent. Republika, GΕ‚ogΓ³w, Poland. Available online: https://ekonomi.republika.co.id/berita/rnp7n1383/ihsg-sepanjang-2022-menguat-409-persen (accessed on August 2024). [4] Kristanto, M. E., & Idris. (2016). Analysis of the Influence of Inflation, Exchange Rates, and Interest Rates on the Joint Movement of JCI Stock Returns and Trading Volume for the January 2006–December 2015 Period. Diponegoro Journal of Management, 3(2337–3806), 1–15. [5] Pardede, N., Hidayat, R. R., & Sulasmiyati, S. (2016). The Influence of World Crude Oil Prices, Inflation, Central Bank Rate, and Exchange Rate on the Stock Price Index of the Mining Sector in ASEAN (Study in Indonesia, Singapore, and Thailand for the Period of July 2013–December 2015. [6] Geetha, C., Chong, V., Mohidin, R., & Chandran, V. V. (2011). The Relationship Between Inflation and Stock Market: Evidence from Malaysia, United States and China. International Journal of Economics and Management Sciences, 1(2), 1–16. [7] Akua Miyanti, G. A. D., & Wiagustini, L. P. (2018). The Effect of Fed Interest Rates, Oil Prices and Inflation on the Composite Stock Price Index (IHSG) on the Indonesia Stock Exchange. E-Journal of Economics and Business, Udayana University, 5, 1261. doi:10.24843/eeb.2018.v07.i05.p02. [8] Utama, E. Y. (2016). The Effect of SBI Interest Rates, Inflation, and Money Supply on the Composite Stock Price Index (IHSG) on the Indonesia Stock Exchange (BEI). Faculty of Economics, Universitas Negeri Yogyakarta, 4(1), 11–13. https://finance.yahoo.com/ https://www.bi.go.id/id/statistik/indikator/data-inflasi.aspx https://ekonomi.republika.co.id/berita/rnp7n1383/ihsg-sepanjang-2022-menguat-409-persen HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 757 [9] Suriyani, N. K., & Sudiarta, G. (2018). The Effect of Interest Rates, Inflation and Exchange Rates on Stock Returns on the Indonesian Stock Exchange. E-Journal of Management of Udayana University, 7(6), 3172–3200. doi:10.24843/EJMUNUD.2018.v07.i06.p12. [10] Pratikno, D. (2009). The Influence of Rupiah Exchange Rate, Inflation, SBI, and Dow Jones Index on the Movement of the Composite Stock Price Index (IHSG) in the Indonesia Stock Exchange. Journal of Economics, Universitas Sumatra Utara Medan, 16–99. [11] Tampubolon, K. M., H. M. Roy Sembel, & Sidharta, J. (2021). The Effect of Inflation, Interest Rate, Exchange Rate and Gross Domestic Product Growth on Composite Stock Price Index (Study at Indonesia Stock Exchange Period of 2010-2019). Fundamental Management Journal, 6(1), 66–90. doi:10.33541/fjm.v6i1.2826. [12] Hasanudin. (2021). The Effect of Inflation, Exchange Rate, BI Rate and Dow Jones Index to Indonesia Composite Index in Indonesia Stock Exchange on the Year 2013-2018. International Journal of Science and Society, 3(3), 50–60. doi:10.54783/ijsoc.v3i3.353. [13] Amin, M. Z. (2012). The Effect of Inflation Rate, SBI Rate, Dollar Value (USD/IDR), and Dow Jones Index (DJIA) to Movement of Composite Stock Price Index In Indonesia Stock Exchange (IDX)(Periode 2008–2011). Thesis Journal FEB UB, 1–17. [14] Christopher, A., Sediono, E. A., & Suliyanto, Dan M. F. (2021). Application of the ARIMAX-GARCH Model in Modeling and Forecasting the Volume of Electronic Money Transactions in Indonesia. Journal of Mathematics Education, Science and Technology, 2(6), 241–256. [15] Hamjah, M. A. (2014). Forecasting Major Fruit Crops Productions in Bangladesh using Box-Jenkins ARIMA Model. Journal of Economics and Sustainable Development, 5(7), 96–107. [16] Chyon, F. A., Suman, M. N. H., Fahim, M. R. I., & Ahmmed, M. S. (2022). Time series analysis and predicting COVID-19 affected patients by ARIMA model using machine learning. Journal of Virological Methods, 301, 114433. doi:10.1016/j.jviromet.2021.114433. [17] Abonazel, M. R., & Abd-Elftah, A. I. (2019). Forecasting Egyptian GDP using ARIMA models. Reports on Economics and Finance, 5(1), 35–47. doi:10.12988/ref.2019.81023. [18] Wei, W. (2006). The Oxford Handbook of Quantitative Methods in Psychology Vol2: Statistical Analysis. Oxford University Press, Oxford, United Kingdom. [19] Kawakita, S., & Takahashi, H. (2022). Time-series analysis of population dynamics of the common cutworm, Spodoptera litura (Lepidoptera: Noctuidae), using an ARIMAX model. Pest Management Science, 78(6), 2423–2433. doi:10.1002/ps.6873. [20] Paul, R. K. (2015). ARIMAX-GARCH-WAVELET model for forecasting volatile data. Model Assisted Statistics and Applications, 10(3), 243–252. doi:10.3233/MAS-150328. [21] Simms, L. E., Engebretson, M. J., & Reeves, G. D. (2022). Removing Diurnal Signals and Longer-Term Trends from Electron Flux and ULF Correlations: A Comparison of Spectral Subtraction, Simple Differencing, and ARIMAX Models. Journal of Geophysical Research: Space Physics, 127(2), 2021 030021. doi:10.1029/2021JA030021. [22] Alharbi, F. R., & Csala, D. (2022). A Seasonal Autoregressive Integrated Moving Average with Exogenous Factors (SARIMAX) Forecasting Model-Based Time Series Approach. Inventions, 7(4), 94. doi:10.3390/inventions7040094. [23] Bari, S. H., Rahman, M. T., Hussain, M. M., & Ray, S. (2015). Forecasting Monthly Precipitation in Sylhet City Using ARIMA Model. Civil and Environmental Research, 7(1), 69–77. [24] Adesi, G. B., Engle, R. F., & Mancini, L. (2014). A GARCH option pricing model with filtered historical simulation. Simulating Security Returns: A Filtered Historical Simulation Approach, 21(3), 66–108. doi:10.1057/9781137465559.0009. [25] Endri, E., Aipama, W., Razak, A., Sari, L., & Septiano, R. (2021). Stock price volatility during the COVID-19 pandemic: The GARCH model. Investment Management and Financial Innovations, 18(4), 12–20. doi:10.21511/imfi.18(4).2021.02. [26] Kim, S., & Kim, H. (2016). A new metric of absolute percentage error for intermittent demand forecasts. International Journal of Forecasting, 32(3), 669–679. doi:10.1016/j.ijforecast.2015.12.003. [27] Yao, K. C., Hsueh, H. W., Huang, M. H., & Wu, T. C. (2022). The Role of GARCH Effect on the Prediction of Air Pollution. Sustainability (Switzerland), 14(8), 4459. doi:10.3390/su14084459. [28] Şentaş, A., Tashiev, Δ°., KΓΌΓ§ΓΌkayvaz, F., Kul, S., Eken, S., Sayar, A., & Becerikli, Y. (2020). Performance evaluation of support vector machine and convolutional neural network algorithms in real-time vehicle type and color classification. Evolutionary Intelligence, 13, 83-91. doi:10.1007/s12065-018-0167-z. [29] Ghozali, I. (2016). Multivariate Analysis Application Using IBM SPSS 20. Semarang: Badan Penerbit Universitas Diponegoro, Indonesia. HighTech and Innovation Journal Vol. 5, No. 3, September, 2024 758 [30] Bagwell, C. B., Hill, B. L., Herbert, D. J., Bray, C. M., & Hunsberger, B. C. (2016). Sometimes simpler is better: VLog, a general but easy-to-implement log-like transform for cytometry. Cytometry Part A, 89(12), 1097–1105. doi:10.1002/cyto.a.23017. [31] Xue, D., & Hua, Z. (2016). ARIMA Based Time Series Forecasting Model. Recent Advances in Electrical & Electronic Engineering (Formerly Recent Patents on Electrical & Electronic Engineering), 9(2), 93–98. doi:10.2174/2352096509999160819164242. [32] Garima Jain, E., & Mallick, B. (2017). A Study of Time Series Models ARIMA and ETS. In International Journal of Modern Education and Computer Science, 9(4), 57-63. doi:10.5815/ijmecs.2017.04.07. [33] Abunofal, M., Poshiya, N., Qussous, R., & Weidlich, A. (2021). Comparative Analysis of Electricity Market Prices Based on Different Forecasting Methods. 2021 IEEE Madrid PowerTech, PowerTech 2021 - Conference Proceedings, 1–6. doi:10.1109/PowerTech46648.2021.9495034. [34] Shadkam, A. (2020). Using SARIMAX to forecast electricity demand and consumption in university buildings. Doctoral dissertation, University of British Columbia, British Columbia, Canada. [35] Krueger, D., Mitman, K., & Perri, F. (2016). Macroeconomics and Household Heterogeneity. Handbook of Macroeconomics, Elsevier, Volume 2, 843–921. doi:10.1016/bs.hesmac.2016.04.003. [36] Moffat, I. U., Akpan, E. A., & Abasiekwere, U. A. (2017). A time series evaluation of the asymmetric nature of heteroscedasticity: an EGARCH approach. International Journal of Statistics and Applied Mathematics, 2(6), 111–117. [37] Aue, A., HorvΓ‘th, L., & Pellatt, D. F. (2017). Functional Generalized Autoregressive Conditional Heteroskedasticity. Journal of Time Series Analysis, 38(1), 3–21. doi:10.1111/jtsa.12192. [38] Dritsaki, C. (2017). An Empirical Evaluation in GARCH Volatility Modeling: Evidence from the Stockholm Stock Exchange. Journal of Mathematical Finance, 07(02), 366–390. doi:10.4236/jmf.2017.72020. [39] Rosyadi, M. I., Widowati, A., & Jamil, P. C. (2020). Analysis of Correlation Between Indonesian Composite Index and Economic Growth in Indonesia. European Journal of Business and Management, 12(8), 42–46. doi:10.7176/ejbm/12-8-07. [40] Lim, K., & McNelis, T. (2013). The Impact of Economic Structure on Stock Market Performance: A Comparative Study of Indonesia and Malaysia. Pacific-Basin Finance Journal, 24(2), 256–271. [41] Sutrisno, B. (2024). Macroeconomic Variables and Stock Market Volatility: An Empirical Analysis of Indonesia. Asian Economic and Financial Review, 14(1), 24-38. [42] Devianto, D., Permathasari, P., Yollanda, M., & Wirahadi Ahmad, A. (2020). The Model of Artificial Neural Network and Nonparametric MARS Regression for Indonesian Composite Index. IOP Conference Series: Materials Science and Engineering, 846(1), 12007. doi:10.1088/1757-899X/846/1/012007. [43] Wijaya, S., & Handoko, F. (2023). Forecasting the Composite Index with SARIMAX-GARCH Models in the Context of the COVID-19 Pandemic. Indonesian Journal of Economics and Business, 32(4), 89–105. [44] Pratama, R., & Susilo, E. (2024). Examining the Relationship between Inflation, Exchange Rates, and Stock Market Indices: Evidence from Indonesia. Global Journal of Economics, 19(1), 75–90.