Cross-frequency coordination in a hippocampo–cortical circuit during probabilistic reversal learning

This repository contains the complete set of custom MATLAB scripts and functions used for processing behavior, neural data, and generating the findings presented in the manuscript: "Cross-frequency routing in a hippocampo-cortical circuit during probabilistic reversal learning" (Aguilera et al., Communications Biology). The scripts are designed to analyze simultaneous local field potential (LFP) and multiunit spiking data recorded across a distributed circuit encompassing the dorsal hippocampus (CA1d), lateral entorhinal cortex (LEC), and medial prefrontal cortex (PFC) in rats performing a two-armed bandit task. Repository Structure & File Contents The repository is organized into four main compressed directories to streamline navigation and pipeline replication: 1. Calculations.zip (Core Analysis Functions) Contains all foundational custom functions developed to process raw time series, handle data segmentation across behavioral epochs, and execute specialized neurophysiological computations. 2. MainFigs.zip (Main Text Figure Generation) Contains the master execution scripts dedicated to reproducing the main display items of the publication. 3. SupplementaryFigs.zip (Extended Data Figure Generation) Includes the scripts required to generate all supplementary materials (e.g., individual behavioral metrics, control covariates, validation plots of surrogate metrics, and histology cross-references). 4. Tables.zip (Statistical Modeling Output Scripts) Contains scripts setting up the Generalized Linear Mixed-Effects Models (GLMMs) used throughout the text (binomial and linear formulations including animal random intercepts and slopes) to evaluate performance markers while controlling for behavioral strategy and training. System Requirements & Software Variables Primary Platform: MATLAB (Code verified and maintained on version R2018a or newer). Toolbox Dependencies: Signal Processing Toolbox, Statistics and Machine Learning Toolbox, and Curve Fitting Toolbox. Specialized secondary computations utilize components from the FieldTrip architecture and the CircStat (Circular Statistics) toolbox. Detailed information on parameter initializations, specific variables, and temporal padding limits used during filtering are documented as comments within individual master routines.