Replace empirical densities with analytical densities

[lnalborczyk]

I am currently simulating RT data using rwiener(n = 1e3, alpha = alpha, beta = beta, delta = delta, tau = tau). From there I compute the number of upper and lower boundaries hits and the densities using density() for each boundary. An alternative would be to directly plot the density for each boundary using dwiener(q, alpha, tau, beta, delta, resp = "upper") and dwiener(q, alpha, tau, beta, delta, resp = "lower"), which would be more robust to weird parameter values. However, this solution does not permit to compute the number of lower and upper hits (or does it?)...