Short-Run Forecasting of Commodity Prices: An Application of Autoregressive Moving Average Models

KE-YOUNG CHU *

Abstract

KE-YOUNG CHU *

KE-YOUNG CHU *

Many time series—whether they are weekly prices, monthly interest rates, or the annual transaction volumes of commodities—have certain regularities, partially obscured by noises. Decomposing a time series into its systematic part and the random part is useful for understanding the nature of the series as well as for predicting its future values. The time-series technique extensively discussed in Box and Jenkins (1970) deals with the identification and estimation of autoregressive integrated moving average (ARIMA) models, which can be used for predicting the systematic parts of time series.

ARIMA models are now used not only for predicting time series but also for analytical purposes. 1 In comparison with forecasting based on reduced forms derived from simultaneous equations models, ARIMA equations are very simple to use and require much less information for forecasting. Furthermore, ARIMA models are not entirely without an economic justification. It was shown by Tinbergen (1939), Zellner and Palm (1974), Palm (1977), and others that the “final equation” in an ARIMA form can be derived for each of the endogenous variables in a simultaneous equations model on a simple set of assumptions.

This paper reviews the relative merits of using ARIMA models for forecasting, compared with reduced forms derived from a simultaneous equations model. The paper also attempts to apply the method of ARIMA models to a short-run simulation of commodity prices. The results suggest that ARIMA models can yield quite satisfactory short-run forecasts for the monthly prices of some commodities and that they can be useful even for forecasting the yearly averages of monthly prices. 2

The paper is divided into three sections. Section I illustrates how the final equation for price in an ARIMA form can be derived from a structural system; it also reviews the relative merits of the two methods. Section II outlines the results obtained from applying the ARIMA scheme to short-run forecasting of the international prices of ten selected commodities (foods and raw materials). Section III summarizes the findings.

I. Two Representations of an Economic System—Reduced Forms and Final Equations

In this section, a simple system of equations containing the most rudimentary elements of a full-fledged commodity model is introduced to illustrate the conceptual relationship between reduced forms and final equations. It is shown that the final equation for price has an ARIMA form. The discussion begins with a deterministic system, which is later made stochastic. Since the model is introduced for an illustration, no attempt is made to discuss appropriateness of the explanatory variables or the lag structure of the model. Suppose that a commodity market is represented by the following system of equations:

Qdt=α1Pt+α2YtQst=β1PtSt=γ1St1+γ2PtQstQdt=StSt1(1)

where

article image
article image

reduced forms and final equations

The four nonstochastic equations in (1) constitute a structural system. Each equation may include current and lagged endogenous and exogenous variables as regressors. From the structural system, a reduced form equation is derived for each of the endogenous variables. In a reduced form equation, only lagged endogenous variables and the current and lagged exogenous variables are used as regressors. The reduced form equation for price derived from the system in (1) is as follows:

Pt=1α1β1+γ2[(1γ1)St1α2Yt](2)

The reduced form equations are determined if the structural system is fully specified. If the coefficient values of the structural system are known, the values of the reduced form coefficients are also determined. Otherwise, the reduced form parameters can be either estimated directly or computed on the basis of estimated structural coefficients by deriving restricted reduced form equations. 3 In any case, forecasts can be obtained if the coefficient values, the future realizations of exogenous variables, and initial conditions are given.

The concept “final equation” was used in econometric analysis as early as 1939 by Tinbergen. For example, Tinbergen derives final equations in the form of autoregressive (AR) processes by eliminating endogenous variables from the right-hand sides of the equations in a structural system. 4 The final equation for price is derived from the structural system in (1) as follows: 5 Suppose that the income series Yt(t = 1, 2, …,) in the system (1) is viewed as generated by a linear stochastic process. Specifically, let the process be specified as the following autoregressive integrated moving average (ARIMA) process:

ΔYt=θ1ΔYt1+ϵt(3)

A combination of the four structural equations in (1) and the equation in (3) gives a final equation for each of the endogenous variables in the system. The final equation for price is as follows:

(1π1L)(1π2L)ΔPt=(1γ1L)ut(4)

Where

L = lag operator, such that LPt = Pt–1 and (1–L)Pt = ΔPt

π1=γ1(β1α1)γ2β1α1γ2π2=θ1ut=α2β1α1γ2ϵt

The final equation in (4) is an ARIMA (2, 1, 1) process. In a final equation, a current endogenous variable is expressed as a function of its own lagged values and the errors. The coefficients of the final equation in the form of (4) can be either estimated directly or identified from the estimates of structural parameters.

The final equation for price can still be derived in a similar way even if a random error term is added to the supply equation. Let the structural system be as follows:

Qdt=α1Pt+α2YtQst=β1Pt+ηtSt=γ1St1+γ2PtQstQdt=StSt1(5)

where ηt has a property such that

lt = ηt – ηt-1 = (1 – L)ηt

and ιt is an identically and independently distributed (iid) random error independent of ϵt. The final equation in this case is as follows:

[1γ1(β1α1)γ2β1α1γ2L](1θ1L)(1L)Pt=α2β1α1γ2(1γ1L)ϵt1β1α1γ2(1γ1L)(1θ1L)ιt(6)

or

(1π1L)(1π2L)ΔPt=(1γ1L)ut+(1γ1L)(1θ1L)wt(6)

where

wt=1β1α1γ2ιt

The final equation in (6’) has an autoregressive (AR) part identical with that in (4) and a moving average (MA) part that is the sum of two independent MA processes. The two MA processes with ut and wt are independent of each other because ϵt and ιt are mutually independent. The second term in

(1 – γ1L)ut + (1 – γ1L)(1 – θ1L)wt

comes from the disturbance term added to the supply function. As Anderson (1975) and Rose (1977) show, the sum of two or more independent and invertible MA processes can be represented by a single invertible MA process that has a different random error. The MA part in (6’) can therefore be written as

(1π1L)(1π2L)ΔPt=(1Φ1L)(1Φ2L)vt(7)

where the MA process (1 – Φ1L)(1 – Φ2L)vt is derived from the multiple MA process in (6′).

relative merits in forecasting

A successful prediction based on a model depends on a number of conditions: The model should be correctly specified; the coefficient values should be either known or consistently estimated with reasonably high asymptotic efficiency; and the system should be stable throughout the sample and prediction periods. Even if all these conditions are met, there can still be prediction errors if the equations have random error terms.

Although both can be derived from a given economic model, a reduced form and a final (ARIMA) equation have different capabilities in forecasting, since they are derived from different sets of assumptions and require different information bases for forecasting.

Reduced forms have an advantage in forecasting because a correct structural analysis can prevent the misspecification of a model. This is particularly true if the reduced form equations are derived from estimated structural equations. An economic and statistical evaluation can be made on whether specification and estimation are appropriate before forecasting is attempted. Furthermore, reduced forms facilitate an evaluation of alternative economic policies. This cannot be done with final equations, since final equations do not include exogenous variables as explanatory variables. Another merit of reduced forms in forecasting results from the fact that they permit an effective use of information on future realizations of the exogenous variables. As is explained later, the information is not so effectively utilized if final equations are used for forecasting.

To use reduced form equations for forecasting even one endogenous variable, however, generally the whole system must be estimated. Therefore, if only yearly data are available for some variables while quarterly series are available for the others, a quarterly model is difficult to use. Furthermore, separate forecasting should be made for the future values of each of the exogenous variables. Some scheme has to be found for predicting the future values of the exogenous variables in the system.

On the other hand, a foremost merit of a final (ARIMA) equation is simplicity. To forecast an endogenous variable, it is not necessary to deal with the entire model if a decision is made to “estimate” the final equation for the variable. Because no variables other than the endogenous variable that is explained enter the equation, the necessary predictions can be made on the basis of the one equation alone without separate predictions for the exogenous variables in the system. Because the whole system does not necessarily have to be estimated, final equations can be used for forecasting a number of endogenous variables even if observations on some variables in a model are not available.

It is obvious from the preceding discussion that final equations cannot be used for policy evaluations. Furthermore, the information on future exogenous variables cannot be effectively utilized if final equations are used for forecasting. Another weakness of final equations results from the insufficient state of present knowledge on their error terms: The error term in a final equation is a function of a number of the error terms in its associated structural system. In (4), for example, it is clear that the final equation error ut is nothing more than α2β1α1γ2ϵt, that is, the income equation error multiplied by a constant that is determined by structural coefficients. In many cases, however, this clear-cut relationship between structural and final equation errors is difficult to establish. In the second example in (7), for instance, the functional form that explicitly transforms the structural errors ϵt and ηt into the final equation error vt cannot be established. Unlike forecasting based on reduced form equations, it is therefore difficult to effectively utilize extraneous information on future production, consumption, and speculation if final equations are used for forecasting.

In spite of the limitations just discussed, final (ARIMA) equations are still useful especially for a short-run prediction of commodity prices. For example, a simultaneous equations model is difficult to estimate for a short-run prediction of monthly prices of many commodities because monthly data are generally hard to obtain for other relevant variables. Furthermore, forecasts obtained through final equations can be modified by incorporating judgments into the forecasts. Several aspects of this problem are discussed in the following subsection.

error terms in final equations

A useful illustration can be made of the relationship between the error terms in a structural equation and a final (ARIMA) equation. The illustration will further clarify some of the merits and limitations of final equations and will suggest possible improvements in their use.

Let the problem of estimation be disregarded. That is, let the coefficient values be known. This assumption will simplify the discussion. The reduced form equation for price is as follows:

Pt=1γ1α1β1+γ2St1α2α1β1+γ2Yt+1α1β1+γ2ηt(8)

On the other hand, the associated final equation with a multiple MA part has been shown to be

[1γ1(β1α1)γ2β1α1γ2L](1θ1L)(1L)Pt=α2β1α1γ2(1γ1L)ϵt1β1α1γ2(1γ1L)(1θ1L)ιt(6)

which can be written as

(1π1L)(1π2L)ΔPt=(1γ1L)ut+(1γ1L)(1θ1L)wt(6)

or

(1π1L)(1π2L)ΔPt=(1Φ1L)(1Φ2L)vt(7)

if the MA part in (6’) is made single. 6 The final equation is derived by substituting the third structural equation for St– 1 and the ARIMA process for Yt in the reduced form equation. If the realizations of future ηt and Yt are not known but reliable information is available for their prediction, using the reduced form equation in (2) or the final equation in the form of equation (6’) will give better results of prediction. On the other hand, the information cannot be effectively utilized for predictions if the final equation in the form of equation (7) is used, since the functional form that transforms ut and wt in (6′) into vt in (7) is not known.

As explained further in the following discussion, the information can be effectively utilized if reduced forms are used for forecasting. A justification for adding an error term to a structural equation can be found in those minor random incidents that affect the dependent variable. Without any information on the likelihood and magnitude of their occurrence, such incidents are expected to cancel each other out, resulting in a zero combined impact on the dependent variable. In other words, it is assumed that the unconditional mean of the error term vanishes. Once those incidents have occurred, however, it is possible to explain the realizations of the error term on the basis of the incidents that are observed in many cases. For example, it is possible to attribute a large positive residual of demand for fuel in a certain year to an unusually low temperature, or to explain a large positive residual of a hoarding equation for gold on the basis of an anticipation of a social or political turmoil. Similarly, it cannot but be maintained that the effects of those incidents on a dependent variable will cancel each other out during any prediction period if no information is available on the likelihood of their occurrence. In many cases, however, especially in short-run predictions, useful information can be obtained on possible magnitudes of those incidents that may occur during a prediction period. For example, a revised weather forecast may be obtained for a prediction period after the estimation of a model has been finished but immediately before a prediction is made. A rational forecaster in such a case will use the anticipated values of the error terms conditional on the new weather forecast. The information can be effectively utilized if either reduced forms derived from structural equations or “restricted” final equations are used for forecasting. Although the use of the information may not be equally effective, it can still improve forecasts obtained through (unrestricted) final equations if the information is utilized for a correct evaluation of those incidents’ net impact on the final equation error.

II. The Results of Estimation and Simulation for Ten Commodities

This section outlines the estimation and simulation results of the ARIMA equations for the international prices of ten commodities. Of these ten commodities, six are foods (bananas, beef, copra, maize, cocoa, and tea), and the others are raw materials (hides, jute, rubber, and copper). For each commodity, two ARIMA equations—one based on monthly price series and the other on annual data—are estimated. The price series used in this study are the spot prices in leading commodity markets.

estimation

In estimating the equations, a decision has to be made on proper prices to be used for the estimation. If a prediction of nominal prices is the main purpose of the estimation, using the nominal price has an obvious advantage of simplicity, because the deflator itself has to be predicted if relative prices are used for the estimation. During a period of inflation, nominal prices will show nonstationarity, whereas relative prices could be stationary. The nonstationarity problem can be dealt with, however, by successive differencing of the series in some cases.

In this study ARIMA equations are fitted to nominal prices. 7 Because most commodity prices rose sharply during the period of commodity price inflation (1973–75), the wisdom of using stationary models may be questioned. Solid evidence does not seem to exist, however, for maintaining that the commodity price inflation will be sustained in the future. Furthermore, graphical analyses of the monthly prices of the sample commodities do not reveal any strong indication of nonstationarity—at least until 1973, the initial year of the high inflationary period. The estimated autocorrelation functions for the ten monthly commodity prices also suggest that differencing operations may not be necessary. Stationary models are therefore used for the estimation, except for the annual models for the prices of bananas, copra, and maize. If the patterns of price movements are obviously different during a certain period, the months or the years (usually beginning or ending part of the period from January 1957 to September 1976) are excluded from the sample for estimation of the monthly models. 8 The sample period for the annual models is uniformly from 1965 to 1974.

Table 1 shows the orders of the lag polynomials used in the estimated equations, and Tables 2 and 3 present the results of estimation. A feature of the monthly models reported in Table 2 is that the models do not explicitly describe the cyclical price movements for durations longer than one year. Some of the models are designed to explain seasonal price fluctuations (3-month, 6-month, or 12-month). The models, although including prices lagged by, at most, 24 months as explanatory variables, could explain longer-run swings of the prices. The equations are selected, however, not for their capability to predict longer-run fluctuations but for their power to predict month-to-month changes in prices. The sample autocorrelation functions of the monthly prices of beef, cocoa, hides, rubber, and copper suggest longer cycles of varying duration. 9 The existence of such cycles is not explicitly reflected in the chosen functional forms, however, either because the evidence is not believed to be sufficiently strong or because the sample period is not sufficiently long. 10

Table 1.

Ten Commodities: Selected Functional Forms1

article image

figures indicate the orders of the lag polynomials included in the equations.

Annual models are based on the first differences of prices.

Table 2.

Ten Commodities: Estimated ARIMA Equations for Monthly Prices1

article image
article image

Multiplicative seasonal models are used for some commodities. The lag operator is indicated by L. The figures accompanying the estimated coefficients are the ratios between the estimated coefficients and the associated standard errors. The ratios would behave asymptotically like normal random variables. (See footnote 8 in the text for the sample periods.) The coefficient of determination (R2) is computed by using the formula 1V^(at)/V^(Pt).

Table 3.

Ten Commodities: Estimated ARIMA Equations for Annual Prices1

article image

No satisfactory equations are obtained for bananas, copra, and maize. The sample period is uniformly from 1965 to 1974. The lag operator is denoted by L. The figures under the coefficients are the ratios between the estimated coefficients and the associated standard errors.

As for seasonal fluctuations, the autocorrelation functions of the levels of the monthly prices reveal that there are conspicuous 12-month and 6-month cycles in the tea price and a strong 12-month cycle in the banana price (Chart 1), while the autocorrelation functions of the first differences of the prices indicate possible 6-month cycles in the prices of beef, cocoa, and hides and 3-month cycles in the prices of bananas and tea. Multiplicative seasonal models are therefore used for bananas, beef, cocoa, tea, and hides. The results show that most of the parameters of the seasonal lag polynomials are statistically significant. For example, the estimated equation for the monthly price of bananas includes two significant seasonal autoregressive coefficients (for Pt–3 and Pt–12), while the equation for the monthly price of tea has two highly significant seasonal (12-month) coefficients (one autoregressive and the other moving average). The equation for the monthly price of hides also includes highly significant seasonal (6-month) autoregressive and moving average parameters.

Chart 1.
Chart 1.

Autocorrelation Functions of Monthly Prices of Bananas and Tea

Citation: IMF Staff Papers 1978, 001; 10.5089/9781451956498.024.A005

The moduli of all the roots of the estimated lag polynomials lie outside the unit circle, indicating that all the prices are stationary during the sample periods and that the moving average parts are invertible. 11 For the monthly models, the statistic R2, defined as

1V^(at)V^(Pt)

where V^(at) and V^(Pt) denote estimated variances of at and Pt, ranges from 0.644 (bananas) to 0.970 (beef). 12 At least 88 per cent of the variations in the monthly price is explained by a second-order AR lag polynomial for copra, maize, jute, and rubber, while a fourth-order AR lag polynomial is required to explain 90 per cent of the variations in the price of copper. The monthly prices of two commodities (tea and hides) are better explained by equations with MA parts. Both of the equations have MA parts with 6-month seasonal lag polynomials, and the estimated parameters of the polynomials are highly significant.

The estimated monthly models may be used for a short-run prediction of the commodity prices: they can be used for a prediction of the monthly prices; they can also be used for a prediction of the quarterly, half-yearly, or yearly averages of the monthly prices. If a prediction of the yearly averages of the monthly prices is all that is needed, however, annual—rather than monthly—models may be used. For a comparison of the predictive powers of monthly and annual models, simple autoregressive or moving average models are estimated on the basis of the yearly price series of the ten commodities. Since the yearly price series are not sufficiently long, the identification process utilized for the monthly models is bypassed. First-order and second-order autoregressive and moving average models are tried, and a model that is considered to be most reasonable is selected for each commodity. The results are presented in Table 3. 13 It is interesting to note that no satisfactory annual models are found for three commodities (bananas, maize, and copra). Therefore, the annual prices of these commodities are simply assumed to be generated by random walk processes. The coefficients of determination for these prices reported in Table 3 are the coefficients of determination between the current prices and the prices lagged by one year. The coefficients are 0.704 and 0.853, respectively, for copra and maize. For bananas, the coefficient is zero; this implies that using the yearly price series for forecasting the price of bananas is not feasible for the sample period. The yearly prices of the remaining seven commodities yield models with coefficients of determination ranging from 0.159 (tea) to 0.735 (cocoa). The coefficients of determination associated with the monthly models are generally higher than those associated with the annual models, since the former are based on monthly predictions while the latter are based on yearly predictions. The difference between the coefficients of determination associated with the monthly model and the annual model is more pronounced in the seasonal commodities than in the others.

simulation

The procedure described in Box and Jenkins (1970) is used for two sets of dynamic simulations for the monthly models—a series of 12-month and a series of 6-month simulations. 14 For example, a simulated path of the monthly maize price is computed as a result of a series of 12-month dynamic simulations beginning in January each year. The simulations are dynamic in the sense that within any one year (or within any 6-month period for 6-month dynamic simulations) some or all of Pt–j+ls (j> 0), the lagged prices that are used as explanatory variables, are replaced by their simulated values. These dynamic simulations are carried out 12 months (or 6 months) at a time. For maize, the equation that is used for simulation is

P^t+e=0.078+1.338P^t1+e0.395P^t2+e(9)

with = 1, 2, …, 12, for 12-month simulations. By nature of the forecasting scheme and because the models are stationary, the simulated prices tend to converge to the series’ mean values. For example, the time path of the simulated maize price will converge to its mean value if goes to infinity

limP^t+=0.07811.338+0.395(10)

Therefore, if the price of a commodity is higher (lower) than the mean during the months immediately preceding a simulation period, the simulated path of the price tends to fall (rise) ultimately.

An evaluation of the scheme’s predictive power of monthly prices is made by comparing the actual path with the simulated path of the price. First, the correlation coefficient (R6) is calculated between the actual monthly price and the simulated monthly price resulting from 6-month dynamic simulations for each commodity for the entire simulation period (1965–74). Second, the correlation coefficient (R12) is calculated between the actual monthly price and the simulated monthly price resulting from 12-month dynamic simulations. The correlation coefficient R12 indicates the average degree of accuracy of the monthly model’s short-run (one-year) predictions of monthly prices, whereas the correlation coefficient R6 indicates that of 6-month forecasts of monthly prices. In other words, R12 indicates the degree of accuracy of predictions made at the end of each year for the following 12 monthly prices, while R6 indicates that of predictions made at the end of June and December of each year for the following 6 months. If predictions are made at the end of each year for the following 12 months and if the forecasts for the second half of each year are revised on the basis of the predictions made at the end of June, the degree of accuracy of the revised monthly forecasts will be indicated by R6, which should generally be higher than R12. Also, the simulation period is a part of the sample period for each commodity. For an evaluation of the scheme’s predictive power of the yearly averages of monthly prices, the simulated monthly prices obtained from the 12-month simulations and the 6-month simulations are converted into yearly averages. The correlation coefficients (R6a and R12a) are computed from the actual yearly price series and the yearly averages of the simulated monthly prices based on the 12-month and 6-month dynamic simulations.

Table 4 shows the four correlation coefficients (R6, R12, R6a, and R12a) and the correlation coefficient (Ra) between the actual yearly price and the price statically simulated by the annual model for each commodity. A few relevant facts should be mentioned: First, the figures in the table are correlation coefficients—not coefficients of determination; second, they are based on short-run (at most 12-month) simulations. It is not surprising therefore that the coefficients are relatively high in many cases and as high as 0.985 (R12a, maize) in one case. Longer-run simulations would inevitably produce lower coefficients. Also, the simulations conducted in this study are based on no information other than the price of the commodity in consideration. The results suggest that, even if a prediction of the yearly average of the monthly prices of a commodity is the only purpose of the model, working with monthly—rather than yearly—data can give better results. For example, the monthly models for bananas, copra, and maize yield, respectively, 0.536, 0.903, and 0.985 as coefficients, although no satisfactory annual models can be found for the prices. 15 The correlation coefficients improve when the monthly forecasts are converted into yearly forecasts, probably because the monthly prediction errors are canceled when yearly averages are derived.

Table 4.

Ten Commodities: Goodness of Fit of Simulated Prices

article image

Simulated paths of two monthly prices (for bananas and tea) are compared with the actual monthly prices in Charts 2 and 3 for an illustration. The two price series show pronounced seasonality, which seems to become less regular for tea during the second half of the simulation period. Also, the monthly model for tea persistently underpredicts the price during the second half of the simulation period. This may be due to the changes in the seasonal pattern of the price that began to emerge during the period.

Chart 2.
Chart 2.

Monthly Price of Bananas: Series of 12-Month Dynamic Simulations, Beginning in January, 1965–74

Citation: IMF Staff Papers 1978, 001; 10.5089/9781451956498.024.A005

Chart 3.
Chart 3.

Monthly Price of Tea: Series of 12-Month Dynamic Simulations, Beginning in January, 1965–74

Citation: IMF Staff Papers 1978, 001; 10.5089/9781451956498.024.A005

possible improvements

The limitations of the ARIMA model in forecasting commodity prices and other economic time series have already been discussed. The most important of these is probably the fact that the model in the form of (7) cannot effectively utilize the information that a forecaster may have on the exogenous variables. Even if accurate forecasts of Yt are available, the forecaster who uses the equation in (7) would have to make predictions as if he did not have the information. This is an important waste of information. This is also true for the information that the forecaster may have on the future realizations of the structural errors.

Although no attempt to do so is made in this paper, judgmental forecasts based on economic analysis of various factors in commodity markets can be used to revise forecasts made through ARIMA models. An autoregressive final form (or a transfer function) may also be used for prediction if one has accurate information on future realizations of the exogenous variables. For example, one may derive the autoregressive final form for price from the structural system in (1) as

(1ψ1L)ΔPt=ψ2(1γ1L)ΔYt+(1γ1L)nt(11)

where

ψ1=γ1(α1β1+γ2)+(1γ1)α1β1+γ2ψ2=α2α1β1+γ2nt=1α1β1+γ2ι2

by combining the stock equation in (5) with the reduced form in (8). The equation can be estimated and used for prediction. An estimation of the whole system is not necessary for an effective utilization of an equation like (11) for forecasting.

III. Summary and Conclusions

The relative merits of final (ARIMA) equations in forecasting economic time series are reviewed in this paper. The foremost advantage is simplicity. No series other than the one that is predicted is required in forecasting based on final equations. This simplicity facilitates effective use of existing data, which may have to be underutilized if forecasting is based on reduced forms. Thus, a monthly series can be fully utilized through a final equation, while it may have to be converted into quarterly, half-yearly, or yearly series if reduced forms are used for forecasting and if other necessary series are not collected on a monthly basis. On the other hand, reduced forms are essential in policy analysis and more effective in utilizing the relevant information on the future realizations of exogenous variables and structural errors.

The relationship between the reduced form and the final equation suggests that there is some economic justification for using ARIMA models for forecasting: ARIMA models are not mere statistical constructs. An ARIMA equation compatible with a correctly specified structural model could be estimated if the estimation of the equation is constrained by the structural relationships among relevant variables. Even if the ARIMA equation is estimated without any constraints, the estimated equation should be expected to be consistent with the correct structural system if identification and estimation are appropriate.

The results of an application of the scheme to short-run simulations of the international prices of ten commodities are also described in the paper. Annual and monthly models were estimated and used for a simulation of the yearly averages of the monthly prices for the period 1965–74. The annual model gave slightly better results of simulation for beef, tea, and jute, while the monthly model performed better for the other commodities. For bananas, copra, and maize, no satisfactory annual model could be found for the sample period, while the monthly models gave quite accurate short-run predictions for copra and maize.

The time-series technique discussed in this study is a simple and economical device that can be used for short-run forecasting of commodity prices. The scheme has, however, only limited capability to predict unusual movements in prices. For example, the scheme can neither predict the sudden increase in the price of tea nor forecast the wide fluctuations in the price of bananas during the period 1972–74. The predictive power can possibly be improved by incorporating into the scheme judgments based on an economic analysis of the market or by introducing relevant exogenous variables into the equation.

BIBLIOGRAPHY

  • Anderson, O. D., “On a Lemma Associated with Box, Jenkins and Granger,” Journal of Econometrics, Vol. 3 (May 1975), pp. 15156.

  • Box, George E. P., and Gwilym M. Jenkins, Time Series Analysis: Forecasting and Control (San Francisco, 1970).

  • Brillembourg, Arturo, “Purchasing Power Parity Tests of Causality and Equilibrium” (unpublished, International Monetary Fund, December 27, 1976).

    • Search Google Scholar
    • Export Citation
  • den Butter, F. A. G., “The Use of Monthly and Quarterly Data in an ARMA Model,” Journal of Econometrics, Vol. 4 (November 1976), pp. 31124.

    • Search Google Scholar
    • Export Citation
  • Dhrymes, Phoebus J., “Restricted and Unrestricted Reduced Forms: Asymptotic Distribution and Relative Efficiency,” Econometrica, Vol. 41 (January 1973), pp. 11934.

    • Search Google Scholar
    • Export Citation
  • Feige, Edgar L., and Douglas K. Pearce, “Economically Rational Expectations: Are Innovations in the Rate of Inflation Independent of Innovations in Measures of Monetary and Fiscal Policy?” Journal of Political Economy, Vol. 84 (June 1976), pp. 499522.

    • Search Google Scholar
    • Export Citation
  • Frenkel, Jacob A., “The Forward Exchange Rate, Expectations, and the Demand for Money: The German Hyperinflation,” American Economic Review, Vol. 67 (September 1977), pp. 65370.

    • Search Google Scholar
    • Export Citation
  • Naylor, Thomas H., Terry G. Seaks, and D. W. Wichern, “Box-Jenkins Methods: An Alternative to Econometric Models,” International Statistical Review, Vol. 40 (August 1972), pp. 12337.

    • Search Google Scholar
    • Export Citation
  • Nelson, Charles R. (1972), “The Prediction Performance of the FRB-MIT-PENN Model of the U. S. Economy,” American Economic Review, Vol. 62 (December 1972), pp. 90217.

    • Search Google Scholar
    • Export Citation
  • Nelson, Charles R. (1976), “The Interpretation of R2 in Autoregressive-Moving Average Time Series Models,” American Statistician, Vol. 30 (November 1976), pp. 17580.

    • Search Google Scholar
    • Export Citation
  • Palm, Franz, “On Univariate Time Series Methods and Simultaneous Equation Econometric Models,” Journal of Econometrics, Vol. 5 (May 1977), pp. 37988.

    • Search Google Scholar
    • Export Citation
  • Pindyck, Robert S., and Daniel L. Rubinfeld, Econometric Models and Economic Forecasts (New York, 1976).

  • Rose, David E., “Forecasting Aggregates of Independent ARIMA Processes,” Journal of Econometrics, Vol. 5 (May 1977), pp. 32345.

  • Saboia, João Luiz Maurity, “Autoregressive Integrated Moving Average (ARIMA) Models for Birth Forecasting,” Journal of the American Statistical Association, Vol. 72 (June 1977), pp. 26470.

    • Search Google Scholar
    • Export Citation
  • Thompson, Howard E., and George C. Tiao, “Analysis of Telephone Data: A Case Study of Forecasting Seasonal Time Series,” Bell Journal of Economics and Management Science, Vol. 2 (Autumn 1971), pp. 51541.

    • Search Google Scholar
    • Export Citation
  • Tinbergen, Jan, Business Cycles in the United States of America, 1919–1932 (Geneva, 1939).

  • Zellner, Arnold, and Franz Palm, “Time Series Analysis and Simultaneous Equation Econometric Models,” Journal of Econometrics, Vol. 2 (May 1974), pp. 1754.

    • Search Google Scholar
    • Export Citation
*

Mr. Chu, economist in the Commodities Division of the Research Department, is a graduate of Kyung Hee University (Seoul) and Columbia University. Before joining the Fund, he was an instructor at Columbia University.

1

There have been studies on applications of time-series models for forecasting macroeconomic variables (Naylor, Seaks, and Wichern, 1972, and Nelson, 1972), for predicting birth rates (Saboia, 1977), and for analyzing telephone data (Thompson and Tiao, 1971); time-series models have also been used for specifying an econometric model (Zellner and Palm, 1974) and for evaluating economic causality (Brillembourg, 1976, Feige and Pearce, 1976, and Frenkel, 1977).

2

All the equations identified and estimated in this study exhibit stationarity: all of them are in autoregressive integrated moving average (ARIMA; p, d, q) forms with d = 0, or autoregressive moving average (ARMA; p, q1) forms, where p, d, and q indicate the orders of the AR part, differencing, and the MA part, respectively.

4

The nonstochastic second-order AR equations in Tinbergen (1939, p. 131) are essentially ARIMA (2, 0, 0) or ARMA (2, 0) processes.

5

For derivation of final equations from a structural system with a more general moving average part, see Zellner and Palm (1974) or Palm (1977).

6

The equation in (6) may be called a restricted final equation if the parameter estimates are derived from estimated structural coefficients.

7

The series used in preparing this paper are similar to but not identical with those published in International Monetary Fund, International Financial Statistics.

8

The sample periods for the estimation of the monthly models are not uniform. Bananas: January 1957–September 1976; beef: January 1957–November 1976; copra: January 1957–December 1972; maize: January 1957–December 1972; cocoa: January 1957–December 1973; tea: January 1957–September 1976; hides: January 1957–September 1976; jute: January 1957–September 1976; rubber: January 1957–September 1976; copper: January 1957–November 1976.

9

Beef: 2.5 years; cocoa: 3 years; hides: 3 years; rubber: 2.5 years; copper: longer than 3.5 years.

10

The sample autocorrelation function (ρ^k) attains a local maximum at k = 24 (months) for the price of hides and at k = 30 (months) for the price of rubber. The maximum values of the correlation functions are 0.38 and 0.26, respectively. For copper, where the cycle is believed to be longer, the sample period is not long enough to make a clear indication of the exact duration of any cycle.

11

The third-order regular autoregressive lag polynomials for bananas and beef have two imaginary roots and a real root each, and the autoregressive lag polynomials for bananas, tea (both second-order), and copper (fourth-order) have all imaginary roots, whereas all the other roots are real. The moduli of the imaginary roots range from 1.705 (tea, regular) to 2.975 (bananas, 3-month seasonal), and those of the real roots of the second-order, or higher-order, lag polynomials range from 1.006 (hides, regular) to 4.306 (bananas, 12-month seasonal). The real roots of the second-order or higher-order lag polynomials for bananas and hides are negative.

12

See Nelson (1976) for an interpretation of the statistic. Since the statistic is a measure of the degree of static simulation and is comparable with R2 in standard multiple regressions, the results of dynamic simulation may be different.

13

No attempt was made to find an annual model consistent with the estimated monthly model for any commodity. The annual models therefore may not be consistent with the monthly models in the way in which equation (6) is consistent with equation (3) in den Butter (1976).

14

The simulations are carried out for the ten-year period 1965–74, during which most prices showed both regular and unusual movements. This period is a subperiod of the sample period used for estimating the equations for all commodity prices.

15

The random walk model gives 0.839 and 0.923 for copra and maize, but the coefficients are lower than those associated with the monthly models. (See the Ra column and the Rl2a column of Table 4.) The difference between Rl2a and Ra (0.536–0.000 for bananas, 0.903–0.839 for copra, and 0.985–0.923 for maize) indicates the increase in the goodness of fit made possible by the use of the monthly model. Random walk models based on yearly price series would generate the simulated prices no better than those generated by the estimated annual AR or MA models in Table 3, since the estimated models maximize the coefficients of determination between the actual and the predicted prices.

IMF Staff papers: Volume 25 No. 1
Author: International Monetary Fund. Research Dept.