1. Introduction
Since early 2020, COVID-19 has been widely spreading around the globe. It has caused people to suffer in many aspects of life. Besides, it forced people to stop performing their activities normally. Amidst the recovery and containment, the world economic system is characterized as experiencing significant, broad uncertainty. Economic forecasts and consensus among macroeconomics experts show significant disagreement on the overall extent, long-term effects, and projected recovery (Khan et al., 2020). Hanoatubun (2020) found that this pandemic made it difficult for Indonesians to obtain their basic needs; it also reduced their consumption, thus decreasing economic growth. Consequently, governments worldwide have attempted to change the way people live by strongly urging them to stay at home to stop the spread of the virus. This policy has also been implemented in Indonesia being one of the countries greatly affected by the virus, especially in the area of banking.
The big four state-owned banks listed on the Indonesia Stock Exchange (IDX) are Bank Rakyat Indonesia (BBRI), Bank Negara Indonesia (BBNI), Bank Mandiri (BMRI) and Bank Tabungan Negara (BBTN). Together, they account for 43% of the total assets of all banks in Indonesia (Otoritas Jasa Keuangan, 2020). They also suffered as the pandemic was announced by the World Health Organization (WHO). This condition has been reflected in their respective semi-annual financial report. In Semester 1 2020, all four banks recorded a significant loss on a year-to-year basis. BBNI had the most significant drop in earning; it suffered a loss of 41.54% compared with that in Semester 1 2019, followed by BBTN, whose net profit significantly decreased by 40%. Furthermore, BBRI as the most capitalized state-owned bank included in the IDX also experienced a net profit loss of 36.88%. BMRI also recorded a decrease in earnings of 23.94%, which was similar to that in Semester 1 2019. In the credit measurement, the non-performing loan (NPL) increased to more than 3% for the state-owned banks, except for BBTN, whose NPL decreased from 2.42% to 2.40% in June 2020. Although the ratio is still at a normal rate, it indicates that the shocks might affect the lack of ability of creditors to repay their debts.
Investors might perceive this phenomenon as an investment risk, as the shock is likely to influence their returns either positively or negatively at a specific case when investing in financial sectors. Investment risk can be defined as the probability or likelihood of occurrence of losses relative to the expected return on any particular investment. Some studies have proven that financial data is highly volatile as market uncertainty is quite high (Hendrawaty et al., 2021). Therefore, the measurement of investment risk is important to minimize the loss (Rahman et al., 2020) risk, and cost-inefficiency for a sample of 30 commercial banks in Bangladesh from 2006 to 2018. To conduct the analysis, we used the Generalized Methods of Moments (GMM. Value at Risk (VaR) a statistical tool to measure and quantify financial risk within a firm or portfolio over a specific time frame. This metric is often used by Banks to determine the extent and probability of occurrence of a potential loss on the advances.
To some extent, in the measurement of market risk, Emenogu et al. (2020) defined volatility as the statistical measurement of the return distribution on a provided security that can be evaluated by using standard deviation or variance among the returns. In this case, the application of Generalized Autoregressive Conditional Heteroscedasticity (GARCH) statistically provides a more accurate model of variance. Bollerslev (1986) first introduced the GARCH model to avoid the high order of autoregressive conditional heteroscedasticity (ARCH), which was previously available to model variance (Engle, 1982).
Akhmadi et al. (2019) estimated VaR using the extreme value theory (EVT) and generalized Pareto distribution (GPD) for state-owned bank enterprises in Indonesia after the occurrence of the global financial crisis (GFC) in 2008. They suggested that the GPD method outperforms the EVT and is very close to the banks’ capital adequacy ratio (CAR); moreover, the two methods have a higher value compared with the other methods. Meanwhile, for extreme data in the financial sector, the empirical study conducted by Budiarti (2019) demonstrated that the AR-GARCH- copula Tawn approach is the best fit for modeling the joint distribution of a portfolio that can be used as the basis for the calculation of the VaR in extreme cases. Vo et al. (2019) return and portfolio diversification at the industry level in four member countries of the ASEAN for which required data are available: Vietnam, Thailand, Malaysia, and Singapore. Market indices are examined for 10 industries from 2007 to 2016, comprising different economic periods, including 2007–2009 (crisis used conditional value at risk (CVaR) to measure extreme risk. Markowitz’s risk-return framework is utilized to determine the optimal weight of industries in the portfolio.
Furthermore, some previous studies have employed the econometric model of GARCH to estimate the volatility of variance to calculate the percentage of the maximum loss from the return of the given portfolio. In Indonesia, Sitorus (2018) stated that one consideration for investors to make an investment decision in stock instruments by monitoring the daily volatility movement and the trend. From the risk perspective, measuring the maximum probability of risk level in the future can be a consideration in making stock investment decisions that have the same movements and characteristics. Risk projection based on historical data with a certain time period can be calculated with Value at Risk (VaR) with a certain level of confidence. Thus the investment decision will be optimal. The Value at Risk model will calculate the expected losses. This research will show that volatility that looks the same but has a different level of risk. Mutia et al. (2018) used the ARMA-GARCH model to determine the VaR of the portfolio of combined stock prices from five companies included in the IDX. She noticed that BBNI has the riskiest assets in the portfolio. Conversely, Sukono et al. (2019) used the ARMA-FIGARCH model to measure the VaR of some IDX stocks and found that BBRI has the highest VaR compared with the other companies included in the sample. Furthermore, in the study conducted by Denkowska and Wanat (2020), the C-DCC- GARCH model was used to measure conditional variances to determine the systematic risks of European insurances during the crisis in 2020, while Le and Tran (2021) used that model to measure the effects of US stock market on both Vietnamese and Philippines Stock Markets.
Therefore, this study aims to measure the risk returns by computing the respective VaRs of the daily stock price of state-owned banks in Indonesia. The GARCH model was used to estimate the parameters of means and variances during the economic shock period.
2. Research Methods
In this study, the daily stock prices of state-owned banks in Indonesia from the COVID-19 pandemic in early 2020 until the fourth quarter of 2020 were observed. Prior to computing the risk-return using the VaR, the return volatility of each bank was measured using the GARCH model to obtain a more accurate measurement in the ES approach. The stages of the GARCH model are as follows.
2.1. Stationarity of Data
In the GARCH model, the time-series data should be stationary, or the means and variances should fluctuate around zero. A stationary test could be conducted visually by reading the plotting data and statistically using the augmented Dickey-Fuller (ADF) test, which was initially introduced by Dickey and Fuller (1979), Y2, ..., Yn be generated by the model Rao (1961). Mathematically, the equation of the ADF test is as follows:
\(\mathrm{DF}_{\tau}=\frac{\gamma_{i}}{\widehat{S e}_{\gamma_{1}}}\) (1)
The above formula was used to construct the following hypothesis:
H0: DFτ > 2.57 = Non-stationary
H1: DFτ < 2.57 = Stationary (Brockwell & Davis, 2002)
Tsay (2010) developed the measurement of stationarity data in which he tested stationary time-series data by computing its autocorrelation function (ACF) and partial autocorrelation function (PACF). The indication of stationary time-series data can be studied from the data movement of the ACF and PACF graphs, in which the slow movement of data indicates non-stationarity.
However, Ambya et al. (2020) FNG price is categorised as a time series data with volatility in both variance and mean, as well as most likely in some cases having heteroscedasticity problem. To come up with an estimated prediction model, some analysis tools, such as autoregressive integrated moving average (ARIMA proved that financial data is classified as non-stationary data. Hence, in this stage, non-stationary data need to be transformed into stationary data. One of the ways to do this is by employed differencing approach. Such an approach was initially conducted by Granger and Joyeux (1980). The aim of this process was to stabilize the means and variances of the time-series data, particularly financial data, by doing changes in the time-series data (Hyndman & Athanasopoulos, 2018). Once the data is already stationary, the next step will be performed.
2.2. ARCH Effect Test
One concern in modeling financial data is heteroscedasticity (Engle, 1982), which renders the model prediction less accurate. To solve this issue, Tsay (2010) suggested conducting a test on the ARCH effect. To determine whether the ARCH effect is included in the model, the Lagrange Multiplier (LM) test should be conducted (Lee & King, 1993). If at any lag of order the p-value obtained from the LM test is significant (<0.001), then it can be concluded that heteroscedasticity exists. Therefore, we need longer lags and generalize it into a proper model (Wong & Li, 1995), which is the GARCH model.
2.3. Mean Model of AR (p) and Variance Model of GARCH (p, q)
The GARCH model is computed by squaring the past residuals to estimate the variances as it has a long memory process and involves heteroscedasticity (Tsay, 2010). The mean model of AR (p) is defined as order lag p, and the variance model of GARCH (p, q) is defined as conditional variances and squared residuals (Tsay, 2010). The equation for the AR (p)-GARCH (p, q) model is as follows:
\(\mathrm{SP}_{x t}=\vartheta+\sum_{i=1}^{p} \omega_{i} \mathrm{SP}_{x t-i}+\varepsilon_{t}\) (2)
\(\sigma_{t}^{2}=c+\sum_{i=1}^{q} \gamma_{i} \varepsilon_{t-i}^{2}+\sum_{j=1}^{p} \theta_{j} \sigma_{t-j}^{2}\) (3)
where SPxt denotes each state-owned bank’s daily stock prices at time t; ϑ, constant; and ωi, coefficient regression of AR (p) with i = 1, 2, 3, …, p; c is constant, and γi is the coefficient regression of MA (q) with i = 1, 2, 3, …, q.
2.4. Value at Risk (VaR) Measurement
Value at risk is a measure of the risk of loss for investments. It estimates how much a set of investments might lose, given normal market conditions, in a set time period such as a day. VaR has become a common tool for risk measurement in financial agencies. The formal definition of VaR can be as a measurement of portfolio value that might be considered a loss for any given probability level. Thus, measuring the downside risk is important to control the internal risk and to make financial regulations (Meng & Taylor, 2020). The study by Akhmadi et al. (2019) demonstrated that measuring risk level by applying the VaR approach can be defined as the estimated maximum losses for portfolio investments at a certain period and probability level. Tsay (2010) mathematically demonstrates the measurement of VaR with a certain time horizon and confidence level as follows:
\(\operatorname{VaR}_{(1-\alpha)}(t)=W_{0} \times(\mu-R) \sqrt{t}\) (4)
where W0 denotes the initial investment value or portfolio; R, the quantile value αth from the distribution of stock prices; σ, volatility; and t, time horizon.
3. Results
3.1. Data Description
The data obtained for this study includes the daily stock prices from four state-owned banks in Indonesia from the start of the COVID-19 pandemic in January 2020 to the fourth quarter of 2020. The distribution of the daily stock prices for each bank is described in Figure 1.
Figure 1: Distribution of the Daily Stock Prices for (a) Bank Rakyat Indonesia (Code: BBRI), (b) Bank Negara Indonesia (Code: BBNI), (c) Bank Tabungan Negara (Code: BBTN) and (d) Bank Mandiri (Code: BMRI)
Overall, Figure 1 demonstrates that the daily stock prices of all state-owned banks in Indonesia during the study period have fluctuated. Since the COVID-19 pandemic has been announced in Indonesia in March 2020, all stock prices have significantly decreased. However, in April, the returns have increased again due to the adjustments made by the government for the Indonesian financial institutions. The implementation of the ‘New Normal’ policy in the early third quarter of 2020 has made the movements of stock prices less volatile for all banks. In the fourth quarter of 2020, it can be noticed that all stock prices gradually increased due to market makers’ beliefs.
The graphs also demonstrate that the means and variances of the stocks are not maintained at around zero, indicating that all the data sets are visually non-stationary. It can be statistically proven by employing the ADF unit-root test.
Statistically, the results of the ADF test presented in Table 1 demonstrate that the probability values for all returns are more than 0.0001. This means that we do not reject the null hypothesis. It implies that the time-series data are not yet stationary in terms of the means and variances. Therefore, it is important to perform the next step of differencing. The non-stationary time series can also be proven using the graph of the ACF and PACF for each bank, which can be illustrated by the following graphs.
Table 1: Augmented Dickey-Fuller Test for all State-Owned Banks in Indonesia
Figure 2 demonstrates that all the ACF graphs exhibit a slow movement, indicating that after lag 1, the data set distributions are not on the circle of zero, of which it is confirmed non-stationary data set. Besides, the data sets on the PACF graphs at lag 1 are not at around zero which assigns time series as non-stationary data set.
Figure 2: ACF and PACF graphs for (a) Bank Rakyat Indonesia (Code: BBRI), (b) Bank Negara Indonesia (Code: BBNI), (c) Bank Tabungan Negara (Code: BBTN) and (d) Bank Mandiri (Code: BMRI)
3.2. Stationarity Transformation
Since the data sets have been statically indicated as non-stationary data, the next step is to transform them into stationary data by employing the differencing method. Table 2 demonstrates that the ADF test once differencing lag 1 (d = 1) was conducted for all data set to obtain stationary data.
Table 2: Augmented Dickey-Fuller Test for all State-Owned Banks in Indonesia After Differencing 1 (d = 1)
To verify this transformation, Figure 3 presents the graphs for the trend and correlation analysis of all the data sets. It can be observed from the figures that differencing 1 (d = 1) makes the data set have the behavioral of residual distributions at around zero mean. It can also be seen that the ACF and PACF graphs exhibit a delayed movement, indicating stationarity.
Figure 3: Trend and Correlation Analysis Graphs For (a) Bank Rakyat Indonesia (Code: BBRI), (b) Bank Negara Indonesia (Code: BBNI), (c) Bank Tabungan Negara (Code: BBTN) and (d) Bank Mandiri (Code: BMRI)
3.3. ARCH Effect Identification
To successfully create an accurate GARCH model, it is important to determine whether or not the stationary data set of the daily stock prices of all four banks have issues with heteroscedasticity. Table 3 presents the Portmanteau Q and LM test with all data set having a significant value of <0.00001 that conclude in having the issue. Therefore, the outcomes suggest that residual data sets can be applied to the GARCH(p, q) model to estimate their volatility.
Table 3: Tests for ARCH Disturbances Based on Residuals for Each Bank
3.4. GARCH Model
Conditional heteroscedasticity is important to create a good-fit forecasting model. The AR (p)-GARCH (p, q) model can then be run to do so, where AR (p) is conditional to have a mean model and GARCH (p, q) is to have a model of variances and squared residuals for each data set.
Table 4 presents the parameter estimates for the construction of the AR (1)-GARCH (1, 1) model for each bank. The models are believed to have a good-fit measurement to make forecasting as probability values for all banks stated in Table 4 are less than the confidence interval of 5%. This indicates that AR (1) and GARCH (1, 1) estimate the means and variances, respectively, of the four banks. For each bank, the AR (1)-GARCH (1, 1) model is applied as follows:
Table 4: Estimation of the Parameters of AR (1)-GARCH (1, 1) For Each Bank
a. BBRI mean and variance models
AR (1): BBRIt = 4240 – 0.9911 BBRIt–1
GARCH (1, 1):
\(\sigma_{t}^{2}=856.162+0.2055 \varepsilon_{t-1}^{2}+0.711 \sigma_{t-1}^{2}\)
b. BBNI mean and variance models
AR (1): BBNIt = 6395 – 0.9958 BBNIt–1
GARCH (1, 1):
\(\sigma_{t}^{2}=4209+0.1751 \varepsilon_{t-1}^{2}+0.6575 \sigma_{t-1}^{2}\)
c. BBTN model
AR (1): BBTNt = 1715 – 0.9937 BBTNt–1
GARCH (1, 1):
\(\sigma_{t}^{2}=1008+0.2605 \varepsilon_{t-1}^{2}+0.3623 \sigma_{t-1}^{2}\)
d. BMRI model
AR (1): BMRIt = 6540 – 0.9898 BMRIt–1
GARCH (1, 1):
\(\sigma_{t}^{2}=2029+0.2107 \varepsilon_{t-1}^{2}+0.7301 \sigma_{t-1}^{2}\)
3.5. Measurement of VaR
The AR (1)-GARCH (1, 1) models are then applied as the basis for the computation of VaR, where AR (1) is used to calculate the means and GARCH (1, 1) to determine the variances. From Appendix 1, the t–1 data of each bank can be obtained: BBRI242 = 4180; BBNI242 = 6305; BBTN242 = 1725; and BMRI242 = 6325. Thus, the mean values for t = 243 are as follows:
\(\begin{aligned} \text { a. } \quad \mathrm{BBRI}_{243}&=4240-0.9911 \mathrm{BBRI}_{242} \\ \mathrm{BBRI}_{243}&=4240-0.9911(4180) \\ \mathrm{BBRI}_{243}&=97.202 \\ \text { b. } \quad \mathrm{BBNI}_{243}&=6395-0.9958 \mathrm{BBNI}_{242} \\ \mathrm{BBNI}_{243}&=6395-0.9958(6305) \\ \mathrm{BBNI}_{243}&=116.48 \\ \text { c. } \quad \mathrm{BBTN}_{243} &=1715-0.9937 \mathrm{BBTN}_{242} \\ \mathrm{BBTN}_{243} &=1715-0.9937(1725) \\ \mathrm{BBTN}_{243} &=0.8675 \\ \text { d. } \quad \mathrm{BMRI}_{243} &=6540-0.9898 \mathrm{BMRI}_{242} \\ \mathrm{BMRI}_{243}&=6540-0.9898(6325) \\ \mathrm{BMRI}_{243}&=279.511 \end{aligned}\)
The volatility values are as follows:
a. Volatility value for BBRI
\(\begin{aligned} \sigma_{243}^{2} &=856.162+0.2055 \varepsilon_{242}^{2}+0.711 \sigma_{242}^{2} . \\ \sigma_{243}^{2} &=856.162+0.2055(4180)+0.711(10253) \\ \sigma_{243}^{2} &=9067.94 \\ \sigma_{243} &=95.22 \end{aligned}\)
b. Volatility value for BBNI
\(\begin{aligned} \sigma_{243}^{2} &=4209+0.1751 \varepsilon_{t-1}^{2}+0.6575 \sigma_{t-1}^{2} \\ \sigma_{243}^{2} &=4209+0.1751(6305)+0.6575(25247.2) \\ \sigma_{243}^{2} &=21913.04 \\ \sigma_{243} &=148.03 \end{aligned}\)
c. Volatility value for BBTN
\(\begin{aligned} \sigma_{243}^{2} &=1008+0.2605 \varepsilon_{t-1}^{2}+0.3623 \sigma_{t-1}^{2} \\ \sigma_{243}^{2} &=1008+0.2605(1725)+0.3623(2672.55) \\ \sigma_{243}^{2} &=2425.62 \\ \sigma_{243} &=49.25 \end{aligned}\)
d. Volatility value for BMRI
\(\begin{aligned} \sigma_{243}^{2} &=2029+0.2107 \varepsilon_{t-1}^{2}+0.7301 \sigma_{t-1}^{2} \\ \sigma_{243}^{2} &=2029+0.2107(6325)+0.7301(34284.7) \\ \sigma_{243}^{2} &=33748.3 \\ \sigma_{243} &=183.70 \end{aligned}\)
The above calculations enable the construction of VaR for each selected stock price with a 5% confidence interval (a 1.65 standard deviation) and time horizon of 15 days, as shown in Table 5.
Table 5: VaR Computation for 15 days Ahead with 95% Confidence Interval
4. Discussion
Table 5 presents the VaR of the stock prices for the next 15 days. With a confidence interval of 95%, all stock prices are believed to decrease by a maximum of Rp232.03 for BBRI, Rp494.84 for BBNI, Rp311.36 for BBTN, and Rp91.36 for BMRI. The decrease in the stock prices in the banking sector is caused by the unstable economic condition around the globe, most particularly in Indonesia. The VaR results confirmed what the banks published in their annual report, i.e. in 2020, the COVID-19 pandemic has slowed down economic activities, causing many creditors to fail in their obligations to the banks.
VaR measurement can be a consideration for investors when putting their funds on the stock market. Table 5 demonstrates that the stock prices have maintained a downward trend for 15 days. Thus, it can be a recommendation for investors with risk-taker characteristics to go short, or keep that investment and wait until some expecting increasing trends occurring for risk-averse investors. Therefore, it can be deduced that the calculation of VaR in stock prices for state-owned banks in Indonesia is to help investors making decisions either to go long or short for their portfolios to minimize the risk of risk.
5. Conclusion
The COVID-19 pandemic has shocked the people worldwide, forcing them to adjust the way they live or implement the ‘New Normal’ policy. This is the government’s way to restabilize economic activities in Indonesia, and it is relatively working to some extent. The COVID-19 pandemic has affected financial institutions. With the unstable economy in Indonesia, all state-owned banks suffered from a significant profit drop, as people have less activity in the banking industries. Thus, it is important to study the potential loss when investing in the Indonesian banking sector.
The VaR is used to measure the maximum loss that might occur in the next couple of days with a certain confidence interval using the GARCH model to estimate the means and variances. AR (1, 1)-GARCH (1, 1) is believed to be a good-fit model for constructing the means and variances, which are accurate elements for the measurement of the VaR. The finding suggests that the VaR varies among banks’ stock prices, with the largest maximum potential drop occurring in BBNI (494.84); BMRI has the smallest value of VaR, i.e. 91.36. The results of the study suggest that investors who have funds in state-owned banks should reconsider their investments.
References
- Akhmadi, Y., Mustofa, I., Rika, H. M., & Hanggraeni, D. (2019). Value at risk assessment uses the extreme value theory approach and generalized Pareto distribution case studies of state-owned banks in Indonesia in the period 2008-2018. Managament Insight: Jurnal Ilmiah Manajemen, 13(1), 63-72. https://doi.org/10.33369/insight.14.1.63-72
- Ambya, G. T., Hendrawaty, E., Kesumah, F. S. D., & Wisnu, F. K. (2020). Future natural gas price forecasting model and its policy implication. International Journal of Energy Economics and Policy, 10(5), 58-63. https://doi.org/10.32479/ijeep.9676
- Bollerslev, T. (1986). Generalized autoregressive conditional heteroskedasticity. Journal of Econometrics, 31, 307-327. https://doi.org/10.1109/TNN.2007.902962
- Brockwell, P., & Davis, R. (2002). Introduction to time series and forecasting. Springer-Verlag.
- Budiarti, R. (2019). Estimated portfolio value-at-risk for extreme data. E-Prosiding Nasional Seminar Nasional Statistika, 8(1), 111-129. http://prosiding.statistics.unpad.ac.id/index.php/prosiding/article/view/447
- Denkowska, A., & Wanat, S. (2020). Dependencies and systemic risk in the European insurance sector: New evidence based on copula-DCC-Garch model and selected clustering methods. Entrepreneurial Business and Economics Review, 8(4), 7-27. https://doi.org/10.15678/EBER.2020.080401
- Dickey, D. A., & Fuller, W. A. (1979). Distribution of the estimators for autoregressive time series with a unit root. Journal of the American Statistical Association, 74(366), 427. https://doi.org/10.2307/2286348
- Emenogu, N. G., Adenomon, M. O., & Nweze, N. O. (2020). On the volatility of daily stock returns of total Nigeria Plc: Evidence from GARCH models, value-at-risk and backtesting. Financial Innovation, 6(1), 1-25. https://doi.org/10.1186/s40854-020-00178-1
- Engle, R. F. (1982). Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation. Econometrica, 50(4), 987. https://doi.org/10.2307/1912773
- Granger, C. W., & Joyeux, R. (1980). An introduction to long-memory time series models and fractional differencing. Journal of Time Series Analysis, 1(1), 15-29. https://doi.org/10.1111/j.1467-9892.1980.tb00297.x
- Hanoatubun, S. (2020). Effect of profitability, leverage, liquidity and company size on stock returns of pharmaceutical companies in Bei. Journal of Education, Psychology and Counseling, 2(1), 46-69. https://ojs.unud.ac.id/index.php/Manajemen/article/view/23852
- Hendrawaty, E., Azhar, R., Kesumah, F. S. D., Sembiring, S. I. O., & Metalia, M. (2021). Modeling and forecasting crude oil prices during the Covid-19 pandemic. International Journal of Energy Economics and Policy, 11(2), 149-154. https://doi.org/10.32479/ijeep.10578
- Hyndman, R., & Athanasopoulos, G. (2018). Forecasting: principles and practice (2nd ed.). Victoria, Australia: O'Texts Publishers.
- Khan, K., Zhao, H., Zhang, H., Yang, H., Shah, M. H., & Jahanger, A. (2020). The impact of COVID-19 pandemic on stock markets: An empirical analysis of world major stock indices. Journal of Asian Finance, Economics, and Business, 7(7), 463-474. https://doi.org/10.13106/jafeb.2020.vol7.no7.463
- Le, T. P. T. D., & Tran, H. L. M. (2021). The contagion effect from U.S. stock market to the Vietnamese and the Philippine stock markets: The evidence of DCC - GARCH model. Journal of Asian Finance, Economics. and Business, 8(2), 759-770. https://doi.org/10.13106/jafeb.2021.vol8.no2.0759
- Lee, J. H., & King, M. . (1993). A locally most mean powerful-based score test for ARCH and GARCH regression disturbances. Journal of Business and Economic Statistics, 11(1), 17-27. https://doi.org/10.1080/07350015.1993.10509930
- Meng, X., & Taylor, J. W. (2020). Estimating value-at-risk and expected shortfall using the intraday low and range data. European Journal of Operational Research, 280(1), 191-202. https://doi.org/10.1016/j.ejor.2019.07.011
- Mutia, E., Rahmawaty, R., & Afrianandra, C. (2018). Value at risk of Sukuk Ijarah and Mudharabah in Indonesia. Journal of Accounting Research, Organization and Economics, 1(1), 65-73. https://doi.org/10.24815/jaroe.v1i1.10751
- Otoritas Jasa Keuangan. (2020). Indonesian banking statistics. https://www.ojk.go.id/id/kanal/perbankan/data-dan-statistik/statistik-perbankan-indonesia/Documents/Pages/Statistik-Perbankan-Indonesia---Desember-2020/Statistik%20Perbankan%20Indonesia%20Des%202020.pdf
- Rahman, M. M., Chowdhury, A. A., & Moudud-Ul-Huq, S. (2020). How do the banks determine regulatory capital, risk, and cost inefficiency in Bangladesh? Journal of Asian Finance, Economics, and Business, 7(12), 211-222. https://doi.org/10.13106/jafeb.2020.vol7.no12.211
- Sitorus, S. (2018). Investment decision-making based on the value at risk (VaR) analysis for stocks of state own bank in Indonesia. Journal of Economics and Business, 2(2), 128-141. https://doi.org/10.25170/jebi.v2i2.42
- Sukono, F., Lesmana, E., Susanti, D., Napitupulu, H., & Hidayat, Y. (2019). Estimation of value-at-risk adjusted under the capital asset pricing model based on the ARMAX-GARCH approach. Jurnal Matematika Integratif, 15(1), 29. https://doi.org/10.24198/jmi.v15i1.20931
- Tsay, R. S. (2010). Analysis of financial time series (3rd ed.). Hoboken, NJ: Wiley & Sons. https://doi.org/10.1002/9780470644560
- Vo, D. H., Pham, T. N., Pham, T. T. V., Truong, L. M., & Cong Nguyen, T. (2019). Risk, return, and portfolio optimization for various industries in the ASEAN region. Borsa Istanbul Review, 19(2), 132-138. https://doi.org/10.1016/j.bir.2018.09.003
- Wong, H., & Li, W. K. (1995). Portmanteau test for conditional heteroscedasticity, using ranks of squared residuals. Journal of Applied Statistics, 22(1), 121-134. https://doi.org/10.1080/757584402