Model of the Mammalian ET cell

Computational model and files to recreate the results from the paper "A Computational Model of the Mammalian External Tufted Cell" https://doi.org/10.1016/j.jtbi.2018.10.003

Author: Ryan Viertel

Usage:

data = ET(input,sampling_rate);

input: input vector, if no input then just a vector of zeros

sampling_rate: rate at which the input vector should be sampled. 1000 for milisecond

returns the following struct:

data.T - time vector data.X - ODE variables at each time step

data.events - list of spike events

data.which - event type

data.current - system currents

example

create the input vector

input = zeros(1,5000);

run the model

data = ET(input,1000);

plot the voltage trace

plot(data.T,data.X(:,1))

ME-PCM

The code used to sample the model throughout parameter space to determine stability and investigate the effect of model parameters on model output is found in the ME-PCM directory

xpp

The ODE file to recreate the bifurcation diagram is found in the xpp directory