Learn R Programming

⚠️There's a newer version (3.5.0) of this package.Take me there.

EMC2: Extended Models of Choice 2:

The R package EMC2 provides tools to perform Bayesian hierarchical analyses of the following cognitive models: Diffusion Decision Model (DDM), Linear Ballistic Accumulator Model (LBA), Racing Diffusion Model (RDM), and Lognormal Racing Model (LNR). Specifically, the package provides functionality for specifying individual model designs, estimating the models, examining convergence as well as model fit through posterior prediction methods. It also includes various plotting functions and relative model comparison methods such as Bayes factors. In addition, users can specify their own likelihood function and perform non-hierarchical estimation. The package uses particle metropolis Markov chain Monte Carlo sampling. For hierarchical models, it uses efficient Gibbs sampling at the population level and supports a variety of covariance structures, extending the work of Gunawan and colleagues (2020).

Installation

To install the R package, and its dependencies you can use

install.packages("EMC2")

Or for the development version:

remotes::install_github("ampl-psych/EMC2",dependencies=TRUE)

Workflow Overview

Pictured below are the four phases of an EMC2cognitive model analysis with associated functions:.

 

 

For details, please see:

Stevenson, N., Donzallaz, M. C., Innes, R. J., Forstmann, B., Matzke, D., & Heathcote, A. (2024, January 30). EMC2: An R Package for cognitive models of choice. https://doi.org/10.31234/osf.io/2e4dq

Bug Reports, Contributing, and Feature Requests

If you come across any bugs, or have ideas for extensions of EMC2, you can add them as an issue here. If you would like to contribute to the package's code, please submit a pull request.

References

Stevenson, N., Donzallaz, M. C., Innes, R. J., Forstmann, B., Matzke, D., & Heathcote, A. (2024, January 30). EMC2: An R Package for cognitive models of choice. https://doi.org/10.31234/osf.io/2e4dq

Gunawan, D., Hawkins, G. E., Tran, M. N., Kohn, R., & Brown, S. D. (2020). New estimation approaches for the hierarchical Linear Ballistic Accumulator model. Journal of Mathematical Psychology, 96, 102368. https://doi.org/10.1016/j.jmp.2020.102368

Copy Link

Version

Install

install.packages('EMC2')

Monthly Downloads

604

Version

3.4.1

License

GPL (>= 3)

Issues

Pull Requests

Stars

Forks

Maintainer

Niek Stevenson

Last Published

January 12th, 2026

Functions in EMC2 (3.4.1)

auto_thin.emc

Automatically Thin an emc Object
credint.emc.prior

Posterior Quantiles
fit.emc

Model Estimation in EMC2
fix_custom_kernel_pointers

Reset pointers of custom C++ trend kernels to an emc object
factor_diagram

Factor diagram plot #Makes a factor diagram plot. Heavily based on the fa.diagram function of the psych package.
ess_summary.emc

Effective Sample Size
get_group_design.emc.prior

Get Group Design
credible.emc

Posterior Credible Interval Tests
convolve_design_matrix

Convolve Events with HRF to Construct Design Matrices
get_pars

Filter/Manipulate Parameters from emc Object
get_data.emc

Get Data
cut_factors

Cut Factors Based on Credible Loadings
get_design.emc.prior

Get Design
get_power_spectra

Compute Power Spectra With Optional Subject-Level Aggregation
get_prior.emc

Get Prior
get_BayesFactor

Bayes Factors
get_custom_kernel_pointers

Extract pointers of custom C++ trend kernels from trend list or emc object
design_fmri

Create fMRI Design for EMC2 Sampling
make_sem_structure

Define Structural Equation Model (SEM) Matrices
make_random_effects

Generate Subject-Level Parameters
design

Specify a Design and Model
gd_summary.emc

Gelman-Rubin Statistic
forstmann

Forstmann et al.'s Data
group_design

Create Group-Level Design Matrices
get_trend_pnames

Get parameter types from trend object
make_SEM_diagram

Make SEM Diagram
init_chains

Initialize Chains
make_emc

Make an emc Object
mapped_pars

Parameter Mapping Back to the Design Factors
make_data

Simulate Data
model_averaging

Model Averaging
make_trend

Create a trend specification for model parameters
merge_chains

Merge Samples
pairs_posterior

Plot Within-Chain Correlations
high_pass_filter

Apply High-Pass Filtering to fMRI Data
hypothesis.emc

Within-Model Hypothesis Testing
parameters.emc.prior

Return Data Frame of Parameters
plot_relations

Plot Group-Level Relations
plot_pars

Plots Density for Parameters
plot_sbc_ecdf

Plot the ECDF Difference in SBC Ranks
plot.emc.design

Plot method for emc.design objects
plot.emc

Plot Function for emc Objects
plot_sbc_hist

Plot the Histogram of the Observed Rank Statistics of SBC
prior_help

Prior Specification Information
prior

Specify Priors for the Chosen Model
plot_cdf

Plot Defective Cumulative Distribution Functions
plot_stat

Plot Statistics on Data
plot_spectrum

Plot Empirical and Posterior Predictive Power Spectra
plot_delta

Plot Difference of Cumulative Distribution Functions
recovery.emc

Recovery Plots
plot_design_fmri

Plot fMRI Design Matrix
plot_fmri

Plot fMRI peri-stimulus time courses
profile_plot

Likelihood Profile Plots
register_trend

Register a custom C++ trend kernel
reshape_events

Reshape events data for fMRI analysis
plot_trend

Plots trends over time
sampled_pars

Get Model Parameters from a Design
run_sbc

Simulation-Based Calibration
plot.emc.prior

Plot a prior
update2version

Update EMC Objects to the Current Version
trend_help

Get help information for trend kernels and bases
set_custom_kernel_pointers

(Re-)Set pointers of custom C++ trend kernels to an emc object
samples_LNR

LNR Model of Forstmann Data (First 3 Subjects)
summary.emc

Summary Statistics for emc Objects
summary.emc.design

Summary method for emc.design objects
plot_caf

Plot conditional accuracy functions
run_hyper

Run a Group-level Model.
run_emc

Fine-Tuned Model Estimation
predict.emc.prior

Generate Posterior/Prior Predictives
plot_design.emc.design

Plot Design
rotate_loadings

Rotate loadings based on posterior median
run_bridge_sampling

Estimating Marginal Likelihoods Using WARP-III Bridge Sampling
plot_density

Plot Defective Densities
summary.emc.group_design

Summary method for emc.group_design objects
split_timeseries

Split fMRI Timeseries Data by ROI Columns
subset.emc

Shorten an emc Object
summary.emc.prior

Summary method for emc.prior objects
DDM

The Diffusion Decision Model
LNR

The Log-Normal Race Model
MRI

GLM model for fMRI data
RDM

The Racing Diffusion Model
MRI_AR1

Create an AR(1) GLM model for fMRI data
SDT

Gaussian Signal Detection Theory Model for Binary Responses
EMC2-package

EMC2: Bayesian Hierarchical Analysis of Cognitive Models of Choice
align_loadings

Reorder MCMC Samples of Factor Loadings
LBA

The Linear Ballistic Accumulator model
DDMGNG

The GNG (go/nogo) Diffusion Decision Model
contr.increasing

Contrast Enforcing Increasing Estimates
contr.bayes

Contrast Enforcing Equal Prior Variance on each Level
compare

Information Criteria and Marginal Likelihoods
apply_kernel

Apply a kernel implied in an emc object
contr.decreasing

Contrast Enforcing Decreasing Estimates
chain_n

MCMC Chain Iterations
contr.anova

Anova Style Contrast Matrix
compare_subject

Information Criteria For Each Participant
check.emc

Convergence Checks for an emc Object