The measurement of large scale anisotropies in cosmic ray arrival directions
is generally performed through harmonic analyses of the right ascension
distribution as a function of energy. These measurements are challenging due to
the small expected anisotropies and meanwhile the relatively large modulations
of observed counting rates due to experimental effects. In this paper, we
present a procedure based on the shuffling technique to carry out these
measurements, applicable to any cosmic ray detector without any additional
corrections for the observed counting rates.Comment: 22 pages, 10 figures, to appear in Astroparticle Physic