RobPer : An R Package to Calculate Periodograms for Light Curves Based on Robust Regression

• Published in: Journal of statistical software Vol. 69; no. 9; pp. 1 - 37
• Main Authors: Thieler, Anita M., Fried, Roland, Rathjens, Jonathan
• Format: Journal Article
• Language: English
• Published: Foundation for Open Access Statistics 01-03-2016
• Subjects: astroparticle physics; beta distribution; Cramer-von-Mises distance minimization; irregular sampling; light curves; measurement accuracies; outlier detection; period detection; periodogram; regression model; robust regression; time series analysis; weighted regression
• ISSN: 1548-7660; 1548-7660
• DOI: 10.18637/jss.v069.i09

An important task in astroparticle physics is the detection of periodicities in irregularly sampled time series, called light curves. The classic Fourier periodogram cannot deal with irregular sampling and with the measurement accuracies that are typically given for each observation of a light curve. Hence, methods to fit periodic functions using weighted regression were developed in the past to calculate periodograms. We present the R package RobPer which allows to combine different periodic functions and regression techniques to calculate periodograms. Possible regression techniques are least squares, least absolute deviations, least trimmed squares, M-, S- and τ -regression. Measurement accuracies can be taken into account including weights. Our periodogram function covers most of the approaches that have been tried earlier and provides new model-regression-combinations that have not been used before. To detect valid periods, RobPer applies an outlier search on the periodogram instead of using fixed critical values that are theoretically only justified in case of least squares regression, independent periodogram bars and a null hypothesis allowing only normal white noise. Finally, the package also includes a generator to generate artificial light curves.