Code for generating Figures 3, 4, 5, & 7 of the following paper:

Osinski & Kay (2016), Granule cell excitability regulates gamma and
beta oscillations in a model of the olfactory bulb dendrodendritic
microcircuit, J Neurophysiol.
doi: 10.1152/jn.00988.2015

Code is designed to run in MATLAB R2015b

Implementation is by Boleslaw Osinski , to whom questions should be
addressed boleszek@uchicago.edu

---------------------------------------------------------------------
The zipped folder includes the following .m files:

InitNetwork_GCE.m - This function initializes network parameters
OB_network_GCE.m  - This is the numerical simulation of the
                    Mitral-Granule network
ILFP_GCE .m       - This is a wrapper function for
                    OB_netwrok_GCE.m. All simulation products as well
                    as the simulated LFP are generated by this
                    function
ParamSweep_GCE.m  - This is a wrapper function around ILFP_GCE.m. It
                    evaluates the model for every combination of the
                    two parameter inputs
Rasterplot.m      - Generates raster plot
---------------------------------------------------------------------


Usage:
1. Unzip GCE_ModelDB.zip
2. Open MATLAB. Set MATLAB path to unzipped folder's directory. This
can be done with MATLAB prompt commands like (where additional folders
in the cd command may be necessary to find the new GCE_ModelDB folder)

cd GCE_ModelDB
addpath(pwd)

3. Open Fig3, Fig4, Fig5, or Fig7 directory and run the associated
script (i.e. Fig3script.m) for producing the figures. Each folder
includes it's own parameter file (i.e. OB_params_GCE_Fig3.txt)

Notes:

- Fig3script.m runs quickly and produces figures like those in the
paper. I remade the figures for the Methods section (including Figure
3) at the very beginning of the revision process, and there were
subsequent minor changes to the final model (such as adjusting
changing the proportionality between Ca concentration and N-type
current and [Ca] threshold for GABA release). The slightly different
trajectories in Figure 3 resulting from these changes do not change
the network behavior of the model.:

Fig 3 A


Fig 3 B


Fig 3 C


Fig 3 D


- Fig4script.m runs for ~ 4min
Fig4script.m uses the tightfig function which can be found at the
MATLAB file exchange at
http://www.mathworks.com/matlabcentral/fileexchange/34055-tightfig

Fig 4 A.i


- Fig5script.m runs for ~ 3.5hr
- Fig5script.m uses the shadedErrorBar function which can be found at
the MATLAB file exchange at
http://www.mathworks.com/matlabcentral/fileexchange/26311-shadederrorbar
and the coherencycpt function from the Chronux toolbox (Version 2.11)
which can be downloaded at http://chronux.org

- Fig7script.m runs for ~ 1.5hr