Forecasting the Covolatility of Coffee Arabica and Crude Oil Prices: A Multivariate GARCH Approach with High-Frequency Data

Forecasting the covolatility of asset return series is becoming the subject of extensive research among academics, practitioners, and portfolio managers. +is paper estimates a variety of multivariate GARCH models using weekly closing price (in USD/barrel) of Brent crude oil and weekly closing prices (in USD/pound) of Coffee Arabica and compares the forecasting performance of these models based on high-frequency intraday data which allows for a more precise realized volatility measurement. +e study used weekly price data to explicitly model covolatility and employed high-frequency intraday data to assess model forecasting performance. +e analysis points to the conclusion that the varying conditional correlation (VCC) model with Student’s t distributed innovation terms is the most accurate volatility forecasting model in the context of our empirical setting. We recommend and encourage future researchers studying the forecasting performance of MGARCH models to pay particular attention to the measurement of realized volatility and employ high-frequency data whenever feasible.


Introduction
Forecasting the covolatility of asset return series is becoming the subject of extensive research among academics, practitioners, and portfolio managers [1]. is has been used in risk management, derivative pricing and hedging, portfolio selection, and policy making. Similarly, the analysis of volatility spillovers between commodity and asset prices has a profound implication for risk management and portfolio maximization by the government and investors [2]. In view of the current bearish behavior of oil price in the international markets, it is arguably of special interest to study the relationship between oil and coffee prices [2]. Recently, a number of papers have studied the comovement of oil prices with equities, agricultural commodities, and precious metals prices.
Prior studies provide evidence on the connectedness between oil and one or more markets. However, the bulk of these studies has focused on price connectedness. Relatively little research has been devoted to volatility comovement.
is is surprising since price volatilities are the main driving forces in options and commodity future markets. is paper seeks to bridge this gap in the literature. In particular, it undertakes a systematic analysis of the volatility interconnectedness between two important commodity markets-oil and coffee markets-using multivariate GARCH (MGARCH) family of models with the view of providing some rigorous price covolatility forecasting models.
Naturally, different papers in the literature also tend to focus on different commodity pairs when studying volatility interconnectedness. e focus of this work is on the covolatility between Brent crude oil and Coffee Arabica prices in the commodity futures markets which are markets where one can buy specific quantities of a commodity at a specified price with delivery set at a specified time in the future. e choice of these commodities is motivated by the fact that these are the most actively traded commodities in the world. Existing work appears to focus mostly on in-sample modeling of covolatility of commodity prices giving less attention to model forecasting performance. is work will also seek to bridge this gap in the literature.
A number of recent papers have investigated the covolatility of energy and commodity price mainly due to increased interest for understanding the drivers and dynamics of such covolatilities in these highly volatile times. Alternative estimation approaches are used in the literature, but thus far there is a dearth of work undertaking a thorough "MGARCH" analysis-a statistical framework which is particularly suited for modeling asset prices covolatility. Moreover, existing work appears to focus mostly on insample modeling of covolatility of commodity prices giving less attention to model forecasting performance. Accordingly, this work will seek to contribute to the literature by undertaking a systematic analysis of volatility forecasting performance.
Most of the existing research on volatility spillovers employ statistical models in order to estimate realized volatilities which turned out to be oftentimes poor approximations of true volatilities. An attractive alternative to model-based statistical volatility is to compute realized volatility based on high-frequency intraday or "tick" data. Realized volatility has been found to be more accurate than model-based volatilities in predicting latent volatility [3].
us, the results provided by previous studies can be potentially misleading in that they may have underestimated or overestimated the extent of the true volatility.
is work avoids this potential pitfall by using realized volatility based on 30-minute intraday data. We argue that doing so represents an important contribution to the literature of covolatility measurement within commodity markets. e work makes the following contributions to the literature: (1) From an applied econometrics stand point, the work provides a robust and systematic investigation of the statistical and econometric properties of the covolatility of the two important commodities. (2) e direction and magnitude of correlation between the prices of coffee and crude oil which are documented in this work would help the design of a more informed policy making on the part of financial authorities as well as other stakeholders and market participants interested in foreseeing the volatility of these two important commodities.

Multivariate GARCH Models.
A large part of the literature deals with univariate model price series. But markets interact, and therefore a generalization from the univariate model to a multivariate one is important. Multivariate GARCH models can be categorized into four types [4].
(i) Models of the conditional covariance matrix: in these models, the conditional covariance is computed in a direct way. e vector error correction (VEC) and Baba-Engle-Kraft-Kroner (BEKK) models are among this type of models. (ii) Factor models: the return process is assumed to consist of a small number of unobservable heteroskedastic factors. is approach benefits from that the dimensionality of the problem reduces when the number of factors compared to the dimension of the return vector is small. (iii) Models of conditional variances and correlations: at first, the univariate conditional variances and correlations are computed and then used to get the conditional covariance matrix. Some models are, for example, the constant conditional correlation (CCC) model and the dynamic conditional correlation (DCC) model. (iv) Nonparametric and semiparametric approaches: models in this class form an alternative to parametric estimation of the conditional covariance structure. e advantage of these models is that they do not impose a particular distribution (that can be misspecified) on the data.
In this paper, conditional variance and correlation models are used. Before detailed discussion of the MGARCH models, some definitions of the series are presented. is study uses two commodity price series, resulting in a bivariate approach.
where r t is the vector of returns which can be decomposed as where r t is the N × 1 vector of returns; µ t is the N × 1 vector of expected return at time t given the available information set I t−1 ; H 1/2 t N × N is the Cholesky factor of time-varying conditional covariance matrix t ; and z t is the N × 1 vector of independent errors with mean zero and variance one. e variance of the conditional unpredictable component can also be defined as In equation (3), the matrix obtained should be symmetric. Furthermore, the matrix t has to be positive definite for all t. e approach of multivariate modeling brings complications. Firstly, it may be hard to ensure that t is positive definite for all t. Secondly, the number of parameters becomes too large to estimate, and at last, it may be difficult to obtain the stationarity condition for t . e various MGARCH models proposed in the literature differ in how they trade off flexibility and parsimony in their specifications for t . Increased flexibility allows a model to capture more complex t . e following sections will describe the theory of the models, where the VEC and DVEC models which are the predecessors to the CCC, VCC, and DCC models are first described.

e Constant Conditional Correlation (CCC) Model.
e constant conditional correlation model was suggested by Bollerslev [5], where the time-varying covariance matrix t at time t is expressed as e right hand side consists of the conditional correlation matrix R that is time invariant, i.e., R t � R.D is a diagonal matrix of (h 1t, h 2t , . . . , h kt ) such as where each h i,t follows a univariate GARCH process. e conditional correlation matrix is given by R � [ρ i,j ], and the nondiagonal elements of t are e desired conditional variances can be expressed in vector form: e first term C is a vector of the intercepts with a dimension of n × 1 and the matrices of the coefficients are × n. e advantage of the CCC model is that the computational procedure is more easily performed since the correlation matrix R is constant. However, this means that the model may be too restrictive [6].

e Dynamic Conditional Correlation (DCC) Model.
e dynamic conditional correlation model is given as where t is the covariance matrix and Rt is n × n matrix of the conditional correlation of the returns which is symmetric by definition. e diagonal matrix Dt is expressed as where each h i,t follows a univariate GARCH process. Furthermore, t has to be positive definite, which is automatically obtained while Rt is a correlation matrix that is symmetric by definition. When this matrix is defined, two requirements are needed. Firstly, t needs to be positive definite, since it is a covariance matrix. Secondly, the parts that belong to Rt need to be less than one. ese requirements are met through decomposition: where where Q is the unconditional covariance and will be finite as the model contains finite parameters and the vector ε t has finite variance. Additionally, the parameters π 1 and π 2 are scalars and diag (Q) is used to rescale the parts of Qt in order to fulfill that |ρ Moreover, the scalars π 1 and π 2 must be larger than zero, but the sum has to be less than one. One may note that these are conditions of the univariate GARCH to be stationary, but which is applied in the DCC model [7].

e Varying Conditional Correlation (VCC) Model.
e varying conditional correlation (VCC) is a multivariate generalized autoregressive conditionally heteroskedastic (MGARCH) model in which the conditional variances are modeled as univariate GARCH models and the conditional covariance is modeled as nonlinear functions of the conditional variances. e conditional correlation parameters that weight the nonlinear combinations of the conditional variance follow the GARCH-like process specified by Tse and Tsui [6]. e VCC MGARCH model is about as flexible as the closely related dynamic conditional correlation MGARCH model and more flexible than the constant conditional correlation model and more parsimonious than the diagonal vector error conditionally heteroskedastic (DVECH) models.
MGARCH models differ in the parsimony and flexibility of their specifications for a time-varying conditional covariance matrix of the disturbances, denoted by t . In the conditional correlation family of MGARCH models, the diagonal elements of t are modeled as univariate GARCH models, whereas the off-diagonal elements are modeled as nonlinear functions of the diagonal terms. In the VCC MGARCH model,

Journal of Probability and Statistics
where the diagonal elements hii, t and hjj, t follow univariate GARCH processes and ρij, t follows the dynamic process specified in Tse and Tsui [6] and discussed below. Because ρij, t varies with time, this model is known as the VCC GARCH model. e VCC-GARCH(1, 1) model proposed by Tse and Tsui [6] can be written as where x t is a k × 1 vector of independent variables, which may contain lags of r t ; θ is an m × k matrix of parameters; r t is an m × 1 vector of dependent variables; H 1/2 t is the Cholesky factor of the time-varying conditional covariance matrix t ; z t is an m × 1 vector of independent and identically distributed innovations; and D t is a diagonal matrix of conditional standard deviation.
R is the matrix of means to which the dynamic process in equation (14) reverts; π 1 and π 2 are parameters that govern the dynamics of conditional correlations, and they are nonnegative and satisfy 0 ≤ π 1 + π 2 < 1.

Determining the Conditional Distribution.
When fitting a GARCH model based on financial data, the conditional distribution of the returns has to be defined. Studies, for example Bollerslev [8], illustrate that returns are not normally distributed. Instead, the Student t distribution captures the observed kurtosis in empirical returns in a more sufficient way than the normal distribution. Returns have excess kurtosis and fatter tails than the normal distribution. erefore, the Student t distribution is more suitable [9]. ere are three assumptions about the conditional distribution of the error term commonly employed when working with GARCH models: normal (Gaussian) distribution, Student's t distribution, and the generalized error distribution (GED).

Estimation Evaluation.
e evaluation is done with several methods. Firstly, the estimations are compared using the method of three measures: MAE, RMSE, and the explanatory power from an ordinary least squares regression. Secondly, an attempt to determine the goodness of fit of the residuals is done through univariate approach consisting of an analysis with the Ljung-Box test.

Defining a Proxy.
In evaluating volatility forecasts, the usual proxy for "true" volatility is "ex post" squared returns or the squared errors. However, as noted by Andersen and Bollerslev [5], although the use of squared returns or the squared errors is justifiable as an unbiased estimate of the volatility process, it provides a very noisy measure due to a large idiosyncratic term. Specifically, the returns innovation can be written as ε t � σ t z t with z t an independent zero mean and unit variance stochastic process and σ t is the volatility process. If the model for σ 2 t is correctly specified, then the conditional expectation , and the squared error is an unbiased estimate of the volatility process; however, it still contains the noisy idiosyncratic term, z 2 t . is typically resulted in a poor performance, which instigated a discussion of the practical relevance of volatility models.
Hansen et al. [10] provided another important argument for using the realized variance rather than the squared return. ey showed that substituting the squared returns for the conditional variance can severely distort the comparison, in the sense that the empirical ranking of models may be inconsistent for the true (population) ranking. So, an evaluation that is based on squared returns may select an inferior model as the "best" with a probability that converges to one as the sample size increases. For this reason, our evaluation is based on the realized variance. Building upon this line of research, Andersen et al. [11] defined the socalled "realized volatility" on day t as where N is the number of equally spaced intervals within a day and r t,i is a logarithmic return on day t at time interval i with i � 1, 2. . ., N and t � 1,. . .,T. us, realized volatility is the sum of squared intraday returns. In principle, letting N tend to large, i.e., continuous time sampling, the measure approaches to the true integrated volatility of the underlying continuous time process and is theoretically free from measurement error. Further, this measure allows a market participant to essentially treat volatility as an observed variable and to allow direct estimation.

e MZ Regression.
A popular way to evaluate volatility models' out-of-sample forecasting performance is in terms of R 2 from a Mincer-Zarnowitz (MZ) regression, that is, squared returns are regressed on the model forecasts of σ 2 t and a constant. Here, σ 2 t is ex post volatility (e.g., realized volatility) at time t, h 2 t is estimated (in-sample) or forecasted (out-of sample) volatility at time t, and v t is the error term which is independent and identically distributed; v i ∼N (0,1). a and b are parameters to be estimated.
If the model for conditional variance is well specified, we should have a � 0, b � 1. According to specific features of financial data series, the value of R 2 is usually low (even less than 5%) [11]. 4 Journal of Probability and Statistics

Mean Absolute Error (MAE).
Another way of determining the goodness of the estimations and forecasts is calculating the MAE. e approach is to measure how the received conditional covariance is close to its corresponding realized value. e formula is where σ 2 t is used as a proxy and h 2 t is the estimated covariance. By comparing the MAE between the estimated models, it can give an indication of which model that makes the best estimations.

Root Mean Square Error (RMSE).
e third measure is the root mean square error (RMSE), which is defined as Using these methods, the estimated models can be compared, and using the same measurements for estimations and forecasts, one can determine if the relatively best estimation model also makes the best forecast.

Data.
e data considered in this paper were the weekly time series of Brent crude oil and Coffee Arabica futures market closing price, given in US dollars per barrel and US dollars per pound, respectively. Both series contain data spanning between first week of January 2005 and last week of October 2016 extracted from Bloomberg database. e full sample is split in to two parts: in-sample data, in order to estimate the parameters of models, and out-of-sample data, in order to make forecasts. e in-sample period spans from first week of January 2005 up to last week of December 2015 and the out-of-sample period spans from first week of January 2016 through last week of October 2016. For the outof-sample period, we have also extracted 30-minute intraday data for realized volatility measuring purpose.

Test of Stationarity and Features of Log Return Series.
Our data consist of two commodity prices: crude oil and Arabica coffee. e two commodity prices present the same trend and direction during the entire study period. Prices for both of the commodities have risen from 2005 up to 2008 and fallen from 2009 up to 2010 with a similar pattern. Another impression is that both commodity level series are nonstationary (see Figure 1).
To visualize the returns series for these two markets, we depict the return time series plots in Figure 2. e weekly return series display volatility persistence properties, indicating that large changes tend to be followed by large changes of both signs and small changes tend to be followed by small changes, and the differenced series suggests stationarity. e plots further show that the volatilities of these two commodity returns not only have a volatility clustering phenomenon during the selected sample period but also have certain relation on their return volatility processes. at is, when the fluctuation of the Arabica coffee price grew larger, the volatility of crude oil return also became larger for most of the times. is is the main motive for discussing the relationships of coffee price returns and crude oil price returns. Table 1 summarizes the DF test results of the two commodity level prices, respectively. As can be seen in the table, p-value is greater than 0.05. Hence, the null hypothesis of unit root would not be rejected, that is, there is a unit root problem in each of the series.
If time series data are nonstationary, it is necessary to look for possible transformations that might bring stationarity. In practice, econometricians usually transform financial prices in to return forms. is is because often return series are found to be stationary such that analysis is possible. e log return series is obtained by where r t is the log return series of the real price multiplied by 52 (which is simply a scaling factor), to annualize as we are using weekly data and each year contains 52 weeks, and p t is the originalprice series and log is the natural logarithm in base ten. Table 2 summarizes the DF test of the log return series for each of the commodities. e table shows that both p-values are too small (p-value <0.05) at 5% level of significance. ese indicate that the null hypothesis of unit root would be rejected in both cases. Hence, the log return series are stationary for each of the items as shown by the time series plots.
Looking at Table 3, the average prices for Arabica coffee (per Pound) and crude oil (per Barrel) were $145.8668 and $80.9327, respectively, over the entire period. e average log returns were then computed, and they are 0.0458 and 0.0119 for coffee and oil, respectively. e average returns for crude oil price were smaller than those of Arabica coffee price return in the sample period considered. One more impression here is that the two commodities' return is positive, which is an implication of the increasing trend in the two commodity price. e estimate of the standard deviation is 2.229 for coffee and 2.3639 for oil which is nearly similar. Hence, we can say that they exhibit approximately similar variations in their prices.

Volatility Modeling.
To build a volatility model for the log return series, the first step is specifying the mean equation. Once we specify the mean equation, we have to test for ARCH effects using the residuals of the mean equation. If ARCH effects are statistically significant, specifying a volatility model and carrying out a joint estimation of the mean and volatility equations are necessary. e conditional mean specification is, in general, arbitrary for GARCH models of the conditional volatility.

Test for ARCH Effects.
Considering a financial return model and on fitting this model, if there is no volatility clustering in each of the returns, the random disturbance term ε t should be a white noise process. e standardized residual plot can be an initial insight to judge the heteroskedastic characteristics of the error term. e standardized residual plot from each commodity return series is shown in Figure 3. e figure depicts the residual plots of the two commodity return series generated from the mean equation.
We see from the figure that for both of the series, there is a prolonged period of low volatility and prolonged period of high volatility. For example, for Arabica coffee return, there is a long period of low volatility from the first week of Crude oil return also exhibits a prolonged period of low volatility than the coffee series, which is from the first week of 2005 to the first week of 2009. In other words, periods of high volatility are followed by periods of high volatility and periods of low volatility tend to be followed by periods of low volatility, which is known as volatility clustering.
is suggests that the residuals or error terms are conditionally heteroskedastic and can be represented by GARCH models.  Year in weeks    Table 4 shows the results of ARCH LM test for the two commodity returns. e last column of the table includes the p-values that indicate rejection of the null hypothesis that "there is no ARCH effect" up to the fourth lag at 5% level of significance. e results indicate that the two commodities' price log return series are volatile and need to be modeled using GARCH models.

Multivariate GARCH Modeling.
At the beginning of reviewing different formulations of MGARCH models, one should consider what specification of a MGARCH model should be imposed in contrast to the univariate case. On the one hand, it should be flexible enough to state the dynamics of the conditional variances and covariance. On the other hand, as the number of parameters in a MGARCH model increases rapidly along with the dimension of the model, the specification should be parsimonious to simplify the model estimation and also reach the purpose of easy interpretation of the model parameters. However, parsimony may reduce the number of parameters, in which situation the relevant dynamics in the covariance matrix cannot be captured.
A number of alternative MGARCH formulations for modeling long-run conditional heteroskedasticity are extensively used in relevant literature. e most widely used are the time-varying conditional correlation (VCC), the constant conditional correlation (CCC), and the dynamic conditional correlation (DCC) models [5]. A more elaborate discussion of these models is presented in section2.1 of chapter two . One complexity of MGARCH models is that the number of parameters grows rapidly and researchers often use a single ARCH and GARCH term in the model specification, and often it is enough in capturing the variation. So, in this paper, we use bivariate-GARCH (1,1) models.
One of the major objectives of this study is choosing an appropriate multivariate GARCH model from the competing models. Using the Mincer-Zarnowitz regression and the forecast error measures, a model with minimum forecast error measure value and maximum MZ R 2 value is selected as an appropriate model. We can chose a single bivariate GARCH model and then estimate the parameters of the selected model. For the distribution of the innovation term, the multivariate normal distribution and the multivariate version of Student's t distribution with eight degree of    Journal of Probability and Statistics 7 freedom are considered. Like the previous analysis, we also consider two volatility measures, that is, the weekly return square and the 30-minute intraday return square. Looking at Tables 5 and 6 for each of the models, the R 2 obtained from Mincer-Zarnowitz regression (R 2 ) value is larger for the high-frequency data. Moreover, the R 2 is larger under Student t distribution assumption of the innovation term. Now combining the high-frequency data and Student t distribution, the time-varying conditional correlation (VCC) model is a better model for the variance covariance estimation of the Arabica coffee return and crude oil returns. e results in Table 7 of the variance equations actually indicated volatility effects in the returns series across the sample period considered. With regard to the estimates of the time-varying correlation (ρ coffee,oil ), the correlation was found to be positive and significant. is implies positive relationship between the conditional volatility series of the two commodity series and indicates that the returns on these commodities rise or fall together. e estimates of the parameters π 1 and π 2 , that is, the parameters dealing with the effects of previous shocks and previous dynamic conditional correlations on the current dynamic conditional correlation, were significantly different from zero, with π 1 � 0.0291 and π 2 � 0.919. e time-varying conditional correlation (VCC) model reduces to constant conditional correlation (CCC) model when π 1 � π 2 � 0 [1]. Hence, employing CCC is too restrictive for these series.
Concerning the estimates of persistence (α i + β i ), the persistence estimates were 0.9698 and 0.9866 at coffee and oil markets, respectively, which is closer to one, implying that volatility is persistent in these markets. e sufficient conditions for the variance matrix to be positive definite are   e aim has been modeling and forecasting the covolatility dynamics of weekly time series of Brent crude oil and Coffee Arabica futures market closing price using bivariate GARCH models. From the preliminary analysis over the time period considered, both of the price series show an increasing trend. To determine whether the series are stationary or not, the augmented Dickey-Fuller (ADF) test was carried out. Often, raw data of commodity prices are nonstationary, which was also the case in this study. For both level time series, the tests indicate to the existence of unit root I (1). e first log difference of each time series was considered as stationary (Table 3).
In the two commodity returns, it appears that the ARCH term coefficients were significant at 5% level of significance. e result is an implication of the presence of volatility clustering, i.e., large changes followed by large changes and small changes followed by small changes.
Noemi and Junior [12] investigated the volatility of two tropical commodities, i.e., world coffee and crude oil, and found that the two commodities were volatile and concluded that fluctuations in coffee and cocoa price follow the oil price fluctuation. Our results shows that Arabica coffee and crude oil volatilities are correlated (in fact, not highly). In addition, the covariance of Arabica coffee and crude oil is found to be similar to the variance of crude oil which indicates that the variation in the Arabica coffee price follows crude oil price, a similar finding with Noemi and Junior.
Overall, the positive relationship between these commodity returns is consistent with previous research studies on energy and commodity covolatility (see for example, [13,14]).

Conclusion and Recommendations
Forecasting the volatility dynamics of asset returns has been the subject of extensive research among academics, practitioners, and portfolio managers.
is thesis estimates a variety of multivariate GARCH models using weekly closing price (in USD/barrel) of Brent crude oil and weekly closing prices (in USD/per pound) of Coffee Arabica and compares the forecasting performance of these models based on highfrequency intraday data which allows for a more precise realized volatility measurement. e analysis points to the conclusion that varying conditional correlation (VCC) model with Student's t distributed innovation terms is the most accurate volatility forecasting model in the context of our empirical setting.
We recommend and encourage future researchers studying the forecasting performance of MGARCH models to pay particular attention to the measurement of realized volatility and employ high-frequency data whenever feasible.
We also recommend that public policy makers interested in foreseeing the price volatility of these two major commodities in the context of the Ethiopian economy consider using the information documented in this study as input in their deliberations given that it is based on some robust econometric work and highly appropriate data. e scope of the analysis in this study has been limited to the covolatility between the two commodities. In order to overcome this limitation and provide a more nuanced analysis, it might be profitable for future researchers to consider incorporating stock, currency, and bond market volatilities into the analysis. e potential complexity of such a research agenda notwithstanding, the results are likely to be rewarding in light of the deeper integration between global financial and commodity markets in recent years, a phenomenon which came to be known as the financialization of commodities.

Data Availability
e data used to support the findings of the study are available from the corresponding author upon request.