21.3 Wiener First Passage Time Distribution
21.3.1 Probability density function
If \(\alpha \in \mathbb{R}^+\), \(\tau \in \mathbb{R}^+\), \(\beta \in [0, 1]\) and \(\delta \in \mathbb{R}\), then for \(y > \tau\), \[ \text{Wiener}(y|\alpha, \tau, \beta, \delta) = \frac{\alpha^3}{(y-\tau)^{3/2}} \exp \! \left(- \delta \alpha \beta - \frac{\delta^2(y-\tau)}{2}\right) \sum_{k = - \infty}^{\infty} (2k + \beta) \phi \! \left(\frac{2k \alpha + \beta}{\sqrt{y - \tau}}\right) \] where \(\phi(x)\) denotes the standard normal density function; see (Feller 1968), (Navarro and Fuss 2009).
21.3.2 Sampling statement
y ~ wiener(alpha, tau, beta, delta)
Increment target log probability density with wiener_lupdf(y | alpha, tau, beta, delta).
Available since 2.7
21.3.3 Stan functions
real wiener_lpdf(reals y | reals alpha, reals tau, reals beta, reals delta)
The log of the Wiener first passage time density of y given boundary
separation alpha, non-decision time tau, a-priori bias beta and drift
rate delta
Available since 2.18
real wiener_lupdf(reals y | reals alpha, reals tau, reals beta, reals delta)
The log of the Wiener first passage time density of y given boundary
separation alpha, non-decision time tau, a-priori bias beta and drift
rate delta dropping constant additive terms
Available since 2.25
21.3.4 boundaries
Stan returns the first passage time of the accumulation process over the upper boundary only. To get the result for the lower boundary, use \[ \text{wiener}(y | \alpha, \tau, 1 - \beta, - \delta) \] For more details, see the appendix of Vandekerckhove and Wabersich (2014).