Package: timedelay 1.0.11

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

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>.

Authors:Hyungsuk Tak, Kaisey Mandel, David A. van Dyk, Vinay L. Kashyap, Xiao-Li Meng, Aneta Siemiginowska, and Zhirui Hu

timedelay_1.0.11.tar.gz
timedelay_1.0.11.zip(r-4.5)timedelay_1.0.11.zip(r-4.4)timedelay_1.0.11.zip(r-4.3)
timedelay_1.0.11.tgz(r-4.4-any)timedelay_1.0.11.tgz(r-4.3-any)
timedelay_1.0.11.tar.gz(r-4.5-noble)timedelay_1.0.11.tar.gz(r-4.4-noble)
timedelay_1.0.11.tgz(r-4.4-emscripten)timedelay_1.0.11.tgz(r-4.3-emscripten)
timedelay.pdf |timedelay.html
timedelay/json (API)

# Install 'timedelay' in R:
install.packages('timedelay', repos = c('https://hyungsuktak.r-universe.dev', 'https://cloud.r-project.org'))

Peer review:

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

On CRAN:

This package does not link to any Github/Gitlab/R-forge repository. No issue tracker or development information is available.

1.00 score 4 scripts 211 downloads 12 exports 2 dependencies

Last updated 4 years agofrom:9a87917b58. Checks:OK: 7. Indexed: yes.

TargetResultDate
Doc / VignettesOKOct 15 2024
R-4.5-winOKOct 15 2024
R-4.5-linuxOKOct 15 2024
R-4.4-winOKOct 15 2024
R-4.4-macOKOct 15 2024
R-4.3-winOKOct 15 2024
R-4.3-macOKOct 15 2024

Exports:bayesianbayesian.multibandentirelogprofilelikelihoodLog.post.deltaLog.post.thetalogpostDeltaorderDatapostBetapostThetapostTheta.multipostXpostX.multi

Dependencies:MASSmvtnorm