There is a newer version of the record available.

Published July 6, 2026 | Version v2

Pretrained Models and Feature Datasets for ncpi

  • 1. ROR icon Universidad de Granada

Description

This dataset contains derived electrophysiological features, trained machine-learning inverse models, test-set predictions, and evaluation metrics generated from large-scale neural circuit simulations.

The resources were produced using the ncpi Python package to support neural circuit parameter inference from simulated electrophysiological signals.

The original raw simulation signals are not included in this deposit. Instead, this record contains derived feature matrices, reusable trained inference models, test-set outputs, and model-evaluation summaries.

Important note

The included features and machine-learning models have not yet been experimentally tested or independently validated. They are provided as research resources for reproducibility, benchmarking, and further methodological development.

Simulation datasets

Results are provided for three simulation datasets.

Hagen_v1

Current-dipole-moment signals generated using the first version of the Hagen neural circuit model. This model uses current-based synapses.

The Hagen_v1 dataset comprises approximately 2 million simulations generated by sampling neural circuit parameters within biologically plausible ranges. These ranges were constrained a priori according to biological principles. Hagen_v1 is processed in streaming batch mode using six simulation batches.

Model description: https://doi.org/10.1371/journal.pcbi.1010353

Hagen_v2

Current-dipole-moment signals generated using the second version of the Hagen neural circuit model. As in Hagen_v1, this model uses current-based synapses.

The Hagen_v2 dataset is a more constrained simulation set, comprising approximately 100,000 simulations. In addition to tighter parameter-range constraints, simulations were further filtered post hoc based on firing rates, spike-based metrics, and power-spectrum properties. Hagen_v2 is loaded by concatenating the first-search and second-search simulation sets.

Model description: https://doi.org/10.1371/journal.pcbi.1010353

Cavallari

Proxy electrophysiological signals generated using the Cavallari neural circuit model. In contrast to the Hagen datasets, the Cavallari model uses conductance-based synapses.

Model description: https://doi.org/10.3389/fncir.2014.00012

All signals were processed assuming a sampling frequency of 1600 Hz.

Inferred neural circuit parameters

Separate inverse models were trained for the Hagen and Cavallari simulation datasets. All inverse models estimate two target variables:

  • E_I

  • J_ext

For the Cavallari dataset, J_ext corresponds to the external-input scaling parameter ext_input_scale.

Hagen model parameters

Each Hagen simulation directly varies seven neural circuit parameters:

Parameter Description Range
J_EE Strength of recurrent excitatory-to-excitatory connections 0.5–4.0
J_IE Strength of excitatory-to-inhibitory connections 0.5–4.0
J_EI Strength of inhibitory-to-excitatory connections −40.0 to −1.0
J_II Strength of recurrent inhibitory-to-inhibitory connections −40.0 to −1.0
tau_syn_E Excitatory synaptic time constant applied to excitatory inputs received by both neuronal populations 0.1–2.0 ms
tau_syn_I Inhibitory synaptic time constant applied to inhibitory inputs received by both neuronal populations 0.1–8.0 ms
J_ext Synaptic weight of the external Poisson input received by the network 10.0–50.0

The four recurrent coupling strengths are combined into the inferred excitatory-to-inhibitory target:

E_I = (J_EE / J_EI) / (J_IE / J_II)

Therefore, each Hagen inverse model estimates:

  • E_I

  • J_ext

The individual recurrent coupling strengths and synaptic time constants are varied during simulation but are not estimated separately by the inverse models.

Cavallari model parameters

Each Cavallari simulation modifies seven parameters relative to the default Cavallari configuration available in the ncpi repository:

Parameter Description Range
g_EE Directly replaces exc_exc_recurrent, the excitatory-to-excitatory recurrent synaptic conductance 0.5–2.0 × default value of 0.178
g_IE Directly replaces exc_inh_recurrent, the excitatory-to-inhibitory recurrent synaptic conductance 0.5–2.0 × default value of 0.233
g_EI Directly replaces inh_exc_recurrent, the inhibitory-to-excitatory recurrent synaptic conductance 0.5–2.0 × default value of −2.01
g_II Directly replaces inh_inh_recurrent, the inhibitory-to-inhibitory recurrent synaptic conductance 0.5–2.0 × default value of −2.70
tau_syn_AMPA_scale Dimensionless scaling factor applied to tau_decay_AMPA in both excitatory and inhibitory neurons 0.5–2.0
tau_syn_GABA_scale Dimensionless scaling factor applied to tau_decay_GABA_A in both excitatory and inhibitory neurons 0.5–2.0
ext_input_scale Dimensionless scaling factor applied to both thalamic external-input conductances, th_exc_external and th_inh_external 0.5–4.0

The default AMPA decay constants are 2.0 ms for excitatory neurons and 1.0 ms for inhibitory neurons. The default GABA-A decay constant is 5.0 ms for both populations. The default values of the thalamic external-input conductances are 0.234 and 0.317, respectively. Cortico-cortical external inputs are not modified.

The four recurrent conductances are combined into the inferred excitatory-to-inhibitory target:

E_I = (g_EE / g_EI) / (g_IE / g_II)

Therefore, each Cavallari inverse model estimates:

  • E_I

  • J_ext, corresponding to ext_input_scale

The individual recurrent conductances, synaptic time-constant scaling factors, and external-input components are varied during simulation but are not estimated separately by the inverse models.

Feature configurations

Three feature configurations are included:

Feature set Description
catch22_22 Complete set of 22 canonical catch22 time-series features
specparam_3 Three spectral-parameterization features: aperiodic slope, dominant peak frequency, and dominant peak power
catch22_specparam_25 Concatenation of the 22 catch22 features and the 3 specparam-derived features

Spectral parameterization was performed over 5–200 Hz, with a goodness-of-fit threshold of R² ≥ 0.9. Invalid or non-finite feature rows were excluded before model training.

Trained inverse models

Four inverse-model types were trained for each dataset and feature configuration:

Model Description
Ridge regression (Ridge) Linear regression model with L2 regularization. The regularization parameter alpha was selected according to feature dimensionality.
Multilayer perceptron regression (MLPRegressor) Feedforward neural-network regressor using ReLU activation, the Adam optimizer, early stopping, adaptive hidden-layer sizes according to feature dimensionality, and a maximum of 1000 iterations.
Random forest regression (RandomForestRegressor) Ensemble tree-based regressor using bootstrap aggregation, with adaptive tree number, depth, leaf-size, feature-subsampling, and sample-subsampling settings according to feature dimensionality.
Neural posterior estimation (NPE) Simulation-based inference model using a neural spline flow estimator, with adaptive hidden features, number of transforms, and batch size according to feature dimensionality. Models were trained for up to 100 epochs.

This produces up to 36 trained dataset–feature–model combinations:

3 datasets × 3 feature configurations × 4 model types

Models were trained using a reproducible random 85% training / 15% testing split with random seed 0.

For ridge regression, multilayer perceptron regression, and random forest regression, test-set predictions and mean squared error metrics were computed on the held-out test set.

For neural posterior estimation, the posterior estimator and related inference assets were trained and saved, but held-out test-set prediction and evaluation were skipped by design to avoid expensive posterior sampling.

Included files

The archive contains two main directories:

new_features
new_models

new_features

This directory includes:

  • Computed feature matrices in NumPy format

  • Valid-row masks

  • Reproducible training and testing indices

  • Cached feature matrices generated by the training workflow

  • Precomputed per-batch Hagen_v1 feature and target files produced in streaming mode

new_models

This directory includes:

  • Serialized trained models and feature scalers

  • NPE inference, density-estimator, and posterior assets

  • Model configuration files

  • Test target arrays, y_test.npy, for evaluated scikit-learn models

  • Model predictions, y_pred.npy, for evaluated scikit-learn models

  • Per-parameter and mean squared error metrics, metrics.json

  • NPE metrics files indicating that test evaluation was skipped

  • Dataset-level result summaries, summary.json

Software

The processing and training workflow was implemented using ncpi and the script:

run_massive_training_v3.py

The ncpi source code is available at:

https://github.com/necolab-ugr/ncpi

The script used to generate these results is also included in this deposit.

Recommended use

This dataset is intended for:

  • Reproducibility of the ncpi inverse-model training workflow

  • Benchmarking feature-based and simulation-based inference methods

  • Reuse of trained inverse models for methodological development

  • Evaluation of derived electrophysiological feature representations

  • Further validation against experimental EEG, MEG, or related electrophysiological datasets

Limitations

The raw simulation signals are not included in this record. The deposit contains derived features, trained models, predictions, metrics, and workflow assets.

The included inverse models should be considered research tools. They have not yet been experimentally validated and should not be interpreted as clinically or biologically validated estimators without further independent testing.

Files

Files (770.7 MB)

Name Size
md5:539e2eacbbc1dd914c3ee0ad8a1f5d4e
770.6 MB Download
md5:de6436be0c724d3714e2fc543cbeb90b
103 Bytes Download
md5:07ba19078f2a97db4a6308fdf678b881
64.5 kB Download

Additional details

Funding

Agencia Estatal de Investigación
PID2022-139055OA-I00 (MCIN/AEI/10.13039/501100011033, ERDF)
Agencia Estatal de Investigación
PID2022-137461NB-C31 (MCIN/AEI/10.13039/501100011033, ERDF)
Agencia Estatal de Investigación
RYC2024-049595-I (MCIN/AEI/10.13039/501100011033, FSE+)

Software