Time-Series Estimation of Structural Import Demand Equations: A Cross-Country Analysis

Contributor Notes

Author’s E-Mail Address: asenhadji@imf.org

This paper derives a structural import demand equation and estimates it for a large number of countries, using recent time series techniques that address the problem of nonstationarity. Because the statistical properties of the different estimators have been derived only asymptotically, econometric theory does not offer any guidance when it comes to comparing different estimators in small samples. Consequently, the paper derives the small-sample properties of both the ordinary-least-squares (OLS) and the fully-modified (FM) estimators using Monte Carlo methods. It is shown that FM dominates OLS for both the short- and long-run elasticities.

Abstract

This paper derives a structural import demand equation and estimates it for a large number of countries, using recent time series techniques that address the problem of nonstationarity. Because the statistical properties of the different estimators have been derived only asymptotically, econometric theory does not offer any guidance when it comes to comparing different estimators in small samples. Consequently, the paper derives the small-sample properties of both the ordinary-least-squares (OLS) and the fully-modified (FM) estimators using Monte Carlo methods. It is shown that FM dominates OLS for both the short- and long-run elasticities.

I. INTRODUCTION

The empirical investigation of the import demand function has been one of the most active research areas in international economics. This is evidenced by the many surveys on this topic, although most focus on industrial countries.2 Perhaps one of the main reasons for its popularity is its application to a wide range of important macroeconomic policy issues such as the international transmission of domestic disturbances, where these elasticities are a crucial link between economies; the impact of expenditure-switching through exchange rate management and commercial policy on a country’s trade balance; and the degree to which the external balance affects a country’s growth.

The traditional import demand function is specified as a log-linear function of the relative price of imports and real income. Because of data constraints and the empirical success of this specification, it has dominated the empirical literature for more than a quarter century. But, questions about its microeconomic foundation arise since it has not been derived from utility maximization. Another issue that has been largely ignored in the literature is the problem of nonstationarity, which is found present in most macroeconomic variables and which invalidates classical statistical inference. Thus, if the variables that enter the import demand equation contain a unit root, ignoring nonstationarity in these variables may cause serious inference problems.

The objective of this paper is twofold. First, the paper seeks to address the two problems discussed above by deriving an empirically tractable import demand equation that can be estimated for a large number of countries, using recent time series techniques that address the issue of nonstationarity present in the data. Second, because the statistical properties of the different estimators have been derived only asymptotically, econometric theory does not offer any guidance when it comes to comparing the performance of different estimators in small samples. Consequently, the paper derives the small sample properties of both the Ordinary Least Squares (OLS) and the Fully-Modified (FM) estimators of the short- and long-run elasticities, using Monte Carlo methods. It is shown that the FM estimators dominates the OLS estimators, even in small samples.

The derived aggregate import demand equation is log-linear in the relative price of imports and an activity variable defined as GDP minus exports.3 An important insight from the explicit derivation of the aggregate import demand equation is that the definition of the activity variable depends on the aggregation level.4 The model predicts a unique cointegrating vector among imports, the relative price of imports and the activity variable. This prediction is not rejected by the data, and the cointegrating vector is estimated efficiently by the Phillips-Hansen FM estimator. Two recent papers follow a similar methodology. Clarida (1994) derives a similar import demand function for U.S. nondurable consumption goods from explicit intertemporal optimization, carefully taking into account data nonstationarity. Similarly, Reinhart (1995) estimates both structural import and export demand functions for twelve developing countries using Johansen’s cointegration approach.

The results underscore the presence of nonstationarity in the data and the adverse consequences of neglecting it. Both price and income elasticities generally have the expected sign and are precisely estimated. The average price elasticity is close to zero in the short run but is slightly higher than one in the long run. It takes five years for the average price elasticity to achieve 90 percent of its long-run level. A similar pattern holds for income elasticities in the sense that imports react relatively slowly to changes in domestic income. The short-run income elasticities are on average less then 0.5, while the long-run income elasticities are close to 1.5. Industrial countries have both higher income and lower price elasticities than do developing countries. On average, these estimates are relatively close to Reinhart’s.5

Empirical researchers are generally interested in two statistical properties of their estimates of import elasticities. First, they are interested in the magnitude of these elasticities. A relevant question, then, is how close the estimates are to their true value. The systematic deviation of estimates from their true value is measured by the bias of the estimates. Second, they are interested in inference, that is, hypotheses testing, about these estimates. For example, are the price and income elasticities significantly different from one? Testing such hypotheses requires knowing the distribution of the t-statistic (defined as the coefficient estimate divided by its standard deviation). The asymptotic distribution of this statistic is unknown for the long-run elasticities because these elasticities are nonlinear transformations of the import demand coefficients. In addition, the definition of the long-run elasticities includes the lagged dependent variable whose t-statistic follows a nonstandard distribution in the nonstationary case. In light of this, using the critical values of the t-distribution for hypothesis testing may be misleading. Consequently, the small sample distribution of the t-statistic for both the short- and long-run elasticities are computed using Monte Carlo methods.

The analysis shows that the OLS bias is significantly higher than the FM bias for both the short- and long-run elasticity estimates. The FM bias reaches its minimum when the relative price of imports and the activity variable are exogenous. Strong endogeneity of the explanatory variables (that is, high correlation between the import demand innovations and the explanatory variables innovations) may induce substantial bias. But for most countries—being “small” relative to the rest of the world—the relative price of imports and the activity variable are only weakly endogenous, leading to a relatively small bias. The bias of long-run elasticities is generally much lower than the bias of short-run elasticities. For the benchmark case in which both explanatory variables are assumed to be exogenous, the t-statistics of the short-run elasticities are symmetric around zero but are flatter than the asymptotic t-distribution. This implies that an inference based on the usual t-or F-statistic may be misleading. For example, the exact confidence intervals are wider than those based on the t-statistic. The t-statistic distribution of the short-run elasticities become skewed and flatter when the relative price of imports and/or the activity variable is allowed to be endogenous. The stronger the endogeneity, the larger is this departure from the asymptotic t-distribution.

II. THE MODEL

Assume that the import decision in each country is made by an infinitely lived representative agent who decides how much to consume from the domestic endowment (dt) and from the imported good (mt).6 The home good is the numeraire. The intertemporal decision can be formalized by the following problem:

  Max  E0t=0(1+δ)1u(dt,mt){dt,mt}t=0(1)

subject to:

bt+1=(1+r)bt+(etdt)ptmt(2)
et=(1ρ)e¯+ρet1+ξt,ξt~(0,σ2)(3)
limTbT+1t=0T(1+r)1=0,(4)

where δ is the consumer’s subjective discount rate; r is the world interest rate; bt+1 is the next period stock of foreign bonds if positive, and the next period’s debt level if negative; et is the stochastic endowment which follows an AR(1) process with unconditional mean e¯ and an unconditional variance σ2/(l-ρ2), where σ2 is the variance of the iid innovation ξt and ρ determines the degree of persistence of the endowment shocks; and pt, is the relative price of the foreign good, that is, the inverse of the usual definition of terms of trade. In this two-good economy, pt also represents the relative price of imports.Equation (2) is the current account equation, equation (3) is the stochastic process driving the endowment shock, and equation (4) is the transversality condition that rules out Ponzi games. The first order conditions of this problem are:

utd=λt(5)
utm=λtpt(6)
λt=(1+δ)1(1+r)Etλt+1,(7)

where λt is the Lagrange multiplier on the current account equation. From equation (5), λt is the marginal utility of the domestic good. Following Clarida (1994, 1996) and Ogaki (1992), it is assumed that the instantaneous utility function u is addilog:7

u(dt,mt)= Atdt1α(1α)1+ Btmt1β(lβ)1α>0, β>0(8)
At=ea0+A,t(9)
Bt=eb0+B,t(10)

where At and Bt are exponential stationary random shocks to preferences, ∈A, t and ∈B, t are stationary shocks and α and β are curvature parameters. Substituting equation (8) into equations (5) and (6) yields:

dt=λt1/αAt1α(11)
mt=λt1βBt1βpt1β.(12)

Substituting equations (9)-(11) into equation (12) and taking logs yields:

m˜t=c1βp˜t+αβd˜t+t,(13)

where c0 = (1/β)(b0a0), and ∈t = (1/β)(∈B, t − ∈A, t). A tilde indicates the log of the corresponding variable. In this model, xt = etdt = GDPtdt, where xt is exports. Consequently, dt = ptxt. Thus, the model yields an equation for import demand that is close to the standard import demand function except that the correct activity variable is GDPtxt rather than GDPt.Equation (13) can be rewritten as:

m˜t=c1βp˜t+αβ(GDPtxt~)+t.(14)

Taking logs of equation (11) yields:

d˜t=1αA˜t1αλ˜t.(15)

Because each of the three variables in the import demand equation (13) can either be trend-stationary (TS) or difference-stationary (DS), four cases need to be considered (Table 1). In the next section, results from unit root tests show that the first case is the most common, with some countries falling into the second category. The prime interest is the estimates of the standard price and income elasticities for import demand, defined (respectively) as the coefficients of the log of the relative price of imports (−1/β) and the log of the activity variable (α/β). Note that m˜t and p˜t are, in general, endogenously determined by import demand and import supply (not modeled here). Therefore, p˜t is likely to be correlated with the error term ∈t in equation (14). Thus OLS would yield biased estimates of the price and income elasticities. The Phillips-Hansen FM estimator corrects for this potential simultaneity bias, as well as for autocorrelation, in the cointegration framework.

Table 1:

The four possible model specifications

article image
article image

p˜t is either exogenous under a perfectly elastic import supply or is endogenously determined by the interaction of the demand and supply of imports. The supply of imports is not modeled explicitly here.

Equation (14) will be estimated in a dynamic form (that is, with the lagged dependent variable included as an explanatory variable) which proved to be more successful in the estimation stage.8 It is obtained by postulating a partial adjustment process of actual imports toward import demand:9

Δm˜ta=ϕ[m˜tm˜t1a],  |ϕ|<1,(16)

where m˜ta and p˜t and m˜ta denote actual and demanded imports, respectively. If ϕ is close to one, it implies that actual imports adjust quickly to import demand. Substituting equation (14) into (16) yields the final import demand equation:

m˜ta=θ0+θ1m˜t1a+θ2p˜t+θ3(gdptx~t)+t,(17)

Where θ0=ϕ c, θ1=1−ϕ, θ2=−ϕ(1/β), and θ2=ϕ(α/β). Note that all the coefficients of the import demand equation can be recovered from equation (17). The discussion above and Table 1 remain valid for equation (17) as long as |ϕ|<1,, which is the case for most countries in our sample (see Table 2).10

Table 2:

Augmented-Dickey-Fuller test for variables entering the import demand equation

article image
article image
Note: Variables are real imports of goods and non-factor services (m), the real exchange rate (p), computed as the ratio of imports deflator to GDP de exports (gdpx). These three variables are tested for the existence of a unit root using the Augmented-Dickey-Fuller (ADF) test. The optimal lag se criterion in the ADF regression is given by k. Critical values are a linear interpolation between the critical values for T=25 and T=50 given in Hamilton 4), where T is the sample size. Significance levels at 1% and 5% are indicated by**and*, respectively.

III. ESTIMATION RESULTS

The import demand equation (17) will be estimated by both the OLS and the FM estimator. The FM estimator is an optimal single-equation method based on the use of OLS with semiparametric corrections for serial correlation and potential endogeneity of the right-hand variables. The method was developed in Phillips and Hansen (1990) and generalized to include deterministic trends by Hansen (1992a). The FM estimator has the same asymptotic behavior as the full systems maximum likelihood estimators.11

The data comes from the World Bank database BESD. The sample includes 77 countries for which the required data are available for a reasonable time span (the list of countries is given in Table 2). In general, the data are available from 1960 to 1993.12 The usual problem is of course the choice of the corresponding proxies for the variables in the model, since the model is usually a crude simplification of reality—which is the case here. Data constraints highly restrict this choice. Total imports and exports of goods and services will be used for mt and xt in equation (17). The relative price of imports pt will be computed as the ratio of the import deflator to the GDP deflator.13 The activity variable will be computed as the difference between GDP and exports.14

A. Unit root test

The Fully-Modified procedure assumes that some of the variables entering the cointegrating equation (17) have a unit root and that there exists a stationary linear combination of these variables. This section tests for the existence of a unit root in all three variables in the import demand equation (17), namely real imports of goods and services (m), the relative price of imports (p) and the activity variable GDP minus exports (gdpx). The unit root hypothesis is tested using the augmented-Dickey-Fuller (ADF) test. The lag length (k) in the ADF regression is selected using the Schwarz criterion (Table 2). For mt, only 4 out of the 77 countries reject the unit root at 5 percent or less (Australia at 1 percent, Nicaragua, Peru and Philippines at 5 percent). Similarly, the null of a unit root in pt is rejected only for 3 countries (China at 1 percent, Papua New Guinea and Uruguay at 5 percent). Finally, as far as gdpxt is concerned, the unit root is rejected for 10 countries (Burundi, Central African Republic, Iceland, Switzerland, Trinidad and Tobago at 1 percent; Korea, Rwanda, Togo, Tunisia and Zaire at 5 percent). Thus for most of the countries, the unit root hypothesis cannot be rejected at conventional significance levels. This finding, of course, may reflect to a certain extent the low power of the ADF.

B. Import demand equations

The results underscore the presence of nonstationarity in the data and the adverse consequences of neglecting it. Table 2 shows that most countries—60 of the 77—fall into the first case of Table 1 (the unit root hypothesis cannot be rejected for all three variables in the import demand equation) and the remaining countries—17 of the 77—into the second case (the unit root hypothesis can be rejected for only one of the three variables). In the first case, the model predicts a cointegrating relationship between the three I(1) variables, and in the second case between the two I(1) variables. No country belongs to either the third or fourth cases.

Table 3 shows both the OLS and FM estimates of the import demand equation. Only countries with the right sign for the price and income elasticities are reported (66 of the 77 countries). The columns of Table 3 labeled x−1,p and gdpx give, respectively, the coefficient estimates of the lagged dependent variable (log of imports of goods and nonfactor services), the short-term price elasticity (that is, the coefficient of the log of relative price of imports) and the short-term income elasticity (the coefficient of the log of gdpx). The long-run price and income elasticities are defined as the short-run price and income elasticities divided by one minus the coefficient estimate of the lagged dependent variable. They are given by Ep and Ey for the FM estimates.15 The column labeled ser reports the standard error of the regression. Finally, the column labeled AC gives Durbin’s autocorrelation test. It amounts to estimating an AR(1) process on the estimated residuals of the import equation. Durbin’s test is simply a significance test of the AR(1) coefficient using the usual t-test. For the OLS regressions, AR(1) autocorrelation is detected (at 10 percent or less) for 17 of the 66 countries.

Table 3:

Import demand equations

article image
article image
article image
article image
Note: The dependent variable is real imports of goods and non-factor services (m). The explanatory variables are the lagged dependent variable (m.j), the real exchange rate (p) computed as the ratio of imports deflator to gdp deflator, and gdp minus exports (gdpx). The import demand equation is estimated using both OLS and the Phillips-Hansen Fully Modified estimator. The long-run price and income elasticities are given by Ep and Ey, respectively. Epc and Eyc give the long-run price and income elasticities corrected for bias (see Table 4). For each country, the estimated coefficients and their t-stat (below the coefficient estimates) are provided. The following statistics are also provided: Durbin’s test for autocorrelation (AC), R2, standard error of the regression (ser), and the number of observations for each country (nobs). Cointegration between the three variables in the import demand equation is tested using the Phillips-Ouliaris residual test given in column P-O. Finally, the columns labeled Ep=−1 and Ey=l report the two-tailed t-test for unit-price and unit-income elasticities, respectively. The asymptotic critical values for the -Phillips-Ouliaris test at 10%, 5% and 1% are, respectively, −3.84, −4.16 and −4.64. The letters a, b and c indicate significance at 1%, 5% and 10%. Exact critical values (from Table 8) are used to compute the significance level of Ep, Ey, Ep=−1 and Ey=1