One of the most unusual aspects of the X-ray binary LSI +61 303 is that at each orbit (P1=26.4960 \pm 0.0028 d) one radio outburst occurs whose amplitude is modulated with Plong, a long-term period of more than 4 yr. It is still not clear whether the compact object of the system or the companion Be star is responsible for the long-term modulation. We study here the stability of Plong. Such a stability is expected if Plong is due to periodic (P2) Doppler boosting of periodic (P1) ejections from the accreting compact object of the system. On the contrary it is not expected if Plong is related to variations in the mass loss of the companion Be star. We built a database of 36.8 yr of radio observations of LSI +61 303 covering more than 8 long-term cycles. We performed timing and correlation analysis. In addition to the two dominant features at P1 and P2, the timing analysis gives a feature at Plong=1628 \pm 48 days. The determined value of Plong agrees with the beat of the two dominant features, i.e. Pbeat=1/(\nu1 -\nu2)=1626 \pm 68 d. The correlation coefficient of the radio data oscillates at multiples of Pbeat. Cycles in varying Be stars change in length and disappear after 2-3 cycles following the well-studied case of the binary system zeta Tau. On the contrary, in LSI +61 303 the long-term period is quite stable and repeats itself over the available 8 cycles. The long-term modulation in LSI +61 303 accurately reflects the beat of periodical Doppler boosting (induced by precession) with the periodicity of the ejecta. The peak of the long-term modulation occurs at the coincidence of the maximum number of ejected particles with the maximum Doppler boosting of their emission; this coincidence occurs every 1/(\nu1 - \nu2) and creates the long-term modulation observed in LSI +61 303.