Next Article in Journal
Refined Measures of Dynamic Connectedness based on Time-Varying Parameter Vector Autoregressions
Previous Article in Journal
Autoencoder-Based Three-Factor Model for the Yield Curve of Japanese Government Bonds and a Trading Strategy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Transmissions and Volatility Spillovers between Global Price and U.S. Producer Price in Agricultural Markets

Department of Economics, Setsunan University 17-8 Ikedanaka, Neyagawa 572-8508, Japan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
J. Risk Financial Manag. 2020, 13(4), 83; https://0-doi-org.brum.beds.ac.uk/10.3390/jrfm13040083
Submission received: 20 March 2020 / Revised: 14 April 2020 / Accepted: 22 April 2020 / Published: 23 April 2020
(This article belongs to the Section Applied Economics and Finance)

Abstract

:
A considerable number of studies have examined the relationship between global prices and local prices in food-importing nations, but the linkages between international prices and the producer prices of large agricultural exporters have been largely ignored. This paper analyzes the connections between world prices and U.S. producer prices in the wheat, soybeans, and corn markets using a vector error correction generalized autoregressive conditional heteroscedastic model with a multivariate Baba-Engle-Kraft Kroner specification (VECM-GARCH-BEKK) and cross-correlation function (CCF). Our findings indicate firstly that a long-run equilibrium relationship exists between international and U.S. producer prices for the three agricultural crops. It also finds a significant bidirectional causality-in-mean and causality-in-variance between international and U.S. producer prices for these crops. Finally, the empirical results suggest that international wheat and corn prices play a leading role in U.S. local markets in return transmissions and that U.S. wheat price can be considered to be a leading indicator of the global wheat price in volatility transmissions.

1. Introduction

A vast body of literature has sought to delineate the price relationships between miscellaneous markets. For instance, Kouyaté and von Cramon-Taubadel (2016) showed that approximately 500 recently published papers could be found with the search phrase “price transmission” on AgEcon-Search. Many of these extant research projects analyze market connectivity within developing nations (e.g., Van Campenhout 2007; Myers 2008; Lutz et al. 2006; Rashid 2004; Abdulai 2000; Baulch 1997), while relatively fewer articles scrutinize price co-movements from global to local markets (e.g., Minot 2011; Mundlak and Larson 1992; Quiroz and Soto 1995; Conforti 2004; Robles and Torero 2010; Guo and Tanaka 2019), and most of these studies employed error correction models (ECMs). Ceballos et al. (2017) produced the first work on international price volatility passthroughs in the cereal industry and used a bivariate T-GARCH-BEKK model to investigate price volatility transmissions from the world market to grain-importing regional markets. Thus, past studies on this issue have highlighted food-importing nations rather than food suppliers, despite the fact that food suppliers are a vital component in pricing as well. For example, when U.S. soybean exports to China fell by 75% from September 2018 to May 2019 compared with the same nine-month period in the previous year due to the escalation of the U.S.-China trade war, some farming operations in South Dakota were destroyed by the lower producer price transmitted via the international price (Plume 2019; Newburger 2019). Such regulations in international or domestic markets can worsen market efficiency, generating extra margins between international, import or export, wholesale, and retail prices.
Copious literature inspects the market value chain and the vertical price transmission of agricultural commodities, namely the relationships between farmgate, wholesale, retail, and international prices. Any exporting firm needs to bear a fixed cost for building a distributional network overseas, and wholesalers connect farming operators to importers (Akerman 2010). Thus, the value chain structure takes several steps to deliver products to consumers, and the price connectivity of each stage has been studied by many articles with ECMs, many of which concentrate on the symmetry of price behavior (e.g., Rezitis and Tsionas 2019; Zingbagba et al. 2019; McLaren 2015; Usman and Haile 2017; Abdulai 2002; Worako et al. 2008).
Nonetheless, these studies do not particularize the impacts of the producer price of a large exporter on the international price or vice versa. Although Ghoshray (2008) analyzes the interconnectivity of rice export prices between two great exporting countries of rice, Thailand, and Vietnam, they do not examine the relationships between international and domestic markets. Wheat, corn, and soybean are vitally important commodities, but there is no information available on whether the producer prices of the U.S., one of the largest exporters, lead international prices or whether the causal relationships involving these crops are uni- or bi-directional with regards to the markets.
There is a growing demand for investigations of producer prices and relevant factors resulting from climate change, which exacerbates the fluctuations in agricultural crop production and jeopardizes global food security. Climate change in the form of extreme weather patterns has resulted in greater yield variability in recent years (Ünal et al. 2018). For example, the poor harvests in Australia and Ukraine in 2007 may have worsened the 2008 food crisis, and Russia’s poor harvest in 2010 unsettled the global wheat market (Tanaka et al. 2012; Welton 2011). To make matters worse, future climate change will be aggravated by the migration of crop production to other areas (Reimer and Li 2009). However, past studies have utterly underestimated the importance of analyzing the interconnection between producer prices of countries producing large amounts of exports and global prices in the agricultural sector, particularly in the light of the fact that an ample number of published articles have found price spillover effects from international to local food-importing markets. Moreover, it has not been fully elucidated as to whether or not international price is actually transmitted to producer price in countries exporting large amounts of agricultural products despite the fact that the opposite directional causality, namely from domestic prices in exporting countries to global prices. For instance, it is understood intuitively from historical events such as the poor Russian wheat harvest in 2010.
This article contributes to establishing the relationship between producer prices in the U.S., one of the greatest agricultural exporters, and international crop prices by employing rigorous econometric methods. In particular, we identify the magnitude and speed of price transmissions, both from global to regional markets and from local to world markets. There is a remarkably limited number of articles examining international–producer price linkages, despite the vast literature on price transmissions from global to food-importing markets. The findings from our analysis can be utilized to calm tempestuous global food markets if U.S. farmers could obtain early and precise climate or weather forecasting information to stabilize their crop production. In addition, to this point, it is not fully understood whether exporters’ producer prices are uni-directionally or bi-directionally connected to international prices.
The outline of the paper is as follows. Section 2 provides detailed information on the data used for our experiments. Then the methodology applied is introduced in Section 3. Section 4 explains the estimation results and Section 5 presents the results of the robustness tests. Finally, Section 6 concludes the article.

2. Data and Sample Statistics

We used monthly producer price indices for wheat, soybean, and corn from the U.S. Bureau of Labor Statistics and quoted monthly international prices for these crops from the International Monetary Fund (IMF) Commodity Prices spanning from January 1980 to December 2017, as the IMF Commodity Prices are broadly utilized as global commodity prices.1 The producer prices of wheat, corn, and soybean are indices standardized to 1982 and have been seasonally adjusted. In particular, No. 1 hard red winter wheat from Kansas City, No. 2 yellow corn from the Gulf of Mexico, and No. 2 yellow soybean future contracts from Chicago are used as the international prices of wheat, corn, and soybean, respectively.
Logarithmic transformations of each price series are used in the empirical analysis. More specifically, continuously compounded price returns are computed as ln ( P t / P t 1 ) = ln P t ln P t 1 , where ln P t is the monthly logarithmic price for wheat, soybean, or corn for the international and U.S. producer markets. The developments in continuously compounded returns for each price series are plotted in Figure 1. Examining Figure 1, we make some interesting observations. First, all the price returns exhibit time-varying volatility, i.e., the volatilities are large in one period but small in another period. Second, the volatility tends to present a clustering behavior, i.e., periods of high (low) volatility tend to be followed by periods of high (low) volatility. Third, it is also interesting that some price returns (e.g., wheat and soybean) fluctuated dramatically during the 2007–2008 food crisis.
Table 1 reports the descriptive statistics for the returns on international and U.S. producer prices for wheat, soybean, and corn. The data in this table reveal that most of the mean and median returns are positive during the study period, suggesting a rise in these prices either in the international market or the U.S. local market. Moreover, it is evident that the standard deviation of the U.S. producer price for corn is relatively higher than those for the other commodities, suggesting that extreme changes tend to occur more frequently for this price return. Further, prior to empirical analysis, it is important to examine the stationarity of each price series to ensure their appropriateness for model estimation. The augmented Dickey-Fuller (ADF)2 and Kwiatkowski-Phillips-Schmidt-Shin (KPSS)3 unit root tests are applied to all the price returns. The results of the unit root tests presented in Table 1 indicate that the ADF test rejects the null hypothesis of nonstationary, while the KPSS test cannot reject the null hypothesis of a stationary series for all commodities. These results suggest that all variables are stationary in their first log-differenced forms4 and integrated order 1 (I [1]). Additionally, Van Dijk et al. (2005) suggest pre-testing for structural breaks before examining causality. Following this, Bai and Perron’s (2003) structural change test, which allows for the simultaneous estimation of unknown multiple structural breaks, is applied to identify the structural breakpoints in each wheat price series.5 The results of the structural breaks test in Table 1 show that there are no structural breaks in all the price returns.

3. Methodology

We employed a vector error correction6 generalized autoregressive conditional heteroscedastic model with a multivariate Baba-Engle-Kraft-Kroner specification (VECM-GARCH-BEKK) to scrutinize the associations between the global prices and the U.S. producer prices of major agricultural products. In addition, we used the cross-correlation function as a robustness test that enables us to detect not only directional relationships but also the number of lead-lag relationships as well. First, in order to capture volatility clustering and analyze the volatility spillover effect between the international prices and the U.S. producer prices for wheat, soybean, and corn, the GARCH model with a multivariate BEKK specification was applied. Second, based upon the above results for the unit root tests and cointegration tests, we combined the VECM and asymmetric GARCH-BEKK model. This model has the advantage of analyzing conditional means and variances to detect dynamic interactions between variables and identifying the direction of causality and spillover effects. More specifically, the VECM can explain price-level behavior by explicitly allowing cointegrated relationships based on the assumptions that the price variances are constant over time. On the other hand, the GARCH-BEKK model considers volatility clustering and captures the patterns of volatility transmissions between pairs of prices. Further, following Grier et al. (2004), our analysis was based upon an asymmetric version of the GARCH-BEKK model, which introduces an asymmetric effect on the model.
Since cointegrating relationships are identified between the pairs of prices, we used the VECM to express the conditional mean equation. Then the representation of a cointegrated system is as follows:
p t = π + γ E C t 1 + i = 1 k Φ i Δ p t i + υ t   ,   υ t | Ω t 1 ~ ( 0 , H t )   ,
where p t is a 2 × 1 vector of monthly returns at time t of the international prices and the U.S. producer prices for each of the three agricultural commodities, and π is a 2 × 1 constant vector. E C t 1 is the lagged error correction term of the cointegration relationship, and γ is a 2 × 1 vector of parameters that denote the speed of the price adjustment of each equation to long-run equilibrium. Next, Φ is a 2 × 2 matrix of parameters associated with the lagged returns implying individual and cross-causality in means, and k is the order of the lag length. The term υ t is defined as υ t = ( υ 1 , t ,   υ 2 , t ) and is conditionally normally distributed with mean vector 0, conditional on past information Ω t 1 and a 2 × 2 conditional variance–covariance matrix H t .
Considering the price volatility transmission to positive and negative shocks, we used the asymmetric BEKK model, which has the property that the conditional variance–covariance matrix H t in Equation (1) is positive definite by structure and can be modeled as:
H t = C C + A υ t 1 υ t 1 A + B H t 1 B + D ω t 1 ω t 1 D ,
where C is a 2 × 2 lower triangular matrix, and A is the 2 × 2 matrix of parameters, which measures the extent to which conditional variances are correlated with past squared errors. B is also a 2 × 2 matrix of parameters, this time indicating how the past conditional variance H t 1 affects the current levels of the conditional variance. The last item D ω t 1 ω t 1 D on the right-hand side introduces the asymmetric effect on the model. Next, ω t 1 = I [ u t 1 < 0 ] u t 1 , where I [ . ] is an indicator function equal to 1 if u t 1 < 0 and 0 otherwise, and “ ” is the Hadamard product. Finally, D is a 2 × 2 parameter matric that measures the potential asymmetric response. Thus, if we expanded Equation (2), the conditional variance equation of the bivariate GARCH (1, 1) model can be represented by
[ h 11 , t h 12 , t h 21 , t h 22 , t ] = [ c 11 0 c 21 c 22 ] [ c 11 0 c 21 c 22 ] + [ a 11 a 12 a 21 a 22 ] [ υ 1 , t 1 2 υ 1 , t 1 υ 2 , t 1 υ 2 , t 1 υ 1 , t 1 υ 2 , t 1 2 ] [ a 11 a 12 a 21 a 22 ]   + [ b 11 b 12 b 21 b 22 ] [ h 11 , t 1 h 12 , t 1 h 21 , t 1 h 22 , t 1 ] [ b 11 b 12 b 21 b 22 ]   + [ d 11 d 12 d 21 d 22 ] [ ω 1 , t 1 2 ω 1 , t 1 ω 2 , t 1 ω 2 , t 1 υ 1 , t 1 ω 2 , t 1 2 ] [ d 11 d 12 d 21 d 22 ] ,
where h i i , t and h j j , t for i, j = 1, 2 represent the conditional variances of the returns of each price i and price j, respectively, at time t. Next, h i j , t is the conditional covariance between price i and price j at time t, and a i i and a j j indicate the degree of impact of the current conditional variance (volatility) of price i and price j, respectively, from their own past squared errors. Here, b i i and b j j represent the degree of impact of the current conditional variance of price i and price j, respectively, from their own past conditional variances. Further, a i j measures the lagged shock transmission from price i to price j, and b i j measures the volatility spillover effect from price i to price j. Next, d i i and d j j represent the degree of impact of the current conditional variance of price i and price j, respectively, from their own past asymmetric shocks (positive or negative), and d i j represents the previous asymmetric shocks transmission from price i to price j. If any element in matrix D is positive and significant, then an asymmetric effect exists and negative shocks (bad news) will cause larger volatility in markets than positive shocks (good news) will. Conversely, a significant negative coefficient in the D matrix implies that bad news may bring less volatility in another market.
Moreover, we used the conditional variances and covariances matrix H t to generate dynamic correlation coefficients so as to capture the time-varying correlations between global prices and U.S. producer prices in the cereal market. The dynamic correlation coefficient σ i j , t between price i and j at time t can be represented as:
σ i j , t = h i j , t h i i , t h j j , t   f o r   i , j = 1 , 2 .
We estimate the parameters of the VECM-GARCH-BEKK model by employing the maximum likelihood estimation method optimized with the Broyden, Fletcher, Goldfarb, and Shanno (BFGS) algorithm. The conditional log-likelihood function L ( θ ) can be expressed as:
L ( θ ) = T ln ( 2 π ) 1 2 t 1 T ln | H t | 1 2 t 1 T υ t H t 1 υ t .
Furthermore, the structures of the causal relationships between variables are investigated through the Granger causality approach. Specifically, to test for causality in means between international and U.S. producer prices, we used the likelihood ratio statistic to check whether the elements ϕ i j of matric Φ in Equation (1) are significantly different from zero. As for the causality in the variance test, the non-diagonal elements a i j , b i j , and d i j of matrices A , B , and D , respectively, in Equation (2) were analyzed via the joint Wald test.

4. Empirical Results

First of all, the results of the cointegration test strongly suggest the existence of a long-run equilibrium relationship between each pair of agricultural commodity prices,7 implying that the three prices move together over time. According to these results, this paper employs the VECM instead of a standard VAR model.
Table 2 reports the coefficient estimates for the conditional mean return equation and the conditional variance–covariance matrix of the VECM-GARCH-BEKK model for wheat. First, in terms of the conditional mean estimations, we can observe that both the parameters ϕ 1 , 1 and ϕ 2 , 2 are statistically significant, which indicates that the returns of the international and U.S. wheat prices are influenced by their own past returns. In contrast, the coefficients ϕ 1 , 2 and ϕ 2 , 1 are statistically significant, which indicates a cross-market (global wheat market and U.S. domestic wheat market) return spillover between lagged and current wheat price returns. In addition, by comparing the absolute values of these coefficients, we can determine that the impacts of the own-market return spillovers are larger than those of the cross-market return spillovers ( | ϕ 1 , 1 | > | ϕ 1 , 2 | and | ϕ 2 , 2 | > | ϕ 2 , 1 | ). Second, from the estimation results for the variance equation, all the coefficients a 1 , 2 , a 2 , 1 , b 1 , 2 , and b 2 , 1 are found to be statistically significant, revealing cross effects. These results suggest that the contemporaneous return and volatility spillovers are bidirectional between international and U.S. producer prices. In other words, lagged shocks and volatility in the international wheat price returns can be exploited to predict future U.S. wheat producer price volatility. Similarly, past changes in the U.S. wheat market could provide valuable information to the global wheat market.
Third, turning to the estimation results for the asymmetric terms d 1 , 2 and d 2 , 1 , we can see that d 1 , 2 shows a significant negative value, suggesting that the volatilities of the U.S. producer price are more sensitive to positive shocks (good news) than negative shocks (bad news) in the global wheat market. On the other hand, d 2 , 1 has a positive significant value, which provides evidence that the arrival of negative shocks in the U.S. wheat market intensifies the global wheat market volatility more than positive shocks of similar magnitude. Finally, the Ljung-Box (Ljung and Box 1978) tests provide evidence that there is no serial correlation in the squared standardized residuals, and the Mcleod-Li test and Lagrange multiplier (LR) test indicate that there is no autoregressive conditional heteroscedastic (ARCH) effect in the wheat model. These diagnostic tests generally support the adequacy of the model specification considered.
Moreover, in order to examine the mean and variance causality between the international and U.S. prices for wheat, we perform the joint Wald tests and present the results in Table 3. The results indicate that all the null hypotheses are rejected at the 1% significance level, suggesting a significant bi-directional causality in mean and causality in variance between international wheat price and U.S. wheat producer price. Additionally, the likelihood ratio tests for the null hypothesis that the parameters of the D matrix in the GARCH specification are zero are conducted to test for the relevance of the asymmetric effects. The results of the tests for asymmetric effects show that the null d 1 , 1 = d 1 , 2 = d 2 , 1 = d 2 , 2 = 0 is clearly rejected at the 1% level, implying strong evidence of asymmetric effects between the global wheat market and the U.S. domestic wheat market.
Next, we examine the results for soybean presented in Table 4. First, using a comparison with the findings for the wheat price, the observed mean spillover effect between global and local soybean prices can be reconfirmed. However, the results for soybean prices indicate that the returns for the U.S. soybean price are influenced more significantly past global soybean prices than by its own price returns ( | ϕ 2 , 1 | > | ϕ 2 , 2 | ). Second, we also find a uni-directional return spillover between the two prices. Specifically, past shocks in the international soybean price affect the present volatility of the U.S. producer price ( a 1 , 2 is significant). On the other hand, lagged shocks in the U.S. producer price do not have an impact on the volatility of international soybean prices ( a 2 , 1 is insignificant). Meanwhile, we can observe that a bi-directional volatility spillover between international soybean price and the U.S. producer price ( b 1 , 2 and b 2 , 1 ) are both significant. Third, a cross-market asymmetric volatility transmission can be identified. Specifically, the positive significant value of d 1 , 2 provides evidence that, compared to positive shocks, negative shocks in the global soybean price cause more volatility in the U.S. producer price. Meanwhile, the negative significant value d 2 , 1 suggests that positive shocks to the U.S. producer price have a larger impact on the volatility of the global soybean price. These results are the opposite of the results for wheat. Moreover, the results of the Ljung-Box test, the Mcleod-Li test, and the LM test indicate that the estimated models adequately fit the data.
Our findings for soybean prices were rechecked with the joint Wald tests reported in Table 5. From Table 5, we can conclude that there exists a bi-directional Granger causality in the means and variances between the international soybean price and U.S. soybean producer price. Further, asymmetry effects are also verified.
Finally, we focus on the empirical results for corn prices presented in Table 6. First, it is interesting to observe that the parameter ϕ 1 , 2 is statistically insignificant in the conditional mean equation, while its counterpart, ϕ 2 , 1 , is statistically significant. These results indicate a one-way asymmetric return spillover from international corn lagged price to current U.S. corn producer price. Similar to the results for soybean prices, past global soybean price returns have a larger impact on the returns for the U.S. corn price than do its own price returns ( | ϕ 2 , 1 | > | ϕ 2 , 2 | ). Second, Table 6 also shows that the coefficients a 1 , 2 , a 2 , 1 , b 1 , 2 , and b 2 , 1 are all statistically significant, which provides evidence of a bi-directional return and volatility spillovers between global corn price and U.S. corn producer price. These results are similar to the estimation results for wheat prices. However, the results concerning the asymmetry effects are different from those for wheat and soybean prices. We can see that the asymmetry term d 1 , 2 is positive and significant, which suggests that negative shocks from the international corn price cause more variations than positive shocks in the U.S. corn producer price. In contrast, the estimated value d 2 , 1 is not significant, indicating that there is no evidence of asymmetric effect and the shocks from the U.S. corn producer price do not have an impact on the international corn price. Further, all the diagnostic tests show the robustness of our results.
In addition, Table 7 reports on the Granger causality tests for corn prices. At the 1% significance level, the results show that bi-directional Granger causality in means and variances is established between the international corn price and the U.S. corn producer price. Similarly, the Wald tests confirm the asymmetric effects between these two price series.
Figure 2 plots the time-varying dynamic correlations between international prices and U.S. producer prices for wheat, soybean, and corn. There are several features worth noting. First, it is interesting to ascertain that all conditional correlations display positive values throughout the entire sample period. This evidence visually confirms that strong co-movements exist between the global prices and the U.S. producer prices in the cereal market. Second, the correlations tend to exhibit extreme variability in some specific time periods. For instance, a substantial decrease in correlation between the global and U.S. soybean markets is apparent in the middle of 1988. Third, the soybean price has the highest average value (0.878) for its dynamic correlation coefficients compared to the wheat (0.751) and corn (0.836) prices. On the other hand, the wheat price has the highest standard deviation (0.071) compared to the soybean (0.039) and corn (0.047) prices.

5. Robustness Check

In the previous section, we found a significant bi-directional causality in mean and causality in variance between international and U.S. producer prices for three crops. This section presents some results on robustness.
Some of the recent studies use multivariate GARCH (MGARCH) models with dynamic conditional correlation (DCC) to analyze volatility spill-over mechanisms.8 Huang et al. (2010) compare the performances of GARCH-BEKK and GARCH-DCC models and document the superiority of the BEKK model over the DCC model due to the large number of parameters treated by the BEKK model. In contrast, Bauwens et al. (2006) argue that the DCC model outperforms the BEKK model based on a comparison of the goodness of fit. Since the results vary across different datasets and the debate still lacks consensus, we first check the robustness of our empirical results by employing the GARCH-DCC model. The econometric framework of the model was formulated as follows:
p t = μ + i = 1 k φ i p t i + u t = μ + i = 1 k φ i p t i + H t z t ,
u t | Ω t 1 ~ N ( 0 , H t )
H t = D t R t D t .
Equation (6) is an autoregression (AR) model k process for p t conditional on the information set Ω t 1 . Specifically, μ = ( μ 1 , μ 2 ) is the vector of conditional means, φ i is the parameter vector, k is the lag lengths of the mean equations, u t = ( u 1 , t , u 2 , t ) is the vector of innovations, H t is a 2 × 2 conditional variance–covariance matrix, z t is a 2 × 1 i.i.d vector of standardized residuals, D t is the diagonal matrix containing the conditional standard deviations on the diagonal, and R t is the conditional correlation matrix given by:
R t = Q t * 1 Q t Q t * 1 ,
where Q t is the conditional correlation matrix of standardized residuals. Q t * 1 = d i a g ( q 11 , t 1 / 2 , q 22 , t 1 / 2 ) and q is the elements of matrix Q t . Moreover, the matrix D t can be obtained by estimating a univariate GARCH (p, q) and EGARCH9 (p, q) model, with h i , t   ( i = 1 , 2 ) on the ith diagonal as follows:
h i , t = ϖ i + p = 1 p i α i , p u i , t p 2 + q = 1 q i β i , q h i , t q ,  
ln ( h i , t ) = ϖ i + p = 1 p i ( α i , p | u i , t p h i , t p | + κ i , p u i , t p h i , t p ) + q = 1 q i β i , q ln ( h i , t q ) ,
where h i , t is a 2 × 1 conditional variance vector of the price series and ϖ i is a 2 × 1 constant vector, the lag lengths of variance equations are represented as p and q. α and β are the parameters of the GARCH and ARCH terms, respectively.
Furthermore, Engle’s (2002) DCC model and Cappiello et al.’s (2006) asymmetric DCC (henceforth, A-DCC)10 model were used to determine the volatility spill-over between international and U.S. producer prices for each pair of crops. According to Engle (2002), the dynamic correlation structure is given as:
Q t = ( 1 ψ ζ ) P ¯ + ψ ν t 1 ν t 1 + ζ Q t 1 ,
where Q t is a symmetric positive definite matrix in Equation (8) and P ¯ is the 2 × 2 unconditional correlation matrix of the standardized residuals ν t . The parameters ψ and ζ are non-negative, with a sum of less than unity. Cappiello et al. (2006) modified the correlation evolution equations in the following expression:
Q t = ( 1 ψ ζ ) P ¯ δ N ¯ + ψ ( ν t 1 ν t 1 ) + ζ Q t 1 + δ ( η t 1 η t 1 ) .
Equation (13) is a standard A-DCC in which asymmetric terms are included. δ is the coefficient of the asymmetry term. N ¯ represents the unconditional matrices of η t = I [ ν t < 0 ] ν t , I [ . ] is an indicator function equal to 1 if ν t < 0 and 0 otherwise, and “ ” is the Hadamard product. Finally, the parameters11 of the DCC and A-DCC models were estimated by employing the Gaussian quasi-maximum likelihood estimation (QMLE)12 with the BFGS13 optimization algorithm.
First, in order to choose the most suitable GARCH-based models, an extensive specification testing procedure is conducted to ensure that the conditional means and variances processed fit the price series well. In this study, each model is estimated by using the maximum likelihood method, and we determine the lag length of the mean and variance equations based on the BIC. Based on the residual diagnostics and BIC, the AR (1)-EGARCH (1, 1) model was chosen for the international prices for wheat and soybean, and the AR (1)-GARCH (1, 2) model was chosen for the international corn price. On the other hand, for the U.S. producer prices, we chose the AR (1)-EGARCH (1, 2) model, the AR (1)-EGARCH (1, 1) model, and the AR (1)-GARCH (1, 1) model for wheat, soybean, and corn producer prices, respectively.
The estimation results for the GARCH-based models for the international prices and the U.S. producer prices for wheat, soybean, and corn are summarized in Table 8 and Table 9, respectively. Each table reports the parameter estimates and their corresponding standard errors. It is noticeable that almost all coefficients of the ARCH and GARCH terms are statistically significant. Meanwhile, the asymmetric terms ( κ i ) capturing the leverage effects are statistically significant for either international price or U.S. producer price for wheat and soybean. These results imply that in these price returns, a negative shock increases volatility more than a positive shock of equal magnitude. Furthermore, the generalized error distribution (GED) is applied in all the models because the GED outperforms the normal distribution. Finally, the results of the diagnostics14 suggest that the selected GARCH-based model specification explains the data well.
Next, we compare the goodness of fit between the GARCH-BEKK and the GARCH-DCC models using the log-likelihood statistics. Table 10 presents the log-likelihoods of three model specifications for wheat, soybean, and corn. We can observe that the BEKK model outperforms the DCC and A-DCC models in terms of higher log-likelihood values. As such, it is reasonable to select the BEKK model as our empirical methodology from the previous section.
Finally, we employ the approach proposed by Cheung and Ng (1996) and Hong (2001) to check the robustness of the causality results in mean and variance reported in the previous section. This procedure is based on the residual cross-correlation function (CCF) and is robust to distributional assumptions. The advantage of using the CCF approach is that it detects not only the direction of causality in means and variances but also the number of lead-lag relationships. Specifically, from the estimated GARCH-based model, we calculate the standardized residuals and the standardized squared residuals for international prices and for U.S. producer prices, respectively. Following Cheung and Ng (1996),15 the sample cross-correlation coefficient at lag k, r ^ ε ξ ( k ) and r ^ u v ( k ) are calculated from the consistent estimates of the conditional means and variances. Under the condition of regularity, we detect the null hypothesis that there is no causality in mean using the following CCF statistic:
T ( r ^ ε ξ ( k 1 )   ,     ,   r ^ ε ξ ( k m ) ) L N ( 0 ,   I m ) ,
and there is no causality in variance using the test statistic, which is given by:
T ( r ^ u v ( k 1 )   ,     ,   r ^ u v ( k m ) ) L N ( 0 ,   I m ) ,
where k 1   ,     ,   k m are m different integers, and L indicates the convergence in distribution. To test a causal relationship at a specified lag k, we compare T r ^ ε ξ ( k ) and T r ^ u v ( k ) with the standard normal distribution. If the test statistic is larger than the critical value for the normal distribution, then we reject the null hypothesis.
The empirical results for wheat are given in Table 11. Lags are measured in months, which range from 1 to 12. From Table 11, we can observe that the test statistic for the null hypothesis that global wheat price does not Granger-cause U.S. producer price in the mean is rejected at lag 1, which indicates that the global price leads the U.S. producer price by approximately one month. Meanwhile, the results display that the causality in mean from the U.S. wheat producer price to the international wheat price exists in lags 3, 10, and 12, respectively. Therefore, there is evidence of a feedback effect in mean between these price pairs. In addition, it is worth noting that mean spillovers from the U.S. producer price to the international price take a longer time than going in the opposite direction in the cross-market information transmission. Turning to causality in variance, the results reveal that the global wheat price Granger-causes the U.S. producer price in variance (in lag 5), with a volatility feedback effect (in lags 2 and 12). These findings imply that the volatility transmission from the U.S. producer price to the global price takes less time than going in the opposite direction.
As for the causality test for soybean prices, the results are presented in Table 12. The estimations simultaneously provide strong evidence that there exists either causality in mean or causality in variance between the international soybean price and the U.S. soybean producer price. In particular, the results indicate significant feedback effects in mean and variance that differ in terms of the lag orders. More specifically, we can observe a bi-directional causality in mean with a one-month lag, while the bi-directional causality in variance has about an 11-month lag.
Finally, Table 13 reports the results for the causality test for corn prices. According to Table 13, first, we can identify that there is a bi-directional causality in mean and variance between the international corn price and the U.S. corn producer price. Second, similar to the results for wheat, the mean price transmission from the U.S. to international corn price displayed a longer delayed reaction (a five-month lag) than going in the opposite direction (a one-month lag). In addition, we found that the bi-directional volatility transmission between the international and U.S. producer corn prices exhibited the same lag length in both directions (approximately nine months).
Overall, these empirical results based on the CCF approach are qualitatively similar to those of our VECM-GARCH-BEKK analysis in the previous section. Therefore, our empirical results can be considered robust.

6. Conclusions

This paper has explored the Granger causal relationships between international and U.S. producer prices for wheat, soybean, and corn using a VECM-GARCH-BEKK model. Our main findings are as follows. First, a long-run equilibrium relationship exists between the international and U.S. producer prices for wheat, soybean, and corn. Second, significant bi-directional causality in mean and causality in variance are found between the international and U.S. producer prices for the three cereal crops. Third, the results generated in the main analysis using the aforementioned model were endorsed by the CCF sensitivity tests. Finally, the results of the CCF suggest that international wheat and corn prices play a leading role in U.S. local markets in return transmission. Meanwhile, U.S. wheat price can be considered a leading indicator of the global wheat price in volatility transmissions.
It would be meaningful to compare the results from Guo and Tanaka (2019) with ours to contribute to the literature because they also estimated price volatility transmissions between world and local markets in the wheat sector with similar methods. One of the elements to be discussed is causal direction. They found uni-directionality between global and domestic markets, while our results indicated bidirectionality, which is likely to have occurred because Guo and Tanaka (2019) targeted wheat importing regions, whereas we concentrated on a large exporter. This is partly endorsed by An et al. (2016), who investigate price volatility transmissions between domestic wheat and wheat flour prices in Ukraine under export restrictions with a GARCH-BEKK model. The work does not directly analyze the linkage from international to domestic prices, but presents the efficacy of export control policies to insulate domestic from global markets, suggesting that one of the largest exporters in the world may also suffer from tempestuous external markets, which, namely, implies a bi-directional relationship, assuming that Ukraine’s export is great enough to influence the world’s wheat markets.
Another point of view regards transmission speed. While global wheat price led local prices by one month in our analysis, the counterpart in Guo and Tanaka (2019) estimated it to be five months. This could be caused by the difference between the farm gate price of wheat and the retail price of wheat flour. To deliver wheat flour to consumers, imported wheat needs to be processed at milling factories after trading wheat in wholesale markets, and then to retail stores through wholesale markets of wheat flour. Therefore, the results from Guo and Tanaka (2019) may take five times as long as ours (spillover from international to producer prices).
We found that mean causality from the U.S. producer price to the global price for wheat, corn, and soybean takes more months compared with the opposite direction, which takes just one month. This finding could imply that the U.S. government could have spare time to impose export controls if the country faced a severe poor harvest of those crops to calm domestic markets with increased supply as well as to prevent price transmissions to world markets, which would be beneficial for importing regions of the crops.
It was found that global price also affects prices in significant exporting countries. This finding implies that large exporters need to protect their own local markets from stormy global markets, which is particularly important for low-income households even in advanced economies such as that of the U.S. Our experiments focused on producer price, but consumer or retail price are correlated with producer price movements. As a matter of fact, it is proved that international price volatility of wheat can be transmitted to retail wheat flour prices in exporting countries (Tanaka and Guo 2020). Hence, large crop-exporting governments such as the U.S. government also need to make policy interventions to segregate retail from foreign markets.
We uncovered that U.S. producer price volatility leads to international price volatility in the wheat, soybean, and corn markets, implying that the steadier the U.S. producer prices, the steadier the world prices. The level of productivity (good or bad crops) is the most significant factor leading to changes in producer prices. Therefore, the stability of agricultural production in large exporters is crucial in calming international markets. Hence, investments in weather or climate information sectors with regards to farming could alleviate or lessen price volatility in global markets, consequently improving food security in food-importing countries.
Despite highlighting our analysis on the interlinkages between producer and global prices for three crops, there exists a complex supply chain structure in those sectors. Haile et al. (2017) scrutinized the relationships between the international, wholesale, and retail prices of wheat, flour, and bread in Ethiopia. By encompassing additional explanatory variables such as wholesale price in a model, more detailed connections between markets can be delineated, which is left for future research.

Author Contributions

J.G. analyzed and interpreted the data regarding price volatility transmission from international to local markets in U.S. agricultural markets. T.T. performed the background of the paper and discussed the policy implications of the empirical results. All authors read and approved the final manuscript.

Funding

The research of the second author is in part supported by a Grant-in-Aid from the Japan Society for the Promotion of Science (Grant Number [A] 18K14533).

Acknowledgments

We deeply appreciate anonymous referees for their insightful comments.

Conflicts of Interest

The authors declare that they have no competing interests.

Data Availability

The datasets generated and analyzed in this study are available in the Global Information and Early Warning System repository: http://www.fao.org/giews/en/. The data that support the findings of this study can be obtained from the corresponding author upon request.

Appendix A

Table A1. Johansen’s cointegration test for wheat.
Table A1. Johansen’s cointegration test for wheat.
Null HypothesisMaximum Eigenvalue TestTrace Test
Statistic5% Critical Valuep-ValueStatistic5% Critical Valuep-Value
None20.50315.8920.00925.30520.2620.009
At most 14.8029.1650.3064.8029.1650.306
Notes: The lag length of the model is selected as 2 based on the BIC. Test statistics are based on the specifications provided by the constant term. The restriction test on the cointegrating vector is done via the likelihood ratio test. The critical values are computed using Osterwald-Lenum (1992). The p-values are based on MacKinnon-Haug-Michelis’s p-values (Mackinnon et al. 1999).
Table A2. Johansen’s cointegration test for soybean.
Table A2. Johansen’s cointegration test for soybean.
Null HypothesisMaximum Eigenvalue TestTrace Test
Statistic5% Critical Valuep-ValueStatistic5% Critical Valuep-Value
None65.68515.8920.00071.16920.2620.000
At most 15.4859.1650.2345.4859.1650.234
Notes: The lag length of the model is selected as 2 based on the BIC. Test statistics are based on the specifications provided by the constant term. The restriction test on the cointegrating vector is done via the likelihood ratio test. The critical values are computed using Osterwald-Lenum (1992). The p-values are based on MacKinnon-Haug-Michelis’s p-values (Mackinnon et al. 1999).
Table A3. Johansen’s cointegration test for corn.
Table A3. Johansen’s cointegration test for corn.
Null HypothesisMaximum Eigenvalue TestTrace Test
Statistic5% Critical Valuep-ValueStatistic5% Critical Valuep-Value
None65.68515.8920.00071.16920.2620.000
At most 15.4859.1650.2345.4859.1650.234
Notes: The lag length of the model is selected as 2 based on the BIC. Test statistics are based on the specifications provided by the constant term. The restriction test on the cointegrating vector is done via the likelihood ratio test. The critical values are computed using Osterwald-Lenum (1992). The p-values are based on MacKinnon-Haug-Michelis’s p-values (Mackinnon et al. 1999).

References

  1. Abdulai, Awudu. 2000. Spatial Price Transmission and Asymmetry in the Ghanaian Maize Market. Journal of Development Economics 63: 327–49. [Google Scholar] [CrossRef]
  2. Abdulai, Awudu. 2002. Using Threshold Cointegration to Estimate Asymmetric Price Transmission in the Swiss Pork Market. Applied Economics 34: 679–87. [Google Scholar] [CrossRef]
  3. Akerman, Anders. 2010. A Theory on the Role of Wholesalers in International Trade based on Economies of Scope. Working Paper. Stockholm: Department of Economics, Stockholm University. [Google Scholar]
  4. An, Henry, Feng Qiu, and Yanan Zheng. 2016. How Do Export Controls Affect Price Transmission and Volatility Spillovers in the Ukrainian Wheat and Flour Market? Food Policy 62: 142–50. [Google Scholar] [CrossRef]
  5. Bai, Jushan, and Pierre Perron. 2003. Computation and Analysis of Multiple Structural Change Models. Journal of Applied Econometrics 18: 1–22. [Google Scholar] [CrossRef] [Green Version]
  6. Basher, Syed Abul, and Perry Sadorsky. 2016. Hedging Emerging Market Stock Prices with Oil, Gold, VIX and Bonds: A Comparison between DCC ADCC and GO-GARCH. Energy Economics 54: 235–47. [Google Scholar] [CrossRef] [Green Version]
  7. Baulch, Bob. 1997. Transfer Costs, Spatial Arbitrage and Testing for Food Market Integration. American Journal of Agricultural Economics 79: 477–87. [Google Scholar] [CrossRef] [Green Version]
  8. Bauwens, Luc, Sébastien Laurent, and Jeroen VK Rombouts. 2006. Multivariate GARCH Models: A Survey. Journal of Applied Econometrics 21: 79–109. [Google Scholar] [CrossRef] [Green Version]
  9. Bollerslev, Tim. 1986. Generalized Autoregressive Conditional Heteroskedasticity. Journal of Econometrics 31: 307–27. [Google Scholar] [CrossRef] [Green Version]
  10. Bollerslev, Tim, and Jeffrey M. Wooldridge. 1992. Quasi-Maximum Likelihood Estimation and Inference in Dynamic Models with Time-Varying Covariances. Econometric Reviews 11: 143–72. [Google Scholar] [CrossRef]
  11. Cappiello, Lorenzo, Robert F. Engle, and Kevin Sheppard. 2006. Asymmetric Correlations in the Global Equity and Bond Returns. Journal of Financial Econometrics 4: 537–72. [Google Scholar] [CrossRef]
  12. Ceballos, Francisco, Manuel A. Hernandez, Nicholas Minot, and Miguel Robles. 2017. Grain Price and Volatility Transmission from International to Domestic Markets in Developing Countries. World Development 94: 305–20. [Google Scholar] [CrossRef]
  13. Cheung, Yin-Wong, and Lilian K. Ng. 1996. A Causality in Variance Test and its Application to Financial Market Prices. Journal of Econometrics 72: 33–48. [Google Scholar] [CrossRef]
  14. Conforti, Piero. 2004. Price Transmission in Selected Agricultural Markets. Commodity and Trade Policy Research Working Paper No 7. Rome: Food and Agriculture Organization. [Google Scholar]
  15. Dickey, David A., and Wayne A. Fuller. 1979. Distribution of the Estimators for Autoregressive Time Series with a Unit Root. Journal of the American Statistical Association 74: 427–31. [Google Scholar]
  16. Engle, Robert. 2002. Dynamic Conditional Correlation: A Simple Class of Multivariate Generalized Autoregressive Conditional Heteroscedasticity Models. Journal of Business & Economic Statistics 20: 339–50. [Google Scholar]
  17. Ghoshray, Atanu. 2008. Asymmetric Adjustment of Rice Export Prices: The Case of Thailand and Vietnam. International Journal of Applied Economics 5: 80–91. [Google Scholar]
  18. Grier, Kevin B., Ólan T. Henry, Nilss Olekalns, and Kalvinder Shields. 2004. The Asymmetric Effects of Uncertainty on Inflation and Output Growth. Journal of Applied Econometrics 19: 551–65. [Google Scholar] [CrossRef]
  19. Guo, Jin. 2018. Co-Movement of International Copper Prices, China’s Economic Activity and Stock Returns: Structural Breaks and Volatility Dynamics. Global Finance Journal 36: 62–77. [Google Scholar] [CrossRef]
  20. Guo, Jin, and Tetsuji Tanaka. 2019. Determinants of International Price Volatility Transmissions: The Role of Self–Sufficiency Rates in Wheat–Importing Countries. Palgrave Communications 5: 124. [Google Scholar] [CrossRef]
  21. Haile, Mekbib G., Matthias Kalkuhl, Bernardina Algieri, and Samuel Gebreselassie. 2017. Price Shock Transmission: Evidence from the Wheat-Bread Market Value Chain in Ethiopia. Agricultural Economics 48: 769–80. [Google Scholar] [CrossRef]
  22. Hong, Yongmiao. 2001. A Test for Volatility Spillover with Applications to Exchange Rates. Journal of Econometrics 103: 183–224. [Google Scholar] [CrossRef]
  23. Huang, Yiyu, Wenjing Su, and Xiang Li. 2010. Comparison of GARCH-BEKK and GARCH-DCC models: An empirical study. Advanced Data Mining and Applications. ADMA 2010. In Lecture Notes in Computer Science. Berlin/Heidelberg: Springer, pp. 99–110. [Google Scholar]
  24. Johansen, Søren, and Katarina Juselius. 1990. Maximum Likelihood Estimation and Inference on Cointegration––With Applications to the Demand for Money. Oxford Bulletin of Economics and Statistics 52: 169–210. [Google Scholar] [CrossRef]
  25. Kouyaté, Carolin, and Stephan von Cramon-Taubadel. 2016. Distance and Border Effects on Price Transmission: A Meta-Analysis. Journal of Agricultural Economics 67: 255–71. [Google Scholar] [CrossRef] [Green Version]
  26. Kwiatkowski, Denis, Peter CB Phillips, Peter Schmidt, and Yongcheol Shin. 1992. Testing the Null Hypothesis of Stationarity Against the Alternative of a Unit Root. Journal of Econometrics 54: 159–78. [Google Scholar] [CrossRef]
  27. Ljung, Greta M., and George E. P. Box. 1978. On Measuring of Lag of Fit in Time Series Models. Biometrics 67: 297–303. [Google Scholar] [CrossRef]
  28. Lutz, Clemens, W. Erno Kuiper, and Aad van Tilburg. 2006. Maize Market Liberalisation in Benin: A Case of Hysteresis. Journal of African Economies 16: 102–33. [Google Scholar] [CrossRef]
  29. Mackinnon, Andrew, Anthony F. Jorm, Helen Christensen, Ailsa E. Korten, Patricia A. Jacomb, and Bryan Rodgers. 1999. A Short Form of the Positive and Negative Affect Schedule: Evaluation of Factorial Validity and Invariance Across Demographic Variables in a Community Sample. Personality and Individual Differences 27: 405–16. [Google Scholar] [CrossRef]
  30. McLaren, Alain. 2015. Asymmetry in Price Transmission in Agricultural Markets. Review of Development Economics 19: 415–33. [Google Scholar] [CrossRef]
  31. Minot, Nicholas. 2011. Transmission of World Food Price Changes to Markets in Sub–Saharan Africa. IFPRI Discussion Paper 1059. Washington, DC: International Food Policy Research Institute. [Google Scholar]
  32. Mundlak, Yair, and Donald F. Larson. 1992. On the Transmission of World Agricultural Prices. The World Bank Economic Review 6: 399–422. [Google Scholar] [CrossRef]
  33. Myers, Robert J. 2008. Evaluating the Efficiency of Inter-Regional Trade and Storage in Malawi Maize Markets. Report for the World Bank. East Lansing: Michigan State University. [Google Scholar]
  34. Nelson, Daniel B. 1991. Conditional Heteroskedasticity in Asset Returns: A New Approach. Econometrica 59: 347–70. [Google Scholar] [CrossRef]
  35. Newburger, Emma. 2019. “Trump is ruining our markets”: Struggling farmers are losing a huge customer to the trade war—China. CNBC. Available online: https://www.cnbc.com/2019/08/10/trump-is-ruining-our-markets-farmers-lose-a-huge-customer-to-trade-war----china.html (accessed on 28 August 2019).
  36. Newey, Whitney K., and Kenneth D. West. 1994. Automatic Lag Selection in Covariance Matrix Estimation. Review of Economic Studies 61: 631–53. [Google Scholar] [CrossRef]
  37. Osterwald-Lenum, Michael. 1992. A Note with Quantiles of the Asymptotic Distribution of the Maximum Likelihood Cointegration Rank Test Statistics. Oxford Bulletin of Economics and Statistics 54: 461–72. [Google Scholar] [CrossRef]
  38. Phillips, Peter CB, and Pierre Perron. 1988. Testing for a Unit Root in Time Series Regression. Biometrika 75: 335–46. [Google Scholar] [CrossRef]
  39. Plume, Karl. 2019. On the front lines: Trade war sinks North Dakota soybean farmers. Reuters. Available online: https://www.reuters.com/article/us-usa-trade-china-soybean/on-the-front-lines-trade-war-sinks-north-dakota-soybean-farmers-idUSKCN1VC0ZX (accessed on 28 August 2019).
  40. Quiroz, Jorge, and Raimundo Soto. 1995. International Price Signals in Agricultural Prices: Do Governments Care? Documento de Investigacion 88. Santiago: ILADES Postgraduate Economics Program, Georgetown University. [Google Scholar]
  41. Rashid, Shahidur. 2004. Spatial Integration of Maize Markets in Post–Liberalised Uganda. Journal of African Economies 13: 102–33. [Google Scholar] [CrossRef]
  42. Reimer, Jeffrey J., and Man Li. 2009. Yield Variability and Agricultural Trade. Agricultural and Resource Economics Review 38: 258–70. [Google Scholar] [CrossRef]
  43. Rezitis, Anthony N., and Mike Tsionas. 2019. Modeling Asymmetric Price Transmission in the European Food Market. Economic Modeling 76: 216–30. [Google Scholar] [CrossRef] [Green Version]
  44. Robles, Miguel, and Máximo Torero. 2010. Understanding the Impact of High Food Prices in Latin America. Economia 10: 117–64. [Google Scholar]
  45. Tanaka, Tetsuji, and Jin Guo. 2020. How Does the Self-Sufficiency Rate Affect International Price Volatility Transmissions in the Wheat Sector: Evidence from Wheat-Exporting Countries. Palgrave Communications. forthcoming. [Google Scholar]
  46. Tanaka, Tetsuji, Nobuhiro Hosoe, and Huanguang Qiu. 2012. Risk Assessment of Food Supply: A Computable General Equilibrium Approach. Newcastle upon Tyne: Cambridge Scholars Publishing. [Google Scholar]
  47. Ünal, Gözde Erhan, Baris Karapinar, and Tetsuji Tanaka. 2018. Welfare-at-Risk and Extreme Dependency of Regional Wheat Yields: Implications of a Stochastic CGE Model. Journal of Agricultural Economics 69: 18–34. [Google Scholar] [CrossRef] [Green Version]
  48. Usman, Muhammed A., and Mekbib G. Haile. 2017. Producer to Retailer Price Transmission in Cereal Markets of Ethiopia. Food Security 9: 815–29. [Google Scholar] [CrossRef]
  49. Van Campenhout, Bjorn. 2007. Modeling Trends in Food Market Integration: Method and an Application to Tanzanian Maize Markets. Food Policy 32: 112–27. [Google Scholar] [CrossRef] [Green Version]
  50. Van Dijk, Dick, Denise R. Osborn, and Marianne Sensier. 2005. Testing for causality in variance in the presence of breaks. Economics Letters 89: 193–99. [Google Scholar] [CrossRef] [Green Version]
  51. Welton, George. 2011. The Impact of Russia’s 2010 Grain Export Ban. Oxfam Research Reports. Available online: http://pantheon-live.oxfam.org.nz/sites/default/files/reports/Impact-of-russia-grain-2010-280611-en.pdf (accessed on 22 November 2019).
  52. Worako, Tadesse Kumma, Herman D. van Schalkwyk, Zerihun Gudeta Alemu, and Gezahegne Ayele. 2008. Producer Price and Price Transmission in a Deregulated Ethiopian Coffee Market. Agrekon 47: 492–508. [Google Scholar] [CrossRef]
  53. Zingbagba, Mark, Rubens Nunes, and Muriel Fadairo. 2019. The Impact of Diesel Price on Upstream and Downstream Food Prices: Evidence from São Paulo. Energy Economics 85: 1–45. [Google Scholar] [CrossRef]
1
The Federal Reserve Economic Data also quotes the dataset as world prices.
2
3
4
The ADF and KPSS unit root tests indicate that all the variables have unit root processes in their levels. We do not report these results for the sake of brevity. The results can be obtained from the authors upon request.
5
6
The characteristic of price series allows for the possibility that there is a station ary long-run equilibrium relationship (cointegrating relationship) between individual price series. If there are no cointegrating relationships among the variables, after n (n is equal to the order of integrated variables) differences, then the standard vector autoregressive (VAR) model is employed. In contrast, if cointegrating relationships are identified, the vector error-correction model (VECM) can be used for the empirical analysis. Given that all the price series are nonstationary I (1) series, we will use the Johansen-Juselius procedure (Johansen and Juselius 1990) to examine the cointegrating relationship between international prices and U.S. producer prices for wheat, soybean, and corn.
7
8
9
The GARCH model was developed by Bollerslev (1986) and the exponential GARCH (EGARCH) model by Nelson (1991).
10
The A-DCC model modified the original DCC model by including asymmetries in the correlation dynamics. See Cappiello et al. (2006) for an extensive analysis of these models’ advantages.
11
The estimated results of the DCC coefficients were not reported for the sake of brevity. The results can be obtained from the authors upon request.
12
13
BFGS (Broyden, Fletcher, Goldfarb, and Shanno) is a quasi-Newton optimization method that uses information about the gradient of the function at the current point to calculate where to find a better point.
14
The Ljung-Box and ARCH tests show that there is no autocorrelation up to order 10 for the standard residuals and squared standard residuals and no further ARCH effect in all of the models.
15
We did not choose Hong’s (2001) method since one of our paper’s aims is to detect the number of lead-lag relationships.
Figure 1. Development in returns for wheat, soybean, and corn prices.
Figure 1. Development in returns for wheat, soybean, and corn prices.
Jrfm 13 00083 g001
Figure 2. Conditional correlations estimated with the VECM-GARCH-BEKK model.
Figure 2. Conditional correlations estimated with the VECM-GARCH-BEKK model.
Jrfm 13 00083 g002
Table 1. Summary statistics and unit root tests for price returns.
Table 1. Summary statistics and unit root tests for price returns.
MeanMedianStd. Dev.ADF TestKPSS TestStructural Break Tests
International wheat price0.000−0.0020.060−13.718 *** (1)0.087 (3)No break
U.S. producer price of wheat0.0000.0000.050−11.296 *** (0)0.048 (10)No break
International soybean price0.0010.0010.057−15.497 *** (0)0.041 (2)No break
U.S. producer price of soybean0.0010.0060.069−12.185 *** (2)0.040 (10)No break
International corn price0.0010.0010.058−15.910 *** (0)0.042 (2)No break
U.S. producer price of corn0.0010.0010.075−16.719 *** (0)0.034 (7)No break
Notes: *** denotes rejection of the null hypothesis at the 1% significance level. Numbers in brackets are the lag length and bandwidth. Lag length selection is based on the Bayesian-Schwarz information criterion (henceforth, BIC) in the augmented Dickey-Fuller (ADF) tests. The bandwidth for the KPSS test is determined using the Newey–West bandwidth selection algorithm (Newey and West 1994). We implement all the unit root tests with intercept and trend terms. We used Bai-Perron’s sequential test for the hypothesis of k breaks versus k + 1 breaks, employing the F-statistics. Lag length selection was based on the BIC in the test.
Table 2. Empirical results of the asymmetric VECM-GARCH-BEKK model for wheat.
Table 2. Empirical results of the asymmetric VECM-GARCH-BEKK model for wheat.
ParameterInternational Price (i = 1)U.S. Producer Price (i = 2)
EstimateSEEstimateSE
π i 0.124 ***0.024−0.130 ***0.034
γ i 0.003 ***0.001−0.003 ***0.001
ϕ i , 1 0.213 ***0.0730.217 ***0.062
ϕ i , 2 0.129 *0.0720.255 ***0.062
c i , 1 0.011 ***0.0030.0170.003
c i , 2 0.005 *0.003
a i , 1 0.492 ***0.075−0.276 ***0.096
a i , 2 0.422 ***0.077−0.0740.129
b i , 1 0.978 ***0.031−0.192 ***0.065
b i , 2 0.086 **0.0350.659 ***0.065
d i , 1 0.0160.1200.232 *0.122
d i , 2 −0.351 ***0.1200.721 ***0.118
Q ( 10 ) 6.71215.053
Mcleod-Li (10)11.93911.370
LM test1.3451.077
Notes: SE: standard error; *, **, and ***: statistical significance at 10%, 5% and 1% levels.
Table 3. Joint Wald tests for causality in mean and variance: wheat.
Table 3. Joint Wald tests for causality in mean and variance: wheat.
Causality in MeanCausality in Variance
Null HypothesisChi-Squaredp-ValueNull HypothesisChi-Squaredp-Value
ϕ 1 , 2 = ϕ 2 , 1 = γ 1 = γ 2 = 0
(GP LP)
122.570 (4)0.000 a 1 , 2 = a 2 , 1 = b 1 , 2 = b 2 , 1 = d 1 , 2 = d 1 , 2 = 0
(GP LP)
57.277 (6)0.000
ϕ 2 , 1 = γ 2 = 0
(GP LP)
28.141 (2)0.000 a 1 , 2 = b 1 , 2 = d 1 , 2 = 0
(GP LP)
49.630 (3)0.000
ϕ 1 , 2 = γ 1 = 0
(LP GP)
31.411 (2)0.000 a 2 , 1 = b 2 , 1 = d 2 , 1 = 0
(LP GP)
24.549 (3)0.000
Joint Wald tests for asymmetry effects
d 1 , 1 = d 1 , 2 = d 2 , 1 = d 2 , 2 = 0 65.677 (4)0.000
Notes: GP: Global wheat price; LP: Local (U.S.) wheat producer price. The arrow indicates the direction of Granger causality. The likelihood ratio test statistic for causality testing is a χ 2 ( k ) statistic, where the k within the parentheses is the number of degrees of freedom.
Table 4. Empirical results of the asymmetric VECM-GARCH-BEKK model for soybean.
Table 4. Empirical results of the asymmetric VECM-GARCH-BEKK model for soybean.
ParameterInternational Price (i = 1)U.S. Producer Price (i = 2)
EstimateSEEstimateSE
π i −0.139 ***0.048−0.474 ***0.051
γ i −0.004 ***0.001−0.014 ***0.001
ϕ i , 1 0.778 ***0.0991.065 ***0.103
ϕ i , 2 −0.503 ***0.071−0.760 ***0.079
c i , 1 0.033 ***0.0040.031 ***0.004
c i , 2 −0.0000.003
a i , 1 0.708 ***0.193−0.2340.209
a i , 2 0.477 **0.204−0.1400.218
b i , 1 0.329 **0.1290.211 **0.094
b i , 2 −0.621 ***0.1210.959 ***0.087
d i , 1 0.717 ***0.224−0.599 ***0.213
d i , 2 0.835 ***0.218−0.717 ***0.186
Q ( 10 ) 9.14011.548
Mcleod-Li (10)3.4914.724
LM test0.3780.426
Notes: SE: standard error; ** and ***: statistical significance at 5% and 1% levels.
Table 5. Joint Wald tests for causality in mean and variance: soybean.
Table 5. Joint Wald tests for causality in mean and variance: soybean.
Causality in MeanCausality in Variance
Null HypothesisChi-Squaredp-ValueNull HypothesisChi-Squaredp-Value
ϕ 1 , 2 = ϕ 2 , 1 = γ 1 = γ 2 = 0
(GP LP)
383.601 (4)0.000 a 1 , 2 = a 2 , 1 = b 1 , 2 = b 2 , 1 = d 1 , 2 = d 1 , 2 = 0
(GP LP)
45.141 (6)0.000
ϕ 2 , 1 = γ 2 = 0
(GP LP)
323.960 (2)0.000 a 1 , 2 = b 1 , 2 = d 1 , 2 = 0
(GP LP)
39.103 (3)0.000
ϕ 1 , 2 = γ 1 = 0
(LP GP)
88.770 (2)0.000 a 2 , 1 = b 2 , 1 = d 2 , 1 = 0
(LP GP)
15.950 (3)0.001
Joint Wald tests for asymmetry effects
d 1 , 1 = d 1 , 2 = d 2 , 1 = d 2 , 2 = 0 20.490 (4)0.000
Notes: GP: Global soybean price; LP: Local (U.S.) wheat producer price. The arrow indicates the direction of Granger causality. The likelihood ratio test statistic for causality testing is a χ 2 ( k ) statistic, where the k within parentheses is the degrees of freedom.
Table 6. Empirical results of the asymmetric VECM-GARCH-BEKK model for corn.
Table 6. Empirical results of the asymmetric VECM-GARCH-BEKK model for corn.
ParameterInternational Price (i = 1)U.S. Producer Price (i = 2)
EstimateSEEstimateSE
π i −0.025 ***0.005−0.037 ***0.006
γ i −0.015 ***0.002−0.023 ***0.003
ϕ i , 1 0.159 *0.0830.518 ***0.100
ϕ i , 2 0.0860.064−0.1090.077
c i , 1 0.010 ***0.0030.018 ***0.005
c i , 2 0.008 ***0.002
a i , 1 −0.147 *0.0830.155 **0.064
a i , 2 −0.355 ***0.1260.475 ***0.098
b i , 1 0.928 ***0.005−0.161 ***0.009
b i , 2 0.293 ***0.0400.672 ***0.063
d i , 1 0.215 **0.107−0.0760.084
d i , 2 0.390 **0.1860.0080.170
Q ( 10 ) 3.62010.207
Mcleod-Li (10)4.0408.591
LM test0.4050.881
Notes: SE: standard error; *, **, and ***: statistical significance at 10%, 5% and 1% levels.
Table 7. Joint Wald tests for causality in mean and variance: corn.
Table 7. Joint Wald tests for causality in mean and variance: corn.
Causality in MeanCausality in Variance
Null HypothesisChi-Squaredp-ValueNull HypothesisChi-Squaredp-Value
ϕ 1 , 2 = ϕ 2 , 1 = γ 1 = γ 2 = 0
(GP LP)
190.032 (4)0.000 a 1 , 2 = a 2 , 1 = b 1 , 2 = b 2 , 1 = d 1 , 2 = d 1 , 2 = 0
(GP LP)
397.089 (6)0.000
ϕ 2 , 1 = γ 2 = 0
(GP LP)
132.134 (2)0.000 a 1 , 2 = b 1 , 2 = d 1 , 2 = 0
(GP LP)
66.346 (3)0.000
ϕ 1 , 2 = γ 1 = 0
(LP GP)
37.197 (2)0.000 a 2 , 1 = b 2 , 1 = d 2 , 1 = 0
(LP GP)
369.807 (3)0.000
Joint Wald tests for asymmetry effects
d 1 , 1 = d 1 , 2 = d 2 , 1 = d 2 , 2 = 0 20.250 (4)0.000
Notes: GP: Global corn price; LP: Local corn price. The arrow indicates the direction of Granger causality. The likelihood ratio test statistic for causality testing is a χ 2 ( k ) statistic, where the k within parentheses is the degrees of freedom.
Table 8. Empirical results for international price.
Table 8. Empirical results for international price.
Wheat
AR (1)-EGARCH (1, 1)
Soybean
AR (1)-EGARCH (1, 1)
Corn
AR (1)-GARCH (1, 2)
θ 0.0000.0020.001
φ 1 0.258 ***0.264 ***0.233 ***
ϖ −0.099 *−0.245 *0.001 ***
α 1 0.0540.0360.013 ***
κ 1 0.078 ***0.140 ***-
β 1 0.990 ***0.964 ***1.655 ***
β 2 --−0.969 ***
GED parameter1.354 ***1.195 ***1.156 ***
Q ( 10 ) 9.0878.2709.707
Q 2 ( 10 ) 13.5212.0285.763
ARCH test (10)8.3670.2250.552
BIC−2.964−3.085−3.038
Notes: * and ***: statistical significance at the 10% and 1% levels. The standard error follows Bollerslev and Wooldridge’s (1992) robust standard error. Q ( 10 ) and Q 2 ( 10 ) : test statistics (Ljung and Box 1978) for the null hypotheses of no autocorrelation up to order 10 for standardized residuals and standardized residuals squared. The ARCH test: the Lagrange multiplier test statistic is used to check ARCH effects in residuals (it is distributed as chi-square).
Table 9. Empirical results for the U.S. producer price.
Table 9. Empirical results for the U.S. producer price.
Wheat
AR (1)-EGARCH (1, 2)
Soybean
AR (1)-EGARCH (1, 1)
Corn
AR (1)-GARCH (1, 1)
θ −0.0010.0040.003
φ 1 0.431 ***−0.0010.197 ***
ϖ −2.119 ***−0.153 *0.003 **
α 1 0.299 ***0.0400.157 *
κ 1 −0.076 **0.114 ***
β 1 1.529 ***0.977 ***0.324
β 2 −0.829 ***
GED parameter1.155 ***1.324 ***1.189 ***
Q ( 10 ) 14.56115.32011.486
Q 2 ( 10 ) 13.5218.47914.163
ARCH test (10)1.2940.8821.416
BIC−3.477−2.575−2.414
Notes: *, **, and ***: statistical significance at the 10%, 5%, and 1% levels. The standard error follows Bollerslev and Wooldridge’s (1992) robust standard error. Q ( 10 ) and Q 2 ( 10 ) : test statistics (Ljung and Box 1978) for the null hypotheses of no autocorrelation up to order 10 for standardized residuals and standardized residuals squared. The ARCH test: the Lagrange multiplier test statistic is used to check ARCH effects in residuals (it is distributed as chi-square).
Table 10. Log-likelihood values for the DCC, A-DCC, and BEKK estimated models.
Table 10. Log-likelihood values for the DCC, A-DCC, and BEKK estimated models.
WheatSoybeanCorn
ModelLog-likelihoodLog-likelihoodLog-likelihood
DCC1572.2891571.5071471.103
A-DCC1573.0731573.5011480.206
BEKK1651.099 *1743.960 *1559.138 *
Notes: * denotes the largest value of log-likelihood.
Table 11. Results for the causality test for wheat.
Table 11. Results for the causality test for wheat.
Lag kMean CausalityVariance Causality
International Price → U.S. PriceU.S. Price → International PriceInternational Price → U.S. PriceU.S. Price → International Price
12.394 **1.2930.019−0.093
2−0.292−0.3301.4221.956 *
3−0.730−1.763 *0.3450.093
40.2561.048−0.9590.226
51.1830.3531.649 *0.066
6−1.217−0.5400.794−0.288
70.631−1.041−0.5740.796
8−0.411−1.1621.2150.453
9−0.800−0.2101.448−0.506
10−0.0590.3640.749−1.035
111.1791.911 *0.4091.484
120.745−1.674 *0.4783.141 ***
Notes: *, **, and ***: statistical significance at the 10%, 5%, and 1% levels. The arrow indicates the direction of Granger causality.
Table 12. Results of the causality test for soybean.
Table 12. Results of the causality test for soybean.
Lag kMean CausalityVariance Causality
International Price → U.S. PriceU.S. Price → International PriceInternational Price → U.S. PriceU.S. Price → International Price
15.776 ***−2.821 ***1.3250.157
21.890 *0.7730.7940.677
3−1.441−1.867 *−0.150−0.148
4−0.4000.210−0.764−0.809
5−0.938−1.812 *−0.0360.436
6−0.2560.764−0.002−0.506
7−0.127−1.0350.3771.609
80.222−0.133−0.557−1.549
9−0.282−0.779−0.525-0.203
10−2.222 **−0.6160.3560.984
111.1180.4022.620 ***1.651 *
120.019−0.0341.5641.854 *
Notes: *, **, and ***: statistical significance at the 10%, 5%, and 1% levels. The arrow indicates the direction of Granger causality.
Table 13. Results of the causality test for corn.
Table 13. Results of the causality test for corn.
Lag kMean CausalityVariance Causality
International Price → U.S. PriceU.S. Price → International PriceInternational Price → U.S. PriceU.S. Price → International Price
13.348 ***0.237−0.694−0.182
21.441−0.510−0.3770.392
30.4360.2750.275−0.787
4−1.090−0.0321.4560.135
5−0.607−2.982 ***−0.428−0.519
6−1.816−0.946−0.099−0.176
70.3921.035−1.2360.436
8−2.170 **−1.966 **−0.351−0.707
90.1520.3242.038 **2.951 ***
10−0.6600.533−1.012−0.345
112.277 **1.5770.8812.570 **
120.106−1.3140.6180.070
Notes: ** and ***: statistical significance at the 5%, and 1% levels. The arrow indicates the direction of Granger causality.

Share and Cite

MDPI and ACS Style

Guo, J.; Tanaka, T. Dynamic Transmissions and Volatility Spillovers between Global Price and U.S. Producer Price in Agricultural Markets. J. Risk Financial Manag. 2020, 13, 83. https://0-doi-org.brum.beds.ac.uk/10.3390/jrfm13040083

AMA Style

Guo J, Tanaka T. Dynamic Transmissions and Volatility Spillovers between Global Price and U.S. Producer Price in Agricultural Markets. Journal of Risk and Financial Management. 2020; 13(4):83. https://0-doi-org.brum.beds.ac.uk/10.3390/jrfm13040083

Chicago/Turabian Style

Guo, Jin, and Tetsuji Tanaka. 2020. "Dynamic Transmissions and Volatility Spillovers between Global Price and U.S. Producer Price in Agricultural Markets" Journal of Risk and Financial Management 13, no. 4: 83. https://0-doi-org.brum.beds.ac.uk/10.3390/jrfm13040083

Article Metrics

Back to TopTop