Estimation methods of bivariate fractional cointegration models are numerous. In most cases they have non-equivalent asymptotic and finite sample properties, implying diffculties in determining an optimal estimation strategy. In this paper, we address this issue by means of simulations and provide useful guidance to practitioners. Our Monte Carlo study reveals the superiority of techniques that estimate jointly all parameters of interest, over those operating in two steps. In some cases, it also shows that estimators originally designed for the stationary cointegration, have good finite sample properties in non-stationary regions of the parameter space.
