A multi-frequency whole-brain neural mass model with homeostatic feedback inhibition.
The 17 matches
- [1] § 4. Methods › 4.5 Hemodynamic model ↔ JR ReDLat + Sleep/BOLDModel.py, lines 25–43 · score 0.93 · oxygen extraction, signal decay, stiffness constant, blood volume, deoxyhemoglobin content, resistance
- [2] § 2. Results ↔ JR Model Tutorial/JansenRitModelMulti.py, lines 1–84 · score 0.90 · Jansen Rit model, neural activity, modified version, cortical columns, pyramidal neurons, synaptic plasticity
- [3] § 2. Results ↔ JR model and running example/JansenRitModelMulti.py, lines 1–84 · score 0.90 · Jansen Rit model, neural activity, modified version, cortical columns, pyramidal neurons, synaptic plasticity
- [4] § 4. Methods › 4.2 Functional connectivity estimation ↔ JR Model Tutorial/Whole-Brain Part 1.ipynb, lines 206–267 · score 0.81 · 0.5–4 Hz, 13–30 Hz, 8–13 Hz, EEG frequency bands, 4–8 Hz, Bessel
- [5] § 2. Results › 2.3 Fitting to EEG source functional connectivity data ↔ JR Model Tutorial/Whole-Brain Part 1.ipynb, lines 206–267 · score 0.77 · 30–40 Hz, 0.5–4 Hz, 13–30 Hz, 8–13 Hz, frequency bands, 4–8 Hz
- [6] § 4. Methods › 4.4 Whole-brain model ↔ JR Model Tutorial/JansenRitModelMulti.py, lines 1–84 · score 0.76 · neural mass, modified version, pyramidal neurons, Jansen Rit, brain region, excitatory
- [7] § 2. Results › 2.1 Analysis of Jansen-Rit single-node dynamics ↔ JR ReDLat + Sleep/Single_Node.py, lines 232–272 · score 0.75 · single node, Peak frequency, Noise free, feedback inhibition, steady state, target firing rate
- [8] § 4. Methods › 4.2 Functional connectivity estimation ↔ JR Model Tutorial/runJR.py, lines 74–82 · score 0.72 · 0.5–4 Hz, 13–30 Hz, 8–13 Hz, 4–8 Hz, Bessel, filtered
- [9] § 4. Methods › 4.4 Whole-brain model ↔ JR ReDLat + Sleep/JansenRitModelMulti.py, lines 1–79 · score 0.72 · Jansen Rit neural, neural masses, modified version, mass model, pyramidal neurons, excitatory
- [10] § 2. Results › 2.3 Fitting to EEG source functional connectivity data ↔ JR Model Tutorial/runJR.py, lines 74–82 · score 0.71 · 30–40 Hz, 0.5–4 Hz, 13–30 Hz, 8–13 Hz, 4–8 Hz, filtered
- [11] § 4. Methods › 4.4 Whole-brain model ↔ JR model and running example/JansenRitModelMulti.py, lines 1–84 · score 0.70 · sigmoid function, pyramidal neurons, target firing rate, bounding, min, inverse
- [12] § 2. Results › 2.1 Analysis of Jansen-Rit single-node dynamics ↔ JR ReDLat + Sleep/Single_Node.py, lines 1–43 · score 0.66 · single node simulations, feedback inhibition, power spectrum, target firing rates, external, Jansen
- [13] § 4. Methods › 4.6 Simulations ↔ JR Model Tutorial/Whole-Brain Part 1.ipynb, lines 43–142 · score 0.66 · reduce memory consumption, spectral density, Power Spectrum, Welch, downsampled, plasticity
- [14] § 4. Methods › 4.5 Hemodynamic model ↔ JR ReDLat + Sleep/BOLDModel.py, lines 54–83 · score 0.65 · blood inflow, deoxyhemoglobin content, vasodilatory, firing rate, signals, simulated
- [15] § 4. Methods › 4.5 Hemodynamic model ↔ JR ReDLat + Sleep/BOLDModel.py, lines 25–43 · score 0.56 · venous blood, deoxygenated, fraction, V0, model
- [16] § 4. Methods › 4.2 Functional connectivity estimation ↔ JR ReDLat + Sleep/Sleep_EEG_fMRI.py, lines 80–150 · score 0.53 · 0.01–0.08 Hz, FC matrix, Pearson, correlation, 0.01 Hz, filtered
- [17] § 2. Results › 2.1 Analysis of Jansen-Rit single-node dynamics ↔ JR ReDLat + Sleep/Single_Node.py, lines 1–43 · score 0.51 · fast limit cycles, Jansen Rit, external, node, dynamics, model
Paper
Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
Python · 157 lines · 4.5 KB · no license · 3 matches
- #!/usr/bin/env python3
- # -*- coding: utf-8 -*-
- """
- Created on Thu Aug 8 15:18:02 2019
- @author: Carlos Coronel
- Generalized Hemodynamic Model to reproduce fMRI BOLD-like signals.
- [1] Stephan, K. E., Weiskopf, N., Drysdale, P. M., Robinson, P. A., & Friston, K. J.
- (2007). Comparing hemodynamic models with DCM. Neuroimage, 38(3), 387-401.
- [2] Deco, Gustavo, et al. "Whole-brain multimodal neuroimaging model using serotonin
- receptor maps explains non-linear functional effects of LSD." Current Biology
- 28.19 (2018): 3065-3074.
- """
- import numpy as np
- from numba import jit,float64
- from numba.core.errors import NumbaPerformanceWarning
- import warnings
- warnings.simplefilter('ignore', category=NumbaPerformanceWarning)
- #PARAMETERS
- taus = 0.65 #time constant for signal decay
- tauf = 0.41 #time constant for feedback regulation
- tauo = 0.98 #time constant for volume and deoxyhemoglobin content change
- #itauX = inverse of tauX
- itaus = 1 / taus
- itauf = 1 / tauf
- itauo = 1 / tauo
- nu = 40.3 #frequency offset at the outer surface of the magnetized
- #vessel for fully deoxygenated blood at 1.5 Tesla (s^-1)
- r0 = 25 #slope of the relation between the intravascular relaxation
- #rate and oxygen saturation (s^-1)
- alpha = 0.32 #resistance of the veins; stiffness constant
- ialpha = 1 / alpha
- epsilon = 0.5 #ratio of intra and extravascular signal
- E0 = 0.4 #resting oxygen extraction fraction
- TE = 0.04 #echo time (!!determined by the experiment)
- V0 = 0.04 #resting venous blood volume fraction
- #Kinetics constants
- k1 = 4.3 * nu * E0 * TE
- k2 = epsilon * r0 * E0 * TE
- k3 = 1 - epsilon
- def update():
- BOLD_response.recompile()
- BOLD_signal.recompile()
- @jit(float64[:,:](float64[:,:],float64[:],float64), nopython = True)
- def BOLD_response(y, rE, t):
- """
- This function generates a BOLD response using the firing rates rE.
- Parameters
- ----------
- y : numpy array.
- Contains the following variables:
- s: vasodilatory signal. If it increases, the blood vessel experiments
- vasodilatation.
- f: blood inflow. Increases with the vasodilatation.
- v: blood volumen. Increases with blood inflow.
- q: deoxyhemoglobin content.
- rE: numpy array.
- Firing rates of neural populations/neurons.
- t : float.
- Current simulation time point.
- Returns
- -------
- Numpy array with s, f, v and q derivatives at time t.
- """
- s, f, v, q = y
- s_dot = 1 * rE + 0 - itaus * s - itauf * (f - 1)
- f_dot = s
- v_dot = (f - v ** ialpha) * itauo
- q_dot = (f * (1 - (1 - E0) ** (1 / f)) / E0 - q * v ** ialpha / v) * itauo
- return(np.vstack((s_dot, f_dot, v_dot, q_dot)))
- @jit(float64[:,:](float64[:,:],float64[:,:]), nopython = True)
- def BOLD_signal(q, v):
- """
- This function returns the BOLD signal using deoxyhemoglobin content and
- blood volumen as inputs.
- Parameters
- ----------
- q: numpy array.
- deoxyhemoglobin content over time.
- v: numpy array.
- blood volumen over time.
- Returns
- -------
- Numpy array with BOLD-like signals.
- """
- return(V0 * (k1 * (1 - q) + k2 * (1 - q / v) + k3 * (1 - v)))
- def Sim(rE, nnodes, dt):
- """
- Simulate the BOLD-like signals (raw non-filtered) with the current parameter values.
- Note that the time unit in this model is seconds.
- Parameters
- ----------
- rE : numpy array
- time x nodes matrix which values contains the firing rates of each node.
- nnodes : integer
- number of nodes.
- dt : float.
- actual integration step (the inverse of the sampling rate)
- Raises
- ------
- ValueError
- An error raises if the number of nodes of rE did not match with the given
- number by nnodes
- Returns
- -------
- y : numpy array
- Raw BOLD-like signals for each node
- """
- Ntotal = rE.shape[0]
- ic_BOLD = np.ones((1, nnodes)) * np.array([0.1, 1, 1, 1])[:, None] #initial conditions
- BOLD_vars = np.zeros((Ntotal,4,nnodes)) #matrix for storing the values
- BOLD_vars[0,:,:] = ic_BOLD
- #Solve the ODEs with Euler
- for i in range(1,Ntotal):
- BOLD_vars[i,:,:] = BOLD_vars[i - 1,:,:] + dt * BOLD_response(BOLD_vars[i - 1,:,:], rE[i - 1,:], i - 1)
- y = BOLD_signal(BOLD_vars[:,3,:], BOLD_vars[:,2,:])
- return(y)
- def ParamsBOLD():
- pardict={}
- for var in ('taus','tauf','tauo','nu','r0','alpha','epsilon','E0',
- 'V0','TE','k1','k2','k3'):
- pardict[var]=eval(var)
- return pardict
BOLDModel.py at commit 13fab02, no license · at the source
Overview
15 affiliations
- Latin American Brain Health Institute (BrainLat), Universidad Adolfo Ibañez, Santiago, Chile
- Trinity College Dublin, The University of Dublin, Dublin, Ireland
- Global Brain Health Institute (GBHI), Trinity College Dublin, Dublin, Ireland
- Advanced Center for Electrical and Electronic Engineering, Universidad Técnica Federico Santa María, Valparaíso, Chile
- Department of Psychology, University of the Balearic Islands, Palma de Mallorca, Spain
- Institut du Cerveau - Paris Brain Institute - ICM, Sorbonne Université, Inserm, CNRS, Paris, France
- Network Science Institute, School of Medicine, Northeastern University London, London, United Kingdom
- Precision Imaging, School of Medicine, University of Nottingham, Nottingham, United Kingdom
- Faculty of Psychology, SWPS University of Social Sciences and Humanities, Warsaw, Poland
- Faculty of Biological Sciences, Pontifical Catholic University of Chile, Santiago, Chile
- Escuela de Fonoaudiología, Facultad de Ciencias de la Rehabilitación y Calidad de Vida, Universidad San Sebastián, Santiago, Chile
- Centro Interdisciplinario de Neurociencia de Valparaíso (CINV), Universidad de Valparaíso, Valparaíso, Chile
- Instituto de Neurociencia, Facultad de Ciencias, Universidad de Valparaíso, Playa Ancha, Valparaíso, Chile
- Barcelonaβeta Brain Research Center (BBRC), Pasqual Maragall Foundation, Barcelona, Spain
- Department of Biophysics, School of Medicine, Istanbul Medipol University, Istanbul, Türkiye
Abstract
Whole-brain models are valuable tools for understanding brain dynamics in health and disease by enabling the testing of causal mechanisms and identification of therapeutic targets through dynamic simulations. Among these models, biophysically inspired neural mass models have been widely used to simulate electrophysiological recordings, such as MEG and EEG. However, traditional models face limitations, including susceptibility to over-saturation of the sigmoid function by model hyperexcitability, which constrains their ability to capture the full richness of neural dynamics. Here, we thoroughly characterize a previously introduced multi-frequency Jansen-Rit neural mass model with inhibitory synaptic plasticity (ISP) aimed at overcoming these limitations. The ISP adjusts inhibitory feedback onto pyramidal neurons to clamp their firing rates around a target value. This mechanism allows for fine control of neuronal firing rates, preventing over-saturation in whole-brain simulations. In this model, we analyzed how different model parameters modulate oscillatory frequency and connectivity. As a demonstration, we considered simultaneously fitting EEG and fMRI recordings during NREM sleep. Bifurcation analysis showed that ISP widened the range of parameters in which the model exhibited sustained oscillations; the target firing rate can modulate oscillatory dynamics, producing different oscillatory regimes, from slower (δ, θ and α) to faster (β and γ) oscillations. High-frequency activity emerged from low global coupling, high firing rates, and a high proportion of γ versus α subpopulations. The ISP was necessary in the multi-frequency model to successfully fit EEG functional connectivity across frequency bands. Finally, ISP-controlled reductions in excitability reproduced both the slow-wave activity and the reduced connectivity in NREM sleep. Altogether, our model is compatible with biological evidence of the effects of excitability on modulating brain rhythms and connectivity, as observed in sleep, neurodegeneration, and chemical neuromodulation. This biophysical model with ISP provides a springboard for realistic brain simulations in health and disease.
Reproduced under the paper's license (CC BY), from the paper cited above.
Repositories
Its files are read in the Code ↔ Paper reader above, with 17 matches between paragraphs and lines of code.
carlosmig/EEG-Dementias
13fab025b3f04b6faf47c5704c6f9b893ee60cb0, 2 September 2026Availability: 1 check, the latest on 28 September 2026: the link answers
- 28 September 2026: the link answers
23 files
- JR Model Tutorial/
JansenRitModelMulti.py , Python, 276 lines, 2 matches - JR Model Tutorial/
Whole-Brain Part 1.ipynb , Jupyter, 385 lines, 3 matches - JR Model Tutorial/
Whole-Brain Part 2.ipynb , Jupyter, 508 lines - JR Model Tutorial/
Whole-Brain Part 3.ipynb , Jupyter, 239 lines - JR Model Tutorial/
runJR.py , Python, 143 lines, 2 matches - JR ReDLat + Sleep/
BOLDModel.py , Python, 157 lines, 3 matches - JR ReDLat + Sleep/
EEG_fMRI_showcase.py , Python, 149 lines - JR ReDLat + Sleep/
EEG_traces.py , Python, 169 lines - JR ReDLat + Sleep/
ISP_slopes.py , Python, 140 lines - JR ReDLat + Sleep/
JansenRitModelMulti.py , Python, 337 lines, 1 match - JR ReDLat + Sleep/
ReDLat_EEG_FCs.py , Python, 191 lines - JR ReDLat + Sleep/
Single_Node.py , Python, 393 lines, 3 matches - JR ReDLat + Sleep/
Sleep_EEG_fMRI.py , Python, 243 lines, 1 match - JR ReDLat + Sleep/
Tau_Convergence.py , Python, 93 lines - JR Single Subject Fitting/
JR_fitting.py , Python, 192 lines - JR Single Subject Fitting/
JansenRitModelMulti.py , Python, 276 lines - JR model and running example/
JansenRitModelMulti.py , Python, 276 lines, 2 matches - JR model and running example/
runJR.py , Python, 143 lines - Metaconnectivity/
metaconnectivity.py , Python, 158 lines - SVM and FC matrices/
SVM_main.py , Python, 414 lines - SVM and FC matrices/
SVM_reg_toy_data.py , Python, 277 lines - SVM and FC matrices/
data_augmentation.py , Python, 81 lines - README.md, Text, 38 lines
Zenodo 18841259
Availability: 1 check, the latest on 28 September 2026: the link answers (HTTP 200)
- 28 September 2026: the link answers (HTTP 200)
The paper's code and data availability statement is in the Data section.
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 2 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 22 scripts, each with its path and the digest of its content;
- 17 matches between paragraphs of the paper and lines of the code (method lexical-v1);
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
Datasets cited
- zenodo:16755776, at Zenodo; found in “Data Availability”
Data Availability
All the codes and basic data used to perform the simulations are freely available in GitHub at https://
Reproduced under the paper's license (CC BY), from the paper cited above.
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 1, 28 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 14 authors, 12 MeSH terms, 10 funders, 81 references.
Cite
This paper
Coronel-Oliveros, C., Lehue, F., Herzog, R., Mindlin, I., Gatica, M., Kowalczyk-Grębska, N., Medel, V., Cruzat, J., Gonzalez-Gomez, R., Hernandez, H., Tagliazucchi, E., Prado, P., Orio, P., & Ibáñez, A. (2026). A multi-frequency whole-brain neural mass model with homeostatic feedback inhibition. PLoS computational biology, 22(5), e1013463. https://
BibTeX
@article{coroneloliveros
author = {Coronel-Oliveros, Carlos and Lehue, Fernando and Herzog, Rubén and Mindlin, Iván and Gatica, Marilyn and Kowalczyk-Grębska, Natalia and Medel, Vicente and Cruzat, Josephine and Gonzalez-Gomez, Raul and Hernandez, Hernán and Tagliazucchi, Enzo and Prado, Pavel and Orio, Patricio and Ibáñez, Agustín},
title = {{A multi-frequency whole-brain neural mass model with homeostatic feedback inhibition}},
journal = {PLoS computational biology},
year = {2026},
month = may,
volume = {22},
number = {5},
pages = {e1013463},
publisher = {PLOS},
issn = {1553-734X},
doi = {10.1371/
url = {https://
pmid = {42127149},
pmcid = {PMC13183287}
}
RIS
TY - JOUR
AU - Coronel-Oliveros, Carlos
AU - Lehue, Fernando
AU - Herzog, Rubén
AU - Mindlin, Iván
AU - Gatica, Marilyn
AU - Kowalczyk-Grębska, Natalia
AU - Medel, Vicente
AU - Cruzat, Josephine
AU - Gonzalez-Gomez, Raul
AU - Hernandez, Hernán
AU - Tagliazucchi, Enzo
AU - Prado, Pavel
AU - Orio, Patricio
AU - Ibáñez, Agustín
TI - A multi-frequency whole-brain neural mass model with homeostatic feedback inhibition
T2 - PLoS computational biology
J2 - PLoS Comput Biol
PY - 2026
DA - 2026/
VL - 22
IS - 5
SP - e1013463
SN - 1553-734X
PB - PLOS
DO - 10.1371/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1371/
"type": "article-journal",
"title": "A multi-frequency whole-brain neural mass model with homeostatic feedback inhibition",
"container-title": "PLoS computational biology",
"author": [
{
"family": "Coronel-Oliveros",
"given": "Carlos"
},
{
"family": "Lehue",
"given": "Fernando"
},
{
"family": "Herzog",
"given": "Rubén"
},
{
"family": "Mindlin",
"given": "Iván"
},
{
"family": "Gatica",
"given": "Marilyn"
},
{
"family": "Kowalczyk-Grębska",
"given": "Natalia"
},
{
"family": "Medel",
"given": "Vicente"
},
{
"family": "Cruzat",
"given": "Josephine"
},
{
"family": "Gonzalez-Gomez",
"given": "Raul"
},
{
"family": "Hernandez",
"given": "Hernán"
},
{
"family": "Tagliazucchi",
"given": "Enzo"
},
{
"family": "Prado",
"given": "Pavel"
},
{
"family": "Orio",
"given": "Patricio"
},
{
"family": "Ibáñez",
"given": "Agustín"
}
],
"container-title-short":
"volume": "22",
"issue": "5",
"page": "e1013463",
"DOI": "10.1371/
"PMID": "42127149",
"PMCID": "PMC13183287",
"ISSN": "1553-734X",
"publisher": "PLOS",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
5,
13
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.1038/s42003-026-10205-z [code]
- Source-space EEG alpha activity reveals brain age gaps due to neurodegeneration and disparity.Journal: Communications biologyIn common: Nilearn, scikit-learn, SciPy, 2 other tools, EEG, 3 references, author Agustin Ibanez
- [2] doi:10.1073/pnas.2532072123 [code]
- Spatially structured heterogeneity shapes large-scale cortical dynamics in a model of the human cortex.Journal: Proceedings of the National Academy of Sciences of the United States of AmericaIn common: SciPy, Matplotlib, NumPy, 8 references
- [3] doi:10.1002/alz.71232
- Neural embedding of frailty in cognitively unimpaired aging and dementia across Latin America.Journal: Alzheimer's & dementia : the journal of the Alzheimer's AssociationIn common: 3 references, 2 authors
- [4] doi:10.1162/imag.a.1276 [code]
- High-resolution whole-brain magnetic resonance spectroscopic imaging in youth at risk for psychosis.Journal: Imaging neuroscience (Cambridge, Mass.)In common: Brain Connectivity Toolbox, Numba, Nilearn, 7 other tools
- [5] doi:10.1038/s41467-026-74466-2 [code]
- Neuromorphic hierarchical modular reservoirs.Journal: Nature communicationsIn common: Brain Connectivity Toolbox, Numba, Nilearn, 6 other tools, computational
- [6] doi:10.1038/s42003-026-10011-7 [code]
- Learning brain dynamics across distinct scaling regimes reveals psychiatric signatures.Journal: Communications biologyIn common: Numba, Nilearn, NetworkX, 6 other tools, 1 reference
- [7] doi:10.1162/imag.a.1278 [code]
- Young and old adult brains experience opposite effects of acute sleep restriction on the functional connectivity network.Journal: Imaging neuroscience (Cambridge, Mass.)In common: Brain Connectivity Toolbox, Nilearn, NetworkX, 4 other tools, 2 references
- [8] doi:10.1371/journal.pbio.3003916 [code]
- Arousal-driven critical roaming reproduces human functional connectivity dynamics.Journal: PLoS biologyIn common: SciPy, Matplotlib, NumPy, 6 references
- [9] doi:10.1002/ana.78203 [code]
- AI-Driven Mapping of Seizure Spread Patterns.Journal: Annals of neurologyIn common: Brain Connectivity Toolbox, NetworkX, scikit-image, 5 other tools, EEG, 1 reference
- [10] doi:10.1212/wnl.0000000000218472 [code]
- Lesion-Level Subtypes of White Matter Hyperintensity Evolution Beyond Spatial Location.Journal: NeurologyIn common: Nilearn, NiBabel, scikit-learn, 3 other tools, author Raul Gonzalez-Gomez
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 2 repositories of the authors' code, each at its verified commit and with its license, 22 scripts, and 17 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:b836d6343bc8f648…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[, paste the snippet at the top, then “Commit changes…” and, to review it first, “Create a new branch and start a pull request”. You open the pull request; OSCR asks for no permission.
Request its removal
To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).
Discussion, reproductions, activity
Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.
Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
