Looking for a Stata Expert to optimize highly advanced postgraduate International Economics code

Hello, we are estimating the effect of currency unions on bilateral trade using the gravity equation and both OLS and PPML, but we get too high estimates compared to the literature. We use the traditional CEPII gravity dataset and control variables. We checked for collinearity, outliers, etc. and all seems to be in order. We only need to find what is wrong with my code - or to adjust it properly, so that We get much "normal" estimates. This is why we are looking for someone familiar with International Economics, the gravity equation, OLS, and PPML. Please note that our code is already finished. Sincerely.