Neural transformations on spike timing information (Tripp and Eliasmith 2007)

 Download zip file 
Help downloading and running models
" ... Here we employ computational methods to show that an ensemble of neurons firing at a constant mean rate can induce arbitrarily chosen temporal current patterns in postsynaptic cells. ..."
1 . Tripp B, Eliasmith C (2007) Neural populations can induce reliable postsynaptic currents without observable spike rate changes or precise spike timing. Cereb Cortex 17:1830-40 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network;
Brain Region(s)/Organism:
Cell Type(s): Abstract integrate-and-fire leaky neuron;
Gap Junctions:
Simulation Environment: MATLAB;
Model Concept(s): Activity Patterns; Temporal Pattern Generation;
Implementer(s): Tripp, Bryan [bryan.tripp at]; Eliasmith, Chris [celiasmith at];
% function plotCorr(datafile, jitterIndex) plots a summary of MSE vs
% correlation. 
% dataFile: filename of saved data from a exp_correlation.m
% jitterIndex: index of jitter level in this file 

function plotCorr(datafile, jitterIndex)
    data = load(datafile); 
    hold on
    plotData(data, data.epochRateIndices, 'k-', jitterIndex);
    plotData(data, data.thresholdIndices10, 'k:', jitterIndex);
    plotData(data, data.thresholdIndices22, 'k:', jitterIndex);
    plotData(data, data.thresholdIndices55, 'k:', jitterIndex);
    xlabel('Pairwise correlation')
    addLabel(data, data.thresholdIndices10(1), jitterIndex, 'a');
    addLabel(data, data.thresholdIndices22(1), jitterIndex, 'b');
    addLabel(data, data.thresholdIndices55(1), jitterIndex, 'g');

function addLabel(data, caseIndex, jitterIndex, labelText)
    h = text(data.meanPeakCorrelation(caseIndex)+.005, data.meanE(caseIndex,jitterIndex), labelText)
    set(h, 'FontName', 'Symbol')  

function plotData(data, caseIndices, symbol, jitterIndex)
    x = data.meanPeakCorrelation(caseIndices);
    y = squeeze(data.meanE(caseIndices,jitterIndex,:));
    meanY = mean(y, 2);
    sdY = std(y, 0, 2);
 	plot(x, meanY, symbol);
	hold on    
	for i = 1:length(x)
        plot([x(i) x(i)], [meanY(i)-sdY(i) meanY(i)+sdY(i)], symbol) 

Loading data, please wait...