Package: ream
Title: Density, Distribution, and Sampling Functions for Evidence
        Accumulation Models
Version: 1.0-12
Authors@R: 
    c(person("Raphael", "Hartmann", email = "raphael.hartmann@protonmail.com", role = c("aut", "cre"), 
             comment = c(ORCID = "0000-0003-4686-9329")),
      person("Matthew", "Murrow", email = "matthew.a.murrow@vanderbilt.edu", role = c("aut")))
Description: Calculate the probability density functions (PDFs) for two threshold evidence 
    accumulation models (EAMs). These are defined using the following Stochastic 
    Differential Equation (SDE), dx(t) = v(x(t),t)*dt+D(x(t),t)*dW, where x(t) is
    the accumulated evidence at time t, v(x(t),t) is the drift rate, D(x(t),t) is 
    the noise scale, and W is the standard Wiener process. The boundary conditions 
    of this process are the upper and lower decision thresholds, represented by b_u(t) 
    and b_l(t), respectively. Upper threshold b_u(t) > 0, while lower threshold b_l(t) < 0. 
    The initial condition of this process x(0) = z where b_l(t) < z < b_u(t). We 
    represent this as the relative start point w = z/(b_u(0)-b_l(0)), defined as 
    a ratio of the initial threshold location. This package generates the PDF using
    the same approach as the 'python' package it is based upon, 'PyBEAM' by Murrow and Holmes
    (2023) <doi:10.3758/s13428-023-02162-w>. First, it converts the SDE model into the 
    forwards Fokker-Planck equation dp(x,t)/dt = d(v(x,t)*p(x,t))/dt-0.5*d^2(D(x,t)^2*p(x,t))/dx^2,
    then solves this equation using the Crank-Nicolson method to determine p(x,t). Finally, 
    it calculates the flux at the decision thresholds, f_i(t) = 0.5*d(D(x,t)^2*p(x,t))/dx 
    evaluated at x = b_i(t), where i is the relevant decision threshold, either upper (i = u) 
    or lower (i = l). The flux at each thresholds f_i(t) is the PDF for each threshold, 
    specifically its PDF. We discuss further details of this approach in this package 
    and 'PyBEAM' publications. Additionally, one can calculate the cumulative distribution 
    functions of and sampling from the EAMs.
Suggests: knitr, rmarkdown
NeedsCompilation: yes
SystemRequirements: C++17
VignetteBuilder: knitr
License: GPL (>= 2)
Encoding: UTF-8
RoxygenNote: 7.3.3
URL: https://github.com/RaphaelHartmann/ream
BugReports: https://github.com/RaphaelHartmann/ream/issues
Packaged: 2026-06-30 11:14:16 UTC; hartmann
Author: Raphael Hartmann [aut, cre] (ORCID:
    <https://orcid.org/0000-0003-4686-9329>),
  Matthew Murrow [aut]
Maintainer: Raphael Hartmann <raphael.hartmann@protonmail.com>
Repository: CRAN
Date/Publication: 2026-06-30 13:30:09 UTC
Built: R 4.5.2; x86_64-apple-darwin20; 2026-06-30 16:18:03 UTC; unix
Archs: ream.so.dSYM
