It is known that while estimating system of equations by using SUR and OLS estimator we obtain the same results when matrix of independent variables (X) is identical in each equation. I've tried this and results are ... different. Maybe you could figure out why?

I've attached files (with data and script). Could this be effect of cross equation residuals correlations or it's just ... an error?