logLik_dt {CircularDDM}

Log-Likelihood for Circular First Passage Time

Description

Calculate circular log-likelihood of the first passage time, using equation (22) on p 432.

Usage

logLik_dt(x, pVec, k = 141L)

Arguments

x

a matrix storing a first column as RT and a second column of continuous responses/reports/outcomes. Each row is a trial.

pVec

a parameter vector with the order [a, vx, vy, t0, s], a stands for response threshold, vx is the drift rate along x axis, vy is the drift rate along y axis, t0 is the non-decision time, and s is the within-trial standard deviation.

k

a precision for bessel function. The larger the k is, the larger the memory space is required. Default is 141.

Value

a vector

References

Smith, P. L. (2016). Diffusion Theory of Decision Making in Continuous Report, Psychological Review, 123 (4), 425–451.

Examples

x <- cbind(
RT=c(1.2595272, 0.8693937, 0.8009044, 1.0018933, 2.3640007, 1.0521304),
R =c(1.9217430, 1.7844653, 0.2662521, 2.1569724, 1.7277440, 0.8607271)
)
pVec <- c(a=2.45, vx=1.5, vy=1.25, t0=.1, s=1)
den  <- logLik_dt(x, pVec=pVec);
den