Issue |
A&A
Volume 615, July 2018
|
|
---|---|---|
Article Number | A111 | |
Number of page(s) | 12 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/201732524 | |
Published online | 23 July 2018 |
Estimating activity cycles with probabilistic methods
I. Bayesian generalised Lomb-Scargle periodogram with trend
1
ReSoLVE Centre of Excellence, Aalto University, Department of Computer Science,
PO Box 15400,
00076
Aalto, Finland
e-mail: nigul.olspert@aalto.fi
2
Tartu Observatory,
61602
Tõravere, Estonia
3
Max-Planck-Institut für Sonnensystemforschung,
Justus-von-Liebig-Weg 3,
37077
Göttingen, Germany
Received:
21
December
2017
Accepted:
27
March
2018
Context. Period estimation is one of the central topics in astronomical time series analysis, in which data is often unevenly sampled. Studies of stellar magnetic cycles are especially challenging, as the periods expected in those cases are approximately the same length as the datasets themselves. The datasets often contain trends, the origin of which is either a real long-term cycle or an instrumental effect. But these effects cannot be reliably separated, while they can lead to erroneous period determinations if not properly handled.
Aims. In this study we aim at developing a method that can handle the trends properly. By performing an extensive set of testing, we show that this is the optimal procedure when contrasted with methods that do not include the trend directly in the model. The effect of the form of the noise (whether constant or heteroscedastic) on the results is also investigated.
Methods. We introduced a Bayesian generalised Lomb-Scargle periodogram with trend (BGLST), which is a probabilistic linear regression model using Gaussian priors for the coefficients of the fit and a uniform prior for the frequency parameter.
Results. We show, using synthetic data, that when there is no prior information on whether and to what extent the true model of the data contains a linear trend, the introduced BGLST method is preferable to the methods that either detrend the data or opt not to detrend the data before fitting the periodic model. Whether to use noise with other than constant variance in the model depends on the density of the data sampling and on the true noise type of the process.
Key words: methods: statistical / methods: numerical / stars: activity
© ESO 2018
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.