timedelay: Time Delay Estimation for Stochastic Time Series of Gravitationally Lensed Quasars
Version 1.0.7

We provide a toolbox to estimate the time delay between the brightness time series of gravitationally lensed quasar images via Bayesian and profile likelihood approaches. The model is based on a state-space representation for irregularly observed time series data generated from a latent continuous-time Ornstein-Uhlenbeck process. Our Bayesian method adopts scientifically motivated hyper-prior distributions and a Metropolis-Hastings within Gibbs sampler, producing posterior samples of the model parameters that include the time delay. A profile likelihood of the time delay is a simple approximation to the marginal posterior distribution of the time delay. Both Bayesian and profile likelihood approaches complement each other, producing almost identical results; the Bayesian way is more principled but the profile likelihood is easier to implement.

Package details

AuthorHyungsuk Tak, Kaisey Mandel, David A. van Dyk, Vinay L. Kashyap, Xiao-Li Meng, and Aneta Siemiginowska
Date of publication2017-05-27 21:23:42 UTC
MaintainerHyungsuk Tak <[email protected]>
LicenseGPL-2
Version1.0.7
Package repositoryView on CRAN
Installation Install the latest version of this package by entering the following in R:
install.packages("timedelay")

Try the timedelay package in your browser

Any scripts or data that you put into this service are public.

timedelay documentation built on May 29, 2017, 9:37 a.m.