Modeling Stochastic Variability in Multi-Band Time Series Data

Abstract

In preparation for the era of the time-domain astronomy with upcoming large-scale surveys, we propose a state-space representation of a multivariate damped random walk process as a tool to analyze irregularly-spaced multi-filter light curves with heteroscedastic measurement errors. We adopt a computationally efficient and scalable Kalman-filtering approach to evaluate the likelihood function, leading to maximum $O(k^{3}n)$ complexity, where $k$ is the number of available bands and $n$ is the number of unique observation times across the $k$ bands. This is a significant computational advantage over a commonly used univariate Gaussian process that can stack up all multi-band light curves in one vector with maximum $O(k^3{n}^3)$ complexity. Using such efficient likelihood computation, we provide both maximum likelihood estimates and Bayesian posterior samples of the model parameters. Three numerical illustrations are presented; (i) analyzing simulated five-band light curves for a comparison with independent single-band fits; (ii) analyzing five-band light curves of a quasar obtained from the Sloan Digital Sky Survey (SDSS) Stripe~82 to estimate the short-term variability and timescale; (iii) analyzing gravitationally lensed $g$- and $r$-band light curves of Q0957+561 to infer the time delay. Two R packages, Rdrw and timedelay, are publicly available to fit the proposed models.

References in corpus (34)

• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Are the Variations in Quasar Optical Flux Driven by Thermal Fluctuations?
• Modeling the Time Variability of SDSS Stripe 82 Quasars as a Damped Random Walk
• Quantifying Quasar Variability As Part of a General Approach To Classifying Continuously Varying Sources
• H0LiCOW I. $H_0$ Lenses in COSMOGRAIL's Wellspring: Program Overview
• An Alternative Approach To Measuring Reverberation Lags in Active Galactic Nuclei
• Is Quasar Optical Variability a Damped Random Walk?
• A Description of Quasar Variability Measured Using Repeated SDSS and POSS Imaging
• Space Telescope and Optical Reverberation Mapping Project. II. Swift and HST Reverberation Mapping of the Accretion Disk of NGC 5548
• Flexible and Scalable Methods for Quantifying Stochastic Variability in the Era of Massive Time-Domain Astronomical Data Sets
• Kepler Observations of Rapid Optical Variability in Active Galactic Nuclei
• Time Delay Cosmography
• Strong Lens Time Delay Challenge: II. Results of TDC1
• Are the Variability Properties of the Kepler AGN Light Curves Consistent with a Damped Random Walk?
• COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XI. Techniques for time delay measurement in presence of microlensing
• Limitations on the recovery of the true AGN variability parameters using Damped Random Walk modeling
• Application of Stochastic Modeling to Analysis of Photometric Reverberation Mapping Data
• A novel variability-based method for quasar selection: evidence for a rest frame ~54 day characteristic timescale
• AGN Variability Analysis Handbook
• Strong Lens Time Delay Challenge: I. Experimental Design
• Assessment of stochastic and deterministic models of 6304 quasar lightcurves from SDSS Stripe 82
• Disentangling Time Series Spectra with Gaussian Processes: Applications to Radial Velocity Analysis
• A Degeneracy in DRW Modelling of AGN Light Curves
• Robust Strong Lensing Time Delay Estimation
• A Repelling-Attracting Metropolis Algorithm for Multimodality
• A Refined QSO Selection Method Using Diagnostics Tests: 663 QSO Candidates in the LMC
• Fast Bayesian inference of the multivariate Ornstein-Uhlenbeck process
• Multivariate CARMA processes, continuous-time state space models and complete regularity of the innovations of the sampled processes
• How proper are Bayesian models in the astronomical literature?
• Bayesian Estimates of Astronomical Time Delays between Gravitationally Lensed Stochastic Light Curves
• Robust and Accurate Inference via a Mixture of Gaussian and Student's t Errors
• A 5.5-year robotic optical monitoring of Q0957+561: substructure in a non-local cD galaxy
• Rgbp: An R Package for Gaussian, Poisson, and Binomial Random Effects Models with Frequency Coverage Evaluations
• Data transforming augmentation for heteroscedastic models

Cited by in corpus (18)

• Optical Variability of Quasars with 20-Year Photometric Light Curves
• Examining AGN UV/optical Variability Beyond the Simple Damped Random Walk
• Toward Extremely Precise Radial Velocities: II. A Tool For Using Multivariate Gaussian Processes to Model Stellar Activity
• Quasars with Periodic Variability: Capabilities and Limitations of Bayesian Searches for Supermassive Black Hole Binaries in Time-Domain Surveys
• The LSST era of supermassive black holes accretion-disk reverberation mapping
• Non-detection of ${}^6$Li in Spite plateau stars with ESPRESSO
• LSST Cadence Strategy Evaluations for AGN Time-series Data in Wide-Fast-Deep Field
• Incorporating Measurement Error in Astronomical Object Classification
• How can the optical variation properties of active galactic nuclei be unbiasedly measured?
• Detecting the periodicity of highly irregularly sampled light-curves with Gaussian processes: the case of SDSSJ025214.67-002813.7
• Latent Stochastic Differential Equations for Modeling Quasar Variability and Inferring Black Hole Properties
• Novel bivariate autoregressive model for predicting and forecasting irregularly observed time series
• TD-CARMA: Painless, accurate, and scalable estimates of gravitational-lens time delays with flexible CARMA processes
• Tracking hidden magnetospheric fluctuations in accretion-powered pulsars with a Kalman filter
• GausSN: Bayesian Time-Delay Estimation for Strongly Lensed Supernovae
• Six Maxims of Statistical Acumen for Astronomical Data Analysis
• Scalable Gaussian Processes for Integrated and Overlapping Measurements Via Augmented State Space Models
• A Robust Bayesian Meta-Analysis for Estimating the Hubble Constant via Time Delay Cosmography