Slowly activating outward membrane currents generate input-output sub-harmonic cross frequency coupling in neurons
Description
Data and code files for manuscript published in the Journal of Theoretical Biology (https://www.sciencedirect.com/science/article/abs/pii/S0022519320303647)
Files
Steps to reproduce
The set of code files comprises of two parts: 1) The julia code files (.jl) which perform the simulation and store the .mat files in the 'Na_K_M_x' folders for each step (30%) of increase in M-channel conductance (raw data). Julia packages required for running 'main.jl' are listed in the top of the file in the 'include' command. The code was written in Julia 1.3.1. The file 'HH.jl' contains the neuron model. and 'Data_analysis.jl' contains the code to extract spike train and related data from the output of the julia SDE simulations. 2) The matlab file 'my_plotter.m' has code sub-segments that generates the output results into the folder 'input_vs_spikes_HH'. The first segment of the code titled 'extract data' gets the data from the raw data ('Na_K_M_xx' folders). The second segment, titled 'EPSP convolution' generates the continuous signal for coherence analysis. The third segment 'NMC' generates the coherence plots and data. The final 'Plots' segment generates the plots for the manuscript. For these files, we have already extracted the raw data into 'NI.mat'. Hence, please load this file and run the segments from 'EPSP_convolution' segment onwards of 'my_plotter.m' file. We have not included the raw data in the 'Na_K_M_xx' folders because of space constraints (the total size is around 18GB). One can generate all the raw data from julia files as described above. Please contact the authors to get a link to the raw data in case you do not have the resources to run the Julia codes.