Code to calc. spike-trig. ave (STA) conduct. from Vm (Pospischil et al. 2007, Rudolph et al. 2007)

 Download zip file 
Help downloading and running models
Accession:116086
PYTHON code to calculate spike-triggered average (STA) conductances from intracellular recordings, according to the method published by Pospischil et al., J Neurophysiol, 2007. The method consists of a maximum likelihood estimate of the conductance STA, from the voltage STA (which is calculated from the data). The method was tested using models and dynamic-clamp experiments; for details, see the original publication (Pospischil et al., 2007). The first application of this method to experimental data was from intracellular recordings in awake cat cerebral cortex (Rudolph et al., 2007).
References:
1 . Pospischil M, Piwkowska Z, Rudolph M, Bal T, Destexhe A (2007) Calculating event-triggered average synaptic conductances from the membrane potential. J Neurophysiol 97:2544-52 [PubMed]
2 . Rudolph M, Pospischil M, Timofeev I, Destexhe A (2007) Inhibition determines membrane potential dynamics and controls action potential generation in awake and sleeping cat cortex. J Neurosci 27:5280-90 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Synapse;
Brain Region(s)/Organism:
Cell Type(s): Abstract integrate-and-fire leaky neuron;
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: Python;
Model Concept(s): Methods; Sleep; Activity Patterns;
Implementer(s): Destexhe, Alain [Destexhe at iaf.cnrs-gif.fr]; Pospischil, Martin ;
  
/
demo_STA
readme.txt
ge_sta.txt
gi_sta.txt
header.py
main.py
methods.py
vm_sta.txt
                            

############################
#       parameters         #
############################

vmFile      = 'vm_sta.txt'  # file containing the voltage STA

dt          = 0.048         # time step in ms
t_cut       = 1.0           # time cut before spike in ms
n_cut       = int(t_cut/dt) 
n_smooth    = 0             # SD for smoothing in timesteps

Iext    = -0.2              # DC in nA
ge      = 0.032             # mean exc. cond. in uS
gi      = 0.096             # mean inh. cond. in uS
se      = 0.012             # exc. cond. fluctuations in uS
si      = 0.0072            # inh. cond. fluctuations in uS
gl      = 0.015             # leak cond. in uS
vl      = -92.66            # leak reversal potential in mV
cap     = 0.26              # capacitance in nF
te      = 2.728             # exc. correlation time in ms
ti      = 10.49             # inh. correlation time in ms
ve      = 0.                # exc. reversal potential in mV
vi      = -75.              # inh. reversal potential in mV

Loading data, please wait...