1. Introduction
This study aims to examine the Baltic Dry Index (BDI), the Investor Sentiment Index (VIX), and the MSCI index, which are considered important indicators regarding providing early signals to predict down-turns and crises in the global economy. Among the three indices, the VIX and MSCI are the more commonly known indices in contrast to the BDI. Though the prediction of future fluctuations in markets is highly relevant, the modeling of indices that are considered leading indicators is decidedly challenging. To overwhelm the uncertainty of the future, some researchers debated indicators such as the VIX and MSCI in addition to the BDI and questioned whether the so-called leading indicators provide significant insight, especially for the indices reported in daily frequency.
The VIX is based on future contracts and options to achieve information regarding future behavior in financial markets. The MSCI is calculated through a weighting scheme through the use of developed and emerging stock markets. However, though relatively less attended, the BDI is based on worldwide shipment costs of goods and materials through major shipment routes, and the index is directly related to supply and demand conditions, a characteristic that relates the index directly to economic production. On the other hand, all three indices have a common characteristic: their highly volatile and nonlinear nature, which leads to a loss in efficiency in empirical models which aim at predicting and forecasting these series to utilize them as leading indicators. To this end, the paper aims at the utilization of Logistic Smooth Transition Autoregressive (LSTAR) augmented variants of traditional GARCH (LSTAR-GARCH) in addition to LSTAR augmented Asymmetric Power GARCH (LSTAR-APGARCH) models to capture the regime-dependent nonlinear and volatility dynamics to achieve important tools to predict and forecast the movements in the BDI, VIX, and MSCI indices. The general framework of these indices is evaluated, which will also provide important insights.
Though the index is not well explored in the literature, the BDI has a long historical background dating back to 1744. The BDI was originated by the first committee of 23 transport merchants. This development led to the establishment and development of a company in 1857. The Baltic Dry was developed in 1985 as the Baltic Dry Index (BDI) [
1]. The BDI index is calculated by the Baltic Exchange, which is a London-based organization that provides vital information on maritime trade and shipment rates. In its calculation, Baltic Exchange utilizes major routes and shipment rates for various different types of vessel sizes. In the calculation, the shipment costs for the four major vessel sizes, Capesize, Panamax, Supramax, and Handysize, are added to create the BDI composite index. In the weighting scheme, the typical number of ships in each vessel type and their relative daily shipping rates are used in the calculation of the BDI. The Handysize vessels are the smallest and cheapest, whilst the Capesize vessels are the biggest and most expensive. The weighted rates are then determined by the average number of ships in each category, taking into account the size of each of the four vessels. The BDI is measured in points, with higher values corresponding to higher shipping costs and, thus, increased demand for raw materials. Investors, analysts, and traders can use the index to track international trade activities and evaluate the state of the world economy. To provide a more complete view of economic trends, the BDI should be used in conjunction with other data sources since it is only one indication of economic activity. With these characteristics, the BDI has been considered an indicator that effectively represents shipping costs in addition to being taken as keeping and determining the pulse of trade and production activities around the Globe [
2].
Among the three indices subject to the study, the BDI index is a relatively less explored index in the literature though the BDI is directly linked to global economic conditions. Since the index is based on the movement of sea transportation, it is also an important indicator for crisis prediction. Nowadays, in addition to being an important source of information for shipping costs of commercial goods, the BDI is also a tool used in economic forecasts.
The BDI is a different index relative to the other indices subject to our research. The BDI is an important index based on shipment costs, and the index is highly linked to global economic demand and supply conditions in the goods markets. While the VIX is based on future options and contracts in financial markets, and the MSCI is directly linked to global stock market conditions, the BDI, as an index linked to the conditions in the international shipments’ costs, reacts differently to economic fluctuations. The BDI is under the influence of port and dock costs and placement fees, in addition to transportation costs and is also linked to fluctuations in oil prices. Though the BDI is dependent on these by its structure, these factors also make it exposed to global demand and supply of manufactured products. Due to various shipment routes and different cost structures, the BDI is formed as a daily collection of average prices of 23 major shipment routes worldwide. In this respect, the BDI is not influenced by economic indicators such as the unemployment rate, inflation and oil prices which could be affected or manipulated by states and speculators. However, the BDI is difficult to manipulate or influence. The conduct of the BDI being affected by factors such as supply and demand relationships can be shown as a reason for this. Another reason is that the number of ships available on earth that can affect or manipulate demands is limited. Building more ships will result in more costs [
3].
The BDI is a vehicle and an indicator with five different sub-indices for the investors [
3]. After the economic recession and during economic growth, the increase in demand for raw materials increases investments in this area, and consequently, growth in shipping volumes is seen [
4]. In general, when combined with changes in global demand for raw materials used in production [
5], inflexible vessel supply has strong effects on the conditional average of the BDI index and sharp fluctuations in its variant to cope with the short-term increase in demand due to its lack of flexibility [
6].
The Cboe Volatility Index (VIX), introduced by the Chicago Board Options Exchange (CBOE) in 1993, is the second index to be explored in the study. The VIX is a real-time index that indicates the market’s expectations for the relative strength of the S&P 500 Index’s near-term price fluctuations (SPX). The index generates a 30-day forward estimate of the implied value of option prices in the calculation. Therefore, the index represents implied volatility in the prices of SPX index options with close expiration dates. The VIX is frequently considered a good indicator of market emotions and expectations and, particularly, the level of fear among traders measured by the VIX volatility index. The level of volatility or risk rises in direct proportion to how dramatically the index’s price swings occur. In the mathematical calculation of the VIX index, the implied value of option prices is the probability that the price of a particular stock will increase or decrease by an amount sufficient to reach the strike price, also known as the exercise price, which determines the price of options. Only SPX options with expiration dates that are between 23 and 37 days are taken in the calculation of the index. Further, investors are known to utilize financial derivatives, including ETFs, VIX futures and options on fluctuations in the VIX. Without presenting the mathematical formula here, the formula theoretically operates as follows: the weighted values of numerous SPX puts and calls over a wide range of strike prices are combined to calculate the anticipated volatility of the S&P 500. All such qualified options must have legitimate, non-zero bid and ask prices that reflect the opinion of the market regarding which options’ strike prices will be reached by the underlying equities in the time left before expiration [
7].
The VIX has gained popularity among investors and has been modeled by empirical methods as a subject of a lot of research in literature. The VIX has been commonly known as the investor sentiment or fear index and reacts sharply to changes in trends, especially to high levels of market turmoil. Past experiences have shown that sharp inclines in the VIX effectively point to harsh drops in global financial market indices. An example is the 2008 sharp incline in the VIX before the 2009 Global Crisis. Furthermore, though the VIX is among the largest indices in terms of volume and liquidity among all other financially volatile indices, its efficiency in the predictability of financial crises is questioned in the literature [
8]. The third index, MSCI, is an index generally considered an important barometer or indicator of the health of the global economy in addition to being considered as representing the performance in the global stock markets. If the history of the MSCI is investigated, it is observed that the MSCI dates back to 1969, the introduction by Capital International of several stock indices to reflect international markets outside the United States, making it among the first global stock indices. Morgan Stanley then purchased the data license rights of Capital International in 1986 and started to use the MSCI acronym instead. As of 2007, the New York Stock Exchange (NYSE) has become an independent public company. Though the MSCI world index has been assumed to be composed of stocks from different countries since then, our investigation shows that the index, and the included sectors in the MSCI world index (MSCI), is, in fact, an index that constitutes around 85% of total capitalization of stock markets in the world with a major share belonging to USA stocks. Since the MSCI captures large and mid-cap representation across 23 stock markets of developed countries (DC), with 1508 constituents, the index covers approximately 85% of the free float-adjusted market capitalization. As a result, though the included stocks are not from all over the world, the MSCI index is widely accepted as a benchmark for international stock portfolio positions. In addition to the MSCI, Morgan Stanley publishes the MSCI All Country World Index (MSCI ACWI). Comparatively, while the MSCI’s coverage corresponds to 100% of indices from DC countries, the MSCI ACWI is 88% DC and 12% major emerging market stocks of the MSCI Emerging Market Index (MSCI EMI). In terms of market capitalization, large-cap stocks take the largest share in both the MSCI World and MSCI ACWI, and both indices have very close weighting schemes. The shares of sectors in both indices are very close in terms of share and the included sectors. In both indices, included sectors with the largest shares are financial, information tech, consumer durables, industrial and health care sectors with 14.6%, 20.8%, 11.6%, 9.5%, and 12.9% weights in the MSCI ACWI and 13.5%, 21.1%, 11.3%, 9.95%, and 14.02% for the MSCI World. Further, though indices are reported for the world, the majority of stocks are USA oriented. The share of shares of stocks from USA markets is 69.7% in the MSCI World and 61.9% in MSCI ACWI. A major difference is Chinese stocks. While no Asian stocks are included in the MSCI World, the share of Chinese stocks is 3.49% in MSCI ACWI. In both indices, stocks with large capitalizations are included to both indices. The mean cap of stocks is 30 and 18 billion US dollars for MSCI World and MSCI ACWI indices. A third index is the MSCI Emerging Markets Index. However, its coverage is largely for the Asian markets instead of the major emerging markets of the world. China’s share is 33.4%, Taiwan’s is 14.4%, in addition to South Korea with a share of 11.8%, India with 13% and Brazil with 5% only. Hence, the sectoral coverage of the index is largely high-tech and blue-chip stocks. Under these factors, all three indices are indices that are drivers of the world stock market fluctuations.
MSCI indices are considered indicators of investors’ behavior and expectations. Inclines in the MSCI index suggest that the investors are optimistic about the global economy and the good performance of the stock markets around the globe. Conversely, falls in the index are attributed to poor performance in the global economy and stock markets in addition to being linked to investors with negative expectations about future performances. Because all MSCI indices are reviewed every quarter and re-balanced on the basis of the number and reflectivity of shares twice a year, all ETF and investment funds follow these indices as a benchmark in addition to making similar and also identical changes to the shares of assets in their portfolio. For this reason, the MSCI index and related indices have significant influence and power in affecting and changing the financial market, and their influence on diverting the markets cannot be diminished. Moreover, it is common to observe that the recently added stocks to the index increase their share price, whereas the opposite, the removal of a stock from the index, generally reduces the share price. Therefore, it is also subject to discussions and criticism whether the index is open to manipulation or not.
This study also hypothesized that the BDI MSCI and VIX indices are pioneering real and financial indicators, respectively, with good predictive performances. However, the predictability and forecastability of the three indices are challenging. For the reasons stated, the sample covers the period of 1 November 2007–30 May 2022, which includes the Global Crisis and, recently, the COVID-19 economic shutdown. The three indices are explored with a family of non-linear volatility models, which aim at their regime-dependent asymmetric behavior in conditional volatilities in these indices that make the predictive power of econometric models for these series challenging. As will be seen, all series are subject to nonlinear and leptokurtic distribution with heavy tails with smooth changes between distinct regimes. As a result, the BDI, VIX and MSCI indicators are modeled with smooth logistic transition autoregressive—generalized autoregressive conditional heteroskedasticity (LSTAR-GARCH) and its generalization to nonlinear LSTAR type asymmetric power GARCH (LSTAR-APGARCH) model which are shown to produce effective results in forecasting relative to their single-regime variants. With this respect, we propose that to utilize the three indices as leading economic indicators; they should be efficiently forecasted in order to achieve effective investment strategies along with achieving effective economic policy recommendations. Therefore, the paper focuses on the examination of whether the BDI, MSCI and VIX indices retain forward-looking properties. The paper has three central contributions. First is the modeling of the BDI, VIX and MSCI indices concurrently for comparative purposes by integrating the nonlinear characteristics of these indices. The second is to recommend a hybrid modeling method that benefits from the logistic transition auto-regressive model (LSTAR) and the STAR-GARCH model proposed in the literature to be further extended to the asymmetric power GARCH model (APGARCH) to achieve the LSTAR-APGARCH model. The third is in terms of sample size. The sample covers one of the largest in the literature, especially for the BDI, which corresponds to 1 November 2007–30 May 2022 with daily series. The sample has the potential to provide important information given the existence of the recent COVID-19 and 2008 Great Recession periods in addition to recessionary and expansionary business cycle periods. As a result, the period is expected to cover quite diverse data-generating processes, possibly to be modeled with nonlinear regime-dependent models instead of traditional linear time series approaches.
The paper is organized as follows. The literature review is given in Part 2. The data, unit root tests, nonlinearity tests and chaos tests are in Part 3. Methods are presented in Part 4, where the econometric methodology of LSTAR-GARCH and LSTAR-APGARCH models are presented. The results section is given in Part 5, which also includes a discussion section. The conclusion is given in the last section.
2. Literature Review
If an overlook to the empirical papers with econometric approaches is evaluated, among the indices to be analyzed in this study, papers focusing on modeling the BDI is rather limited. Especially lack of studies focusing on the utilization of daily BDI is noteworthy. In contrast, the body of research with empirical models of forecasting the VIX and MSCI is relatively less scarce. In the first stage of this section selected research focusing on the BDI will be provided for their distinguished contributions.
By studying the predictability of the futures of the dry load market on the main roads and by utilizing monthly changes in the BDI, Chen et al. showed that spot prices and trade routes do not have long-term relationships with the size of all three shipment routes [
9]. Leonov and Nikolov used wavelength and nonlinear artificial neural network models to estimate dry load transport fluctuations and transport rates of the Baltic Panamax 2A route and Baltic Panamax 3A route [
10]. Duru et al. exploited the Fuzzy Integrated Logical Forecasting model with an integrated logic estimation model to investigate the performance of the effect of dry cargo transport’s limited vessel rental contract rates over time [
11]. Chen and Wang explored the freight market indices, and their findings with the EGARCH model underlined leverage effects in their volatility, and their findings emphasized the asymmetric volatility response of freight indices amid negative and positive shocks [
12]. Geman and Smith explored different shipping markets in financial literature and put forth the rates of the freight that drive the indicator and their relationship with the BDI in the global economy [
4]. For the BDI, Geman and Smith’s findings emphasized the existence of structural breaks between 1988 and 2003 in mid-2003 with their Constant Elasticity of Variance model; therefore, their findings underlined the necessity of utilization of models before and after the break periods separately [
4]. Oomen’s study reviewed the impact of BDI returns on the futures trading returns and the potential impact of the BDI on the stock market indices of 11 developed countries [
3]. The study was pivotal in displaying the association between the BDI and stock market returns as well as showing the BDI as an important predictor of these markets [
3]. An important finding of Oomen that we should emphasize is the determination of the predictive power of the BDI in modeling financial markets in the world [
3]. Further, the findings also underlined the existence of structural changes that hamper the efficiency of the models, and Oomen’s suggested division of the sample into four sub-periods proved significant results [
3].
The existence of nonlinearity in the BDI is a central characteristic of the BDI. The index is directly related to shipment costs that are under the influence of business cycles in the world. Various studies have hinted at the nonlinear behavior in the BDI. Among these, Papaillias and Thomakos stated that the annual growth of the BDI is subject to cyclical behavior, which became disturbed during the 2007 crisis, and the cyclical component of the BDI series could be used for forecasting with the use of spectral analysis and trigonometric regressions [
13]. The findings of [
13], in this respect, pointed at not only nonlinearity but also structural changes in time in the distributional properties of the BDI. Lin and Sim’s research investigated the effects of international trade on GDP growth in less developed nations, and showed that the BDI is an effective measure of the costs of trade transport [
2]. To investigate the volatility of BDI, ref. [
14] proposed a method based on the EMD approach, which parsed the original freight price series into several independent internal modes to improve the BDI prediction of neural networks models. The proposed method undertakes nonlinearity in these series more effectively than the traditional linear econometric approaches [
14]. Further, [
6] debated the challenges of the BDI’s predictability due to nonlinearity in addition to putting forth the ability of the BDI to model economic growth. Regarding nonlinearity, ref. [
8] put forth the threshold effects in the daily BDI and VIX, and their findings underline regime-specific mean and volatility behavior in daily BDI price changes with TAR-TR-GARCH models. To this end, ref. [
8] underlined the necessity to augment the GARCH models with TAR-TR-GARCH models to capture nonlinearity in the BDI.
A set of literature has focused on multivariate models to investigate the relations of the BDI with various variables. The traditional ARIMA time series, fractional ARIMA and various machine learning methods are evaluated in forecasting BDI and it is shown that no significant improvement could be achieved and results improve only if model forecasts are combined [
15]. Ruan et al. investigated cross-correlations between the BDI and crude oil spot prices using the cross-correlation analysis (MF-DCCA) model, which benefits from multifractal trending [
16]. Through the MF-DCCA, the cross-correlation between each pair of BDI and crude oil prices is demonstrated, which is strong in the short term, but the cross-correlations in the long term are weakly permanent [
16]. By utilizing Markov-Switching VAR, the nonlinear relations between the BDI, economic growth and gold prices are determined in addition to providing the BDI as an important tool to model economic series and precious metals [
17]. The impact of the BDI on commodity futures and foreign exchange and stock markets was investigated using the BEKK-GARCH-X model for a sample covering 1 October 2007 to 31 October 2018, and findings confirmed spillovers between the BDI and financial markets [
18]. Kamal et al. examine short and long-term forecasts of the BDI with the proposed DERN deep learning model [
19]. They showed that forecasts improve with the BDI, and modeling nonlinearity with neural networks provided tools to overcome uncertainty for market participants and shipowners [
19]. Their findings also underlined the need to utilize nonlinear models for short and long-term maritime business decisions and to avoid market risk [
19]. The decline in the trend of the BDI before an economic crisis is an important phenomenon [
20]. The decline in the trend of the BDI became evident during the 2008 global crises and afterward, the commodity and stock markets took a hard hit [
20]. The study also emphasized existence of a similar pattern before and during COVID-19, which underlined the predictive power of the BDI as a leading indicator for markets [
20]. Yang et al. [
21] proposed the necessity of nonlinear modeling of the Baltic Panama Index with support vector machines and wavelet transformations, the former for augmenting forecasts and the latter for the denoisification of data. Han et al. examined the predictive potential of BDI and various sub-indices in forecasting exchange rates, and they emphasized the efficiency of the BDI in in-sample and out-of-sample forecasts in modeling economic variables [
22]. Zhang et al. inspected the BDI’s forecasting potential with various nonlinear models [
23]. They emphasized the importance of neural network models in capturing nonlinearity in addition to suggesting a combination of methods, i.e., dynamic fluctuation nets and AI, to increase the effectiveness in modeling and forecasting the BDI [
23]. Chen et al. investigated the effects of BDI between crude oil and commodity markets with the Copula-VAR-BEKK-GARCH-X model [
24]. Accordingly, dependence and spillover behavior between the BDI, iron, and brent oil are time-varying. Hence the results revealed the necessity of nonlinearity modeling [
24].
A recent set of studies have associated freight transport costs with COVID-19. Michail and Melas explored COVID-19’s impacts on the sea transport industry. Their findings emphasized that freight prices are highly affected caused from conditions for dry load, clean and dirty transfer tankers through VAR models and GARCH models [
25]. The findings of [
25] further confirmed pandemics’ effects on the freight industry and information obtained through volatility modeling could be used as a hedge for risk in stock markets. Ref. [
26] showed that COVID-19 had strong effects on the container transport operations of firms in addition to effects on their financial structure and employment in the sector. Ref. [
27] investigated the effects of COVID-19 on the world stock market, proxied with the World-MSCI, world energy market, proxied with the MSCI-Energy index and freight costs measured with the BDI, by using structural VAR models for 21 January 2020–26 February 2021 period. Though their findings revealed negative effects of COVID-19 on all three variables, among these, the MSCI and MSCI-Energy were most badly hit by COVID-19 relative to a lessened effect on the BDI [
27]. However, the lessened effect could also be due to the unexpectedly sharp decline in oil prices during the COVID-19 shutdown, which counteracted the sharp inclines in shipment costs [
27]. Chen et al. showed the significance of spillovers of volatility between oil, iron-ore prices and BDI during COVID-19 [
24]. Ref. [
28] explored the effects of COVID-19 on the S&P index in the US in addition to global economic activity proxied by the BDI, and the spread among 10-year treasury and federal funds rates with SVAR models. His findings underlined the negative effects of COVID-19 [
28]. The lockdown policies had strong effects on shipping during COVID-19 [
29]. Zhao et al. showed that the global dry carriage was largely affected by quarantine policies during the second month of COVID-19 when the BDI dropped about 35.5% compared with the same month a year before [
29]. Ref. [
30] stress the impact of COVID-19 on BDI volatility during COVID-19. With GARCH-MIDAS models, their findings signified that as the number of infected inclined, so did the BDI volatility in addition to similar results obtained for sub-variables [
30]. The infections also influenced the volatility of shipment cost variables such as crude oil price, freight rate, container idle rate, global port calls, and port congestion levels [
30]. Hence COVID-19 infections significantly influenced the volatility of the BDI and other shipment variables [
30]. In order to achieve the goal of cost or risk control, the shipping industry is suggested to benefit from the development of the epidemic in numerous nations as a reference in addition to the utilization of models aiming at modeling the trend of and volatility of the BDI, which is found to be under the influence of changes in the confirmed Covid cases [
30]. Jeris and Nath, by using the wavelet approach and daily data for 21 January 2020–30 October 2020, showed that COVID-19 inclined volatility of the BDI in addition to its similar effects on economic policy uncertainty, crude oil price volatility and banking stock market volatility in the USA [
31].
The VIX is among one of the indices studied extensively in the literature. Copeland and Copeland (1999) explored the impacts of the VIX on large and small-capitalization stocks and inclines in the daily VIX exhibits effects in the following days on portfolios containing large-capitalization stocks and the VIX was shown to help on investor decisions once used effectively in their portfolios [
32]. Ref. [
33] debated the return and volatility relationship in the implied volatility indices. According to [
33], sharp inclines in the VIX point at oversold markets and future returns are positive (negative) after the periods with very high (low) levels of VIX. Ref. [
34] advocated TAR-type threshold responses of the S&P stock market to low investment fear and high investment fear regimes measured with the VIX index. These results confirmed that the best performance is achieved if the VIX determines the threshold variable compared to the other four candidate indicators [
34]. By using a sample covering 1995–2022, a recent study [
35] demonstrated that the VIX successfully estimates future stock market volatility in addition to the VIX being confirmed as a tool for portfolio insurance prices [
7]. Consequently, due to the fact that the S&P 500 index options market is affected by hedgers that buy index put options, risk inclines hence, leads to investor sentiment increases with anticipations of a potential drop in the market [
7]. By the use of linear VAR models, Ref. [
36] investigated the VIX and stock market ETF relation and emphasized the negative effects of VIX shocks on S&P stock returns in addition to European ETF markets. Nevertheless, in addition to putting forth the correlation between the markets and VIX, in the spirit of the seminal studies of [
32,
33], high and rising levels of VIX were linked with a future rebound in the stock markets for the investors [
36]. The lead-lag relation between the stock markets and VIX was noted to co-move significantly [
36]. Further, if the VIX is in a high-volatility regime, the negative impact of the VIX could be comparatively more severe [
33]. One of the early signs of such a finding was Guo and Wohar, who employed econometric techniques that put forth the mean shifts and, therefore, nonlinearity in the selected implied volatility indices in addition to shifts of the standard deviation of the VXO and VIX indices [
37]. With this respect, their identified subperiods due to shifts are pre-1992, 1992–1997, and post-1997 periods [
37].
Regarding the market risk and VIX relationship, Durand et al. (2011) showed that changes in the VIX affect market risk premium. An additional factor to the VIX is shown as the value premium by the use of the three-factor model within the Fama and French factors approach [
38]. The VIX is also calculated and examined for India [
39], and the study displayed that the VIX is a gauge of investor fear when markets are declining and when stock market volatility is rising [
39]. Furthermore, the expected volatility is influenced by the actual return volatility signifying a direct association between VIX inclines and market turmoil [
39].
A selected set of literature will be highlighted to evaluate the links between the MSCI and the world stock market performance in addition to forecast performance for the MSCI. Recent empirical research [
40] on the MSCI world stock index utilized VAR and EGARCH models and demonstrated the efficiency of the MSCI in capturing global economic crisis by employing daily dataset covering 3 June 2002–22 March 2013 [
40]. Further, Ref. [
41] applied four asymmetric-GARCH models, including NA-GARCH, TGARCH, GJR-GARCH, and AV-GARCH, to model and forecast the MSCI’s volatility. Within a risk exposure framework in the context of value-at-risk, their findings emphasized nonlinearity and asymmetry in the MSCI, in addition to confirming these factors improving the performances of GARCH models in forecasting MSCI [
41]. The MSCI and MSCI-Energy sector indices were modeled for the effects of COVID-19 on these series in addition to the BDI as another leading indicator [
27]. By using SVAR models, their study confirmed the negative impacts of COVID-19 on investigated leading indicators [
27]. The literature above indicates that the modeling of the MSCI, in addition to the BDI and VIX, also necessitates taking nonlinearity and volatility in addition to leverage effects into consideration.
Though considered nonlinear in nature, the traditional Autoregressive Conditional Heteroskedasticity (ARCH) model [
42] and Generalized ARCH (GARCH) model [
43] do not take regime-dependency into account. For this purpose, the models in our study benefit from STAR-GARCH models [
44] and their extensions. The STAR-GARCH model assumes STAR-type nonlinearity in the conditional mean process, and the conditional variance follows the GARCH process and under outliers, robustness of QMLE conditions were shown by [
44]. Our study follows the LSTAR-GARCH and LSTAR-APGARCH modeling of [
45]. STAR-GARCH models are employed to investigateSTAR-type nonlinearity in Bank of Finland’s Banking and Finance Indices [
46]. Various studies of Bildirici and Ersin highlight nonlinear volatility in crude oil prices, which is also an important cost factor for BDI dry cargo freight shipment costs index. Among these, Ref. [
47] showed that modeling and forecasting crude oil prices requires LSTAR-type nonlinearity in the conditional mean and variance processes. The findings of Bildirici and Ersin in [
48] confirmed the integration of fractional integration into the LSTAR-LST-GARCH family of models, including LSTAR-LST-FIGARCH and LSTAR-LST-FIAPGARCH, which benefit from STAR-type nonlinearity [
47]. For forecasting gains in future crude oil prices, Ref. [
49] proposed a new set of STAR-GARCH models with exponential AR components, and the results confirmed better forecasting capabilities of the STAR-GARCH, LSTAR-GARCH and ESTAR-GARCH models over single-regime GARCH models in modeling and forecasting exchange rates of the British pound and Botswana Pula against the US dollar [
49]. Bildirici et al. [
8] examined the regime-dependent volatility of the daily VIX and BDI with TAR-TR-GARCH models. LSTAR-GARCH and LSTAR-APGARCH models will be evaluated in the next section.
5. Results
In this section, the BDI, VIX and MSCI indices are modeled with nonlinear LSTAR-GARCH and LSTAR-APGARCH models. As shown in
Section 4, both models utilize logistic smooth transition functions to capture transition dynamics between sub-regression spaces in which
,
and
series follow regime-dependent and asymmetric conditional volatility processes.
5.1. STAR-GARCH Type Nonlinearity Tests and Model Selection
As an initial step, following the methodology, all three series were modeled with single-regime GARCH and APGARCH models. The estimated models were tested for STAR-type nonlinearity following [
51,
52]’s remaining STAR-type nonlinearity tests. The STAR-type nonlinearity testing was based on Taylor expansions and also helped distinguish between logistic and exponential transition functions in addition to testing linearity (For a discussion for multiple STAR modeling that allows more than two regimes, readers are referred to [
54]. The testing process necessitated testing the single-regime model under the null to be tested for a two-regime model alternative. In case of rejecting the null hypothesis of the single-regime model, a second step involved testing the remaining STAR-type nonlinearity with two regimes under the null to be tested against a three-regime STAR-type nonlinearity which assumes two logistic transition functions. Therefore, the above-mentioned linearity testing methodology not only allowed for testing single-regime GARCH against two-regime STAR-GARCH models, but also allowed the determination of the number of regimes in the STAR-GARCH model. To save space, the test results are not reported but are available upon request. The overall test results favored two-regime LSTAR-GARCH and two-regime LSTAR-APGARCH against their single-regime variants). The test results are given in
Table 5.
Due to model specification, all tests assumed the transition variable as εt−1. Accordingly, F test results favored nonlinearity for all modeled series, and F4, F3, and F2 tests resulted in the selection of logistic transition functions for all models at a 5% significance level. Results also confirmed STAR-GARCH and STAR-APGARCH models against their single-regime counterparts.
In the previous section, in the determination of orders of stationarity, KPSS and KSS tests were conducted to obtain results robust to certain forms of nonlinearity [
59,
60], The existence of nonlinearity of series were further confirmed with the BDS test [
61]. At this section, single-regime GARCH and APGARCH models are tested against their two-regime nonlinear LSTAR-GARCH variants with Lagrange-Multiplier tests (LR) under nuisance parameters [
62]. After estimation, LR tests confirmed the acceptance of two-regime variants in all single-regime models. We also estimated three regime models however the estimation results led to inconsistent results if more than one transition function. Ref. [
44] discussed testing for and modeling of STAR-GARCH models. In the next stage, the models were estimated with QMLE (Quasi-Maximum Likelihood Estimator). QMLE is susceptible to the initial values. According to [
63], the statistical and structural characteristics of STAR-GARCH models generally limit these models to two regimes, though m-regime extensions are possible [
63]. In addition to two-regime models, we estimated three-regime models. However, these led to inconsistent results. As a result, we maintained two-regime models similar to the general tendency in the literature.
5.2. Model Estimation Results
The model estimation results are reported in
Table 6 in three subsections; (a), (b) and (c), which correspond to parameter estimates and diagnostic testsfor BDI, VIX and MSCI daily % returns, respectively. In all estimated models, the stability condition for GARCH models, the sum of ARCH and GARCH coefficients being less than 1, is satisfied. The calculated
t-test statistics show that the majority of parameters are statistically significant with minor exceptions. Among these, interesting results have been found for the BDI and MSCI indices.
5.2.1. Comparison of Regime-Dependent ARCH and GARCH Parameters
There were major changes in the coefficient of ARCH and GARCH between regimes 1 and 2 for the BDI and MSCI. The GARCH coefficients were higher than the ARCH parameters as expected for all regimes and models, with an exception.
For the BDI, the GARCH coefficient for (regime) 1 for the LSTAR-GARCH was too low, estimated as 0.005, close to zero. In the second regime, the coefficient was calculated as 0.547. For the LSTAR-APGARCH, this problem was solved, with GARCH parameters estimated as 0.742 and 0.508 in the respectful regimes. For the MSCI, the GARCH parameter in regime 1 was very close to 0.954, while the second regime had a value of 0.151. ARCH coefficients were inevitably the opposite for these two regimes. For the LSTAR-APGARCH estimated for the MSCI, this problem was solved again since the GARCH coefficients was determined as 0.858 and 0.821 for regimes 1 and 2. The results show that the LSTAR-APGARCH models are better in capturing the expectations for the ARCH and GARCH parameters with these respects.
5.2.2. Comparisons of Asymmetric Power Effects
The asymmetric power terms were estimated as being regime-dependent in all LSTAR-APGARCH models. For the BDI, the asymmetric power parameter was estimated as 0.595 for regime 1 and as 0.683 for regime 2 in the LSTAR-APGARCH model. The LogL statistic was highest in regime 2 for the BDI’s LSTAR-APGARCH model. For the BDI, the LogL values for the LSTAR-GARCH model were relatively closer and in addition, if evaluated together, the fit of LSTAR-GARCH was higher for the BDI in this respect. However, one should also keep in mind that if regime 2 is attained, the overall fit is the highest for the BDI series modeled with LSTAR-APGARCH. LSTAR-APGARCH in regime 2 is more meaningful in modeling the BDI data than the LSTAR-GARCH for the BDI however, if regime 1 is attained, the fit of the model would be worsened relatively.
5.2.3. Evaluation of In-Sample Fit
When the AIC and SIC values were evaluated, similar results were obtained: in both models, fit was relatively higher in regime 2. Interestingly, the APARCH (Gamma1) parameter was statistically no different than zero. This parameter is similar to the leverage parameter of the traditional GJR model and the parameter introduces asymmetry between negative and positive shocks. As a result, in regime 1 of the LSTAR-APGARCH model, the model uses this leverage effect, therefore, reduces to LSTAR-PGARCH model.
Though fit of all models is evaluated in this section, note that these results are for in-sample fit of the models only. The out-of-sample forecast evaluation will be crucial, for which the results will be given in the following section. If an overlook is presented, since the total sample covers 1 November 2007–30 May 2022, the out-of-sample forecast results are obtained for the 10 days ahead, corresponding to the last two weeks of the dataset since the data is for working days. As a result, we decided not to reject models in terms of in-sample results, and we maintained the models to be analyzed for their out-of-sample forecast performances in the next section. So far, the models provided interesting characteristics which distinguish the volatility dynamics in two different regimes for all models estimated for the BDI daily % returns.
In section (b) of
Table 6, the models estimated for the VIX led to important findings. The stability condition is satisfied for both models for the VIX. For in-sample fit, the use of LSTAR-GARCH models for the VIX is clearer. The AIC and SIC and in addition LogL statistics show better fit for the LSTAR-GARCH relative to the LSTAR-APGARCH if two regimes are evaluated for both models. However, again it should be stated that out-of-sample performances would provide decisive information with this respect. The ARCH term in regime 1 of LSTAR-APGARCH is statistically significant at 10% significance level. However, the parameter is significant at 5% in regime 2 in addition to significant GARCH parameters in both regimes.
5.2.4. APARCH Parameters: Comparative Analysis
Similar to the BDI’s LSTAR-APGARCH model, the APARCH(Gamma1) parameter is not significant. Therefore, it is zero in regime 1. As a result, the process that the conditional variance follows in regime 1 reduces to an LSTAR-PGARCH model since the leverage effect is insignificant for this regime due to the insignificant gamma parameter. This interpretation for both the BDI and VIX should be made with caution since the model is still asymmetric as a whole, especially the power terms being statistically different in each regime. As a typical example, the APARCH(delta) is estimated as 0.93 in regime 1 and as 1.305 in regime 2, both significant, leading to higher levels of power terms and more drastic levels of volatility in regime 2. Further, the Q statistics and ARCH tests in diagnostics suggest no autocorrelation and no remaining ARCH effects.
The results for the MSCI index daily % returns are reported in section (c) of
Table 6. For the MSCI, the use of LSTAR-APGARCH models in both regimes is more meaningful considering the AIC and SIC results in addition to LogL statistics. With asymmetric power terms, the explanatory power increases for the MSCI. The stability condition is satisfied for both models. However, the ARCH term is statistically insignificant in regime 1 of the LSTAR-APGARCH model.
Further, similar to the model for the VIX, the APARCH(Gamma1) is statistically insignificant in regime 1. Therefore, the process in this regime becomes an LSTAR-PGARCH. However, again, the asymmetric power terms are still statistically different than zero, and for both regimes of the LSTAR-APGARCH model, the power terms are regime-dependent. Note that though certain minor insignificant results occur for the model parameters, it should also be noted that the models are selected over their single-regime counterparts, and the out-of-sample forecast performances will provide important insights in the next section.
5.3. Out-of-Sample Forecast Results
As noted in the previous section, the daily sample covers 1 November 2007–30 May 2022. The last 10 days are not utilized in the estimation of models and are kept out for forecast evaluation. As a result, this section aims at providing a comparative analysis on forecast performances of the GARCH, APGARCH, LSTAR-GARCH and LSTAR-APGARCH models for 10-working-days-ahead, which corresponds to a two-week period of future forecasts. The models are evaluated for their out-of-sample performances and the results are reported in
Table 7 where the median squared error (MedSE) statistics are reported. The MedSE criteria is selected due to its non-sampling estimation capabilities and the statistic is robust to outliers.
If the single-regime GARCH and APGARCH models given in columns 1 and 2 are evaluated, the APGARCH model is observed to lead to lower MedSE, suggesting better forecast performance for all series. If their nonlinear counterparts are evaluated, there is a significant improvement in forecast performance in terms of the drastic decline in MedSE statistics for LSTAR-GARCH and LSTAR-APGARCH models. The results are obtained by assuming that the conditional volatility process is either in regime 1 or regime 2 instead of reporting a single MedSE for the model. With this approach, we stress that uncertainty is likely, and the analyzed variable could shift to another regime during the forecast period due to an external shock. Though the model assumes a smooth transition between regimes, such a forecast practice could lead to important insights. The overall result shows that, for the LSTAR-GARCH model, regime 2 led to the lowest MedSE in all regimes for all variables forecasted. In addition, whatever the regime is, both regimes of the LSTAR-GARCH perform better than the GARCH and APGARCH, the single-regime models. It is observed that, in the first regime, the MedSE value for the LSTAR-GARCH model for the BDI is 0.00007, and this value becomes 0.00000004083 in regime 2 for the LSTAR-GARCH model. The difference between the MedSE in regimes 1 and 2 should be interpreted as the difference in the volatility dynamics. A similar reduction in MedSE is also observed for all estimated LSTAR-GARCH models including those for the VIX and MSCI. Hence, the results confirm higher uncertainty in the first regimes of the LSTAR-GARCH models. However, even in the high-forecast-uncertainty regime, the MedSE statistic is still lower for the LSTAR-GARCH compared to the single-regime models, signifying forecast performance for the LSTAR-GARCH specification.
For the LSTAR-APGARCH models, forecast improvement is achieved for all series forecasted over the single-regime GARCH and APGARCH models. However, if regime results are evaluated, differences should be taken into consideration. For the BDI, both regimes provided similar forecast performance since MedSE = 0.00382 in both regimes of the LSTAR-APGARCH. The MedSE statistics are lower in both regimes of LSTAR-GARCH for the BDI. Therefore, LSTAR-GARCH performs better in forecasting compared to LSTAR-APGARCH for this series. For the MSCI, a similar finding holds. The LSTAR-APGARCH model provided better forecast accuracy for the MSCI in both regimes relative to single-regime GARCH and APGARCH. The MedSE statistics are calculated as 0.00000383 and 0.00000787 for regimes 1 and 2. Though the model provides an improvement, it falls behind if compared to the LSTAR-GARCH specification for which the MedSE statistics are strikingly less, 0.00000000852 and 0.000000000164, hence, for the MSCI, LSTAR-GARCH provided best forecast performance over all models analyzed. For the VIX, on the other hand, it is clear that the LSTAR-APGARCH model performs best in regime 2 and regime 1. However, looking closely, LSTAR-GARCH regime 2 performs similar to LSTAR-APGARCH regime 1. For the VIX, it is clear that LSTAR-APGARCH specification leads to significant performance improvement for regime 2. Therefore, the overall forecast performances show that both LSTAR-GARCH and LSTAR-APGARCH models provide important forecast improvement over the traditional GARCH and the asymmetric power upgraded APGARCH model. Further results also suggest that regime-dependence is an important phenomenon in forecast practices. Therefore, following the findings given above, policymakers interested in utilizing the BDI, VIX and MSCI indices as economic indicators should keep their regime-dependent characteristics and most importantly, ignoring nonlinearity would lead to drastic reduction in forecast performances.
5.4. Discussion
The study reveals the significant consequences of the predecessor models. The empiric results obtained in the study have important outcomes. The key results are summarized in
Table 8.
The following results are obtained through the models estimated regarding the BDI, VIX and MSCI indices. i. The nonlinearity and regime-dependency in the conditional variance of the daily BDI, MSCI and VIX returns cannot be rejected. ii. To model the BDI, MSCI, and VIX indices, which are the leading indicator candidates examined in the study, regime-dependent leverage effects must be included in the conditional variance regimes. iii. In addition to the regime effects captured by LSTAR-GARCH parameters, modeling regime changes governed with smooth changes provide significant improvement in the accuracy of in-sample predictions for the VIX, MSCI and BDI series investigated. Given that all three series are considered as leading economic indicators, investors and policymakers should consider the utilization of nonlinearity in forecasting practices. iv. Important policy recommendations are generated from the findings of the study. It is necessary to improve prediction and forecast accuracy by implying smooth transition into single-regime models in order to consider the VIX, BDI and MSCI as leading economic indicators. Otherwise, due to nonlinear volatility dynamics, the traditional GARCH and APGARCH models suffer in producing effectiveness in forecasts. Given the nature of the series analyzed, careful investigation of the distributional properties of these series are necessary. Investors and policymakers should consider the asymmetry, nonlinearity and regime-dependent asymmetric power dynamics to hinder success in assessing future economic conditions.
In these respects, in addition to the VIX and MSCI economic indicators variables, the volatility of the Baltic Dry Load Index (BDI), is shown be an important economic indicator that can be considered a leading indicator of financial conditions. Moreover, changes in the calculation methods made from time to time in financial indicators such as the VIX or MSCI, which have volatility in their nature, or changing the weights of the criteria that reweighted the index, such as a certain share or stock or sectors, are sometimes the subject of discussions in the literature. To this end, the forecasting and modeling performances of the VIX and MSCI indices are criticized as not reflecting reality. As a matter of fact, a set of research studies criticize the performance of the MSCI index in terms of its failure to predict the 2008 Global Financial Crisis in the USA.
It is very difficult to manipulate the BDI because the BDI is an index representing shipping costs with more than 20 shipping routes worldwide, and the index provides a measure of raw material costs in the world. The BDI is considered an economic indicator that is independent of political influences and directly impacted by supply/demand economic conditions. aims to explore financial crises with modeling techniques. The fact is that investors and policymakers are interested in practical economic indicators that they can use in policy and investment decisions. The Baltic Dry Load Index (BDI) is considered a candidate index, which can be taken by many researchers as an economic indicator on a global scale. The BDI aims to provide an index based on the cost of raw materials such as iron, coal, cement, and grain worldwide. In addition, the BDI Index includes both the volatility inherited by the rigidity of crude oil prices and the docking fees, such as the port and load download and fill. Therefore, given the cost of the materials used in production and the demand conditions that encourage manufacturers to change production decisions, ship berthing, loading, filling charges and freight volumes are largely affected by global supply and demand conditions. As a result, the BDI is highly responsive to the global demand-supply conditions for raw materials used for manufactured products, in addition to the cost of oil prices and freight prices affected by the volatility. However, although the BDI is an important indicator, leptokurtic distribution and its nature creates difficulties in modeling and predicting the index as a result of the nonlinear state of the BDI time series.
While the literature focused on modeling, estimating and estimating the BDI, VIX and MSCI indices is limited, studies aimed at modeling these variables with various econometric techniques are still in progress, and studies are increasing. With the accuracy motivation to explain the above-mentioned financial conditions, these series are important to review the asymmetric power (APGARCH) amplified conditions of the extended conditional mean and conditional variance with the non-linear threshold autoregressive model (TAR), Smooth Transition Autoregressive (STAR) model family, Logistic Smooth Transition Autoregressive (LSTAR) and the GARCH model family.
6. Conclusions
BDI, MSCI and VIX predictive variable data were analyzed for the periods of 1 November 2007 through 30 May 2022 with LSTAR-GARCH and LSTAR-APGARCH models. According to the results of the models, investors and policy makers should consider the threshold effects and nonlinearity of the MSCI, VIX and BDI series in assessing the three series as pioneering indicators of future economic activity. In addition, as the literature research revealed, both series were successful in the global recession forecast, which began in late 2008. Efficiency in the estimation of the off-sampling path followed by the two series is significant at this time. When assessed in terms of the compatibility of the predicted models with the series, nonlinear models with threshold characteristics have provided high performance for improved in-sample modeling. In addition, it has been observed that the family of STAR type nonlinear volatility models (LSTAR-GARCH and LSTAR-APGARCH) provides significant gains in sampling accuracy, significantly improving the model’s descriptive effectiveness. In addition, given the results of other applied models, the generalized GARCH, LSTAR-GARCH and LSTAR-APGARCH models have been compared with their predictive accuracy for 10 business days in different conditional volatility models with APGARCH architecture, and the out-of-sample estimation of LSTAR-GARCH and LSTAR-APGARCH models has proven to be extremely strong.
Given the successful estimation capabilities of the proposed models, asymmetric power cannot be disallowed by the power of incorporating regime transitions and threshold effects into the analysis at the same time. Therefore, in a non-linear threshold-type environment, when the model is extended with a smooth transient asymmetric force for each regime, it is proof that the prediction challenges in the MSCI and BDI can be overcome and that the MSCI, VIX and BDI indices can be used as an economic indicator. For the BDI, it is recommended to use predictive model with LSTAR-APGARCH method for all fluctuations, where predictive modeling results are meaningful; The MSCI and VIX financial pioneering indicators, which are discussed as more open to manipulations due to the changing calculation method in the literature, have been observed to be interpreted with LSTAR type models, especially in wavy periods, as they cause deviations on basic GARCH and APGARCH models.
As a result, the volatility and non-linearity of the BDI, VIX and MSCI indices must be kept in mind in order to make them more important today, such as the VIX and MSCI indices, and to be used as leading indicators and even considered as benchmarks.
It is necessary to select the disposing models. In addition, in future studies to achieve predictive improvements, focusing on models that aim to model different types of nonlinear GARCH models and non-linear series that allow not only conditional averages to change but also conditional variance processes between different regimes is recommended for policymakers and investors. The BDI is a leading economic indicator and is becoming a useful and important predictive for economic crises. Because the BDI’s volatility is susceptible to volatility in the world’s economies and the BDI Index has the ability to reflect the volatility of the financial conditions of the global economy with more than 20 shipments worldwide, our analysis results suggest that it can predict the direction of the global economy for investors and policymakers in advance.