image thumbnail

burstindex(varargin​)

version 1.0.0.0 (5.27 KB) by Jordan Sorokin
Calculates the burst index of a rebound-bursting neuron from a CCIV

46 Downloads

Updated 15 Feb 2016

View License

% ----------- [BI,Ib,It,Sb,St] = burstindex(varargin) -----------------
%
% Calculates the burst index of a rebound-bursting neuron from a CCIV ranging from
% negative to postivie currents. Assumes that the CCIV protocol sweeps
% through a large enough range of currents to elicit rebound bursts and
% tonic firing.
%
% Burst index is defined as follows:
%
% BI = 1 - (B / (B+T))
% where
% B = abs(Ib)/sqrt(Sb)
% T = abs(It)/sqrt(St)
% and
% Ib = threshold current for eliciting rebound burst
% It = threshold current for eliciting tonic spikes
% Sb = (# spikes in rebound) / (stim length in seconds)
% St = (# of tonic spikes) / (stim length in seconds)
% * dividing by stim length normalizes the spike counts
%
% We take the square root of St and Sb to avoid excessive weighting of spike
% rate in the equation, and instead emphasize the threshold currents.
% This is due to the physiology of rebound bursts, which typically plateau
% at 4-5 spikes despite increasing hyperpolarizing current, while tonic
% spikes plateau at much higher counts, and depend on the stim length
%
%
% >>> INPUTS >>>
% Optional:
% stimon = on time of pulse from CCIV protocol in seconds
% stimoff = off time of pulse from CCIV protocol in seconds
% Spiketimes = spiketimes from CCIV in column format.
% units of spiketimes needs to equal units of stimon/stimoff
% current = current levels used in CCIV protocol
% * if you provide these 4 optional inputs, this function
% will calculate the parametrs Ib, It, Sb, and St for you.
%
% If you already have [Ib,It,Sb,St] parameters calculated,
% leave the 1st 4 inputs blank and proceed with the following four inputs.
%
% Ib = threshold burst current
% It = threshold tonic current
% Sb = # burst spikes / length of pulse (seconds)
% St = # tonic spikes / length of pulse (seconds)
%
% <<< OUTPUTS <<<
% BI = burst index, ranging between [0,1]. Higher values = more "bursty"
% Ib = threshold burst current (optional)
% It = threshold tonic current (optional)
% Sb = normalized # rebound spikes (optional)
% St = normalized # tonic spikes (optional)
%
% By JMS, 02/12/2016
%
%-------------------------------------------------------------------------------------

Cite As

Jordan Sorokin (2021). burstindex(varargin) (https://www.mathworks.com/matlabcentral/fileexchange/55411-burstindex-varargin), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2013a
Compatible with any release
Platform Compatibility
Windows macOS Linux
Tags Add Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!