Next Article in Journal
Economic Fraud and Associated Risks: An Integrated Bibliometric Analysis Approach
Previous Article in Journal
The Impact of Firm Risk and the COVID-19 Crisis on Working Capital Management Strategies: Evidence from a Market Affected by Economic Uncertainty
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation and Prediction of Commodity Returns Using Long Memory Volatility Models

by
Kisswell Basira
1,*,
Lawrence Dhliwayo
1,
Knowledge Chinhamu
2,
Retius Chifurira
2 and
Florence Matarise
1
1
Department of Mathematics and Computational Sciences, University of Zimbabwe, Harare P.O. Box MP167, Zimbabwe
2
School of Mathematics, Statistics and Computer Science, University of KwaZulu Natal, Durban 3630, South Africa
*
Author to whom correspondence should be addressed.
Submission received: 19 February 2024 / Revised: 26 March 2024 / Accepted: 27 March 2024 / Published: 23 April 2024

Abstract

:
Modelling the volatility of commodity prices and creating more reliable models for estimating and forecasting commodity price returns are crucial. The body of research on statistical models that can fully reflect the empirical characteristics of commodity price returns is lacking. The main aim of this research was to develop a modelling framework that could be used to accurately estimate and forecast commodity price returns by combining long memory models with heavy-tailed distributions. This study employed dual hybrid long-memory generalised autoregressive conditionally heteroscedasticity (GARCH) models with heavy-tailed innovations, namely, the Student-t distribution (StD), skewed-Student-t distribution (SStD), and the generalised error distribution (GED). Based on the smallest forecasting metrics values for mean absolute error (MAE) and mean squared error (MSE) values, the best performing LM-GARCH-type model for lithium is the ARFIMA (1, o , 1)-FIAPARCH (1, ξ , 1) with normal innovations. For tobacco, the best model is ARFIMA (1, o , 1)-FIGARCH (1, ξ , 1) with SStD innovations. The robust performing model for gold is the ARFIMA (1, o , 1)-FIGARCH (1, ξ , 1)-GED model. The best performing forecasting model for crude oil and cotton returns are the F I A P A R C H   1 , ξ ,   1 S S t D model and H Y G A R C H   1 , ξ ,   1 S t D model, respectively. The results obtained from this study would be beneficial to those concerned with financial market modelling techniques, such as derivative pricing, risk management, asset allocation, and valuation.

1. Introduction

When examining time series data, a unique trait known as long memory (LM) emerges. Time series exhibiting long memory demonstrate persistence in their patterns as the autocorrelation function (ACF) decays more gradually than with an exponential decay. Conversely, a time series with an ACF that decays rapidly is characterised as having short memory. Significant research has taken effect, especially covering long memory dynamics in financial markets. However, relatively little research has been conducted on commodity markets, particularly emerging agricultural assets and emerging resources such as lithium. Globally, derivative and futures commodity markets are characterised by a lack of consensus. Some researchers have come to varying suppositions with regard to the proficiency of financial commodity markets if, indeed, prices can be regarded as unbiased predictors of their spot prices. According to Moosa and Al-Loughani (1994), the future prices of crude oil are not efficient. On the one hand, market efficiency as objective crude oil price predictions were supported by data discovered by Kumar (1992). Avsar and Goss (2001) note that inefficiencies are likely to be exacerbated in relatively young, and shallow, markets, for example, the electricity market. In such a market, forecast errors may indicate a market still coming to terms with the actual market dynamics. However, Brenner and Kroner (1995) suggested that the inconsistencies observed between futures and spot prices could be the result of carrying costs rather than a failure of the efficient market hypothesis (EMH).
Commodities play an essential role in shaping the global economy (ARDA 2004). Developing countries anchor their economies on commodities such as crude oil, gold, lithium, tobacco, and cotton, to mention just a few. Lithium, for example, has received renewed interest with the impending battery vehicle takeover from the gas-driven ones. Commodities range from agricultural products, tropical beverages, vegetable oil seeds and oils, mineral ores, and metals. Variations in commodity prices continue to affect global economic activity. Commodities remain an integral source of export earnings, and commodity price movements significantly impact overall macroeconomic performance for many developing countries. Understanding the dynamics of commodity prices helps alleviate problems bedevilling countries in the Global South and has the potential to set economic activity on a positive trajectory.
LM behaviour of commodity prices has caught the attention of scholars in recent years, and debate is rife, as seen in research by Cuddington (1992), who found little evidence to support the view that prices of primary commodities were on a downward trend over the long term. This was supported by Cashin and McDermott (2002), although they also found that such trends had alternate high volatility of commodity prices.
The efficient market hypothesis proposed by Fama (1970) greatly influences theoretical and empirical finance. LM characteristics are important pointers for uncovering non-linear dependence in the mean and volatility of commodity prices. This has captured the interest of researchers; among them are Arouri et al. (2012), Diaz (2016), Ranganai and Khubeka (2016), and Chinhamu et al. (2022). Dual LM (DLM) plays an integral part in market performance and pricing modalities of commodities. Furthermore, LM implies that markets do not promptly respond to information flows but do so gradually over a period of time. Contrary to EMH, if commodity returns exhibit dual DLM, the data bear a predictable part; hence, past returns can predict future returns, a feature often abused by speculators and arbitrage traders.
The main aim of this study is to explore the dual LM (DLM) characteristics in the returns and variance of commodities. Autoregressive fractionally integrated moving average (ARFIMA)- and generalised autoregressive conditional heteroskedasticity (GARCH)-type models can provide useful insight into the relationship between the return and volatility of a process displaying the LM property instantaneously. Furthermore, we consider the distributional properties of commodity returns under heavy-tailed and asymmetric distributional assumptions. In particular, five (5) commodities forming part of the International Monetary Fund’s commodities price index are considered. They include crude oil, gold, lithium, cotton, and tobacco. We investigate dual long-memory dynamics inherent in these commodities. Furthermore, the discrepancy of the unpredictable or asymmetric response of the market to information flow will be explored. Recent studies have revealed the existence of spurious LM in the conditional variance when sudden changes or structural breaks arise Ural and Küçüközmen (2011)
The rest of the paper is organised as follows: The literature review is discussed in Section 2. In Section 3, we provide theoretical background on DLM-GARCH-type models, the Student-t distribution (StD), skewed-Student-t distribution (SStD), and generalized error distribution (GED). Section 4 presents the empirical results and offers a discussion. Finally, Section 5 concludes this work.

2. Literature Review

LM duality, which speaks to the presence of LM in both returns and conditional volatility, has been a subject of debate and has generated renewed interest in recent studies. Shi and Yu (2023), Chinhamu et al. (2022), Diaz (2016), Ranganai and Khubeka (2016), and Arouri et al. (2012) have accomplished tremendous work by considering joining GARCH-type LM models. Kasman et al. (2009) and Barkoulas et al. (2000) bemoaned the presence of LM in volatility and the importance of uncertainty or risk in characterising financial assets. This subject has generated tremendous interest towards modelling LM properties in returns and volatility in economics and finance. Recent literature attaches much significance to the estimation and forecasting of the volatility of asset returns and their applications in financial markets, such as derivative pricing, risk management (actuarial services, hedge funds, and general insurance), asset allocation, and valuation (Tsay 2005).
Several techniques have been used to analyse financial data in the literature. Chinhamu et al. (2022) used DLM heavy-tailed asymmetric volatility models for estimating value-at-risk (VaR) in precious metal returns. They evaluated the performances of LM GARCH models to estimate VaR for daily returns for platinum, gold, and silver. Results confirmed the adequacy of such models in risk management valuations and hedging stratagems. ARFIMA-fractionally integrated GARCH (FIGARCH), ARFIMA-hyperbolic GARCH (HYGARCH), and ARFIMA-fractionally integrated asymmetric power autoregressive conditional heteroskedasticity (FIAPARCH) models under normal inverse Gaussian, the variance gamma and the Pearson type-IV innovations, respectively, in addition to VaR estimation, were found to be suitable models for modelling the extreme risk of metal prices.
Krezłek (2012) used what he referred to as stable distributions to the quantity investment hazard rate of selected non-ferrous metals. The results recommended the stable distributions as risk assessment tools. Volatility clustering, heavy tails, and long memory were identified. Arouri et al. (2012), Cochran et al. (2012), and Chkili et al. (2014) also studied the application of non-linear volatility models with long LM in modelling precious metals. Their results confirmed the suitability and capability of FIGARCH models to produce good forecasts, precisely and accurately estimating precious metals’ VaR. Bentes (2015) studied the performances of GARCH, integrated GARCH (IGARCH) and FIGARCH on daily gold price returns. Results compared well and, in some instances, were similar to those of Cochran et al. (2012) and Ranganai and Khubeka (2016). Chkili et al. (2014) and Cochran et al. (2012) confirmed that FIGARCH (1, d, 1) was the best model to describe the variance processes for metal prices.
Kasman et al. (2009) used the ARFIMA, GPH, FIGARCH, and HYGARCH models to investigate the presence of LM in the stock markets of eight central and eastern European countries. The findings suggested that the ARFIMA-FIGARCH model might be used to model LM undercurrents in the returns and volatility. Additionally, they detected LM in volatility. They also found that uncertainty or risk is a significant factor influencing the daily stock data’s behaviour in the Turkish stock market. In their study, Kang and Yoon (2007) investigated the DLM property using two daily Korean stock price indexes and the ARFIMA-FIGARCH model. Their empirical findings show that the combined ARFIMA-FIGARCH model can accurately estimate the LM dynamics in the returns and volatility. Additionally, they suggested that it is preferable to incorporate the tendency of asymmetric leptokurtosis under the premise of an SStD. Tang and Shieh (2006) estimated FIGARCH and HYGARCH models with various distributions to examine the LM characteristics for daily closing prices of three stock index futures markets: the SP500, Nasdaq100, and Dow Jones.
Using a bivariate dual long-memory model, Karanasos and Kartsaklas (2009) investigated the dynamics of range-based volatility, turnover volume, and their corresponding uncertainties in the Korean market from 1995 to 2005. They found that the conditional variance series’ order of integration significantly drops when structural discontinuities are considered. In the Malaysian stock markets, Cheong et al. (2008) investigated the effects of a structural break on the fractionally integrated time-varying volatility model between 1996 and 2006. After accounting for structural breaks in volatility during the Asian crisis, their empirical results showed a significant decrease in the volatility of LM clustering. Forecasting models that can fully capture the empirical characteristics such as long memory, non-normality, asymmetry, and volatility of commodity price returns simultaneously are lacking in the literature. To the best of our knowledge, there are limited studies evaluating the suitability of ARFIMA-GARCH-type models on the five commodities.

3. Methodology

In this section, we present the theoretical background on the long-memory process and the long-memory volatility models. Additionally, we discuss the heavy-tailed distributions that are employed to capture non-normal innovations of the volatility models.

3.1. Long-Memory Process

A stationary process y t is a long-memory process if there is a real number H such that 0 < H < 1 , and ρ ( υ ) , the autocorrelation function, has a hyperbolic decay, lim υ ρ υ = K 2 H 1 , K > 0 , where K is a finite constant, and H is the Hurst exponent.

3.2. Long-Memory Mean and Volatility Models

In this section, we discuss the ARFIMA, FIGARCH, HYGARCH, FIAPARCH, ARFIMA-FIGARCH, ARFIMA-HYGARCH, and ARFIMA-FIAPARCH processes. The notation for the fractional parameters for the mean and volatility models shall be o and ξ , respectively.
The ARFIMA model is characterised by the autocorrelation function, which decays at a slower rate (Kang and Yoon 2007). The ARFIMA ( p , o , q ) can be expressed as follows:
ϕ B 1 B o y t μ = θ B ε t ,
ε t = z t σ t , with z t ~ N ( 0 ,   1 ) , which is independently and identically distributed with a common variance σ 2 . B denotes the lag operator or the backward shift; o is the fractional difference parameter, ϕ B = 1 ϕ 1 B ϕ 2 B 2 ϕ p B p ; and θ B = 1 θ 1 B θ 2 B 2 θ q B q are the AR and MA polynomials with necessary and sufficient stationarity conditions, in that order.
Hosking (1981) postulated that if 0.5 < o < 0.5 ,   y t is stationary and invertible. Therefore, the effect of shocks to ε t on y t depicts LM characteristics. If o = 0 , then the process is stationary, and short memory and the effect of shocks to ε t on y t decays geometrically. For o = 1 , the process is non-stationary. If 0 < o < 0.5 , then y t exhibits positive dependence between distant observations that are LM. If 0.5 < o < 0 , then y t exhibits negative dependence, thus intermediate memory (Kang and Yoon 2007). Baillie et al. (1996) proposed the FIGARCH model to capture LM in volatility. The FIGARCH (p, ξ ,   q ) can be expressed as follows:
ϕ B 1 B ξ ε t 2 = ω + [ 1 θ B ] v t ,
where, ϕ B = 1 ϕ 1 B ϕ 2 B 2 ϕ p B p and θ B = 1 θ 1 B θ 2 B 2 θ q B q are the autoregressive AR and moving-average MA polynomials as given above. An alternative representation for the FIGARCH (p, ξ ,   q ) model was presented by Ural and Küçüközmen 2011) as follows:
1 θ B σ t 2 = ω + 1 θ B ϕ B 1 B ξ ε t 2 .
The conditional of ε t 2 is obtained from:
σ t 2 = ω 1 θ B 1 + 1 ϕ B 1 B ξ 1 θ B 1 ε t 2 ,
which implies that σ t 2 = ω 1 θ B 1 + γ ( B ) ε t 2 , where γ B = 1 γ 1 B γ 2 B 2 γ 3 B 3 . Baillie et al. (1996) argued that the influence of the errors on the conditional variance of the FIGARCH ( p ,   ξ ,   q ) processes diminishes at a hyperbolic rate when 0 < ξ < 1 . Thus, the long-term dynamics of the volatility are considered by the parameter, and the short-term dynamics are modelled through the normal GARCH parameters. Of significance is the covariance stationarity of the FIGARCH. Whether or not it is strictly stationary remains open as argued by Conrad and Haag (2006). Conditions on the parameters have to be imposed on the FIGARCH model to ensure positive conditional variances, based on the work of Conrad and Haag (2006).
The HYGARCH model of Davidson (2004) is an extension of the FIGARCH model, introducing weights in the difference parameter as follows:
σ t 2 = ω 1 θ B 1 + { 1 - 1 θ B 1 ρ [ 1 + α 1 B ξ ] }   ε t 2 ,  
where ω > 0 , 0 < ξ < 1 and θ > 1 .
The HYGARCH allows the existence of the variance at more extreme degrees than the modest IGARCH and FIGARCH models. Conrad (2010) suggested the modification of Equation (5) to take the form:
ϕ B 1 α + α 1 B ξ ε t 2 = ω + θ B ε t 2 σ t 2 .
If α = 0 and α = 1 , the HYGARCH model becomes the stable GARCH and FIGARCH, respectively. When ξ = 1, the parameter α becomes an AR root, and the HYGARCH condenses to either a stationary GARCH ( α < 1 ), an IGARCH ( α = 1 ), or an explosive GARCH ( α > 1 ).
Tse (1998) introduced the FIAPARCH model, which is an extension of FIGARCH to simultaneously capture LM and asymmetric effects in the volatility. The FIAPARCH (p,   ξ , q) model can be presented as:
σ t δ = ω 1 θ B 1 + { 1 1 θ B 1   ϕ B 1 B ξ } ( | ε t | γ ε t ) δ ,
where ω > 0 , δ > 0   and   θ < 1 . γ is the asymmetric parameter with condition 1 < γ < 1 .
If γ > 0 , negative shocks will have more influence on the commodity return volatility than positive shocks of equal size. If 0 < ξ < 1 , the FIAPARCH model captures the LM in the volatility. When γ = 0 and δ = 2 , the FIAPARCH model becomes the FIGARCH model. The parameter δ is a power term in the volatility structure and should be inferred from the data.

3.3. Heavy-Tailed Distributions

To capture the non-normality exhibited by financial returns, the StD, SStD, and GED will be used.

3.3.1. Student-t Distribution (StD)

The probability density function of the univariate Student-t distribution is given (Arfken et al. 2013):
f x = Γ ( v + 1 2 ) ( σ v π Γ v 2 ) 1 1 + x 2 v σ 2 v + 1 2 ,
where for < x < with μ being the location parameter, σ > 0 is understood as the scale parameter, and ν > 0 degrees of freedom. The maximum likelihood is obtained through the application of the numerical optimization of the maximum likelihood function as given in Green (2005).

3.3.2. Skewed-Student-t Distribution (SStD)

Hansen (1994) defined the SStD as follows:
f x ,   λ ,   z = ρ Γ v + 1 2 π v 2 Γ v 2 . ( 1 + ς 2 v 2 ) ( v + 1 ) 2 ,
where ς = ρ x + a 1 λ if x < a ρ and ρ x + a 1 + λ if x a ρ . The constant terms a and ρ are defined as:
a = 4   λ c z 2 z 1   and   ρ = 1 + 3 λ 2 a 2 ,   c = Γ ( v + 1 2 ) π ( v 2 ) Γ ( v 2 ) .
In this probability density function, 2 < ν < defines ν (degrees of freedom parameter), and 1 < λ < 1 is the asymmetry parameter.

3.3.3. Generalized Error Distribution (GED)

The GED is a family of distributions that assumes a range of specific types relying on the value of the parameter ν , and which consists of the normal distribution as an exceptional case. The GED is a more flexible generalization of the normal distribution and is consequently defined by three parameters. The mean ( μ ) determines the peak of the distribution. In the standard normal distribution, the median and mode are equal to the mean ( μ ) and the standard deviation ( σ ) , which determines the dispersion and shape parameter ( β ) , which is referred to as the kurtosis and reveals how much data are in the tails. Given the above background, the definition of the GED is as follows (Giller 2005):
F x =   λ z 2 Γ ( k ) e x p ( λ z x μ 2 ) ,  
for < x < , where Γ ( k ) is a Euler function, z is a shape parameter, λ is a scale parameter, and μ is a location parameter. The new approach to determine the estimates of the GED through MLE estimation was conceived by Purczyński and Bednarz-Okrzyńska (2014).

3.4. Long-Memory Tests

In this study, we utilise the Geweke and Porter-Hudak (GPH) to test for long memory in returns and squared returns.

GPH Long-Memory Test

The Geweke and Porter-Hudak (1983) estimator is widely used to distinguish between long-memory and short-memory effects and is called the spectral regression method. The spectral density of the fractionally integrated process y t is:
f ω = 4 sin 2 ω 2 d f u ω ,
where ω is the Fourier frequency, and f u ( ω ) is the spectral density corresponding to μ t . The difference parameter d can be estimated as:
ln f ω j = β d l n [ 4 s i n 2 ω 2 ] + ε j ,
where, for j = 1 ,   2 ,   ·   ·   ·   ·   ·   · ,   n f T , Geweke and Porter-Hudak showed that the least squares estimate d ^ using regression is normally distributed in large samples if n f   ( T ) = T α , with 0 < α < 1 : with d ^ ~ N [ d ,   π 2 6 j = 1 n f U j U ¯ 2 ] , where U j = 4 sin 2 ω 2 and U ¯ is the sample mean of U j , under the null hypothesis of no memory d = 0 , the test statistic is:
t d = 0 = d ^ [ d ,   [ π 2 6 j = 1 n f U j U ¯ 2 ] ] 0.5 .

4. Empirical Results and Discussion

The data comprised 5722 daily closing prices of five commodities, crude oil, gold, cotton, lithium, and tobacco, running from 2 January 2001 to 26 October 2023. The data set was divided into two samples. The period running from 2001 to 26 October 2018 was regarded as the in-sample data. The period from 27 October 2018 to 26 October 2023 was the out-of-sample data. A daily measure of price volatility was calculated after log transformation of the data. It is also possible to choose from several alternative measures depending on availability and the nature of the pricing regime. The commodity (crude oil, gold, lithium, cotton, and tobacco) prices extracted from the LBMA database were used. The return series for each commodity was generated using first forward differences of the natural logarithm of the commodity price defined as:
r t = ln P t l n ( P t 1 ) ,
where P t is the closing price at day t .

4.1. Summary Statistics

The descriptive statistics, correlation, normality, heteroscedasticity, unit root, and stationary tests are reported in Table 1 for the log-returns of the five commodities.
Table 1 shows descriptive statistics, split into three categories as follows: the measures of central tendency; skewness and kurtosis, independence, and normality; and stationarity tests. A general increasing trend in commodity prices is evidenced by all positive return means. Skewed distributions and volatility clustering are again evident on all commodity return indices. Heavy-tailed distributions are recommended.
The JB test gives significant p-values for all five commodities return series. This implies that the commodity returns are not normality distributed. The Ljung–Box test, at Q(5), Q(10) and Q 2 (5), and Q 2 (10) indicates that autocorrelation is insignificant for returns but highly significant in squared returns indicating persistence in the volatility process of commodity returns. Heavy-tailed asymmetric GARCH-type models are recommended under these circumstances.
The augmented Dickey–Fuller (ADF) and Phillips and Perron (PP) tests rejected the unit root hypothesis while the Kwiatkowski test failed to reject the unit root hypothesis of stationarity. Hence, commodity price returns are largely stationary in the mean.

4.2. Preliminary Analysis

Further analysis was performed using time series plots of the closing prices, time series plots of the log-returns, Q-Q and box plots of the returns, and the ACF plots of both the returns and squared returns. The results are reported in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5.
Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5 show plots of the time series, return series, box plots, Q-Q plots, ACF of return, and ACF of square returns for the five commodities. The returns of all the indices exhibit volatility clustering. All the ACFs and PACFs of the returns series and squared returns series suggest the existence of LM duality. All the QQ plots show that the tails of all the commodity indices’ returns are heavier than the tails of the normal distribution. They indicate the presence of heavy-tailed distributions and asymmetric dispersion of all the returns. Such evidence advocates the use of volatility models such as asymmetric models or heavy-tailed distribution to account for kurtosis and skewness aspects of the data. Adjudged by box plots, commodity returns are skewed and heavy-tailed.

4.3. Testing for LM

A long-memory test was performed using GPH. The GPH’s long-memory test gave plausible results indicating the presence of DLM in three commodities as presented in Table 2.
Based on d -values f o r   d = 0.5 ,   0.6 , and 0.7 , GPH depicts the presence of DLM in both lithium and tobacco return series. Gold at a bandwidth of 0.7 shows LM in returns (mean). It is worth trying to fit a joint ARFIMA-GARCH-type model.

4.4. Summary of Empirical Properties of Commodity Returns

From the data exploration, it can be concluded that all three commodities considered exhibit heavy tails, volatility clustering, and long memory of both returns and squared returns while one commodity exhibited LM in squared returns only. The models are all stationary; hence, the suggested models are long-memory ARFIMA-GARCH-type models with heavy-tailed innovations. This study will explore ARFIRMA-FIGARCH, ARFIMA-HYGARCH, and ARFIRMA-FIAPARCH with heavy-tailed innovations. Distribution such as StD, SStD, and GED will be considered among other distributions. We also found that the assumption of a SStD is commendable for incorporating the tendency of asymmetric leptokurtosis in a return distribution as postulated by Kang and Yoon (2007).

4.5. Empirical Analysis

In this section, we present the parameter estimation results of the several joint GARCH-type models encompassing ARFIMA-FIGARCH, ARFIMA-HYGARCH, and ARFIMA-FIAPARCH models. Several models were experimented with, and the models presented in Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9 and Table 10 are the most plausible ones under heavy-tailed innovation distributions StD, SStD, and GED benchmarked against the normal distribution results. The ARFIMA model was used for modelling returns, and for volatility we used FIGARCH, FIEGARCH, FIAPARCH, and HYGARCH models under the same error distributions.

4.6. Discussion

Heavy-tailed and asymmetric LM volatility models under a variety of conditional distributed innovations were fitted. Because the residuals suffer from excess kurtosis and skewness, the assumption of a normal distribution is not suitable for capturing asymmetry and fat tails in most return series. These included Student distribution, skewed-Student distribution, and generalised error distribution. The models fitted include FIGARCH, HYGARCH, FIAPARCH, ARFIMA-FIGARCH, ARFIMA-HYGARCH, and ARFIMA-FIAPARCH under the above innovations.
As is standard procedure, our model selection was multipronged. Firstly, attention was given to significance levels of parameters. Statistics and their corresponding p-values were considered. Secondly, when the models were fitted, Q-statistics on standardised residuals and squared standardised residuals were extracted for further tests to confirm post-fitting serial correlation status. In addition, Akaike Information Criteria was to be minimised. Table 5, Table 6, Table 7, Table 8, Table 9 and Table 10 show these diagnostic benchmarks in addition to further normality tests for skewness and excess kurtosis, along with the Jarque–Bera test. It has become standard procedure to mine the standardised residuals when fitting a model and expose the extracted residuals to further ARCH LM test. Suffice to say that all the models selected passed all these litmus tests including efficiency and sufficiency tests. Hence, they can be considered adequate.
Table 3 and Table 7 show parameter estimates adjacent to the associated adequacy diagnostic checking mechanisms of the FIGARCH, HYGARCH, and FIAPARCH models for crude oil and cotton returns, respectively. In both cases, the coefficients satisfy the BLUE conditions for the positivity of the conditional variances as referenced by Conrad and Haag (2006). Our gold results compare well with the those obtained by Chinhamu et al. (2022).
The joint ARFIMA-FIGARCH, ARFIMA-HYGARCH, and ARFIRMA-FIAPARCH models for gold, lithium, and tobacco under different error distributions are presented in Table 4, Table 5, Table 6 and Table 8. While lithium and tobacco models are adequate under different conditional distributions, gold is suitable under the normal and GED assumptions. ARFIMA-HYGARCH for gold failed under heavy-tailed and asymmetric error distribution assumptions. However, for gold (normal errors) and tobacco models, the long-range dependence parameter of the ARFIMA model was negative 0.5 < o < 0 indicating intermediate persistence. This illustrates that log returns of both gold and tobacco are mean reverting, and hence, will revert to the mean overtime. Gold has a positive fractional parameter under the heavy-tailed distributions. For volatility, the LM parameter is positive; 0 < ξ < 1 shows strong LM. This confirms the results from Chinhamu et al. (2022) as gold shows high persistence.

4.7. Sensitivity/Forecast Evaluation

We discuss forecasting performance of all the fitted commodity models in the table below. Evaluation of forecasts for models is important as it helps us understand the forecasting accuracy of the models estimated. There are a number of forecast evaluation measures available in the literature. The three measures commonly used are the mean square error (MSE), the mean absolute error (MAE), the Theil Inequality Coefficient (TIC) and R-squared.
The performance of the forecasts is evaluated by the measures MSE, MAE, TIC, and R-squared as shown in Table 11 and Table 12. The MSE for gold gives low prediction errors for all models. Based on the MAE, the ARFIMA-FIGARCH under the ND and StD and the ARFIMA-HYGARCH model under the ND give fewer prediction errors. Overall, based on the TIC, the ARFIMA-FIGARCH under StD gives fewer prediction errors for gold. Hence, the ARFIMA-FIGARCH under the GED performs best, and thus it is the selected model ahead of all other models.
For lithium, the ARFIMA-FIGARCH and the ARFIMA-FIAPARCH under the ND were the only plausible models. All other assumptions did not yield convergence parameters. Based on the MSE, MAE, and TIC, the ARFIMA-FIAPARCH gives fewer prediction errors. The ARFIMA-FIAPARCH model was the best model. Model selection for the other commodities is summarised in Table 13 below.

4.8. Model Selection

Table 13 shows the best models selected for each commodity based on a balance of various criterions as argued above. The superior forecasting model was chosen based on a balance of probabilities considering the AIC and best volatility predictions as well as on information in the table below.

Out-of-Sample Validation of Top-Performing Models

In order to validate the adequacy of the best forecasting model for each commodity price return, we fitted them to the out-of-sample data set. The best forecasting model for crude oil price returns is a FIAPARCH (1,   ξ , 1). The out-of-sample parameter estimates of the hybrid FIAPARCH (1,   ξ , 1)-StD model for crude oil returns are shown in Table 14. Table 15, Table 16, Table 17 and Table 18 show the results of out-of-sample parameter estimates for the best-fitting forecasting models for gold, lithium, cotton, and tobacco price returns, respectively.
From Table 14, Table 15, Table 16, Table 17 and Table 18, it is clear that the ML parameter estimates for the selected best-fitting models for each commodity price return are statistically significant. The Ljung–Box test and ARCH-LM test results confirm that the selected best-fitting forecasting models for each commodity price return are adequate. Thus, the proposed models can be used as alternative models for forecasting commodity price returns.

5. Conclusions and Further Research

This study extended the work of Youssef et al. (2015), Ranganai and Khubeka (2016), and Chinhamu et al. (2022) by determining the best forecasting models with which to capture dual long-memory, non-normality, asymmetry, and volatility-clustering features of commodity price returns. We included commodity returns in the portfolio of financial assets with DLM characteristics. Research findings revealed that commodity returns are characterised by asymmetry, heavy tail, and DLM. We examined the DLM GARCH models under the StD, SStD, and GED assumptions. The conditional variance and LM in mean and volatility were modelled by ARFIMA-FIGARCH, ARFIMA-HYGARCH, and ARFIMA-FIAPARCH models with normal assumptions while the StD, SStD, and GED distributions were applied to capture the fat-tail behaviour for the extracted standardised residuals. Model adequacy checking was performed using p-values, the ARCH LM test, and the Ljung–Box Q statistics of standardised and squared residuals. Results show that ARFIMA-FIGARCH, ARFIMA-HYGARCH, and ARFIMA-FIAPARCH models with StD, SStD, and GED governing the innovations are suitable for depicting crude oil, cotton and gold returns. The ARFIMA-HYGARCH with GED governing the innovations is the overall best model for gold returns. Our results are consistent with the results of Chinhamu et al. (2022), Ranganai and Khubeka (2016), and Youssef et al. (2015). For future study, we recommend obtaining a one-day-ahead return forecast for the proposed models in order to determine the sign for the next day’s return. A robust model should correctly forecast the direction (up or down) of the next day’s return. This will be useful for commodity market practitioners.

Author Contributions

Conceptualization, K.B.; Methodology, K.C. and R.C.; Software, K.B., L.D., K.C. and R.C.; Validation, K.C. and R.C.; Formal analysis, K.B., L.D. and F.M.; Resources, F.M.; Writing—original draft, K.B.; Writing—review & editing, L.D., K.C., R.C. and F.M.; Supervision, L.D. and F.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. ARDA. 2004. Agricultural and Rural Development Authority Stratigic plans. Zimbabwe: The National Development Strategy 1: 2021–25. [Google Scholar]
  2. Arfken, George B., Hans J. Weber, and Frank E. Harris. 2013. Mathematical Methods for Physicists. Cambridge: Academic Press. [Google Scholar]
  3. Arouri, Mohamed El Hedi, Shawkat Hammoudeh, Amine Lahiani, and Duc Khuong Nguyen. 2012. Long memory and struc tural breaks in modelling the return and volatility dynamics of precious metals. The Quarterly Review of Economics and Finance 52: 207–18. [Google Scholar] [CrossRef]
  4. Avsar, Gulay S., and Barry A. Goss. 2001. Forecast Errors and Efficiency in the U.S. Electricity Futures Market. Australian Economic Papers 40: 479–99. [Google Scholar] [CrossRef]
  5. Baillie, Richard T., Tim Bollerslev, and Hans Ole Mikkelsen. 1996. Fractionally Integrated Generalized Autoregressive Conditional Heteroscedasticity. Journal of Econometrics 74: 3–30. [Google Scholar] [CrossRef]
  6. Barkoulas, John T., Christopher F. Baum, and Nickolaos Travlos. 2000. Long Memory in the Greek Stock Market. Applied Financial Economics 10: 177–84. [Google Scholar] [CrossRef]
  7. Bentes, Sonia R. 2015. Forecasting volatility in gold returns under the GARCH, IGARCH and FIGARCH frameworks: New evidence. Physica A 438: 355–64. [Google Scholar] [CrossRef]
  8. Brenner, Robin J., and Kenneth F. Kroner. 1995. Arbitrage, Cointegration, and Testing the Unbiasedness Hypothesis in Financial Markets. The Journal of Financial and Quantitative Analysis 30: 23–42. [Google Scholar] [CrossRef]
  9. Cashin, Paul, and C. John McDermott. 2002. The long-run behavior of commodity prices: Small trends and big variability. IMF Staff Papers 49: 175–99. [Google Scholar] [CrossRef]
  10. Cheong, Chin Wen, Zaidi Isa, and Nor Abu Hassan Shaari Mohd. 2008. Fractionally Integrated Time-varying Volatility under Structural Break: Evidence from Kuala Lumpur Composite Index. Sains Malaysiana 37: 405–11. [Google Scholar]
  11. Chinhamu, Knowledge, Retius Chifurira, and Edmore Ranganai. 2022. Value-at-Risk Estimation of Precious Metal Returns using Long Memory GARCH Models with Heavy-Tailed Distribution. Journal of Statistics Applications & Probability 11: 89–107. [Google Scholar] [CrossRef]
  12. Chkili, Walid, Shawkat Hammoudeh, and Duc Khuong Nguyen. 2014. Volatility forecasting and risk management for commodity markets in the presence of asymmetry and long memory. Energy Economics 41: 1–18. [Google Scholar] [CrossRef]
  13. Cochran, Steven J., Iqbal Mansur, and Babatunde Odusami. 2012. Volatility persistence in metals returns: A FIGARCH approach. Journal of Economics and Business 64: 287–305. [Google Scholar] [CrossRef]
  14. Conrad, Christian, and Berthold R. Haag. 2006. Inequality Constraints in the Fractionally Integrated GARCH Model. Journal of Financial Econometrics 4: 413–49. [Google Scholar] [CrossRef]
  15. Conrad, Christian. 2010. Non-negativity for the hyperbolic GARCH model. Journal of Econometrics 157: 441–57. [Google Scholar] [CrossRef]
  16. Cuddington, John T. 1992. Long-run Trends in 26 Primary Commodity Prices: A Disaggregated Look at the Prebisch-Singer Hypothesis. Journal of Development Economics 9: 207–27. [Google Scholar] [CrossRef]
  17. Davidson, John. 2004. Moment and memory properties of linear conditional heteroscedastic models, and a new model. Journal of Business and Economic Statistics 22: 16–29. [Google Scholar] [CrossRef]
  18. Diaz, John Francis T. 2016. Do scarce precious metals equate to safe habour investments? The case of platinum and palladium. Economics Research International 12: 2361954. [Google Scholar] [CrossRef]
  19. Fama, Eugene F. 1970. Efficient Capital Markets: A Review of Theory and Empirical Work. Journal of Finance 25: 383–417. [Google Scholar] [CrossRef]
  20. Geweke, John, and Susan Porter-Hudak. 1983. The estimation and application of long memory time series models. Journal Time Series Analysis 4: 221–38. [Google Scholar] [CrossRef]
  21. Giller, Graham L. 2005. A Generalized Error Distribution. Available online: https://ssrn.com/abstract=2265027 (accessed on 1 March 2024).
  22. Green, Christopher. 2005. Heavy-Tailed Distributions in Finance: An Empirical Study. Seattle: Department of Statistics, University of Washington. [Google Scholar]
  23. Hansen, Bruce E. 1994. Autoregressive Conditional Density Estimation. International Economic Review 35: 705–30. [Google Scholar] [CrossRef]
  24. Hosking, J. R. M. 1981. Fractional Differencing. Biometrika 68: 165–76. [Google Scholar] [CrossRef]
  25. Kang, Sang Hoon, and Seong-Min Yoon. 2007. Long Memory Properties in Return and Volatility: Evidence from the Korean Stock Market. Physica A 385: 591–600. [Google Scholar] [CrossRef]
  26. Karanasos, Menelaos G., and Aris Kartsaklas. 2009. Dual Long-Memory, Structural Breaks and the Link between Turnover and The Range-Based Volatility. Journal of Empirical Finance 16: 838–51. [Google Scholar] [CrossRef]
  27. Kasman, Adnan, Saadet Kasman, and Erdost Torun. 2009. Dual Long Memory Property in Returns and Volatility: Evidence from the CEE Countries’ Stock Markets. Emerging Markets Review 10: 122–39. [Google Scholar] [CrossRef]
  28. Krezłek, Dominik. 2012. Non-Classical measures of investment Risk on the Market of Precious Non-Ferrous Metals. Dynamic Economic Models 12: 89–103. [Google Scholar] [CrossRef]
  29. Kumar, Manmohan S. 1992. The Forecasting Accuracy of Crude Oil Futures Prices. IMF Economic Review (Staff Papers) 39: 432–61. [Google Scholar] [CrossRef]
  30. Moosa, Imad A., and Nabeel E. Al-Loughani. 1994. Unbiasedness and Time-varying Risk Premia in the Crude Oil Futures Market. Energy Economics 16: 99–105. [Google Scholar] [CrossRef]
  31. Purczyński, Jan, and Kamila Bednarz-Okrzyńska. 2014. Estimation of the shape parameter of GED distribution for a small sample size. Folia Oeconomica Stetinensia 14: 35–46. [Google Scholar] [CrossRef]
  32. Ranganai, Edmore, and Sihle B. Khubeka. 2016. Long memory mean and volatility models of platinum and palladium return series under heavy tailed distributions. Springer Plus 5: 2089. [Google Scholar] [CrossRef]
  33. Shi, Shuping, and Jun Yu. 2023. Volatility puzzle: Long memory or antipersistency. Management Science 69: 3861–83. [Google Scholar] [CrossRef]
  34. Tang, Ta-Lun, and Shwu-Jane Shieh. 2006. Long Memory in Stock Index Futures Markets: A Value-at-Risk Approach. Physica A 366: 437–48. [Google Scholar] [CrossRef]
  35. Tsay, Ruey S. 2005. Analysis of Financial Time Series. Hoboken: John Wiley & Sons, pp. 97–143. [Google Scholar]
  36. Tse, Yiu Kuen. 1998. The Conditional Heteroscedasticity of the Yen-Dollar Exchange Rate. Journal of Applied Econometrics 13: 49–55. [Google Scholar] [CrossRef]
  37. Ural, Mert, and C. Coşkun Küçüközmen. 2011. Analyzing the Dual Long Memory in Stock Market Returns. Ege Academic Review 11: 19–28. [Google Scholar]
  38. Youssef, Manel, Lofti Balkacem, and Khaled Mokni. 2015. Value-at-Risk estimation of energy commodities: A long-memory GARCH–EVT approach. Energy Economics 51: 99–110. [Google Scholar] [CrossRef]
Figure 1. Descriptive plots: Top-left panel—Time series plot of the daily closing price; Top-middle panel—Time series plot of the daily returns; Top-right panel—QQ plot of the daily returns; Bottom-left panel—Box plot of the returns; Bottom-middle panel—ACF plot of the returns; Bottom-right panel—ACF of the squared returns for gold.
Figure 1. Descriptive plots: Top-left panel—Time series plot of the daily closing price; Top-middle panel—Time series plot of the daily returns; Top-right panel—QQ plot of the daily returns; Bottom-left panel—Box plot of the returns; Bottom-middle panel—ACF plot of the returns; Bottom-right panel—ACF of the squared returns for gold.
Risks 12 00073 g001
Figure 2. Descriptive plots: Top-left panel—Time series plot of the daily closing price; Top-middle panel—Time series plot of the daily returns; Top-right panel—QQ plot of the daily returns; Bottom-left panel—Box plot of the returns; Bottom-middle panel—ACF plot of the returns; Bottom-right panel—ACF of the squared returns for crude oil.
Figure 2. Descriptive plots: Top-left panel—Time series plot of the daily closing price; Top-middle panel—Time series plot of the daily returns; Top-right panel—QQ plot of the daily returns; Bottom-left panel—Box plot of the returns; Bottom-middle panel—ACF plot of the returns; Bottom-right panel—ACF of the squared returns for crude oil.
Risks 12 00073 g002
Figure 3. Descriptive plots: Top-left panel—Time series plot of the daily closing price; Top-middle panel—Time series plot of the daily returns; Top-right panel—QQ plot of the daily returns; Bottom-left panel—Box plot of the returns; Bottom-middle panel—ACF plot of the returns; Bottom-right panel—ACF of the squared returns for tobacco.
Figure 3. Descriptive plots: Top-left panel—Time series plot of the daily closing price; Top-middle panel—Time series plot of the daily returns; Top-right panel—QQ plot of the daily returns; Bottom-left panel—Box plot of the returns; Bottom-middle panel—ACF plot of the returns; Bottom-right panel—ACF of the squared returns for tobacco.
Risks 12 00073 g003
Figure 4. Descriptive plots: Top-left panel—Time series plot of the daily closing price; Top-middle panel—Time series plot of the daily returns; Top-right panel—QQ plot of the daily returns Bottom-left panel—Box plot of the returns; Bottom-middle panel—ACF plot of the returns; Bottom-right panel—ACF of the squared returns for cotton.
Figure 4. Descriptive plots: Top-left panel—Time series plot of the daily closing price; Top-middle panel—Time series plot of the daily returns; Top-right panel—QQ plot of the daily returns Bottom-left panel—Box plot of the returns; Bottom-middle panel—ACF plot of the returns; Bottom-right panel—ACF of the squared returns for cotton.
Risks 12 00073 g004
Figure 5. Descriptive plots: Top-left panel—Time series plot of the daily closing price; Top-middle panel—Time series plot of the daily returns Top-right panel—QQ plot of the daily returns; Bottom-left panel—Box plot of the returns; Bottom-middle panel—ACF plot of the returns; Bottom-right panel—ACF of the squared returns for lithium.
Figure 5. Descriptive plots: Top-left panel—Time series plot of the daily closing price; Top-middle panel—Time series plot of the daily returns Top-right panel—QQ plot of the daily returns; Bottom-left panel—Box plot of the returns; Bottom-middle panel—ACF plot of the returns; Bottom-right panel—ACF of the squared returns for lithium.
Risks 12 00073 g005
Table 1. Descriptive statistics, normality, stationarity and ARCH test.
Table 1. Descriptive statistics, normality, stationarity and ARCH test.
Crude OilGoldLithiumCottonTobacco
Descriptive Statistics
Minimum−0.5100−0.0982−0.0597−0.3184−0.5878
Maximum0.31100.08420.07590.17970.7865
Mean0.00020.00040.00020.00010.0003
Standard
Deviation
0.02740.01180.01160.02350.0836
Skewness−1.4653−0.47580.6375−1.42710.0757
Kurtosis39.795421.64839.594626.743814.1570
Testing for correlation (Q-test), normality (JB test), and heteroscedasticity ( Q 2 test) (p-values in brackets)
Q(5)40.5
(<0.0001)
3.64
(<0.0001)
208.7
(<0.0001)
185.5
(<0.0001)
260.0
(<0.0001)
Q(10)45.3
(<0.0001)
17.6
(<0.0001)
3453.7
(<0.0001)
210.4
(<0.0001)
272.8
(<0.0001)
JB test30,000
(<0.0001)
6385
(<0.0001)
31,123
(<0.0001)
10,000
(<0.0001)
28,634
(<0.0001)
Q 2 (5)118.5
(<0.0001)
43.7
(<0.0001)
5295.4
(<0.0001)
666.94
(<0.0001)
674.3
(<0.0001)
Q 2 (10)207.5
(<0.0001)
164.5
(<0.0001)
571.7
(<0.0001)
1352
(<0.0001)
1123
(<0.0001)
Unit root and stationary tests
ADF−17.4
(<0.0001)
−18.3
(<0.0001)
−18.3
(<0.0001)
−32.2
(<0.0001)
−37.9
(<0.0001)
PP test−5764
(<0.0001)
−5467
(<0.0001)
−5466
(<0.0001)
−3983
(<0.0001)
−4725
(<0.0001)
KPSS0.0640.2550.2550.0040.024
(>0.1000)(>0.1000)(>0.1000)(>0.1000)(<0.1000)
Table 2. GPH’s long-memory test for returns and squared returns.
Table 2. GPH’s long-memory test for returns and squared returns.
ReturnsSquared Returns
CommodityBandwidth d m Std Devp-Value d v Std Devp-Value
Crude oil d = 0.5 0.00120.08120.99600.13360.04170.0110
d = 0.6 0.05770.05430.25200.38590.0286<0.0001
d = 0.7 0.00840.03390.80100.35380.0122<0.0001
Gold d = 0.5 −0.02330.80000.80300.44670.0986<0.0001
d = 0.6 −0.07760.05130.14300.54240.0674<0.0001
d = 0.7 0.17820.03370.02300.46900.0559<0.0001
Cotton d = 0.5 0.20780.36870.15060.82700.0932<0.0001
d = 0.6 0.16280.16370.23870.66330.0578<0.0001
d = 0.7 0.51170.26450.17430.15470.0469<0.0001
Tobacco d = 0.5 0.10940.17350.18000.60270.1236<0.0001
d = 0.6 −0.35260.04600.04400.66430.0796<0.0001
d = 0.7 −0.47390.0385<0.00010.72320.0475<0.0001
Lithium d = 0.5 0.52360.0967<0.00010.27430.0893<0.0001
d = 0.6 0.65640.0873<0.00010.39540.0567<0.0001
d = 0.7 0.54280.0574<0.00010.3486 0.0558<0.0001
Table 3. ML parameter estimates for the LM-GARCH-type models with different error distributions and diagnostic tests (crude oil returns).
Table 3. ML parameter estimates for the LM-GARCH-type models with different error distributions and diagnostic tests (crude oil returns).
F I G A R C H ( 1 , ξ , 1 ) HYGARCH   ( 1 ,   ξ , 1) FIAPARCH   ( 1 , ξ , 1)
InnovationsParametersEstimatep-ValueEstimatep-ValueEstimatep-Value
ND ξ 0.5296<0.00010.5340<0.00010.3456<0.0001
α 1 0.2625<0.00010.2632<0.00010.3449<0.0001
β 1 0.7123<0.00010.7488<0.00010.6844<0.0001
γ 1 ----0.51740.0421
δ 1 ----2.1067<0.0001
Q (5)2.12760.86722.25780.87623.35850.6423
Q (10)3.14610.98973.01340.98434.63740.9358
Q (20)8.35941.0008.28451.00009.53330.9743
Q 2 (5)6.26120.34765.46810.28575.23590.1674
Q 2 (10)11.87450.235211.08750.234211.77650.2254
Q 2 (20)20.24360.446318.86470.531219.5420.4456
ARCH(5)7.38560.28646.00510.39856.12650.4143
ARCH(10)11.86770.359510.93280.354311.32870.3368
AIC−4.8463−4.8484−4.8434
StD ξ 0.5323<0.00010.5394<0.00010.3163<0.0001
α 1 0.2699<0.00010.2602<0.00010.8432<0.0001
β 1 0.7148<0.00010.7184<0.00010.7657<0.0001
γ 1 ----0.25180.04580
δ 1 ----1.7294<0.0001
Q (5)2.12770.83122.08760.836913.34440.2647
Q (10)3.03320.98063.00070.98144.25610.9351
Q (20)8.35650.98928.27540.98999.45160.9771
Q 2 (5)6.05880.10885.44900.14175.16020.1665
Q 2 (10)11.51060.174410.44380.235210.97940.2329
Q 2 (20)19.02660.390218.18990.443218.18420.4438
ARCH(5)5.99610.30665.38080.37125.11440.4221
ARCH(10)11.5930.313210.5270.395511.30570.3343
AIC−4.8564−4.8574−4.8748
SStD ξ 0.5484<0.00010.5557<0.00010.31480.0435
α 1 0.3001<0.00010.2534<0.00010.3478<0.0001
β 1 0.71675<0.00010.7145<0.00010.6531<0.0001
γ 1 ----0.56900.4221
δ 1 ----1.6988<0.0001
Q (5)2.34210.83362.10180.83773.47910.6266
Q (10)3.44470.98782.98690.98444.39710.9277
Q (20)7.21450.98328.24111.00009.64530.9741
Q 2 (5)6.61430.23255.42980.27796.21380.1017
Q 2 (10)10.91670.23759.72450.282312.45480.1320
Q 2 (20)18.64230.417915.76310.477919.42810.3659
ARCH(5)6.62370.37384.65850.431866.14110.2927
ARCH(10)11.31000.35829.83780.4639112.7560.2376
AIC−4.8821−4.8871−4.8894
GED ξ 0.5454<0.00010.5668<0.00010.4357<0.0001
α 1 0.3429<0.00010.2667<0.00010.3448<0.0001
β 1 0.7654<0.00010.7432<0.00010.6683<0.0001
γ 1 ----0.4954<0.0001
δ 1 ----1.7528<0.0001
Q (5)2.14010.82942.09030.83653.29870.6575
Q (10)3.03490.98062.98960.98174.21620.9389
Q (20)8.35840.98928.24420.99019.34370.9654
Q 2 (5)4.69970.19514.09830.25104.65360.3267
Q 2 (10)9.93110.26998.733170.36538.77660.3614
Q 2 (20)17.69820.475716.88680.530916.27650.5776
ARCH(5)4.66070.45874.04690.54273.25720.6512
ARCH(10)10.0290.43798.80810.550410.12020.5219
AIC−4.8830−4.8866−4.8876
Table 4. ML parameter estimates for the ARFIMA (1, o , 1)-FIGARCH (1, ξ , 1) with different error distributions and diagnostic tests (gold returns).
Table 4. ML parameter estimates for the ARFIMA (1, o , 1)-FIGARCH (1, ξ , 1) with different error distributions and diagnostic tests (gold returns).
N D StD S S t D GED
ParameterEstimatep-ValueEstimatep-ValueEstimate p-ValueEstimatep-Value
ϕ 0 0.00030.00030.00040.00010.00040.01300.0025<0.0001
o −0.08550.04450.01650.45260.01580.39640.03450.0174
θ 1 0.6938<0.0001−0.06690.0216−0.05540.0158−0.05280.0032
β 0 0.06120.00760.04690.00580.04480.00360.05590.0044
ξ 0.3354<0.00010.4275<0.00010.4156<0.00010.3820<0.0001
α 1 0.2981<0.00010.2568<0.00010.2645<0.00010.2853<0.0001
β 1 0.5924<0.00010.6667<0.00010.6575<0.00010.6277<0.0001
Q (5)2.46590.48565.46850.14115.66120.22615.94110.2036
Q (10)5.75570.67489.42970.31119.87890.361510.29340.3273
Q (20)11.67730.866416.16430.585616.56650.621616.84770.6002
Q 2 (5)3.12310.393818.80940.000619.36490.00069.79280.0554
Q 2 (10)5.35340.723421.86280.03820.22550.008911.93870.1540
Q 2 (20)6.98820.995622.43420.213121.86850.237713.64110.7522
ARCH(5)3.21370.701320.25780.026917.90990.00359.59230.0877
ARCH(10)6.01250.967818.26370.002819.81980.031211.81500.2976
AIC−6.3239−6.4068−6.4084−6.4176
Table 5. ML parameter estimates for the ARFIMA (1, o , 1)-HYGARCH (1, ξ , 1) with different error distributions and diagnostic tests (gold returns).
Table 5. ML parameter estimates for the ARFIMA (1, o , 1)-HYGARCH (1, ξ , 1) with different error distributions and diagnostic tests (gold returns).
N D StD S S t D GED
ParameterEstimatep-ValueEstimatep-ValueEstimate p-ValueEstimatep-Value
ϕ 0 0.0030.00030.0005---0.0005<0.0001
o −0.0764<0.00010.01130.12000.03430.25430.03340.0235
θ 1 0.6985<0.00010.14470.57620.69780.29390.08120.6987
β 0 0.06450.0342−0.21640.4132--−0.15160.4732
ξ 0.3450<0.00010.2912<0.00010.2963<0.00010.2683<0.0001
α 1 0.3015<0.00010.3065<0.00010.3061<0.00010.3209<0.0001
β 1 0.5954<0.00010.6438<0.00010.6376<0.00010.6085<0.0001
Q (5)2.44820.48474.761540.19015.45290.14485.75670.1923
Q (10)5.77780.672510.13970.255410.30580.24428.97540.2700
Q (20)11.67680.863517.75360.472017.70540.4752117.31380.5268
Q 2 (5)3.29420.388118.53930.000318.18500.00049.57620.0735
Q 2 (10)5.37540.716920.20170.009619.84970.010911.51870.1726
Q 2 (20)6.99721.000021.82310.239921.46780.256513.21290.7784
ARCH(5)3.12960.695418.14700.00287.78300.00329.41570.1943
ARCH(10)5.35650.86619.96100.029619.33230.033211.44460.3274
AIC−6.3852−6.4013−6.4063−6.4127
Table 6. ML parameter estimates for the LM-GARCH-type models with normal innovations and diagnostic tests (lithium returns).
Table 6. ML parameter estimates for the LM-GARCH-type models with normal innovations and diagnostic tests (lithium returns).
ARFIMA   ( 1 ,   o , 1)-FIGARCH (1, ξ , 1) ARFIMA   ( 1 ,   o , 1)-FIAPARCH (1, ξ , 1)
ParameterEstimatep-ValueEstimatep-Value
o 0.8156<0.00010.8396<0.0001
θ 1 0.21240.00340.20320.0096
ϕ 1 −0.7690<0.0001−0.7845<0.0001
Ξ 0.3879<0.00010.35280.0015
α 1 0.3772<0.00010.3383<0.0001
β 1 0.6915<0.00010.6082<0.0001
γ 1 --−0.1173<0.0001
δ 1 --2.2703<0.0001
Q (5)2.59450.45971.48760.6851
Q (10)6.35770.63295.92940.6551
Q (20)22.52640.235421.66940.2470
Q 2 (5)1.64480.64670.81650.8455
Q 2 (10)1.72680.99371.00420.9982
Q 2 (20)12.73760.813210.73210.9054
ARCH(5)1.63340.89560.81490.9761
ARCH(10)1.73260.99910.98851.0000
AIC−7.0437−7.0664
Table 7. ML parameter estimates for the LM-GARCH-type models with different error distributions and diagnostic tests (cotton returns).
Table 7. ML parameter estimates for the LM-GARCH-type models with different error distributions and diagnostic tests (cotton returns).
FIGARCH ( 1 , ξ , 1 ) HYGARCH   ( 1 ,   ξ , 1) FIAPARCH   ( 1 , ξ , 1)
InnovationsParametersEstimatep-ValueEstimatep-ValueEstimatep-Value
ND ξ 0.3467<0.00010.04570.89630.2488<0.0001
α 1 −0.9343<0.00010.56340.0190−0.9235<0.0001
β 1 −0.9448<0.00010.32680.04713−0.9383<0.0001
γ 1 ----−0.06770.2540
δ 1 ----2.1374<0.0001
Q (5)3.24990.66753.22970.66863.50560.6261
Q (10)6.67680.65976.95260.73427.16380.7286
Q (20)17.73610.639418.83650.568719.40130.5593
Q 2 (5)0.58670.84450.56900.91460.56240.9317
Q 2 (10)0.86381.00000.77241.00000.66861.0000
Q 2 (20)1.33541.00001.18681.00001.21321.0000
ARCH(5)0.58770.98990.53560.99460.45891.0000
ARCH(10)0.86121.00000.74551.00000.66741.0000
AIC7.0268−7.0270−7.0136
StD ξ 0.5323<0.00010.5397<0.00010.3943<0.0001
α 1 0.2699<0.00010.2612<0.00010.3505<0.0001
β 1 0.7148<0.00010.7187<0.00010.6526<0.0001
γ 1 ----0.51480.0853
δ 1 ----1.7324<0.0001
Q (5)2.12770.83122.16800.83693.34740.6443
Q (10)3.03320.98063.05120.98144.25490.9357
Q (20)8.35650.98928.27590.98999.45130.9774
Q 2 (5)6.05880.10885.44860.14175.16270.1614
Q 2 (10)11.51060.174410.44780.235210.98140.2149
Q 2 (20)19.02660.390218.18680.443218.18550.4430
ARCH(5)5.99610.30665.38180.37125.11350.4271
ARCH(10)11.5930.313210.52700.395511.30620.3351
AIC−7.0268−7.0384−7.0136
SStD ϕ 0 0.0050.00720.00030.15680.07430.0730
ξ 0.3630<0.00010.03320.92450.2775<0.0001
α 1 0.8943<0.00010.56960.03497−0.9268<0.0001
β 1 0.9549<0.00010.33050.0976−0.9369<0.0001
γ 1 ----−0.0487<0.0001
δ 1 ----2.1201<0.0001
Q (5)4.34760.50453.22140.66583.50250.6233
Q (10)9.27570.6116.91710.73567.17240.7183
Q (20)19.96580.460018.3170.568418.46770.5568
Q 2 (5)0.35660.94780.535130.91980.454420.9287
Q 2 (10)0.46451.00000.74611.00000.66481.0000
Q 2 (20)1.37351.00001.17791.00001.21381.0000
ARCH(5)0.64710.89910.53091.00000.45760.8694
ARCH(10)0.46481.00000.74391.00000.66141.0000
AIC−7.0198−7.0267−7.0227
Table 8. ML parameter estimates for the ARFIMA (1, o , 1)-FIGARCH (1, ξ , 1) models with different error distributions and diagnostic tests (tobacco returns).
Table 8. ML parameter estimates for the ARFIMA (1, o , 1)-FIGARCH (1, ξ , 1) models with different error distributions and diagnostic tests (tobacco returns).
N D StD S S t D
ParameterEstimatep-ValueEstimatep-ValueEstimate p-Value
o −0.3964<0.0001−0.1533<0.0001−0.1531<0.0001
θ 1 −0.36450.00100.26360.01220.27740.0083
ϕ 1 0.39670.0010−0.01740.0329−0.0180<0.0001
β 0 5.5548<0.00012.8488<0.00012.8493<0.0001
ξ 0.2467<0.00010.4127<0.00010.4132<0.0001
α 1 0.6624<0.00010.8503<0.00010.8505<0.0001
β 1 0.4576<0.00010.5445<0.00010.5453<0.0001
Q (5)9.43670.463517.47620.07899.82370.1932
Q (10)11.38280.537919.22760.124114.64590.2064
Q (20)22.5140.123428.30180.035425.81240.0639
Q 2 (5)7.87690.654310.54380.361711.75850.3217
Q 2 (10)14.93110.671715.68500.316216.68930.4418
Q 2 (20)16.67980.748717.65920.438719.37640.0935
ARCH(5)0.46660.33880.98740.41631.54170.2235
ARCH(10)0.57830.75650.79430.65321.09450.4387
AIC−2.5982−3.1060−3.5197
Table 9. Forecasting evaluation metrics of the LM-GARCH-type models combined with different error distributions for crude oil.
Table 9. Forecasting evaluation metrics of the LM-GARCH-type models combined with different error distributions for crude oil.
FIGARCH (1, ξ , 1)FIAPARCH (1, ξ , 1)
Forecasting MeasureNDStDSStDGEDNDStDSStDGED
MSE0.00130.00110.00110.00100.00140.00090.00080.0009
MAE0.00780.00770.00810.00820.008110.00760.00740.0084
TIC0.57230.57260.57280.57270.58230.68280.68590.6832
R 2 0.09210.08450.08650.08540.08470.07850.08840.0875
Table 10. Forecasting evaluation metrics of the LM-GARCH-type models combined with different error distributions for gold.
Table 10. Forecasting evaluation metrics of the LM-GARCH-type models combined with different error distributions for gold.
ARFIMA   ( 1 , σ , 1)-FIGARCH (1, ξ , 1) ARFIMA   ( 1 ,   σ , 1)-HYGARCH (1, ξ , 1)
Forecasting MeasureNDStDSStDGEDNDStDSStDGED
MSE0.00070.00060.00060.00050.00590.00600.00610.0061
MAE0.02160.02200.02210.01850.02600.02700.02700.0275
TIC0.68420.67780.67810.69580.63450.63630.64030.6353
R 2 0.07370.04910.04950.04700.10760.09980.09660.0998
Table 11. Forecasting evaluation metrics of the LM-GARCH-type models combined with different error distributions for tobacco and cotton.
Table 11. Forecasting evaluation metrics of the LM-GARCH-type models combined with different error distributions for tobacco and cotton.
TobaccoCotton
ARFIMA   ( 1 ,   σ , 1)-FIGARCH (1, ξ , 1)HYGARCH (1, ξ , 1)
Forecasting MeasureNDStDSStDGEDNDStDSStDGED
MSE0.00030.00020.0002N/A0.00840.00780.0081N/A
MAE0.00670.00690.00650.00580.00610.0063
TIC0.54910.55280.55290.61830.65740.6574
R20.08210.08450.08650.08320.08770.0830
Table 12. Forecasting evaluation metrics of the LM-GARCH-type models combined with different error distributions for lithium.
Table 12. Forecasting evaluation metrics of the LM-GARCH-type models combined with different error distributions for lithium.
ARFIMA   ( 1 ,   σ , 1)-FIGARCH (1, ξ , 1) ARFIMA   ( 1 ,   σ , 1)-FIAPARCH (1, ξ , 1)
Forecasting MeasureNDStDSStDGEDNDStDSStDGED
MSE0.0005N/A0.0014N/A
MAE0.00120.0015
TIC0.57230.5823
R20.09210.0811
Table 13. Model selection based on the smallest values of forecasting metric values.
Table 13. Model selection based on the smallest values of forecasting metric values.
CommodityBest ModelSelection CriteriaDistribution
Crude oil FIAPARCH 1 , ξ , 1 Lowest AIC (−4.8876), lowest MSE and MAE, higher TIC (69%), R 2 (0.05)SStD
GoldARFIMA (1, o , 1)-FIGARCH (1, ξ , 1)Lowest AIC (−6.4176), significance of parameter estimates, lower MSE and MAE, and higher TIC (70%)GED
CottonHYGARCH (1, ξ , 1)Lowest AIC (−7.0384), lowest MSE and MAE, and favourable TIC (66%) and R 2 (0.0877)StD
LithiumARFIMA (1, o , 1)-FIAPARCH (1, ξ , 1)The only plausible models with parameter convergenceND
TobaccoARFIMA (1, o , 1)-FIGARCH (1, ξ , 1)Lowest AIC (−3.5197), lower MSE and MAE, TIC (55%), and R 2 (0.0865)SStD
Table 14. Out-of-sample ML parameter estimates for the FIAPARCH (1,   ξ , 1) with SStD error distributions and diagnostic tests (crude oil returns).
Table 14. Out-of-sample ML parameter estimates for the FIAPARCH (1,   ξ , 1) with SStD error distributions and diagnostic tests (crude oil returns).
ParametersEstimatep-Value
ξ 0.3005<0.0001
α 1 0.8469<0.0001
β 1 0.76175<0.0001
γ 1 0.121670.0258
δ 1 1.5770<0.0001
Q (5)15.45360.3246
Q (10)19.84210.3923
Q (20)21.08221.0000
Q 2 (5)0.00081.0000
Q 2 (10)0.00181.0000
Q 2 (20)0.00321.0000
ARCH(5)9.36440.1058
ARCH(10)9.60090.4762
AIC: −4.8953
Table 15. Out-of-sample ML parameter estimates for the ARFIMA (1, o , 1)-FIGARCH (1, ξ , 1) with GED innovations and diagnostic tests (gold returns).
Table 15. Out-of-sample ML parameter estimates for the ARFIMA (1, o , 1)-FIGARCH (1, ξ , 1) with GED innovations and diagnostic tests (gold returns).
ParameterEstimatep-Value
ϕ 0 0.0020<0.0001
o 0.0327<0.0001
θ 1 −0.05880.0242
β 0 0.0659<0.0001
ξ 0.4013<0.0001
α 1 0.3246<0.0001
β 1 0.6552<0.0001
Q (5)7.18490.3221
Q (10)15.54370.4743
Q (20)9.37210.6714
Q 2 (5)8.29640.1463
Q 2 (10)12.01750.2160
Q 2 (20)16.64110.5874
ARCH(5)6.32850.1877
ARCH(10)13.25170.3428
AIC: −6.4086
Table 16. Out-of-sample ML parameter estimates for the ARFIMA (1, σ , 1)-FIAPARCH (1, ξ , 1) with ND innovations and diagnostic tests (lithium returns).
Table 16. Out-of-sample ML parameter estimates for the ARFIMA (1, σ , 1)-FIAPARCH (1, ξ , 1) with ND innovations and diagnostic tests (lithium returns).
ParameterEstimatep-Value
o 0.7865<0.0001
θ 1 0.3317<0.0001
ϕ 1 −0.5934<0.0001
ξ 0.3882<0.0001
α 1 0.4162<0.0001
β 1 0.6753<0.0001
γ 1 −0.2318<0.0001
δ 1 2.6435<0.0001
Q (5)1.76830.7231
Q (10)6.47450.6879
Q (20)19.49260.3276
Q 2 (5)1.13260.9375
Q 2 (10)1.53720.9052
Q 2 (20)0.93960.9054
ARCH (5)0.88671.0000
ARCH (10)0.98851.0000
AIC: −7.0660
Table 17. Out-of-sample ML parameter estimates for the HYGARCH (1, ξ , 1) with StD innovations and diagnostic tests (cotton).
Table 17. Out-of-sample ML parameter estimates for the HYGARCH (1, ξ , 1) with StD innovations and diagnostic tests (cotton).
ParametersEstimatep-Value
ξ 0.5385<0.0001
α 1 0.3003<0.0001
β 1 0.7410<0.0001
Q (5)2.16700.8766
Q (10)3.09470.98836
Q (20)7.79261.0000
Q 2 (5)5.4630.2367
Q 2 (10)10.50480.2358
Q 2 (20)17.79190.4642
ARCH(5)5.3770.3765
ARCH(10)10.67510.3736
AIC−7.0378
Table 18. Out-of-sample ML parameter estimates for the ARFIMA (1, o , 1)-FIGARCH (1, ξ , 1) models with SStD innovations and diagnostic tests (tobacco returns).
Table 18. Out-of-sample ML parameter estimates for the ARFIMA (1, o , 1)-FIGARCH (1, ξ , 1) models with SStD innovations and diagnostic tests (tobacco returns).
ParameterEstimate p-Value
o −0.1832<0.0001
θ 1 0.3076<0.0001
ϕ 1 −0.0254<0.0001
β 0 2.8441<0.0001
ξ 0.4053<0.0001
α 1 0.8512<0.0001
β 1 0.5487<0.0001
Q (5)10.81630.1954
Q (10)14.71790.2791
Q 20)25.86740.0601
Q 2 (5)12.16380.3263
Q 2 (10)16.62250.4672
Q 2 (20)22.16820.0855
ARCH(5)2.00360.2421
ARCH(10)1.09450.4387
AIC: −3.5197
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Basira, K.; Dhliwayo, L.; Chinhamu, K.; Chifurira, R.; Matarise, F. Estimation and Prediction of Commodity Returns Using Long Memory Volatility Models. Risks 2024, 12, 73. https://0-doi-org.brum.beds.ac.uk/10.3390/risks12050073

AMA Style

Basira K, Dhliwayo L, Chinhamu K, Chifurira R, Matarise F. Estimation and Prediction of Commodity Returns Using Long Memory Volatility Models. Risks. 2024; 12(5):73. https://0-doi-org.brum.beds.ac.uk/10.3390/risks12050073

Chicago/Turabian Style

Basira, Kisswell, Lawrence Dhliwayo, Knowledge Chinhamu, Retius Chifurira, and Florence Matarise. 2024. "Estimation and Prediction of Commodity Returns Using Long Memory Volatility Models" Risks 12, no. 5: 73. https://0-doi-org.brum.beds.ac.uk/10.3390/risks12050073

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop