An adaptive rejection sampler for sampling from the Wiener diffusion model

The Wiener diffusion model with two absorbing boundaries is one of the most frequently applied models for jointly modeling responses and response latencies in psychological research. We consider four methods for sampling from the model with and without variability in drift rate, starting point, and non-decision time: Inverse transform sampling, rejection sampling, and two new methods based on adaptive rejection sampling (ARS). We implement these four methods in an R package, validate the methods, and compare their sampling speed in different settings. All four implemented methods provide samples that follow the intended distributions. The ARS-based methods, however, outperform the other methods in sampling speed as the requested sample size increases. We provide guidelines for when using ARS is more efficient than using traditional methods and vice versa.