Oversampling method to extract excitatory and inhibitory conductances (Bedard et al. 2012)

 Download zip file 
Help downloading and running models
" ... We present here a new method that allows extracting estimates of the full time course of excitatory and inhibitory conductances from single-trial Vm recordings. This method is based on oversampling of the Vm . We test the method numerically using models of increasing complexity. Finally, the method is evaluated using controlled conductance injection in cortical neurons in vitro using the dynamic-clamp technique. ..."
1 . Bédard C, Béhuret S, Deleuze C, Bal T, Destexhe A (2012) Oversampling method to extract excitatory and inhibitory conductances from single-trial membrane potential recordings. J Neurosci Methods 210:3-14 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s):
Gap Junctions:
Simulation Environment: MATLAB;
Model Concept(s): Parameter Fitting; Methods;
Implementer(s): Bedard, Claude [bedard at unic.cnrs-gif.fr];
%  generates the oversampled voltage from known conductances

function [GGG PPP] = signalR(oversampling,gge,ggi,tt,PPP)

                   % transformation of conductances (first step)
                 for n=1:1:length(tt)-1    
                       for m=oversampling-1:-1:0
                            tta(oversampling*n-m) = tt(n) + (oversampling-m-1)*(tt(n+1)-tt(n))/oversampling;
                 % compute preconductances : (second step)

                   condinitiale =  PPP.Initial; %condition initial of tension
                   % computes preconductances normals               
                  ggl(1,1:length(ggi)) = PPP.gl;
                  ggs1 = -(gge+ggi+ggl)/PPP.cap; 
                  ggs2 = (PPP.Ee*gge +PPP.Ei*ggi +PPP.El*ggl)/PPP.cap;
                   % compute voltage normal             
                 [tt vv] = eqdifferentielleR(tt,ggs1,ggs2,condinitiale);
                   % computes preconductances for oversampled voltage
                  ggla(1,1:length(ggia)) = PPP.gl;
                  ggs1a = -(ggea+ggia+ggla)/PPP.cap; 
                  ggs2a = (PPP.Ee*ggea +PPP.Ei*ggia +PPP.El*ggla)/PPP.cap;
                   % compute oversampled voltage
                 [tta vva] = eqdifferentielleR(tta,ggs1a,ggs2a,condinitiale);

                 % structure
                   GGG.tt=tt; GGG.vv= vv;  % voltage, normal
                   GGG.tta=tta; GGG.vva= vva;  % voltage, oversampled
                   GGG.gge=gge; GGG.ggi=ggi;  % inhibitors and excitators, originals
                   GGG.ggea=ggea; GGG.ggia=ggia;  % inhibitors and excitators, with oversampling
                   GGG.ggs1=ggs1; GGG.ggs2=ggs2; % preconductances, normal
                   GGG.ggs1a=ggs1a; GGG.ggs2a=ggs2a; % preconductances, oversampled

Loading data, please wait...