A method for the estimation of the significance of cross-correlations in unevenly sampled red-noise time series

We present a practical implementation of a Monte Carlo method to estimate the significance of cross-correlations in unevenly sampled time series of data, whose statistical properties are modelled with a simple power-law power spectral density. This implementation builds on published methods; we introduce a number of improvements in the normalization of the cross-correlation function estimate and a bootstrap method for estimating the significance of the cross-correlations. A closely related matter is the estimation of a model for the light curves, which is critical for the significance estimates. We present a graphical and quantitative demonstration that uses simulations to show how common it is to get high cross-correlations for unrelated light curves with steep power spectral densities. This demonstration highlights the dangers of interpreting them as signs of a physical connection. We show that by using interpolation and the Hanning sampling window function we are able to reduce the effects of red-noise leakage and to recover steep simple power-law power spectral densities. We also introduce the use of a Neyman construction for the estimation of the errors in the power-law index of the power spectral density. This method provides a consistent way to estimate the significance of cross-correlations in unevenly sampled time series of data.

[1]  J. Chiang,et al.  GAMMA-RAY LIGHT CURVES AND VARIABILITY OF BRIGHT FERMI-DETECTED BLAZARS , 2010, 1004.0348.

[2]  B. M. Peterson,et al.  Central Masses and Broad-Line Region Sizes of Active Galactic Nuclei. II. A Homogeneous Analysis of a Large Reverberation-Mapping Database , 2004, astro-ph/0407299.

[3]  竹安 数博,et al.  Time series analysis and its applications , 2007 .

[4]  M. F. Aller,et al.  The University of Michigan Radio Astronomy Data Base. I. Structure Function Analysis and the Relation between BL Lacertae Objects and Quasi-stellar Objects , 1992 .

[5]  J. Scargle.  Studies in astronomical time series analysis. II - Statistical aspects of spectral analysis of unevenly spaced data , 1982 .

[6]  M. Aller,et al.  Linear Polarization as a Probe of Gamma Ray Flaring Blazar Jets , 2009, 0912.3176.

[7]  A. P. Marscher,et al.  Multiwavelength monitoring of the BL Lacertae object PKS 2155-304. 4: Multiwavelength analysis , 1995 .

[8]  Dan Maoz,et al.  Astronomical Time Series , 1997 .

[9]  Y. Avni,et al.  Energy spectra of X-ray clusters of galaxies , 1976 .

[10]  Gino Tosti,et al.  Fermi Gamma-ray Space Telescope Observations of the Gamma-ray Outburst from 3C454.3 in November 2010 , 2011 .

[11]  D. Emmanoulopoulos,et al.  Generating artificial light curves: revisited and updated , 2013 .

[12]  Bradley M. Peterson,et al.  REVERBERATION MAPPING OF ACTIVE GALACTIC NUCLEI , 1993 .

[13]  W. Welsh,et al.  On the Reliability of Cross‐Correlation Function Lag Determinations in Active Galactic Nuclei , 1999, astro-ph/9911112.

[14]  T. Brown,et al.  Detection of Planetary Transits Across a Sun-like Star , 1999, The Astrophysical journal.

[15]  I. Papadakis,et al.  X-ray variability of active galactic nuclei : a universal power spectrum with luminosity-dependent amplitude , 1993 .

[16]  J. Scargle.  Studies in astronomical time series analysis. III - Fourier transforms, autocorrelation functions, and cross-correlation functions of unevenly spaced data , 1989 .

[17]  Bradley M. Peterson,et al.  On Uncertainties in Cross‐Correlation Lags and the Reality of Wavelength‐dependent Continuum Lags in Active Galactic Nuclei , 1998, astro-ph/9802103.

[18]  B. Peterson,et al.  The Accuracy of Cross-Correlation Estimates of Quasar Emission-Line Region Sizes , 1987 .

[19]  T. Deeming.  Fourier analysis with unequally-spaced data , 1975 .

[20]  I. M. McHardy,et al.  Measuring the broad-band power spectra of active galactic nuclei with RXTE , 2002 .

[21]  Ian M. McHardy,et al.  On the use of structure functions to study blazar variability: caveats and problems , 2010, 1001.2045.

[22]  Alan D. Martin,et al.  Review of Particle Physics (RPP) , 2012 .

[23]  William H. Press,et al.  Numerical Recipes in FORTRAN - The Art of Scientific Computing, 2nd Edition , 1987 .

[24]  R Edelson,et al.  The Discrete Correlation Function: a New Method for Analyzing Unevenly Sampled Variability Data , 1988 .

[25]  Technion,et al.  Correlated X-ray/optical variability in the quasar MR 2251−178 , 2008, 0807.2451.

[26]  Stefan J. Wagner,et al.  Intraday Variability in Quasars and BL LAC Objects , 1995 .

[27]  F. Harris.  On the use of windows for harmonic analysis with the discrete Fourier transform , 1978, Proceedings of the IEEE.

[28]  Jonathan D. Cryer,et al.  Time Series Analysis , 1986 .

[29]  H. R. Miller,et al.  Correlated Multi-Wave Band Variability in the Blazar 3C 279 from 1996 to 2007 , 2008, 0808.2194.

[30]  Bradley M. Peterson,et al.  Correlated Long-Term Optical and X-Ray Variations in NGC 5548 , 2003, astro-ph/0301216.

[31]  J. Bregman,et al.  Optical and radio variability in blazars , 1992 .

[32]  J. Chiang,et al.  THE LARGE AREA TELESCOPE ON THE FERMI GAMMA-RAY SPACE TELESCOPE MISSION , 2009, 0902.1089.

[33]  M. van der Klis,et al.  Fourier Techniques in X-Ray Timing , 1988 .

[34]  K. Grainge,et al.  BLAZARS IN THE FERMI ERA: THE OVRO 40 m TELESCOPE MONITORING PROGRAM , 2010, 1011.3111.