simple.band1

A 50 by 5 matrix containing a data set of a doubly-lensed quasar observed in band 1, which is modified from realistic data.

datasets

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. A new functionality is added in version 1.0.9 for estimating the time delay between doubly-lensed light curves observed in two bands. See also Tak et al. (2017) <doi:10.1214/17-AOAS1027>, Tak et al. (2018) <doi:10.1080/10618600.2017.1415911>, Hu and Tak (2020) <arXiv:2005.08049>.

Hyungsuk Tak

timedelay

Time Delay Estimation for Stochastic Time Series of
Gravitationally Lensed Quasars

simple.band1 function

A simple data set of a doubly-lensed quasar:<dl class="dl-horizontal">
 <dt><code>time</code></dt><dd>observation times</dd>
 <dt><code>x</code></dt><dd>magnitudes of light A observed in band 1</dd>
 <dt><code>se.x</code></dt><dd>measurement errors of x</dd>
 <dt><code>y</code></dt><dd>magnitudes of light B observed in band 1</dd>
 <dt><code>se.y</code></dt><dd>measurement errors of y</dd>
 </dl>

Format

Simulated simple data of a doubly-lensed quasar observed in band 1 — simple.band1

A simple data set of a doubly-lensed quasar:<dl class='dl-horizontal'>
 <dt><code>time</code></dt><dd>observation times</dd>
 <dt><code>x</code></dt><dd>magnitudes of light A observed in band 1</dd>
 <dt><code>se.x</code></dt><dd>measurement errors of x</dd>
 <dt><code>y</code></dt><dd>magnitudes of light B observed in band 1</dd>
 <dt><code>se.y</code></dt><dd>measurement errors of y</dd>
 </dl>

simple.band1: Simulated simple data of a doubly-lensed quasar observed in band 1

Description

Usage

Arguments

Format

Examples