1. a quick "getting started" introduction

Suppose you have run an experiment in which a colorful stimulus was presented inside the receptive field of a visually responsive, isolated neuron. The experiment consisted of 4 different conditions, characterized by the color of the stimulus:

cnd #1: "red"
cnd #2: "green"
cnd #3: "blue"
cnd #4: "yellow"

At the end of the run, the data file collected by CORTEX contains a few hundred trials of randomly interleaved conditions. You now want to know the average response rate of the recorded neuron to the "blue" stimulus in a 200-ms time window starting 50 ms after the stimulus onset. Here's how to do that.

First you must create a @cortex object from the source data file, i.e., the file CORTEX created during the experiment:

>> c = cortex ('datafilename', 1:2);

The range [1:2] tells the function that this data file contains two trains of spike data encoded with the numbers 1 and 2.
Next, you create a @key object that will allow quick selection of an epoch and of trials that "belong together" in the experiment

>> k = key;
>> k.conds = 3;
>> k.start = 50;
>> k.finish = 250;
>> k.cue = 25;

the k @key object will isolate a time window starting and ending, respectively, 50 ms and 250 ms after the occurrence of event #25 in all trials of condition #3. You can now put things together:

>> mfr = meanfiringrate (c, 1, k)

which returns the average firing rate of spike #1 over said epoch for all trials that match the criteria specified by k.

Previous <-> Next