The log of gravity revisited

Abstract This paper evaluates the performance of alternative estimation methods for gravity models with heteroskedasticity and zero trade values. Both problematic issues, recently addressed by Santos Silva and Tenreyro in an influential paper, are re-examined here. We use Monte Carlo simulations to compare the Pseudo Poisson Maximum Likelihood (PPML) estimator recommended by Santos Silva and Tenreyro, a Gamma pseudo-maximum-likelihood (GPML), a Non-Linear Least Squares (NLS) estimator and a Feasible Generalized Least Squares (FGLS) estimator with more traditional techniques. Additionally, estimates of the gravity equation are obtained for three different data sets with the abovementioned methods. The results of the simulation study indicate that, although the PPML estimator is less affected by heteroskedasticity than others are, its performance is similar, in terms of bias and standard errors, to the FGLS estimator performance. GPML presents however the lowest bias and standard errors in the simulations without zero values. The results of the empirical estimations, using three different samples containing real data, indicate that the choice of estimator has to be made for each specific dataset. There is not a general? best? estimator and it is highly recommended to follow a model selection approach using a number of tests to select the more appropriate estimator for any application.