Movatterモバイル変換

DMCfun: Diffusion Model of Conflict (DMC) in Reaction Time Tasks

Description

DMC model simulation detailed in Ulrich, R., Schroeter, H., Leuthold, H., & Birngruber, T. (2015). Automatic and controlled stimulus processing in conflict tasks: Superimposed diffusion processes and delta functions. Cognitive Psychology, 78, 148-174. Ulrich et al. (2015)doi:10.1016/j.cogpsych.2015.02.005. Decision processes within choice reaction-time (CRT) tasks are often modelled using evidence accumulation models (EAMs), a variation of which is the Diffusion Decision Model (DDM, for a review, see Ratcliff & McKoon, 2008). Ulrich et al. (2015) introduced a Diffusion Model for Conflict tasks (DMC). The DMC model combines common features from within standard diffusion models with the addition of superimposed controlled and automatic activation. The DMC model is used to explain distributional reaction time (and error rate) patterns in common behavioural conflict-like tasks (e.g., Flanker task, Simon task). This R-package implements the DMC model and provides functionality to fit the model to observed data. Further details are provided in the following paper: Mackenzie, I.G., & Dudschig, C. (2021). DMCfun: An R package for fitting Diffusion Model of Conflict (DMC) to reaction time and error rate data. Methods in Psychology, 100074.doi:10.1016/j.metip.2021.100074.

Author(s)

Maintainer: Ian G. Mackenzieian.mackenzie@uni-tuebingen.de

Authors:

• Carolin Dudschigcarolin.dudschig@uni-tuebingen.de

See Also

Useful links:

• https://github.com/igmmgi/DMCfun

• https://CRAN.R-project.org/package=DMCfun

• https://www.sciencedirect.com/science/article/pii/S259026012100031X

addDataDF

Description

Add simulated ex-gaussian reaction-time (RT) data andbinary error (Error = 1, Correct = 0) data to an R DataFrame. This functioncan be used to create simulated data sets.

Usage

addDataDF(dat, RT = NULL, Error = NULL)

Arguments

Value

DataFrame with RT (ms) and Error (bool) columns

Examples

# Example 1: default dataframedat <- createDF()dat <- addDataDF(dat)head(dat)hist(dat$RT, 100)table(dat$Error)# Example 2: defined overall RT parametersdat <- createDF(nSubjects = 50, nTrl = 50, design = list("Comp" = c("comp", "incomp")))dat <- addDataDF(dat, RT = c(500, 150, 100))boxplot(dat$RT ~ dat$Comp)table(dat$Comp, dat$Error)# Example 3: defined RT + Error parameters across conditionsdat <- createDF(nSubjects = 50, nTrl = 50, design = list("Comp" = c("comp", "incomp")))dat <- addDataDF(dat,                 RT = list("Comp_comp"   = c(500, 80, 100),                           "Comp_incomp" = c(600, 80, 140)),                 Error = list("Comp_comp"   = 5,                              "Comp_incomp" = 15))boxplot(dat$RT ~ dat$Comp)table(dat$Comp, dat$Error)# Example 4:# create dataframe with defined RT + Error parameters across different conditionsdat <- createDF(nSubjects = 50, nTrl = 50, design = list("Comp" = c("comp", "incomp", "neutral")))dat <- addDataDF(dat,                 RT = list("Comp_comp"      = c(500, 150, 100),                           "Comp_neutral"   = c(550, 150, 100),                           "Comp_incomp"    = c(600, 150, 100)),                 Error = list("Comp_comp"    =  5,                              "Comp_neutral" = 10,                              "Comp_incomp"  = 15))boxplot(dat$RT ~ dat$Comp)table(dat$Comp, dat$Error)# Example 5:# create dataframe with defined RT + Error parameters across different conditionsdat <- createDF(nSubjects = 50, nTrl = 50,                design = list("Hand" = c("left", "right"),                              "Side" = c("left", "right")))dat <- addDataDF(dat,                 RT = list("Hand:Side_left:left"   = c(400, 150, 100),                           "Hand:Side_left:right"  = c(500, 150, 100),                           "Hand:Side_right:left"  = c(500, 150, 100),                           "Hand:Side_right:right" = c(400, 150, 100)),                 Error = list("Hand:Side_left:left"   = c(5,4,2,2,1),                              "Hand:Side_left:right"  = c(15,4,2,2,1),                              "Hand:Side_right:left"  = c(15,7,4,2,1),                              "Hand:Side_right:right" = c(5,8,5,3,1)))boxplot(dat$RT ~ dat$Hand + dat$Side)table(dat$Error, dat$Hand, dat$Side)

addErrorBars: Add errorbars to plot.

Description

Add error bars to current plot (uses base arrows function).

Usage

addErrorBars(xpos, ypos, errorSize, arrowSize = 0.1)

Arguments

Value

Plot (no return value)

Examples

# Example 1plot(c(1, 2), c(450, 500), xlim = c(0.5, 2.5), ylim = c(400, 600), type = "o")addErrorBars(c(1, 2), c(450, 500), errorSize = c(20, 20))# Example 2plot(c(1, 2), c(450, 500), xlim = c(0.5, 2.5), ylim = c(400, 600), type = "o")addErrorBars(c(1, 2), c(450, 500), errorSize = c(20, 40), arrowSize = 0.1)

calculateBinProbabilities

Description

Calculate bin probabilities in observed data

Usage

calculateBinProbabilities(resOb, quantileType = 5)

Arguments

Value

resOb Observed data with additional $probSubject/$prob table

Examples

# Example 1:resOb <- flankerDataresOb <- calculateBinProbabilities(resOb)resOb$prob

calculateCAF

Description

Calculate conditional accuracy function (CAF).The DataFrame should contain columns defining the participant, compatibility condition,RT and error (Default column names: "Subject", "Comp", "RT", "Error"). The "Comp" column shoulddefine compatibility condition (Default: c("comp", "incomp")) and the "Error" column shoulddefine if the trial was an error or not (Default: c(0, 1) ).

Usage

calculateCAF(  dat,  nCAF = 5,  columns = c("Subject", "Comp", "RT", "Error"),  compCoding = c("comp", "incomp"),  errorCoding = c(0, 1))

Arguments

Value

calculateCAF returns a DataFrame with conditional accuracy function (CAF) data (Bin, comp, incomp, effect)

Examples

# Example 1dat <- createDF(nSubjects = 1, nTrl = 10000, design = list("Comp" = c("comp", "incomp")))dat <- addDataDF(dat,                 RT = list("Comp_comp"   = c(500, 80, 100),                           "Comp_incomp" = c(600, 80, 140)),                 Error = list("Comp_comp"   = c( 5, 4, 3, 2, 1),                              "Comp_incomp" = c(20, 8, 6, 4, 2)))caf <- calculateCAF(dat)# Example 2dat <- createDF(nSubjects = 1, nTrl = 10000, design = list("Congruency" = c("cong", "incong")))dat <- addDataDF(dat,                 RT = list("Congruency_cong"   = c(500, 80, 100),                           "Congruency_incong" = c(600, 80, 140)),                 Error = list("Congruency_cong"   = c( 5, 4, 3, 2, 1),                              "Congruency_incong" = c(20, 8, 6, 4, 2)))head(dat)caf <- calculateCAF(dat, columns = c("Subject", "Congruency", "RT", "Error"),                    compCoding = c("cong", "incong"))

calculateCostValueCS

Description

Calculate cost value (fit) using chi-square (CS) from correct and incorrect RT data.

Usage

calculateCostValueCS(resTh, resOb)

Arguments

Value

cost value (CS)

Examples

# Example 1:resTh <- dmcSim()resOb <- flankerDataresOb <- calculateBinProbabilities(resOb)cost  <- calculateCostValueCS(resTh, resOb)

calculateCostValueGS

Description

Calculate cost value (fit) using likelihood-ratio chi-square statistic (GS) from correct and incorrect RT data.

Usage

calculateCostValueGS(resTh, resOb)

Arguments

Value

cost value (GS)

Examples

# Example 1:resTh <- dmcSim()resOb <- flankerDataresOb <- calculateBinProbabilities(resOb)cost  <- calculateCostValueGS(resTh, resOb)

calculateCostValueRMSE

Description

Calculate cost value (fit) using root-mean-square error (RMSE) from a combination of RT and error rate.

Usage

calculateCostValueRMSE(resTh, resOb)

Arguments

Value

cost value (RMSE)

Examples

# Example 1:resTh <- dmcSim()resOb <- dmcSim()cost <- calculateCostValueRMSE(resTh, resOb)# Example 2:resTh <- dmcSim()resOb <- dmcSim(tau = 150)cost <- calculateCostValueRMSE(resTh, resOb)

calculateCostValueSPE

Description

Calculate cost value (fit) using squared percentage errror (SPE) from combination of RT and error rate.

Usage

calculateCostValueSPE(resTh, resOb)

Arguments

Value

cost value (SPE)

Examples

# Example 1:resTh <- dmcSim()resOb <- dmcSim()cost <- calculateCostValueSPE(resTh, resOb)# Example 2:resTh <- dmcSim()resOb <- dmcSim(tau = 150)cost <- calculateCostValueSPE(resTh, resOb)

calculateDelta

Description

Calculate delta plot. Here RTs are split into n bins (Default: 5) for compatible andincompatible trials separately. Mean RT is calculated for each condition in each bin thensubtracted (incompatible - compatible) to give a compatibility effect (delta) at each bin.

Usage

calculateDelta(  dat,  nDelta = 19,  tDelta = 1,  columns = c("Subject", "Comp", "RT"),  compCoding = c("comp", "incomp"),  quantileType = 5)

Arguments

Value

calculateDelta returns a DataFrame with distributional delta analysis data (Bin, comp, incomp, meanBin, Effect)

Examples

# Example 1dat <- createDF(nSubjects = 1, nTrl = 10000, design = list("Comp" = c("comp", "incomp")))dat <- addDataDF(dat,                 RT = list("Comp_comp"   = c(500, 80, 100),                           "Comp_incomp" = c(600, 80, 140)))delta <- calculateDelta(dat)# Example 2dat <- createDF(nSubject = 1, nTrl = 10000, design = list("Congruency" = c("cong", "incong")))dat <- addDataDF(dat,                 RT = list("Congruency_cong"   = c(500, 80, 100),                           "Congruency_incong" = c(600, 80, 140)))head(dat)delta <- calculateDelta(dat, nDelta = 9, columns = c("Subject", "Congruency", "RT"),                        compCoding = c("cong", "incong"))

createDF

Description

Create dataframe (see also addDataDF)

Usage

createDF(  nSubjects = 20,  nTrl = 50,  design = list(A = c("A1", "A2"), B = c("B1", "B2")))

Arguments

Value

DataFrame with Subject, Factor(s) columns

Examples

# Example 1dat <- createDF()# Example 2dat <- createDF(nSubjects = 50, nTrl = 50, design = list("Comp" = c("comp", "incomp")))# Example 3dat <- createDF(nSubjects = 50, nTrl = 50, design = list("Comp" = c("comp", "incomp"),                                                         "Side" = c("left", "right", "middle")))

dmcCombineObservedData

Description

Combine observed datasets

Usage

dmcCombineObservedData(...)

Arguments

Value

dmcCombineObservedData returns a list of objects of class "dmcob"

Examples

# Example 1dat <- dmcCombineObservedData(flankerData, simonData)  # combine flanker/simon dataplot(dat, figType = "delta", xlimDelta = c(200, 700), ylimDelta = c(-20, 80),     cols = c("black", "darkgrey"), legend.parameters = list(x=200, y=80,     legend = c("Flanker Task", "Simon Task")))

dmcCppR

Description

dmcCppR

dmcFit

Description

Fit theoretical data generated from dmcSim to observed data byminimizing the root-mean-square error ("RMSE") between a weighted combinationof the CAF and CDF functions using optim (Nelder-Mead). Alternative cost functionsinclude squared percentage error ("SPE"), and g-squared statistic ("GS").

Usage

dmcFit(  resOb,  nTrl = 1e+05,  startVals = list(),  minVals = list(),  maxVals = list(),  fixedFit = list(),  freeCombined = list(),  fitInitialGrid = TRUE,  fitInitialGridN = 10,  fixedGrid = list(),  nCAF = 5,  nDelta = 19,  pDelta = vector(),  tDelta = 1,  deltaErrors = FALSE,  spDist = 1,  drOnset = 0,  drDist = 0,  drShape = 3,  drLim = c(0.1, 0.7),  rtMax = 5000,  costFunction = "RMSE",  printInputArgs = TRUE,  printResults = FALSE,  optimControl = list(),  numCores = 2)

Arguments

Value

dmcfit returns an object of class "dmcfit" with the following components:

Examples

# Code below can exceed CRAN check time limit, hence donttest# Example 1: Flanker data from Ulrich et al. (2015)fit <- dmcFit(flankerData) # only initial search tauplot(fit, flankerData)summary(fit)# Example 2: Simon data from Ulrich et al. (2015)fit <- dmcFit(simonData) # only initial search tauplot(fit, simonData)summary(fit)# Example 3: Flanker data from Ulrich et al. (2015) with non-default# start vals and some fixed valuesfit <- dmcFit(flankerData,  startVals = list(drc = 0.6, aaShape = 2.5),  fixedFit = list(drc = TRUE, aaShape = TRUE))# Example 4: Simulated Data (+ve going delta function)dat <- createDF(nSubjects = 20, nTrl = 500, design = list("Comp" = c("comp", "incomp")))dat <- addDataDF(dat,  RT = list(    "Comp_comp" = c(510, 100, 100),    "Comp_incomp" = c(540, 130, 85)  ),  Error = list(    "Comp_comp" = c(4, 3, 2, 1, 1),    "Comp_incomp" = c(20, 4, 3, 1, 1)  ))datOb <- dmcObservedData(dat, columns = c("Subject", "Comp", "RT", "Error"))plot(datOb)fit <- dmcFit(datOb, nTrl = 5000)plot(fit, datOb)summary(fit)# Example 5: Fitting 2+ datasets within all common parameters valuesfit <- dmcFit(list(flankerData, simonData), nTrl=1000)plot(fit[[1]], flankerData)plot(fit[[2]], simonData)summary(fit)# Example 6: Fitting 2+ datasets within some parameters values varyingfit <- dmcFit(list(flankerData, simonData), freeCombined=list(amp=TRUE, tau=TRUE), nTrl=1000)summary(fit) # NB. amp/tau values different, other parameter values equal

dmcFitDE

Description

Fit theoretical data generated from dmcSim to observed data byminimizing the root-mean-square error (RMSE) between a weighted combinationof the CAF and CDF functions using the R-package DEoptim. Alternative cost functionsinclude squared percentage error ("SPE"), and g-squared statistic ("GS").

Usage

dmcFitDE(  resOb,  nTrl = 1e+05,  minVals = list(),  maxVals = list(),  fixedFit = list(),  freeCombined = list(),  nCAF = 5,  nDelta = 19,  pDelta = vector(),  tDelta = 1,  deltaErrors = FALSE,  spDist = 1,  drOnset = 0,  drDist = 0,  drShape = 3,  drLim = c(0.1, 0.7),  rtMax = 5000,  costFunction = "RMSE",  deControl = list(),  numCores = 2)

Arguments

Value

dmcfit returns an object of class "dmcfit" with the following components:

Examples

# The code below can exceed CRAN check time limit, hence donttest# NB. The following code when using numCores = 2 (default) takes approx 20 minutes on# a standard desktop, whilst when increasing the number of cores used, (numCores = 12),# the code takes approx 5 minutes.# Example 1: Flanker data from Ulrich et al. (2015)fit <- dmcFitDE(flankerData, nTrl = 1000);plot(fit, flankerData)summary(fit)# Example 2: Simon data from Ulrich et al. (2015)fit <- dmcFitDE(simonData, nTrl = 5000, deControl = list(itermax=30))plot(fit, simonData)summary(fit)

dmcFitSubject

Description

Fit theoretical data generated from dmcSim to observed data byminimizing the root-mean-square error ("RMSE") between a weighted combinationof the CAF and CDF functions using optim (Nelder-Mead). Alternative cost functionsinclude squared percentage error ("SPE"), and g-squared statistic ("GS").

Usage

dmcFitSubject(  resOb,  nTrl = 1e+05,  startVals = list(),  minVals = list(),  maxVals = list(),  fixedFit = list(),  fitInitialGrid = TRUE,  fitInitialGridN = 10,  fixedGrid = list(),  freeCombined = list(),  nCAF = 5,  nDelta = 19,  pDelta = vector(),  tDelta = 1,  deltaErrors = FALSE,  spDist = 1,  drOnset = 0,  drDist = 0,  drShape = 3,  drLim = c(0.1, 0.7),  rtMax = 5000,  costFunction = "RMSE",  subjects = c(),  printInputArgs = TRUE,  printResults = FALSE,  optimControl = list(),  numCores = 2)

Arguments

Value

dmcFitSubject returns a list of objects of class "dmcfit"

Examples

# Code below can exceed CRAN check time limit, hence donttest# Example 1: Flanker data from Ulrich et al. (2015)fit <- dmcFitSubject(flankerData, nTrl = 1000, subjects = c(1, 2));plot(fit, flankerData, subject = 1)plot(fit, flankerData, subject = 2)summary(fit)

dmcFitSubjectDE

Description

Fit theoretical data generated from dmcSim to observed data byminimizing the root-mean-square error (RMSE) between a weighted combinationof the CAF and CDF functions using the R-package DEoptim. Alternative cost functionsinclude squared percentage error ("SPE"), and g-squared statistic ("GS").

Usage

dmcFitSubjectDE(  resOb,  nTrl = 1e+05,  minVals = list(),  maxVals = list(),  fixedFit = list(),  freeCombined = list(),  nCAF = 5,  nDelta = 19,  pDelta = vector(),  tDelta = 1,  deltaErrors = FALSE,  costFunction = "RMSE",  spDist = 1,  drOnset = 0,  drDist = 0,  drShape = 3,  drLim = c(0.1, 0.7),  rtMax = 5000,  subjects = c(),  deControl = list(),  numCores = 2)

Arguments

Value

dmcFitSubjectDE returns a list of objects of class "dmcfit"

Examples

# Code below can exceed CRAN check time limit, hence donttest# Example 1: Flanker data from Ulrich et al. (2015)fit <- dmcFitSubjectDE(flankerData, nTrl = 1000, subjects = c(1, 2), deControl = list(itermax=30))plot(fit, flankerData, subject = 1)plot(fit, flankerData, subject = 2)summary(fit)

dmcObservedData

Description

Basic analysis to create data object required for observed data.Example raw *.txt files are flankerData.txt and simonData.txt. There are four critical columns:

1. column containing subject number

2. column coding for compatible or incompatible

3. column with RT (in ms)

4. column indicating of the response was correct

Usage

dmcObservedData(  dat,  nCAF = 5,  nDelta = 19,  pDelta = vector(),  tDelta = 1,  outlier = c(200, 1200),  columns = c("Subject", "Comp", "RT", "Error"),  compCoding = c("comp", "incomp"),  errorCoding = c(0, 1),  quantileType = 5,  deltaErrors = FALSE,  keepRaw = FALSE,  delim = "\t",  skip = 0)

Arguments

Value

dmcObservedData returns an object of class "dmcob" with the following components:

Examples

# Example 1plot(flankerData)  # flanker data from Ulrich et al. (2015)plot(simonData)    # simon data from Ulrich et al. (2015)# Example 2 (Basic behavioural analysis from Ulrich et al. )flankerDat <- cbind(Task = "flanker", flankerData$summarySubject)simonDat   <- cbind(Task = "simon",   simonData$summarySubject)datAgg     <- rbind(flankerDat, simonDat)datAgg$Subject <- factor(datAgg$Subject)datAgg$Task    <- factor(datAgg$Task)datAgg$Comp    <- factor(datAgg$Comp)aovErr <- aov(perErr ~ Comp*Task + Error(Subject/(Comp*Task)), datAgg)summary(aovErr)model.tables(aovErr, type = "mean")aovRt <- aov(rtCor ~ Comp*Task + Error(Subject/(Comp*Task)), datAgg)summary(aovRt)model.tables(aovRt, type = "mean")# Example 3dat <- createDF(nSubjects = 50, nTrl = 500, design = list("Comp" = c("comp", "incomp")))dat <- addDataDF(dat,                 RT = list("Comp_comp"    = c(500, 75, 120),                           "Comp_incomp"  = c(530, 75, 100)),                 Error = list("Comp_comp" = c(3, 2, 2, 1, 1),                            "Comp_incomp" = c(21, 3, 2, 1, 1)))datOb <- dmcObservedData(dat)plot(datOb)plot(datOb, subject = 1)# Example 4dat <- createDF(nSubjects = 50, nTrl = 500, design = list("Congruency" = c("cong", "incong")))dat <- addDataDF(dat,                 RT = list("Congruency_cong"   = c(500, 75, 100),                           "Congruency_incong" = c(530, 100, 110)),                 Error = list("Congruency_cong"   = c(3, 2, 2, 1, 1),                              "Congruency_incong" = c(21, 3, 2, 1, 1)))datOb <- dmcObservedData(dat, nCAF = 5, nDelta = 9,                         columns = c("Subject", "Congruency", "RT", "Error"),                         compCoding = c("cong", "incong"))plot(datOb, labels = c("Congruent", "Incongruent"))plot(datOb, subject = 1)

dmcSim

Description

DMC model simulation detailed in Ulrich, R., Schroeter, H., Leuthold, H., & Birngruber, T. (2015).Automatic and controlled stimulus processing in conflict tasks: Superimposed diffusion processes and delta functions.Cognitive Psychology, 78, 148-174. This function is essentially a wrapper around the c++ function runDMC

Usage

dmcSim(  amp = 20,  tau = 30,  drc = 0.5,  bnds = 75,  resDist = 1,  resMean = 300,  resSD = 30,  aaShape = 2,  spShape = 3,  sigm = 4,  nTrl = 1e+05,  tmax = 1000,  spDist = 0,  spLim = c(-75, 75),  spBias = 0,  drOnset = 0,  drDist = 0,  drShape = 3,  drLim = c(0.1, 0.7),  rtMax = 5000,  fullData = FALSE,  nTrlData = 5,  nDelta = 9,  pDelta = vector(),  tDelta = 1,  deltaErrors = FALSE,  nCAF = 5,  bndsRate = 0,  bndsSaturation = 0,  printInputArgs = TRUE,  printResults = TRUE,  setSeed = FALSE,  seedValue = 1)

Arguments

Value

dmcSim returns an object of class "dmcsim" with the following components:

Examples

# Example 1dmc <- dmcSim(fullData = TRUE) # fullData only needed for activation/trials (left column plot)plot(dmc)dmc <- dmcSim() # faster!plot(dmc)# Example 2dmc <- dmcSim(tau = 130)plot(dmc)# Example 3dmc <- dmcSim(tau = 90)plot(dmc)# Example 4dmc <- dmcSim(spDist = 1)plot(dmc, "delta")# Example 5dmc <- dmcSim(tau = 130, drDist = 1)plot(dmc, "caf")# Example 6dmc <- dmcSim(nDelta = 10, nCAF = 10)plot(dmc)

dmcSimApp

Description

A shiny app allowing interactive exploration of DMC parameters

Usage

dmcSimApp()

Value

Shiny App

dmcSims: Run multiple dmc simulations

Description

Run dmcSim with range of input parameters.

Usage

dmcSims(params, printInputArgs = FALSE, printResults = FALSE)

Arguments

Value

dmcSims returns a list of objects of class "dmcsim"

Examples

# Example 1params <- list(amp = seq(10, 20, 5), tau = c(50, 100, 150), nTrl = 50000)dmc <- dmcSims(params)plot(dmc[[1]]) # full combination 1plot(dmc) # delta plots for all combinationsplot(dmc[c(1:3)]) # delta plots for specific combinationsplot(dmc[c(1, 3)]) # delta plots for specific combinations# Example 2params <- list(amp = seq(10, 20, 5), tau = seq(20, 40, 20), bnds = seq(50, 100, 25))dmc <- dmcSims(params)plot(dmc[[1]]) # combination 1plot(dmc, ncol = 2) # delta plots for all combinationsplot(dmc[c(1:3)]) # delta plots for specific combinations

errDist

Description

Returns a random vector of 0's (correct) and 1's (incorrect) withdefined proportions (default = 10% errors).

Usage

errDist(n = 10000, proportion = 10)

Arguments

Value

double

Examples

# Example 1x <- errDist(1000, 10)table(x)

A summarised dataset: This is the flanker task data from Ulrich et al. (2015)

Description

• $summary –> Reaction time correct, standard deviation correct, standard error correct,percentage error, standard deviation error, standard error error, reaction time incorrect,standard deviation incorrect, and standard error incorrect trials for both compatibleand incompatible trials

• $caf –> Proportion correct for compatible and incompatible trials across 5 bins

• $delta –> Compatible reactions times, incompatible mean reaction times,mean reaction times, incompatible - compatible reaction times (effect), andstandard deviation + standard error of this effect across 19 bins

• $data –> Raw data from flankerData.txt + additional outlier column

Usage

flankerData

Format

dmcob

mean.dmcfit

Description

Aggregate simulation results from dmcFitSubject/dmcFitSubjectDE.

Usage

## S3 method for class 'dmcfit_subject'mean(x, ...)

Arguments

Value

mean.dmcfit return an object of class "dmcfit" with the following components:

Examples

# Code below can exceed CRAN check time limit, hence donttest# Example 1: Fit individual data then aggregatefitSubjects <- dmcFitSubject(flankerData, nTrl = 1000, subjects = c(1, 2))fitAgg <- mean(fitSubjects)plot(fitAgg, flankerData)

plot.dmcfit: Plot observed + fitted data

Description

Plot the simulation results from the output of dmcFit. The plotcan be an overall summary, or individual plots (activation, trials, pdf, cdf,caf, delta, all). Plot type summary1 contains an activation plot, exampleindividual trials, the probability distribution function (PDF), the cumulativedistribution function (CDF), the conditional accuracy function (CAF) anddelta plots. This required that dmcSim is run with fullData = TRUE. Plot typesummary2 contains only the PDF, CDF, CAF and delta plots and does not requirethat dmcSim is run with fullData = TRUE.

Usage

## S3 method for class 'dmcfit'plot(  x,  y,  figType = "summary",  labels = c("Compatible", "Incompatible", "Observed", "Predicted"),  cols = c("green", "red"),  ylimRt = NULL,  ylimErr = NULL,  xlimCDF = NULL,  ylimCAF = NULL,  cafBinLabels = FALSE,  ylimDelta = NULL,  xlimDelta = NULL,  xlabs = TRUE,  ylabs = TRUE,  xaxts = TRUE,  yaxts = TRUE,  xylabPos = 2,  resetPar = TRUE,  legend = TRUE,  legend.parameters = list(legend = c("Observed", "Predicted")),  ...)

Arguments

Value

Plot (no return value)

Examples

# Example 1resTh <- dmcFit(flankerData, nTrl = 5000)plot(resTh, flankerData)plot(resTh, flankerData, figType = "deltaErrors")# Example 2resTh <- dmcFit(simonData, nTrl = 5000)plot(resTh, simonData)

plot.dmcfit_subject: Plot observed + fitted data

Description

Plot the simulation results from the output of dmcFit. The plotcan be an overall summary, or individual plots (activation, trials, pdf, cdf,caf, delta, all). Plot type summary1 contains an activation plot, exampleindividual trials, the probability distribution function (PDF), the cumulativedistribution function (CDF), the conditional accuracy function (CAF) anddelta plots. This required that dmcSim is run with fullData = TRUE. Plot typesummary2 contains only the PDF, CDF, CAF and delta plots and does not requirethat dmcSim is run with fullData = TRUE.

Usage

## S3 method for class 'dmcfit_subject'plot(  x,  y,  subject = NULL,  figType = "summary",  labels = c("Compatible", "Incompatible", "Observed", "Predicted"),  cols = c("green", "red"),  ylimRt = NULL,  ylimErr = NULL,  xlimCDF = NULL,  ylimCAF = NULL,  cafBinLabels = FALSE,  ylimDelta = NULL,  xlimDelta = NULL,  xlabs = TRUE,  ylabs = TRUE,  xaxts = TRUE,  yaxts = TRUE,  xylabPos = 2,  resetPar = TRUE,  legend = TRUE,  legend.parameters = list(legend = c("Observed", "Predicted")),  ...)

Arguments

Value

Plot (no return value)

Examples

# Example 1resTh <- dmcFitSubject(flankerData, nTrl = 5000, subject = c(1,3))plot(resTh, flankerData, subject = 3)

plot.dmcfits: Plot observed + fitted data

Description

Plot the simulation results from the output of dmcFit. The plotcan be an overall summary, or individual plots (activation, trials, pdf, cdf,caf, delta, all). Plot type summary1 contains an activation plot, exampleindividual trials, the probability distribution function (PDF), the cumulativedistribution function (CDF), the conditional accuracy function (CAF) anddelta plots. This required that dmcSim is run with fullData = TRUE. Plot typesummary2 contains only the PDF, CDF, CAF and delta plots and does not requirethat dmcSim is run with fullData = TRUE.

Usage

## S3 method for class 'dmcfits'plot(  x,  y,  figType = "summary",  labels = c("Compatible", "Incompatible", "Observed", "Predicted"),  cols = c("green", "red"),  ylimRt = NULL,  ylimErr = NULL,  xlimCDF = NULL,  ylimCAF = NULL,  cafBinLabels = FALSE,  ylimDelta = NULL,  xlimDelta = NULL,  xlabs = TRUE,  ylabs = TRUE,  xaxts = TRUE,  yaxts = TRUE,  xylabPos = 2,  resetPar = TRUE,  legend = TRUE,  legend.parameters = list(legend = c("Observed", "Predicted")),  ...)

Arguments

Value

Plot (no return value)

Examples

# Example 1resTh <- dmcFit(flankerData, nTrl = 5000)plot(resTh, flankerData)plot(resTh, flankerData, figType = "deltaErrors")# Example 2resTh <- dmcFit(simonData, nTrl = 5000)plot(resTh, simonData)

plot.dmcfits_subject: Plot observed + fitted data

Description

Plot the simulation results from the output of dmcFit. The plotcan be an overall summary, or individual plots (activation, trials, pdf, cdf,caf, delta, all). Plot type summary1 contains an activation plot, exampleindividual trials, the probability distribution function (PDF), the cumulativedistribution function (CDF), the conditional accuracy function (CAF) anddelta plots. This required that dmcSim is run with fullData = TRUE. Plot typesummary2 contains only the PDF, CDF, CAF and delta plots and does not requirethat dmcSim is run with fullData = TRUE.

Usage

## S3 method for class 'dmcfits_subject'plot(  x,  y,  subject = NULL,  figType = "summary",  labels = c("Compatible", "Incompatible", "Observed", "Predicted"),  cols = c("green", "red"),  ylimRt = NULL,  ylimErr = NULL,  xlimCDF = NULL,  ylimCAF = NULL,  cafBinLabels = FALSE,  ylimDelta = NULL,  xlimDelta = NULL,  xlabs = TRUE,  ylabs = TRUE,  xaxts = TRUE,  yaxts = TRUE,  xylabPos = 2,  resetPar = TRUE,  legend = TRUE,  legend.parameters = list(legend = c("Observed", "Predicted")),  ...)

Arguments

Value

Plot (no return value)

Examples

# Example 1resTh <- dmcFit(flankerData, nTrl = 5000)plot(resTh, flankerData)plot(resTh, flankerData, figType = "deltaErrors")# Example 2resTh <- dmcFit(simonData, nTrl = 5000)plot(resTh, simonData)

plot.dmclist: Plot delta plots from multiple dmc simulations.

Description

Plot delta function from multiple dmc simulations (i.e., dmcSims).

Usage

## S3 method for class 'dmclist'plot(  x,  ylim = NULL,  xlim = NULL,  figType = "delta",  xlab = "Time [ms]",  ylab = expression(paste(Delta, "Time [ms]")),  xylabPos = 2,  col = c("black", "lightgrey"),  lineType = "l",  legend = TRUE,  legend.parameters = list(),  ...)

Arguments

Value

Plot (no return value)

Examples

# Example 1params <- list(amp = seq(20, 30, 2))dmc <- dmcSims(params)plot(dmc, col = c("red", "green"), legend.parameters = list(x = "topright", ncol=2))# Example 2params <- list(amp=c(10, 20), tau = c(20, 40), drc = c(0.2, 0.6), nTrl = 50000)dmc <- dmcSims(params)plot(dmc, col=c("green", "blue"), ylim = c(-10, 120), legend.parameters=list(ncol=2))

plot.dmcob: Plot observed data

Description

Plot results from the output of dmcObservedData. The plotcan be an overall summary, or individual plots (rtCorrect, errorRate,rtErrors, cdf, caf, delta, deltaErrors, all).

Usage

## S3 method for class 'dmcob'plot(  x,  figType = "summary",  subject = NULL,  labels = c("Compatible", "Incompatible"),  cols = c("green", "red"),  errorBars = FALSE,  errorBarType = "sd",  ylimRt = NULL,  ylimErr = NULL,  xlimCDF = NULL,  ylimCAF = NULL,  cafBinLabels = FALSE,  ylimDelta = NULL,  xlimDelta = NULL,  xlabs = TRUE,  ylabs = TRUE,  xaxts = TRUE,  yaxts = TRUE,  xylabPos = 2,  resetPar = TRUE,  legend = TRUE,  ...)

Arguments

Value

Plot (no return value)

Examples

# Example 1 (real dataset)plot(flankerData)plot(flankerData, errorBars = TRUE, errorBarType = "se")plot(flankerData, figType = "delta")plot(flankerData, figType = "caf")# Example 2 (real dataset)plot(simonData)plot(simonData, errorBars = TRUE, errorBarType = "se")plot(simonData, figType = "delta", errorBars = TRUE, errorBarType = "sd")# Example 3 (simulated dataset)dat <- createDF(nSubjects = 50, nTrl = 50,                design = list("Comp" = c("comp", "incomp")))dat <- addDataDF(dat,                 RT = list("Comp_comp"   = c(420, 100, 80),                           "Comp_incomp" = c(470, 100, 95)),                 Error = list("Comp_comp"   = c(5, 3, 2, 1, 2),                              "Comp_incomp" = c(15, 8, 4, 2, 2)))datOb <- dmcObservedData(dat)plot(datOb, errorBars = TRUE, errorBarType = "sd")# Example 4 (simulated dataset)dat <- createDF(nSubjects = 50, nTrl = 50,                design = list("Comp" = c("comp", "incomp")))dat <- addDataDF(dat,                 RT = list("Comp_comp"   = c(420, 100, 150),                           "Comp_incomp" = c(470, 100, 120)),                 Error = list("Comp_comp"   = c(5, 3, 2, 1),                              "Comp_incomp" = c(15, 8, 4, 2)))datOb <- dmcObservedData(dat, nCAF = 4)plot(datOb)

plot.dmcobs: Plot combined observed data

Description

Plot delta results from the output of dmcObservedData. The plotcan be an overall rtCorrect, errorRate, rtErrors, cdf, caf, delta, deltaErrors, deltaER, or allof the previous plots.

Usage

## S3 method for class 'dmcobs'plot(  x,  figType = "all",  subject = NULL,  labels = c("Compatible", "Incompatible"),  cols = c("black", "gray"),  ltys = c(1, 1),  pchs = c(1, 1),  errorBars = FALSE,  errorBarType = "sd",  ylimRt = NULL,  ylimErr = NULL,  xlimCDF = NULL,  ylimCAF = NULL,  cafBinLabels = FALSE,  ylimDelta = NULL,  xlimDelta = NULL,  xlabs = TRUE,  ylabs = TRUE,  xaxts = TRUE,  yaxts = TRUE,  xylabPos = 2,  resetPar = TRUE,  legend = TRUE,  legend.parameters = list(),  ...)

Arguments

Value

Plot (no return value)

Examples

# Example 1dat <- dmcCombineObservedData(flankerData, simonData)  # combine flanker/simon dataplot(dat, figType = "all", xlimDelta = c(200, 700), ylimDelta = c(-20, 80),     cols = c("black", "darkgrey"), pchs = c(1, 2))plot(dat, figType = "delta", xlimDelta = c(200, 700), ylimDelta = c(-20, 80),     cols = c("black", "darkgrey"), pchs = c(1, 2), legend = TRUE,     legend.parameters=list(x="topright", legend=c("Flanker", "Simon")))

plot.dmcsim: Plot dmc simulation

Description

Plot the simulation results from the output of dmcSim. The plotcan be an overall summary, or individual plots (activation, trials, pdf, cdf,caf, delta, all). Plot type summary1 contains an activation plot, exampleindividual trials, the probability distribution function (PDF), the cumulativedistribution function (CDF), the conditional accuracy function (CAF) anddelta plot. This requires that dmcSim is run with fullData = TRUE. Plot typesummary2 contains only the PDF, CDF, CAF and delta plots and does not requirethat dmcSim is run with fullData = TRUE.

Usage

## S3 method for class 'dmcsim'plot(  x,  figType = "summary1",  xlimActivation = NULL,  ylimActivation = NULL,  xlimTrials = NULL,  ylimTrials = NULL,  xlimPDF = NULL,  ylimPDF = NULL,  xlimCDF = NULL,  ylimCAF = NULL,  cafBinLabels = FALSE,  ylimDelta = NULL,  xlimDelta = NULL,  ylimRt = NULL,  ylimErr = NULL,  labels = c("Compatible", "Incompatible"),  cols = c("green", "red"),  errorBars = FALSE,  xlabs = TRUE,  ylabs = TRUE,  xaxts = TRUE,  yaxts = TRUE,  xylabPos = 2,  resetPar = TRUE,  legend = TRUE,  ...)

Arguments

Value

Plot (no return value)

Examples

# Example 1dmc = dmcSim(fullData = TRUE)plot(dmc)# Example 2dmc = dmcSim()plot(dmc)# Example 3dmc = dmcSim(tau = 120)plot(dmc)# Example 4dmc = dmcSim()plot(dmc, figType = "all")

rtDist

Description

Returns value(s) from a distribution appropriate to simulate reaction times.The distribution is a combined exponential and gaussian distribution calledan exponentially modified Gaussian (EMG) distribution or ex-gaussian distribution.

Usage

rtDist(n = 10000, gaussMean = 600, gaussSD = 50, expRate = 200)

Arguments

Value

double

Examples

# Example 1x <- rtDist()hist(x, 100, xlab = "RT [ms]")# Example 2x <- rtDist(n=2000, gaussMean=500, gaussSD=100, expRate=300)hist(x, 100, xlab = "RT [ms]")

A summarised dataset: This is the simon task data from Ulrich et al. (2015)

Description

• $summary –> Reaction time correct, standard deviation correct, standard error correct,percentage error, standard deviation error, standard error error, reaction time incorrect,standard deviation incorrect, and standard error incorrect trials for both compatibleand incompatible trials

• $caf –> Proportion correct for compatible and incompatible trials across5 bins

• $delta –> Compatible reactions times, incompatible mean reaction times,mean reaction times, incompatible - compatible reaction times (effect), andstandard deviation + standard error of this effect across 19 bins

• $data –> Raw data from simonData.txt + additional outlier column

Usage

simonData

Format

dmcob

summary.dmcfit: dmc fit aggregate summary

Description

Summary of the simulation results from dmcFit

Usage

## S3 method for class 'dmcfit'summary(object, digits = 2, ...)

Arguments

Value

DataFrame

Examples

# Example 1fitAgg <- dmcFit(flankerData, nTrl = 1000)summary(fitAgg)

summary.dmcfit_subject: dmcfit individual subject

Description

Summary of the simulation results from dmcFitSubjectX

Usage

## S3 method for class 'dmcfit_subject'summary(object, digits = 2, ...)

Arguments

Value

DataFrame

Examples

# Example 1fitSubject <- dmcFitSubject(flankerData, nTrl = 1000, subjects = c(1:3))summary(fitSubject)

summary.dmcfits: dmc fit aggregate summary (2+ data sets)

Description

Summary of the simulation results from dmcFit

Usage

## S3 method for class 'dmcfits'summary(object, digits = 2, ...)

Arguments

Value

DataFrame

Examples

# Example 1fitAggs <- dmcFit(list(flankerData, simonData), nTrl = 1000)summary(fitAggs)

summary.dmcfits_subject: dmc fit aggregate summary

Description

Summary of the simulation results from dmcFitAgg

Usage

## S3 method for class 'dmcfits_subject'summary(object, digits = 2, ...)

Arguments

Value

DataFrame

Examples

# Example 1fitsSubject <- dmcFitSubject(list(flankerData, simonData), nTrl = 1000, subjects = c(1:3))summary(fitsSubject)

summary.dmcsim: dmc simulation summary

Description

Summary of the overall results from dmcSim

Usage

## S3 method for class 'dmcsim'summary(object, digits = 1, ...)

Arguments

Value

DataFrame

Examples

# Example 1dmc <- dmcSim()summary(dmc)# Example 2dmc <- dmcSim(tau = 90)summary(dmc)