Finally, we estimate the model in its original multiplicative form to tackle the problems of zero trade and heteroskedaticity in the error term. As mentioned in Section 3, we follow the most recent developments in the gravity model literature and estimate a multinomial PML as described in Head and Mayer (forthcoming) and as explained in Section 3 above.