For many decades, ultrahigh energy charged particles of unknown origin that can be observed from the ground have been a puzzle for particle physicists and astrophysicists. As an attempt to discriminate among several possible production scenarios, astrophysicists try to test the statistical isotropy of the directions of arrival of these cosmic rays. At the highest energies, they are supposed to point toward their sources with good accuracy. However, the observations are so rare that testing the distribution of such samples of directional data on the sphere is nontrivial. In this paper, we choose a nonparametric framework that makes weak hypotheses on the alternative distributions and allows in turn to detect various and possibly unexpected forms of anisotropy. We explore two particular procedures. Both are derived from fitting the empirical distribution with wavelet expansions of densities. We use the wavelet frame introduced by [SIAM J. Math. Anal. 38 (2006b) 574-594 (electronic)], the so-called needlets. The expansions are truncated at scale indices no larger than some ${J^{\star}}$, and the $L^p$ distances between those estimates and the null density are computed. One family of tests (called Multiple) is based on the idea of testing the distance from the null for each choice of $J=1,\ldots,{J^{\star}}$, whereas the so-called PlugIn approach is based on the single full ${J^{\star}}$ expansion, but with thresholded wavelet coefficients. We describe the practical implementation of these two procedures and compare them to other methods in the literature. As alternatives to isotropy, we consider both very simple toy models and more realistic nonisotropic models based on Physics-inspired simulations. The Monte Carlo study shows good performance of the Multiple test, even at moderate sample size, for a wide sample of alternative hypotheses and for different choices of the parameter ${J^{\star}}$. On the 69 most energetic events published by the Pierre Auger Collaboration, the needlet-based procedures suggest statistical evidence for anisotropy. Using several values for the parameters of the methods, our procedures yield $p$-values below 1%, but with uncontrolled multiplicity issues. The flexibility of this method and the possibility to modify it to take into account a large variety of extensions of the problem make it an interesting option for future investigation of the origin of ultrahigh energy cosmic rays.

In a recent paper, we analyzed the properties of a new kind of spherical wavelets (called needlets) for statistical inference procedures on spherical random fields; the investigation was mainly motivated by applications to cosmological data. In the present work, we exploit the asymptotic uncorrelation of random needlet coefficients at fixed angular distances to construct subsampling statistics evaluated on Voronoi cells on the sphere. We illustrate how such statistics can be used for isotropy tests and for bootstrap estimation of nuisance parameters, even when a single realization of the spherical random field is observed. The asymptotic theory is developed in detail in the high resolution sense.

We investigate invariant random fields on the sphere using a new type of spherical wavelets, called needlets. These are compactly supported in frequency and enjoy excellent localization properties in real space, with quasi-exponentially decaying tails. We show that, for random fields on the sphere, the needlet coefficients are asymptotically uncorrelated for any fixed angular distance. This property is used to derive CLT and functional CLT convergence results for polynomial functionals of the needlet coefficients: here the asymptotic theory is considered in the high-frequency sense. Our proposals emerge from strong empirical motivations, especially in connection with the analysis of cosmological data sets.

Scalar wavelets have been used extensively in the analysis of Cosmic Microwave Background (CMB) temperature maps. Spin needlets are a new form of (spin) wavelets which were introduced in the mathematical literature by Geller and Marinucci (2008) as a tool for the analysis of spin random fields. Here we adopt the spin needlet approach for the analysis of CMB polarization measurements. The outcome of experiments measuring the polarization of the CMB are maps of the Stokes Q and U parameters which are spin 2 quantities. Here we discuss how to transform these spin 2 maps into spin 2 needlet coefficients and outline briefly how these coefficients can be used in the analysis of CMB polarization data. We review the most important properties of spin needlets, such as localization in pixel and harmonic space and asymptotic uncorrelation. We discuss several statistical applications, including the relation of angular power spectra to the needlet coefficients, testing for non-Gaussianity on polarization data, and reconstruction of the E and B scalar maps.

This paper is concerned with density estimation of directional data on the sphere. We introduce a procedure based on thresholding on a new type of spherical wavelets called {\it needlets}. We establish a minimax result and prove its optimality. We are motivated by astrophysical applications, in particular in connection with the analysis of ultra high energy cosmic rays.

We discuss Spherical Needlets and their properties. Needlets are a form of spherical wavelets which do not rely on any kind of tangent plane approximation and enjoy good localization properties in both pixel and harmonic space; moreover needlets coefficients are asymptotically uncorrelated at any fixed angular distance, which makes their use in statistical procedures very promising. In view of these properties, we believe needlets may turn out to be especially useful in the analysis of Cosmic Microwave Background (CMB) data on the incomplete sky, as well as of other cosmological observations. As a final advantage, we stress that the implementation of needlets is computationally very convenient and may rely completely on standard data analysis packages such as HEALPix.