CRAN
timedelay 1.0.8
Time Delay Estimation for Stochastic Time Series of Gravitationally Lensed Quasars
Released May 16, 2018 by Hyungsuk Tak
Dependencies
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.
Installation
Maven
This package can be included as a dependency from a Java or Scala project by including
the following your project's pom.xml file.
Read more
about embedding Renjin in JVM-based projects.
<dependencies>
<dependency>
<groupId>org.renjin.cran</groupId>
<artifactId>timedelay</artifactId>
<version>1.0.8-b1</version>
</dependency>
</dependencies>
<repositories>
<repository>
<id>bedatadriven</id>
<name>bedatadriven public repo</name>
<url>https://nexus.bedatadriven.com/content/groups/public/</url>
</repository>
</repositories>
Renjin CLI
If you're using Renjin from the command line, you load this library by invoking:
library('org.renjin.cran:timedelay')
Test Results
This package was last tested against Renjin 0.9.2644 on Jun 2, 2018.