This paper stresses the bimodality of the likelihood function of the Mixed causal–noncausal AutoRegressions (MAR), and it is shown that the bimodality issue becomes more salient as the causal root approaches unity from below. The consequences are important as the roots of the local maxima are typically interchanged, attributing the noncausal component to the causal one and vice-versa. This severely changes the interpretation of the results, and the properties of unit root tests of the backward root are adversely affected. To circumvent the bimodality issue, this paper proposes an estimation strategy which (i) increases noticeably the probability of attaining the global MLE; and (ii) selects carefully the maximum used for the unit root test against a MAR stationary alternative.