English

Goodness-of-fit test for noisy directional data

Statistics Theory 2013-11-18 v2 Statistics Theory

Abstract

We consider spherical data XiX_i noised by a random rotation εi\varepsilon_i\in SO(3) so that only the sample Zi=εiXiZ_i=\varepsilon_iX_i, i=1,,Ni=1,\dots, N is observed. We define a nonparametric test procedure to distinguish H0:H_0: ''the density ff of XiX_i is the uniform density f0f_0 on the sphere'' and H1:H_1: ''ff022\CψN\|f-f_0\|_2^2\geq \C\psi_N and ff is in a Sobolev space with smoothness ss''. For a noise density fεf_\varepsilon with smoothness index ν\nu, we show that an adaptive procedure (i.e. ss is not assumed to be known) cannot have a faster rate of separation than ψNad(s)=(N/loglog(N))2s/(2s+2ν+1)\psi_N^{ad}(s)=(N/\sqrt{\log\log(N)})^{-2s/(2s+2\nu+1)} and we provide a procedure which reaches this rate. We also deal with the case of super smooth noise. We illustrate the theory by implementing our test procedure for various kinds of noise on SO(3) and by comparing it to other procedures. Applications to real data in astrophysics and paleomagnetism are provided.

Keywords

Cite

@article{arxiv.1203.2008,
  title  = {Goodness-of-fit test for noisy directional data},
  author = {Claire Lacour and Thanh Mai Pham Ngoc},
  journal= {arXiv preprint arXiv:1203.2008},
  year   = {2013}
}