fedea
Fundación de Estudios de Economía Aplicada energy and climate change
http://www.energiaycambioclimatico.com
Gasoline and Diesel Consumption for Road Transport in Spain: a Dynamic Panel Data Approach
Rosa M. González-Marrero Departamento de Análisis Económico, IUDR (ULL)
Rosa M. Lorenzo-Alegría Departamento de Análisis Económico, IUDR (ULL)
Gustavo A. Marrero
Departamento de Análisis Económico (ULL) Programa de Investigación Energía y Cambio Climático Focus-Abengoa-Fedea
Colección Estudios Económicos 04-2011 Programa de Investigación Energía y Cambio Climático Fedea – Focus Abengoa
ISSN 1988-785X
www.fedea.es
Rosa M. González-Marrero1 Departamento de Análisis Económico, IUDR (ULL) Rosa M. Lorenzo-Alegría Departamento de Análisis Económico, IUDR (ULL) Gustavo A. Marrero2 Departamento de Análisis Económico (ULL) Programa de Investigación Energía y Cambio Climático Focus-Abengoa-Fedea
ABSTRACT:
Using a panel data set for Spanish regions between 1998 and 2006, we study the factors explaining per capita fuel consumption for road transport at the macroeconomic level. The contributions of the article are the following. First, we specify a dynamic panel data (DPD) model for gasoline and diesel consumption. Second, we properly apply estimation techniques based on the system Generalized Methods of Moments (GMM) procedure of Arellano and Bover (1995). We find that more traditional estimation procedures (pooling-OLS, the Within-Group or the first difference GMM), which might generate bias estimate in a DPD framework, produce important differences that may even change policy recommendations. Finally, we find important differences between the results for the gasoline and the diesel model. While the estimated equation correctly fits the gasoline consumption behavior, results emphasizes the need to specify a different model for aggregate diesel consumption, which must include aditional determinants than those traditionally used in fuel consumption models.
JEL: R41, O13, O56
KEY WORDS: Fuel consumption, road transport, Dynamic Panel Data model, GMM estimates.
Correspondence to: Rosa M. González-Marrero, Departamento Análisis Económico, Universidad de La Laguna, Campus de Guajara, La Laguna, 38071 Santa Cruz de Tenerife; Fax: 922 317204, Tfno: 922317113, rmglzmar@ull.es
2 Gustavo A. Marrero acknowledges financial support of FOCUS-ABENGOA through Programa de Investigación de Energía y Cambio Climático Focus-Abengoa-Fedea.
1. Introduction
The transport sector in Spain represented almost 40% of final energy consumption in 2007, the highest between most important European countries,3 and the road transport accounts for almost 80% of the total energy consumed by the sector in 2006 (Spanish Government, 2009). Moreover, it generates about 25% of total emissions, with road transport contributing the most to said emissions. Having such a large road transport sector poses a serious roadblock for Spain to reaching the goals set by the Kyoto Protocol and the recently proposed 20/20/20 plan. Hence, fuel consumption policies in road transport and environmental policies should be intimately related.
There is an extensive literature that characterizes the consumption of energy in road transport. For example, see Baltagi and Griffin (1983, 1997), Schipper et al. (1992) and Johansson and Schipper (1997) for OECD countries, Baltagi and Griffin (1984) for the U.S., Mazzarino (2000) for Italy, Kwon (2005) for the United Kingdom, Polemis (2006) for Greece, Tapio et al. (2007) for the EU-15, Zervas (2006) for Ireland, Alves and Bueno (2003) for Brazil, Samimi (1995) for Australia, Nicol (2003) for Canada, Ramanathan (1999) for India, Koshal et al. (2007) for Japan, and Belhaj (2002) for Morocco, among others. For the case of Spain, the existing literature is scarce and mainly focuses on fuel determinants using micro data. For example, Romero-Jordán et al. (2010) estimates price and income elasticities of fuel demand for passenger transport, while Labeaga and López-Nicolás (1997) and Labandeira and López-Nicolás (2002) focus on analyzing the impact of taxes on overall consumption. In this paper, we focus on macroeconomic determinants.
Using a panel data framework, this article studies the factors explaining per capita fuel consumption for road transport for Spanish regions between 1998 and 2006 at the macroeconomic level. It is important to notice that these results might differ from those factors explaining fuel consumption at the microeconomic level (Puller and Grrening, 1999; Kayser 2000; Nicol, 2003). Following the arguments of Shipper et al. (1993), when data come from household surveys (micro data), the collected information is about automobile fuel consumption (mostly of private-use), and it is not about trucks or busses (mostly of public or freight transport vehicles). However, macro data of gasoline and diesel consumption do not distinguish between private-use vehicles, trucks or busses. The problem is that while most of gasoline consumption is absorbed by private-use vehicles, diesel consumption is more distributed between private-use vehicles, the public transport and freight transport vehicles (Sterner, 2007). This fact may drive that results for the diesel consumption model are very different depending on whether we use micro or macro data.
3 According to the EC(2009) report, only Luxemburg, Malta and Cyprus show a larger share of transport in final consumption.
Another relevant aspect of this article is the distinction made between gasoline and diesel consumption. Polemis (2006), Zervas (2006) and Labandeira and López-Nicolás (2002), among others, have already emphasized the importance of making this distinction. For the period analyzed, this distinction is especially relevant for the case of Spain, since it is probably, along with France, one of the countries in which the dieselization process has been the most significant during this period.
A final contribution of this paper is to write the empirical fuel consumption model in a dynamics panel data (DPD) framework.5 With regards to the estimation procedure, we use the one-step system GMM estimator proposed by Arellano and Bover (1995) and developed by Blundell and Bond (1998), which handles the endogeneity and the weak instruments problems arised in these types of models, as we will discuss in Section 2. In the growth literature, Forbes (2000), Shioji (2001), Levine et al. (2000) and Bond et al. (2001), among others, have usede the one-step system GMM estimator that we consider in the paper. However, there are few exceptions in the energy and transport literature that seriously consider the weakness of traditional methods in estimating DPD models. For example, Halkos (2003), Gang (2004) and Metcalf (2008) address the endogeneity problem and use the first difference GMM estimator, but this method does not consider the weak instruments problem of this procedure when time series are persistent (Blundell and Bond, 1998), which is the case for aggregate fuel consumption. Huang et al. (2008), which revisits the causal relationship between energy consumption and GDP, and Marrero (2010), which studies the relationship between emissions and energy in Europe, are exceptions that properly address both the endogeneity and the weak instruments problems by considering the system GMM approach.6
4 For instance, the diesel to gasoline consumption ratio in road transport in Spain, which was 1.7% in 1998, had risen to 3.8% by 2006.
5 The application of panel data estimation techniques to the study of energy demand was pioneered by Balestra and Nerlove (1966) and followed by many others authors. For a discussion on panel data estimation techniques and applications, see Baltagi (2001) and Baltagi and Baldev (1992), among others.
As explanatory variable, the diesel and the gasoline model both include traditional factors in the transport literature: the real prices of gasoline and diesel, the real regional GDP per capita, the fleet of gasoline and diesel vehicles per capita (the motorization rate by fuel type) and the total fleet divided by the total kilometers of road as a proxy for the saturation of the road network. We find that most explanatory variables are significant – and with the appropriate signs according to economic theory - in explaining the evolution of gasoline consumption. However, we find that most of these factors are not significant at explaining the recent evolution of diesel consumption between Spanish regions, which is consistent with the idea, previously exposed, that diesel consumption is heavily used by the professional transport and, as a consequence, explanatory factors can differ from those of gasoline used in private vehicles.7 With respect to the gasoline model, we find that the alternative estimation procedures produce important differences that may even change policy recommendations, which highlights the need to use appropiate econometric techniques when dealing with DPD models.
6 In order to show the importance of considering this estimation method, we follow Blundell et al. (2000), and compare the system GMM estimates with respect to alternative, more traditional methods – pooling-OLS, the within groups, the first difference GMM of Arellano and Bond (1991).
This paper is structured as follows. Section 2 presents the DPD fuel consumption model, describes the data used in the analysis and briefly comments on the system GMM estimation approach. Section 3 estimates the gasoline and diesel model and discusses the main results. Finally, Section 4 presents the main conclusions.
2. Fuel consumption for Spanish regions: dynamic model, data and econometric method In this section we present a dynamic framework for fuel consumption, show a Dynamic Panel Data (DPD) model for gasoline and diesel consumption, describe the data set used and comment main aspects of the estimation procedure.
2.1. A dynamics framework for fuel consumption
A common starting point in the transport literature is to assume that aggregate fuel can be characterized like any aggregate products demand, that is, as a function of a measure of real income and fuel prices.8 Specifically, using a reduced-form specification, aggregate fuel (gasoline or diesel) desired demand at time , can be assumed to be a lineal logarithm function of real gasoline price, real diesel price, , real income, and an error term, , which is assumed to be identically and independently distributed with zero mean and constant variance:
7 For Spanish regions, diesel and gasoline consumption statistics do not distinguish between passenger and freight transport use. While that could be a problem for diesel consumption, it may not be for gasoline consumption, because gasoline is majority consumed by private vehicles.
8 See Dahl, 1986; Dahl and Sterner, 1991; Sterner, (2007); Basso y Oum (2007), among many others.
\[\ln (E _ {t} ^ {*}) = \delta_ {0} + \delta_ {1} \ln (P G _ {t}) + \delta_ {2} \ln (P D _ {t}) + \delta_ {3} \ln (Y _ {t}) + u _ {t}\tag{1}\]
Fuel consumption and income variables would be expressed in per capital terms. The relationship between the desired demand and current consumption (Et ) is characterized by the following partial adjustment behavior,
\[\ln (E _ {t}) - \ln (E _ {t - 1}) = \theta (\ln (E _ {t} ^ {*}) - \ln (E _ {t - 1}))\tag{2}\]
where is previous demand; is parameter with range as (0,1) and is referred to as the coefficient of adjustment. This model states that current fuel consumption growth is proportional to the difference between the current demand and previous consumption. The coefficient θ indicates the speed of adjustment towards the desired level of demand. The time-lags have several explanations, including habit-persistence, slow turnover of car in the vehicle stock, and the time required for automobile manufacturers to change technologies in the cars they sell (Pollack and Wales, 1981).
Combining expressions (2) and (1), and rearranging terms, we obtain the following dynamic model for aggregate fuel consumption:
\[\ln \left(E _ {t}\right) = \alpha_ {0} + \beta \ln \left(E _ {t - 1}\right) + \alpha_ {1} \ln \left(P G _ {t}\right) + \alpha_ {2} \ln \left(P D _ {t}\right) + \alpha_ {3} \ln \left(Y _ {t}\right) + \varepsilon_ {t},\tag{3}\]
\[\text { where: } \alpha_ {0} = \theta \delta_ {0}, \beta = 1 - \theta , \alpha_ {1} = \theta \delta_ {1}, \alpha_ {2} = \theta \delta_ {2}, \alpha_ {3} = \theta \delta_ {3} a n d \varepsilon_ {t} = \theta u _ {t}.\]
This specification (3) is commonly called a lagged endogenous model, where lagged endogenous variable can be seen as representing the inertia of the system.10
2.2. A Dynamic Panel Data model for gasoline and diesel consumption
Taking (3) as the point of reference, we present and estimate a DPD model for gasoline consumption, GASO (the gasoline model) and another for diesel consumption, DISL (the diesel model),11 extended with some additional explanatory variables, such as the road network and the stock of vehicles12:
9 The partial adjustment model is a ad hoc model and, since the seminal paper of Houthaker et al. (1974), they have been used in many works about fuel demand (Baltagi et al, 2003, Pock, 2009).
10 For an interesting review of the approaches and methods that have been used in automobile fuel demand, see Basso and Oum (2007).
\[\begin{array}{l} \ln (G A S O _ {i, t}) = \alpha_ {i} + \beta \ln (G A S O _ {i, t - 1}) + \lambda_ {1} \ln (Y _ {i, t}) + \\ + \lambda_ {2} \ln (P G _ {i, t}) + \lambda_ {3} \ln (P D _ {i, t}) + \lambda_ {4} \ln (F L E E T G _ {i, t}) + \lambda_ {5} \ln (F L E E T D _ {i, t}) + \lambda_ {6} \ln (S A T _ {i, t}) + \varepsilon_ {i, t}, \end{array}\tag{4}\]
\[\begin{array}{l} \ln (D I S L _ {i, t}) = \alpha_ {i} + \beta \ln (D I S L _ {i, t - 1}) + \lambda_ {1} \ln (Y _ {i, t}) + \\ \lambda_ {2} \ln (P G _ {i, t}) + \lambda_ {3} \ln (P D _ {i, t}) + \lambda_ {4} \ln (F L E E T G _ {i, t}) + \lambda_ {5} \ln (F L E E T D _ {i, t}) + \lambda_ {6} \ln (S A T _ {i, t}) + \varepsilon_ {i, t}, \end{array}\tag{5}\]
where FLEETG and FLEETD are the existing fleet of gasoline and gasoil vehicles per capita (motorization rate), respectively, and is the total number of vehicles divided by total kilometers of road, which is a proxy of the saturation of the road network.13 Fixed factors are time-invariant and inherent to each region, and are not observed or not included in the model, such as geographical, social or local policy regional aspects or initial energy efficiency.14 Dynamic variables and control for conditional convergence across regions in terms of fuel consumption, as it is standard in the dynamic literature. 15 Finally, encompasses effects of a random nature which are not considered in the model, and it is assumed to have a standard error component structure:
11 Given the close relationship between energy and income, the dynamic specification is similar to that used in the convergence-growth literature (Barro and Sala-i-Martin, 1995, among many others]. Álvarez et al. (2005) and Marrero (2010) have adapted this dynamic approach to a economic-pollution model.
12 This is one of the alternatives that can be used for analyzing the demand of fuel. It consists of using an equation of fuel demand where some measure of the stock of vehicles is explicitly incluide (ie, Banfi et. al., 2005; Kayser, 2000; and Puller and Greening, 1999). An alternative procedure is to employ a system of simultaneous equations in which the demand for vehicles and the demand for fuel are analyzed jointly (ie. Belhaj, 2002 and Chandarisi, 2006).
13 In addition to the dynamic term, the other explanatory variables assumed in (4) and (5) are among those traditionally considered as indicators for characterizing the behavior of the road transportation sector [Eltony (1993), Bentzen (1994), Kirby et al. (2000), Alves and Bueno (2003), Polemis (2006)]. Moreover, we also distinguish between the total diesel and gasoline fleet, as proposed by Pock (2009).
14 Fixed effects, such as differences in the initial energy efficiency, would be omitted in a standard OLS pool regression, resulting in bias estimates (i.e., the β estimates is upward bias.) See Anderson and Hsiao (1982) and Hsiao (1986) for more details about this point.
15 Indeed, the interpretation of equation (1) depends on the level of β. A β smaller than one is consistent with the conditional convergence hiphothesis, which means that regions relatively close to their steady-state per capita fuel consumption levels will experience a slowdown in their consumption growth. In this case, αi and all αi explanatory variables affect to the steady-state the fuel consumption of region i is converging to. On the other hand, if β is greater than one, there is no convergence effect and αi and all regressors would measure differences αi in steady-state energy consumption growth rates. Estimated will be lower than one in all cases, hence we will β focus on the conditional convergence interpretation.
\[\mathbf {A 1}: E \left[ \varepsilon_ {i t} \right] = 0; E \left[ \alpha_ {i} \varepsilon_ {i t} \right] = 0; E \left[ \varepsilon_ {i t} \varepsilon_ {i s} \right] = 0, i = 1, \dots , N; t = 1, \dots , N a n d s \neq t.\]
We also consider a common assumption in DPD models [Arellano and Bond (1991)], which is that and are predetermined,
\[\mathbf {A 2}: E \left[ G A S O _ {i 1} \varepsilon_ {i t} \right] = E \left[ D I S L _ {i 1} \varepsilon_ {i t} \right] = 0, f o r i = 1, \dots , N a n d t = 2, \dots , T\]
2.3. Data
Table 1 shows annual growth rates between 1998 and 2006 for each variable used in models (4) and (5) at a regional level.
TABLE 1 ABOUT HERE
Per capita gasoline consumption fell between 1998 and 2006 in all Spanish regions at an average rate of 4.1%, while per capita diesel consumption increased for each region at a higher average rate of 5.7%. During this period, gasoline and diesel real prices increased by 2.8% and 4.8% in Spain, respectively. Although significant differences were noted in their growth rates in the time dimension, differences between regions are very small. The large time-volatility of fuel prices resulted from Spain’s enormous (nearly 100%) dependence on foreign oil, on the important fluctuations in the euro/dollar exchange rate, and on changes in fuel taxes and their repercussions on the final price.16 The data on per capita GDP growth rate showed an annual increase of 2.6% nationally, varying between 1.6% in Valencia and 3.8% in Extremadura. Overall, per-capita GDP growth showed a notable regularity among the different regions, as evidenced by the generalized slowdown between 2001 and 2003-2004 and the subsequent recovery until 2006 for most of the regions.17
16 For example, in 2001 the price in dollars of a barrel of Brent crude fell 14%, while in 2004 and 2005 it rose by 33% and 42%, respectively.
17 Spanish regional GDP is built on the basis of regional physical, economic indicators, such as retail sales, industrial production index, car sales, overnight stay of tourists, consumption of cement, etc, which are highly related to the domestic level of activity in each Region. De la Fuente (2002), among others, has used the Spanish regional GDP dataset to study the source of convergence in Spain. Petrol prices are retail prices and deflated by each regional CPI .For more details about how petrol prices are determined in Spain, see Perdiguero (2006). GASO and DISL are measured in Kilo-tones per inhabitant.
The motorization rates by type of fuel (gasoline fleet pc and diesel fleet pc) show a trend similar to those of their associated fuel consumption series. For most regions, it decreased for gasoline vehicles (2.0% in Spain), while it increased for diesel (9.2% in Spain). As is the case with fuel consumption, the dieselization process has changed the vehicle fleet composition in Spain: the ratio of diesel to gasoline vehicles, which was 43% in 1998, had risen to 105% by 2006. However, the intensity of this substitution process varied greatly depending on the Region. Thus, Madrid is the region that experienced the most significant substitution process from gasoline to diesel vehicles, with a decline in the motorization rate for gasoline vehicles of 4.2% and an increase for diesel of 10.6%. Other regions show a different pattern. For example, in Extremadura, the motorization rate for gasoline vehicles remained relatively stable, while that of diesel experienced the largest increase, within all Spanish regions, of 12% per year. As for the number of vehicles versus road network kilometers (the saturation level of road), it increased in all regions, showing the highest average annual growth rate in Castilla La Mancha (5.4%) and the lowest in Aragon (2.1%).
2.4. Estimation procedure
Traditional procedures for estimating a DPD models like (4) and (5) (i.e., fixed or random effects methods or pooling-OLS) are known to be unsuitable [Anderson and Hsiao (1982); Hsiao (1986)]. Holtz-Eakin et al. (1988) and Arellano and Bond (1991) propose an alternative approach, where first differences in the regression equation are taken to remove unobserved time-invariant country specific effects and then particular moment conditions for lagged variables are exploited to find a set of instruments and construct a GMM-based estimator. Their GMM approach (GMM-DIF) allows us to handle endogeneity problems. However, the GMM-DIF approach shows important bias problems in small sample when variables are hihgly persistent, which is the case of aggregate fuel consumption and other macroeconomic variables. Under these circumstances, the instruments used in the GMM-DIF estimator have proven to be weak and the first difference estimator is poorly behaved, and Arellano and Bover (1995) and Blundell and Bond (1998) proposed an alternative GMM procedure which might overcome the weak instruments problem. This procedure estimates a system of equations in both first-differences and levels, where the instruments in the level equations are lagged first differences of the variables. In this paper we use the one-step system GMM estimator. In contrast to the two-step version, the one-step GMM estimator has standard errors that are asymptotically robust to heteroskedasticity and have been found to be more reliable for finite sample inference.18
3. Gasoline and diesel model results
Our goal of this section is twofold. First, we want to emphasize the importance of considering an appropriate quantitative approach when estimating a dynamic fuel consumption model; second, we point out the differences between the gasoline and the diesel model results.
3.1. Alternative estimation procedures
We first find evidence supporting the good properties of the system GMM estimates. We follow the practical rule proposed by Blundell et al. (2000), and compare the results of alternative estimation methods for the gasoline and the diesel model: the OLS pooling estimates (OLS-POOL), the Within Group estimates (WG), the first-difference GMM approach (GMM-DIF) of Arellano and Bond (1991) and the system GMM method (GMM-SYS). For the GMM-based methods, we employ the one-step GMM estimator, with heteroskedasticity-consistent asymptotic standard errors.19 Tables 2 and 3 show the results for the gasoline and diesel model, respectively, for all these alternative methods. The p-value of the t-significance test associated with each parameter is shown. We also show standard specification tests for each model. Notice that the Haussman test rejects the null hypothesis of random effects at any standard level of significance. For any GMM-based estimates, we show the m1 and the m2 tests and conclude that moment conditions underlying GMM estimates seem to be robustly supported.20
18 See Blundell and Bond (1998), Blundell et al. (2000) or Baltagui (2001), among others. See the technical appendix for more details about this point.
INSERT TABLE 2 AND 3 ABOUT HERE
Based on the results shown in Tables 2 and 3 and following the practical rule proposed by Blundell et al. (2000), OLS-POOL seems to give an upward-biased estimate of the β coefficient (0.853 for the gasoline and 0.998 for the diesel model), while WG appears to give a downward-biased estimate of this coefficient (0.387 for the gasoline and 0.495 for the diesel model). Using GMM-DIF, the β coefficient is barely lower than the WG estimates, suggesting the possibility of important finite sample bias due to the weak instruments problem [Blundell and Bond (1998)]. This comparison also highlights how the estimated coefficients of the remaining regressors, which are our main interest, differ among the alternative procedures. Hence, using a method resulting in bias estimates (the OLS-POOL, WG or the GMM-DIF) might lead to misleading conclusions. For example, the coefficients associated with the network saturation variable in the gasoline and diesel model are not significant under the WG and GMM-DIF estimates, while it is significant under the GMM-SYS procedure; for the diesel model, the per capita GDP variable is significant under the WG, while it is not under the GMM-SYS and OLS-POOL; the magnitude of the estimated price-elasticities are smaller under the GMM-SYS than under the WG estimates in the diesel model.
19 For a given cross-sectional sample size, the use of too many instruments in models with endogenous regressors may result in seriously biased estimates (Álvarez and Arellano, 2003; Arellano and Bond, 1998). Hence, even when computing speed is not an issue, these authors recommend not using the entire series history as instruments. We use instruments up to t-3. Including more lags does not change results significantly.
20 The most frequently used tests to validate the assumptions underlying GMM methods are the m1, m2 and Sargan tests. The m1 and m2 tests are based on the standardized average residuals autocovariance, which are asymptotically N(0,1) distributed under the null hypothesis of no autocorrelation. The Sargan test, in contrast, is distributed chi-squared with degrees of freedom equal to the number of moment restrictions minus the number of parameters, estimated under the null hypothesis that moment conditions are valid. However, the Sargan test is less meaningful since it requires that the error terms be independently and identically distributed, which is not expected in our case. Hence, we will consider primarily the m1 and m2 tests.
In summary, this comparison suggests that the WG estimates are severely biased, that there exists a problem with weak instruments and hence that the GMM-DIF is biased similarly to WG, and that the GMM-SYS approach is a convenient way to overcome the weak instruments problem. This conclusion is an important contribution of the paper, and not always properly considered in the related literature.
3.2. System-GMM estimation results
Focusing on the second goal of the paper, we compare the results of the gasoline and diesel models, and we find important differences in the magnitude and significance of the variables. We will mainly focus our attention on the one-step GMM-SYS estimates from now on. In general, we find that the coefficients associated with the explanatory variables are less significant in the diesel model than in the gasoline model.21 This result could be due to different reasons. Following Sterner (2007), the diesel is used heavily in professional transport (buses, heavy trucks, etc.) with different explanatory factors than those used for private use. Moreover, the dieselization process that has taken place in Spain since 1994, aproximatly, has resulted in diesel consumption being exposed to important regulatory factors, which might have modified expected results (Burguillo et al., 2009). Likewise, since fuel consumption data include both passenger and freight transport, and heavy trucks mainly uses diesel, the diesel model must distinguish between consumption in the passenger and freight transport sector. But, diesel consumption statistics for Spanish regions do not distinguish between consumption in the passenger and the freight transport sectors. Moreover, there are also specific technical characteristics of diesel vehicles, such as a different engineering development, which should also lead to introduce additional dieselspecific variables in the diesel model. Therefore, it would be necessary to further study the characteristics of diesel consumption, which goes beyond the objectives of this work, but they represent a promising extension for future research.
21 For the diesel model, we try with alternative specifications, considering variables such as the percentage of heavy trucks in the total fleet, or the percentage of buses. But these alternative specifications also fail in explaining the recent evolution of diesel consumption in Spain, being in line with the results obtained by Baltagi and Griffin (1983).
For the diesel and gasoline model, we next comment the estimations for each variable. The parameters estimated for the variable and the are positive and lower than one at the 1% level of significance. The estimate is 0.867 for the diesel model and 0.558 for the gasoline model. Hence, the evidence for conditional convergence is significant in both cases, though it is greater for the gasoline case. The estimates indicate that the rate of conditional convergence for the per capita fuel consumption ratio is about 13% for diesel consumption and about 44% for gasoline.
A common result in the gasoline and diesel model is that per capita GDP is not significant in explaining per capita fuel consumption. From Section 2.3, we showed how per capita fuel consumption and GDP experienced important increases from 1998 to 2006. However, it seems that the former evolved independently of per capita GDP when, for a given period, we compare different regions. For example, we find that regions with different per capita GDP levels shared similar fuel consumption patterns. This is the case of the Basque Country, with a large per capita GDP, Castilla and Leon, with intermediate per capita GDP, and Extremadura, with one of the smallest per capita real GDP in Spain. However, they have experienced a similar increase in per capita fuel consumption for the analyzed period (between 3.6% and 3.9%).23
22 A similar result is found in Pock (2009), which estimates a demand function for aggregate gasoline, and finds that, in some models, all coefficients are significant with the exception of the income term. He points out that this result may be due to multicollinearity problems. For the case of diesel demand, Burguillo et al. (2009) find the same result: the coefficient associated to real income was negative, but non-significant. Baltagi and Griffin (1997) obtain that the income elasticities are frequently insignificant in the long run.There exits other articles in the literature which estimates a model for aggregate energy consumption or CO2 emissions whose incomeelasticity estimation is very low, and even negative in some cases [Schmalensee et al. (1998), Judson et al. (1999), Holtz-Eakin and Selden (1995) or Marrero (2010)]. Baltagi and Griffin (1997) obtain the income elasticities in long run ar
Regarding other papers studying the relationship between gasoline demand and income, Dahl and Sterner (1991) showed that short-term income elasticity on gasoline demand varied between 0.30 and 0.52 in the different studies they considered.24 However, notice that our GDP elasticity under the WG estimate was significant and about 0.44 for diesel, which is indeed consistent with the Dahl-Sterner range; this, however, is the result of a bias estimation procedure, as discussed above. With this example, we are not claiming that Dahl-Sterner estimates are wrong. In fact, differences between their estimates and ours may only be due to differences in the data sample used. We are just stressing the importance of considering an appropriate estimation approach for handling fuel consumption models, because, otherwise, results could lead to misleading conclusions.
Regarding the real price of fuel, the GMM-SYS procedure estimate for its elasticity is negative and significant for the case of gasoline, though its magnitude is well below one (- 0.29). This result confirms the evidence that the elasticity of the demand price for gasoline is low in the short term, as verified by, among others, Kayser (2000) with data for the United States and Baltagi and Griffin(1997) with data for OECD countries. Results at the international level place the price-elasticity in the -0.2 and -0.3 range [Dahl and Sterner (1991)], which is consistent with our results. This result indicates that fuel demand is highly inelastic, at least at current price levels. Moreover, it supports a result commonly discussed in the literature: fuel taxes are convenient for increasing fiscal revenues, but they are not effective enough to reduce fuel consumption [Kirby et al. (2002)].
23 We are aware that the short time dimension of the data can influence this result. Also, the transport network is very different across regions as it is the availability of public transport services. The response to income of consumers in predominantly rural areas could be different than in predominantly urban areas, at least because of the substitution possibilities. Once again, we do not have data that distinguish between consumption in rural and urban areas. However, the importance of this result deserves further reseach, which goes beyond the purpose of this paper.
24 More recently, studies such as that by Koshal et al. (2007) gave values of 0.29. For a detailed review, see Graham and Glaister (2002) and Goodwin et al. (2004).
In addition, the real price of diesel is significant in explaining short-term changes in per capita gasoline consumption, although its parameter is small (0.21). This last result is also consistent with the transport literature due to the strictness that exists in substituting types of vehicles in the short term [Polemis (2006)]. We should emphasize that the recent intensive switch from gasoline to diesel vehicles has been basically due to regulatory reasons (dieselization) rather than to a change in the price of the alternative fuel. For the case of the diesel consumption model, neither its own price nor the price of gasoline is significant. As commented above, differences between the markets of the gasoline (most for passenger transport use) and the diesel (passenger and freight transport are important) could explain this result for the diesel model. Once again, this finding also highlights the need to consider different models for gasoline and diesel consumption, and to go further in the research of the determinants of diesel consumption in Spain.
The remaining variables are specific to the road transport sector and include relevant aspects that can affect fuel consumption. The per capita diesel and gasoline fleet variables show the motorization rate for each type of fuel vehicle. For the gasoline model, the coefficient of the per capita gasoline fleet variable is highly positive and significant (0.64), while the coefficient of the per capita diesel fleet is negative but much smaller in magnitude (-0.08). Regarding these variables, results for the diesel model are controversial, as with other variables.25
The coefficients of the measure of the degree of saturation of the road network (the ratio between total fleet and road network) are negative and significant in both models. Moreover, their estimates are similar: -0.048 for diesel and -0.059 for gasoline. The fact that estimated coefficients are similar in both models is a clear indication that road saturation affects both diesel and gasoline vehicles in a similar way. This result suggests that a reduction in road congestion promotes mobility, which may induce an increment in per capita fuel consumption.26
25 Further investigation on this important topic would constitute a prominent extension of this paper.
5. Final Remarks
This paper has analyzed the factors explaining the fuel consumption for road transport in Spain in a dynamic panel data framework in a macroecomic perspective. Two features on this study are the use of a balanced panel using regional data and the distinction between gasoline and diesel. As explanatory variables, we considered real GDP and fuel prices, which are the most commonly used in the related literature, as well as other relevant and different variables, such as the motorization rate and the congestion of the road network.
A novelty in this paper is the use of a the one-step system GMM estimation method of Arellano and Bover (1995) and Blundell and Bond (1998), which has been shown to solve many of the problems that arise in traditional panel data procedures. When compared with the system GMM results, we found that traditional panel data estimation procedures (the withingroup estimates, OLS-pooling or the first difference GMM approach of Arellano and Bond, 1991) might exhibit significantly biased estimates, which might even change policy recommendations.
For the sample used, we found that most explanatory variables are significant in explaining the evolution of gasoline consumption in Spain, while diesel consumption was found to have less dependence on most of these factors. These conspicuous differences between the results for the gasoline and diesel models imply that the conclusions derived from an analysis of overall fuel consumption could be misleading.
26 As noted by Goodwin (1996), improving the infrastructure has an induced effect on the demand for transport. Moreover, Cervero and Hansen (2002) provided empirical evidence of the existence of a direct relationship between investing in roads and the demand for transport, namely that an expansion of infrastructure generates demand for transport, which in turn induces the creation of infrastructure.
Our estimates confirm that the price elasticity of demand for fuel consumption is low – even negligible for diesel - in the short term, which supports the view that the policy of taxing fuel has little effect on reducing fuel consumption. Our results are also consistent with the evidence of small cross price elasticities for gasoline and diesel consumption. Another important finding of this work is the negative and significant relationship between the degree of saturation of the road network and both types of per capita fuel consumption. This result shows that reducing road network saturation – i.e., by increasing the road network -, could promote mobility and a higher transport demand (“induced travel demand”), which can favour higher levels of fuel consumption.
Several are the implications of our research. On the one hand, we emphasize the need to revisit DPD fuel consumption results obtained with traditional procedures, and show the relevance of considering a suitable estimation method; ant the other, it necessary to investigate the determinants of diesel consumption using a different model than that used for gasoline consumption, which is a promising extension of this paper.
References
- • Alvaes D., Bueno R. Short-run, long-run and cross elasticities of gasoline demand in Brazil. Energy Economics 2003, 25, 191-199.
References
- • Álvarez J., Arellano M. The time series and cross-section asymptotics of dynamic panel data estimators. Econometrica 2003, 71(4), 1121-1159.
References
- • Álvarez F., Marrero G. A., Puch L. A. Air pollution and the Macroeconomy across European Countries. FEDEA 2005 Working Paper 2005-26.
References
- • Anderson T.W., Hsiao C. Formulation and estimation of dynamic models using panel data. Journal of Econometrics 1982, 18, 47-82.
References
- • Arellano M., Bover O. Another look at the instrumental-variable estimation of errorcomponents models. Journal of Econometrics 1995, 68, pp. 29–52.
References
- • Arellano M., Bond S. Some Tests of Specification for Panel Data: Monte Carlo Evidence and an Application to Employment Equations. Review of Economic Studies 1991, 58, 277-297.
References
- • Balestra, P and Nerlove M. Pooling cross-section and time-series data in the estimation of a dynamic model: The demand for natural gas, Econometrica 1966, 34, 585-612.
References
- • Baltagi, B. Econometric Analysis of Panel Data, second edition, Chichester: John Wiley&Sons, Ltd., 2001.
References
- • Baltagi, B and Baldev, R. A survey of recent theoretical developments in the econometrics of panel data, Empirical Economics 1992, 17(1), 85-109.
References
- • Baltagi B., Bresson G., Griffin J. and Pirotte A. Homogeneous, Heterogeneous or Shrinkage Estimators? Some Empirical Evidence from French Regional Gasoline Consumption. Empirical Economics 2003, 28, 795-811.
References
- • Baltagi, B. and Griffin, J. Gasoline demand in the OECD: An application of pooling and testing procedures, European Economic Review 1983, 22(2), 117-137.
References
- • Baltagi, B. and Griffin J. U.S. gasoline demand: What next?. The Energy Journal 1984, 5(1), 129-140.
References
- • Baltagi, B. and Griffin J. Pooled estimators vs. heterogenous counterparts in the context of dynamic demand for gasoline. Journal of Econometrics 1997, 77, 303-327.
References
- • Banfi S., Filippini M. and Chunt L. Fuel tourism in border regions: The case os Switzerland. Energy Economics 2005, 25, 5, 689-707.
References
- • Barro R. J., Sala-i-Martin X. Economic Growth, Advanced Series in Economics, McGraw-Hill, 1995.
References
- Basso, L and Oum T. Automobile fuel demand: A Critical Assessment of Empirical Methodologies. Transport Review 2007, 4, 449-484.
References
- • Belhaj M. Vehicle and fuel demand in Morocco. Energy Policy 2002, 30, 1163-1171.
References
- • Bentzen J. An Empirical analysis of gasoline demand in Denmark using cointegration techniques. Energy Economics 1994, 16, 139-143.
References
- Blundell R.W., Bond S.R., Windmeijer F. Estimation in dynamic panel data models: improving on the performance of the standard GMM estimator. In: Baltagi, B., Editor, 2000. Nonstationary Panels, Panel Cointegration, and Dynamic PanelsAdvances in Econometrics Vol. 15, JAI Press, Elsevier Science, Amsterdam, pp. 53–91
References
- Blundell R., Bond S. Initial Conditions and Moment Restrictions in Dynamic Panel Data Models. Journal of Econometrics 1998, 87, 115-143.
References
- • Bond S. Dynamic panel data models: a guide to microdata methods and practice. Portuguese Economic Journal 2002, 1, 141-162.
References
- • Bond S., Hoeffler A., Temple J. GMM Estimation of Empirical Growth Models. Economics Papers 2001-W21, Economics Group, Nuffield College, University of Oxford.
References
- • Burguillo-Cuesta M., Jorge-Garía I. and Romero D. Dieselization of passenger cars in the EU-15 in the nineties: environmental implications for transport policy . 8th International Conference of the European society for Ecological Economics, June 2009, Slovenia.
References
- • Cervero R., Hansen M. Induced travel demand and induced road investment. A Simultaneous equation analysis. Journal of Transport Economic and Policy 2002, 36(3), 469-490.
References
- • Chandrasiri (2006) : Demand for road-fuel in a small developing economy: The case of Sri Lanka Energy Policy 2006, 34, 1833-1840.
References
- • Dahl C., Sterner T. Analyzing gasoline demand elasticities: a survey. Energy Economics 1991, 13, 203-310.
References
- • De la Fuente A. On the source of convergence: a close look at the spanish regions. European Economic Review 2002, 46 (3), 569-599.
References
- • Dorand H. E., Schmidt P. GMM estimator with improved finite sample properties using principal components of the weighting matrix, with an application to the dynamic panel data models. Journal of Econometrics 2006, 133(1), 387-409.
References
- • EC . Energy and Transport in figures. Statistical Pocketbook, 2009.
References
- • Eltony M. The demand for gasoline in the GCC: an aapplication of pooling and testing procedures. Energy Economics 1993, 18 (3), 203-209.
References
- • Forbes K. A reassessment of the relationship between inequality and growth. American Economic Review 2000, 90(4), 869-887.
References
- • Gang L. Estimating Energy Demand Elasticities for OECD Countries. A Dynamic Panel Data Approach. Research Department of Statistics Norway 2004, Dicussion Paper 373.
References
- • Goodwin P. Empirical evidence on induced traffic. Transportation 1996, 23(1), 35-54.
References
- • Goodwin P., Dargay J., Hanley M. Elasticities of road traffic and fuel consumption with respect to price and income: A review. Transport Review 2004, 24 (3), 275-292.
References
- • Grahan D. J., Glaister S. The demand for automobile fuel: a survey of elasticities. Journal of Transport Economics and Policy 2002, 36, 1-26.
References
- • Halkos G.E. Environmental Kuznets Curve for sulphur: evidence using GMM estimation and random coefficient panel data models. Environment and Development Economics 2003, 8, 581-601.
References
- • Holtz-Eakin D., Newey W., Rosen H. S. Estimating vector autoregressions with panel data. Econometrica 1998, 56(6), 1371-1395.
References
- Holtz-Eakin and Selden M. Stoking the Fires? CO2 Emissions and economic groyth. Journal of Public Economics 1995, 57, 85-101.
References
- • Houtthakker H., Verleger P^. and Sheekan D. Dynamic demand for gasoline and residential electricity. American Journal of Agricultural Ecnomics 1974, 56, 412-418.
References
- Hsiao C. Analysis of panel data. Econometric Society monographs 11, Cambridge. University Press, 1986.
References
- • Huang B.N., Hwang M. J., Yang C. W. Causal relationship between energy consumption and GDP growth revisited: A dynamic panel data approach. Ecological Economics 2008, 67 (1), 41-54.
References
- • Johansson O., Schipper L. Measuring long-run automobile fuel demand: separate estimations of vehicle stock, mean fuel intensity, and mean annual driving distance. Journal of Transport Economic and Policy 1996, 31(3), 277-292.
References
- • Judson R., Schmalensee R. and Stoker M. Economic Development and the Structure of the Demand for commercial energy. Energy Journal 1999, 20 (2), 29-57.
References
- • Kayser H. A. Gasoline demand and car choice estimating gasoline demand using household information. Energy Economics 2000, 22, 331-348.
References
- • Kirby H. R., Hutton B., McQuaid R. W., Raeside R., Zhang X. Modelling the effects of transport policy levers on fuel efficiency and national fuel consumption, Transportation Research Part D 2000, 5, 265-282.
References
- • Koshal R. K., Manjulika K., Yuko Y., Sasuke M., Keizo Y. Demand for gasoline in Japan. International Journal of Transport Economics 2007, 34, 351-367.
References
- • Kwon T-H. The determinants of the changes in car fuel efficiency in Great Britain (1978-2000). Energy Policy 2005, 2, 261-275.
References
- • Labeaga J. M., López-Nicolás A. A study of petrol consumption using Spanish panel data. Applied Economics 1997, 29, 795-802.
References
- • Labandeira X., López-Nicolás A. La imposición de los carburantes de automoción en España; algunas observaciones teóricas y empíricas. Hacienda Pública Española 2002, 160-1, 177-210.
References
- • Levine R., Loayza N., Beck T. Financial intermediation and growth: causality and causes. Journal of Monetary Economics 2000, 46, 31-77.
References
- • Marrero G.A. Greenhouse gases emissions, growth and the energy mix in Europe: a dynamic panel data approach. Working Paper 2009-16, FEDEA.
References
- • Mazzarino M. The economics of the greenhouse effect: evaluating the climate change impact due to the transport sector in Italy. Energy Policy 2000, 28, 957-966.
References
- Metcalf G. An empirical analysis of energy intensity and its determinants at the state level. The Energy Journal 2008, 29, 1-26.
References
- • Nicol C. J. Elasticities of demand for gasoline in Canada and the United States. Energy Economics 2003, 25, 201-214.
References
- • Perdiguero G. J. Dinámica de precios en el mercado español de gasolina: un equilibrio de colusión tácita, Documento de Trabajo de Funcas 2006, nº 253.
References
- • Pollack R.A. and Wales J. Demographic variables in demand analysis. Econometrica 1981, 49, 1533-1551.
References
- • Polemis M. L. Empirical assessment of the determinants of road energy demand in Greece. Energy Economics 2006, 28, 385-403.
References
- • Puller S. and Greening L. Hosehold adjustment to gasoline price change: an analysis using 9 years of US suvey data. Energy Economics 1999, 21, 37-52.
References
- • Ramanathan R. Short and long-run elasticities of gasoline demand in India: an empirical analysis using cointegration techniques. Energy Economics 1999, 21, 321- 330.
References
- • Romero-Jordán D, Rí R., Jorge-García M and Burguillo M. Price and income elasticities of demand for passenger transport fuel in Spain. Implications for public policies. Energy Policy 2010, 38, 3898-3909.
References
- • Samimi R. Road transport energy demand in Australia. Energy Economics 1995, 17, 329-339.
References
- • Schipper L., Steiner R., Duerr P., An F., Strom S. Energy use in passenger transport in OCDE countries: Changes since 1970. Transportation 1992, 19, 25-42.
References
- • Schipper L., Figueroa M., Espey M. and Price L. Mind the gap: the vicious circle of measuring automobile fuel use. Energy Policy 1993, 21, 1173-1190.
References
- • Schipper L., C. Marie-Lilliu y L.Fulton, Diesels in Europe. Analysis of Characteristics, usage patterns, energy savings and CO2 emission implications, Journal of Transport Economics and Policy 2002, 36(2), 305-340.
References
- Schmalensee R., Stoker M. and Judson R. Worl Carbon dioxide emissions:1950-2050. Review of Economics and Statistics 1998, 80(1), 15-27.
References
- • Shioji E. Public capital and economic growth: a convergence approach. Journal of Economic Growth 2001, 6, 205-227.
References
- • Spanis Government. Estrategia española de movilidad sostenible. Madrid 2009.
References
- • Sterner T. Fuel taxes: An important instrument for climate policy. Energy Policy 2007, 35, 3194-3202.
References
- • Tapio P., Banister D., Luukkanen J., Vehmas J., Willamo R. Energy and transport in comparison: Immaterialisation, dematerialisation and decarbonisation in the EU15 between 1970 and 2000. Energy Policy 2007, 35, 433-451.
References
- • Windmeijer F. A finite sample correction for the variance of linear efficient two-step GMM estimators. Journal of Econometrics 2005, 126, 1, 25-51.
References
- • Zervas E. benefit from the increasing percentage of diesel passenger cars. Case of Ireland. Energy Policy 2006, 34, 2848-2857.
Gasoline Model TABLE 1: LISTS THE AVERAGE ANNUAL VARIATION RATES FOR ALL VARIABLES FOR THE PERIOD 1998-2006
| REGION | consumption gasoline pc | consumption diesel pc | gasoline real price | diesel real price | per capita GDP | gasoline fleet pc | diesel fleet pc | total fleet/ road network |
| Andalusia | -3.44 | 6.19 | 2.76 | 5.09 | 2.97 | -1.48 | 11.47 | 4.84 |
| Aragón | -3.90 | 6.25 | 2.63 | 4.78 | 2.79 | -1.77 | 9.64 | 2.14 |
| Asturias | -3.44 | 5.51 | 3.01 | 5.18 | 2.99 | -1.36 | 8.23 | 2.39 |
| Cantabria | -3.78 | 6.67 | 2.93 | 5.03 | 3.07 | -1.30 | 9.46 | 4.28 |
| Castilla y León | -3.33 | 5.91 | 2.72 | 4.90 | 3.26 | -1.11 | 9.78 | 2.94 |
| Castilla La Mancha | -4.76 | 5.87 | 3.03 | 5.23 | 2.14 | -1.71 | 10.43 | 5.40 |
| Catalonia | -5.16 | 4.11 | 2.71 | 4.85 | 2.03 | -2.53 | 8.33 | 2.85 |
| Valencia | -4.33 | 5.13 | 3.01 | 5.05 | 1.64 | -2.33 | 7.85 | 3.72 |
| Extremadura | -2.89 | 6.96 | 3.03 | 5.06 | 3.75 | -0.82 | 12.01 | 3.94 |
| Galicia | -3.24 | 3.74 | 3.00 | 4.99 | 2.99 | -1.32 | 7.94 | 3.18 |
| Madrid | -6.72 | 6.42 | 2.91 | 3.86 | 2.13 | -4.15 | 10.63 | 3.19 |
| Murcia | -4.62 | 6.58 | 2.29 | 4.54 | 1.98 | -2.42 | 8.94 | 4.82 |
| Navarre | -3.33 | 6.04 | 2.64 | 4.60 | 2.65 | -2.62 | 7.38 | 2.74 |
| Basque Country | -3.99 | 6.28 | 2.49 | 4.59 | 3.26 | -2.01 | 7.62 | 2.86 |
| La Rioja | -4.18 | 4.07 | 2.36 | 4.43 | 1.72 | -2.66 | 8.08 | 3.26 |
| SPAIN | -4.07 | 5.72 | 2.77 | 4.81 | 2.63 | -1.97 | 9.19 | 3.50 |
Table 2: Estimates of the gasoline DPD model
| Traditional methods | GMM methods | |||||||
| OLS-POOL | WG-Fixed effects | GMM-DIF | GMM-SYS | |||||
| estimates | p-value | estimates | p-value | estimates | p-value | estimates | p-value | |
| Lag of gasoline consumption pc | 0.853 | 0.000 | 0.387 | 0.000 | 0.343 | 0.010 | 0.558 | 0.000 |
| Gasoline real price | -0.417 | 0.000 | -0.375 | 0.000 | -0.377 | 0.000 | -0.292 | 0.000 |
| Diesel real price | 0.175 | 0.021 | 0.181 | 0.007 | 0.186 | 0.000 | 0.212 | 0.000 |
| Real GDP pc | 0.009 | 0.622 | 0.241 | 0.244 | 0.293 | 0.254 | -0.011 | 0.777 |
| Gasoline fleet pc | 0.163 | 0.013 | 0.640 | 0.000 | 0.707 | 0.000 | 0.639 | 0.000 |
| Diesel fleet pc | -0.065 | 0.002 | -0.264 | 0.006 | -0.284 | 0.009 | -0.083 | 0.014 |
| Total fleet / Road Network | -0.023 | 0.002 | 0.006 | 0.957 | 0.005 | 0.964 | -0.059 | 0.000 |
| R2 | 0.959 | -- | 0.946 | -- | -- | -- | -- | -- |
| Hausman, random effect test | -- | -- | 53.68 | 0.000 | -- | -- | -- | -- |
| m1-test | -- | -- | -- | -- | -2.564 | 0.010 | -3.368 | 0.001 |
| m2-test | -- | -- | -- | -- | -1.022 | 0.307 | -0.510 | 0.610 |
Note: ‘WG’ is Within Groups estimation, OLS-POOL is OLS applied to the entire pool of data. For GMM estimates, we take as instruments the lagged levels of y and the endogenous regressors dated t-2 and earlier and the pre-determined regressors dates t-1 and earlier. We use the lagged difference of y and all regressors dated t-1 as additional instruments in the system GMM estimation. For the DIF-GMM and SYS-GMM, we report their one-step estimations. The null of the Haussman test is the existence of random effects. The null of the m1 and m2 test is the absence of first- and second-order serial correlation between
regressors and residuals, respectively. Number of regressors: 8; number of cross sections: 15 (all Spanish regions except Ceuta and Melilla, Balears and Canary islands); number of time periods: 9 (1998-2006); number of time periods adjusted for GMM-DIF and GMM-SYS: 6 (2001-2006).
Table 3: Estimates of the Diesel DPD model
| Traditional methods | GMM methods | |||||||
| OLS-POOL | WG-Fixed effects | GMM-DIF | GMM-SYS | |||||
| estimates | p-value | estimates | p-value | estimates | p-value | estimates | p-value | |
| Lag of gasoil consumption pc | 0.998 | 0.000 | 0.495 | 0.000 | 0.453 | 0.000 | 0.867 | 0.000 |
| Gasoline real price | -0.113 | 0.285 | -0.101 | 0.265 | -0.116 | 0.044 | -0.047 | 0.551 |
| Diesel real price | -0.049 | 0.502 | -0.083 | 0.205 | -0.075 | 0.100 | -0.027 | 0.666 |
| Real GDP pc | -0.004 | 0.852 | 0.445 | 0.030 | 0.483 | 0.196 | 0.044 | 0.431 |
| Gasoline fleet pc | 0.012 | 0.732 | -0.010 | 0.920 | 0.052 | 0.638 | 0.206 | 0.003 |
| Diesel fleet pc | -0.047 | 0.088 | 0.199 | 0.059 | 0.221 | 0.142 | 0.095 | 0.106 |
| Total fleet / Road Network | -0.001 | 0.947 | -0.019 | 0.852 | 0.021 | 0.882 | -0.048 | 0.023 |
| R2 | 0.986 | -- | 0.954 | -- | -- | -- | -- | -- |
| Hausman, random effect test | -- | -- | 37.206 | 0.000 | -- | -- | -- | -- |
| m1-test | -- | -- | -- | -- | -2.724 | 0.007 | -3.780 | 0.000 |
| m2-test | -- | -- | -- | -- | -1.288 | 0.198 | -0.827 | 0.409 |
| Note: See Note on Table 2. | ||||||||
Appendix: System GMM Estimation of DPD Models
The fixed effect treatment leads to the well known within group estimator (WG) [Hsiao (1986)], which has been applied to multiple frameworks. However, the within transformation in a panel dynamic model implies a correlation of order 1/T between the lagged dependent term and the error which leads to biased estimates [Anderson and Hsiao (1981); Hsiao (1986)]. In addition, a fuel consumption equation such as (4) or (5) suffer from endogeneity and, maybe, from measurement errors problems. For instance, gasoline and diesel prices are jointly determined with gasoline and diesel consumption. The WG method neither properly handle these problems.
Holtz-Eakin et al. (1988) and Arellano and Bond (1991), among others, point out these problems and propose a GMM-based estimation approach. The current response of these authors is to first difference the model equation, remove the fixed effect term and then use the following orthogonally conditions, which, under assumptions A1 and A2 (see Section 2 of the paper), are valid for the first difference model:
\[E \left[ y _ {i t - s} \Delta \varepsilon_ {i t} \right] = 0, t = 3, \dots , T a n d 2 \leq s \leq t - 1, f o r i = 1, \dots , N,\tag{3}\]
Regressors in the gasoline and diesel models are either endogenous (prices, GDP and registrations) or pre-determined (the road network and the vehicle fleets ratios). Assuming a similar condition to A2 but for the regressors in X,
\[\mathrm{A} 3: E \left[ x _ {i 1} \varepsilon_ {i t} \right] = 0, f o r i = 1, \dots , N a n d t = 2, \dots , T,\]
we have additional 0.5(T-1)(T-2) moment conditions,
27 In the case of exogenous regressors, additional moment conditions are available. See Arellano and Bond (1991) for more detail about this point.
\[E \left[ x _ {i t - s} \Delta \varepsilon_ {i t} \right] = 0, t = 3, \dots , T a n d 2 \leq s \leq t - 1, f o r i = 1, \dots , N,\tag{4.a}\]
for each endogenous regressor, and another (T-1)(T-2) moment conditions,
\[E \left[ x _ {i t - s} \Delta \varepsilon_ {i t} \right] = 0, t = 3, \dots , T a n d 1 \leq s \leq t - 1, f o r i = 1, \dots , N,\tag{4.b}\]
for each pre-determined regressor.
For the case of K=1 and endogenous regressor,28 we have a total of moment conditions. Conditions in (3) and (4.a) can be written more compactly as
\[E \left[ Z _ {i} ^ {\prime} \Delta \varepsilon_ {i} \right] = 0, i = 1, \dots , N,\tag{5}\]
where is a matrix, given by
\[Z _ {i D I F} = \left( \begin{array}{c c c c c} y _ {i 1} x _ {i 1} & 0 & \dots & & 0 \\ 0 & y _ {i 1} y _ {i 2} x _ {i 1} x _ {i 2} & \dots & & \dots \\ . & \dots & \dots & & \dots \\ . & \dots & \dots & y _ {i 1} \dots y _ {i T - 3} x _ {i 1} \dots x _ {i T - 3} & 0 \\ 0 & 0 & \dots & 0 & y _ {i 1} \dots y _ {i T - 2} x _ {i 1} \dots x _ {i T - 2} \end{array} \right)\tag{6}\]
These are the moment conditions exploited by the standard first-difference GMM estimator (GMM-DIF).
However, the GMM-DIF estimator has been found to have large finite sample bias and poor precision when the set of instruments is weak [Blundell and Bond (1998).], which is the case of our fuel consumption model. To deal with this problem, Arellano and Bover (1995) and Blundell and Bond (1998) assume additional conditions to A1, A2 and A3,
(i.e., K=1)
28 For ease of exposition, we restrict notation to the case of only one endogenous regressor (i.e., K=1). The extension of the general case is straighforward.
\[\mathbf {A 4}: E \left[ \eta_ {i} \Delta y _ {i 2} \right] = 0, i = 1, \dots , N,\]
\[\mathbf {A 5}: E \left[ \eta_ {i} \Delta x _ {i 2} \right] = 0, i = 1, \dots , N\]
which allows the use of other 2·(T-2) moment conditions for a model in levels,
\[E \left[ u _ {i t} \Delta y _ {i t - 1} \right] = 0, t = 3, \dots , T,\tag{7}\]
\[E \left[ u _ {i t} \Delta x _ {i t - 1} \right] = 0, t = 3, \dots , T\tag{8}\]
which stay informative even for high persistent time series. Their proposal consists in a stacked system of all (T-2) equations in first differences and all (T-2) equations in levels for , and combine restrictions (3), (4), (7) and (8) to form a linear system GMM estimator (GMM-SYS) based on the following instrument matrices:
\[Z _ {i} = \left( \begin{array}{c c} Z _ {i D I F} & 0 \\ 0 & Z _ {i S Y S} \end{array} \right),\tag{9}\]
with given by (6) and by
\[Z _ {i S Y S} = \left( \begin{array}{c c c c c} \Delta y _ {i 2} \Delta x _ {i 2} & 0 & \dots & & 0 \\ 0 & \Delta y _ {i 3} \Delta x _ {i 3} & \dots & & \dots \\ . & \dots & \dots & & \dots \\ . & \dots & \dots & \Delta y _ {i T - 2} \Delta x _ {i T - 2} & 0 \\ 0 & 0 & \dots & 0 & \Delta y _ {i T - 1} \Delta x _ {i T - 1} \end{array} \right);\]
Monte Carlo analysis has shown that using GMM-SYS greatly reduces the finite sample bias and improves the precision of the estimator in presence of weak instruments.29 The linear GMM estimator is given by , where, for the GMM-DIF,
\[\tilde {X} = \left[ \begin{array}{c} \Delta X _ {1} \\ ... \\ \Delta X _ {N} \end{array} \right]; \tilde {Y} = \left[ \begin{array}{c} \Delta Y _ {1} \\ ... \\ \Delta Y _ {N} \end{array} \right]; Z = \left[ \begin{array}{c} Z _ {1, D I F} \\ ... \\ Z _ {N, D I F} \end{array} \right],\]
while for the GMM-SYS case,
\[\widetilde {X} = \left[ \begin{array}{c} \Delta X _ {1} \\ \dots \\ \Delta X _ {N} \\ X _ {1} \\ \dots \\ X _ {N} \end{array} \right]; \widetilde {Y} = \left[ \begin{array}{c} \Delta Y _ {1} \\ \dots \\ \Delta Y _ {N} \\ Y _ {1} \\ \dots \\ Y _ {N} \end{array} \right]; Z = \left[ \begin{array}{c} Z _ {1} \\ \dots \\ Z _ {N} \end{array} \right].\]
For each GMM-FIF and GMM-SYS case, two different choices of result in two different GMM estimators. The one-step estimator sets
\[H _ {N, G M M 1} = \left(\frac {1}{N} \sum_ {i = 1} ^ {N} Z _ {i} ^ {\prime} H Z _ {i}\right) ^ {- 1},\]
where the H matrix is a (T-2) square matrix with on the main diagonal, -1 on the first offdiagonals and zeros elsewhere. The two-step GMM estimator uses
\[H _ {N, G M M 2} = \left(\frac {1}{N} \sum_ {i = 1} ^ {N} Z _ {i} ^ {\prime} \Delta \hat {u} _ {i} \Delta \hat {u} _ {i} ^ {\prime} Z _ {i}\right) ^ {- 1},\]
29 Indeed, Blundell and Bond (1998) and Bond et al. (2001) shows that an optimal combination of differenced and level equations allow us to calculate a GMM estimator using the full set of linear moment conditions implied by assumptions A1-A5.
where estimated residuals are from a consistent one-step estimator (i.e., the one-step), which is an asymptotically efficient GMM estimator.
Under spherical disturbances, GMM1 and GMM2 are equivalent in the first-difference model. Otherwise, GMM2 is more efficient. However, Monte Carlo studies have shown that the efficiency gains of the two-step estimator are generally small. It also has the problem of converging to its asymptotic distribution relatively slowly. Hence, in finite samples, its variance-covariance matrix can be seriously biased. Moreover, for the case where the total number of instruments is large relative to the cross-section dimension of the panel, there may be computational problems in calculating the two-step estimates and serious estimation errors may arise [Arellano and Bond (1998); Doran and Schmidt (2006)]. With this in mind, most empirical works with a relatively small cross-section dimension report results of the one-step GMM estimator, which has standard errors that are asymptotically robust to heteroskedasticity and have been found to be more reliable for finite sample inference [Blundell and Bond (1998), Blundell et al. (2000); Windmeijer (2005); Bond (2002)]. This is the strategy considered in this paper.
There exist some tests to validate the assumptions underlying GMM methods. The standard approach for testing the validity of the moment conditions in GMM estimation is the Sargan test of overidentifying restrictions and the m2 second-order serial correlation test [Arellano and Bond (1991)]. Under the null hypothesis that moment conditions are valid, the Sargan test is distributed chi-squared with degrees of freedom equal to the number of moment restrictions minus the number of parameters estimated. The m2 test is normally distributed under the null hypothesis of the absence of second-order serial correlation between regressors and residuals.