This article evaluates the performance of alternative estimation methods for gravity models with heteroscedasticity 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 Nonlinear 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 heteroscedasticity than others are, its performance is similar, in terms of bias and SEs, to the FGLS estimator performance, in particular for small samples. GPML presents however the lowest bias and SEs 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. It is highly recommended to follow a model selection approach using a number of tests to select the more appropriate estimator for any application.