Estimating gravity equation models in the presence of sample selection and heteroscedasticity