Data Analysis Flashcards
What is spike sorting?
After monitoring the activity of a large number of neurons simultaneously:
Spike sorting is knowing which neurons fired when
How do we classify single neurons from other neurons and noise
which neuron fired when
How many neurons does EEG capture?
millions
what is the regional domain of EEG?
3-5cm
Two types of recording for single cell?
- Wide band continuous recording - continuous
* Filtered, spike-triggered recordings - action potential crossing threshold would trigger the recording
what is High Pass Filtering?
allowing high freq to pass through and discarding lower freq
why high pass filter?
- Local field (the group activity from 100 and 1000s of neurons) potential is primarily at low frequencies - so get rid of this noise.
- Spikes are at higher frequencies.
- So use a high pass filter (800-1000 Hz cutoff is good)
• Preferably a non-causal filter
(causal filter changes the signal slightly - might produce phase distortions)
how do we tell what is a spike and what is not?
• Mostly used is a voltage threshold
• Check the quality of detection
– Visual inspection of spikes superimposed (Study interspike interval histogram)
then cut an epoch around the peak - to get individual neuron
Aids to understanding what is and what is not one single neuron?
No spike interval less than 1 ms (spike refractory period)
- this condition is only applicable once the spikes are sorted
- but this indicates that if two come in sucession <1ms - has to be two signals
The choice of threshold is a trade-off between
a) too high threshold (Type II error)
b) too low threshold (Type I error)
how many standard deviations away to distinguish spikes?
5
but you can also apply modified using the median of the spikes … is more realistic
What is Feature Analysis?
- Two clear action potentials that have roughly the same height but are different in shape.
- If the shape could be characterized, we could use this information to classify each spike.
How do we characterize the shape?
extract the features
each time a neuron fires, it has the same features - like people with a distinct voice. So we can charcterise them and separate them.
What are the main features?
Peak amplitude
peak-to-peak amplitude
peak width
energy etc.
What are more advanced methods for spike sorting?
Principal components analysis (cousin of factor analysis)
Wavelet transform
How does Principal components analysis work?
The idea behind PCA is to find an ordered set of orthogonal basis vectors that capture the directions in the data of largest variations
How is Wavelet transform better?
Wavelet transform based feature extraction
- localized properties (i.e. spike shape details) are emphasized - More adaptive than PCA
Which feature is used to discriminate spikes?
Select those features that best separate the different clusters of spike shapes.
A wavelet coefficient -or any other feature- that is good for distinguishing different spike shapes should have a multimodal distribution, unless there is only one cluster.
What is Clustering?
Group spikes with similar features into clusters, corresponding to different neurons
What is manual clustering?
Manual clustering
– Drawing polygons in 2-D projections of the spike features
issue with manual clustering?
– Prone to error, subjective bias, time consuming, not practical for high dimensional data
what is an alternative clustering approach?
Nearest neighbor k-means clustering
features of Nearest neighbor k-means clustering?
– Define the cluster locations as the mean of data within that cluster
– Cluster membership is defined by Euclidean distance
– This defines a set of implicit decision boundaries separating the
clusters
– Works well when the clusters are separated but not otherwise
name a Distribution free approach?
Super paramagnetic clustering
main points for Super paramagnetic clustering?
- Groups the spikes into clusters as a function of a single parameter, ‘the temperature’.
- The approach is based on tuning this temperature
- For low temperature, the data are grouped together to a single cluster, and for high temperature, the data are dispersed to too many clusters
- For a middle range temperature, corresponding to super paramagnetic regime, the data are optimally sorted into a few clusters but with large membership
what is Rate Estimation?
how many times a neuron fire in a second
Most widely used estimate of neural activity = Average firing rate over a time interval
Rate Estimation is usually expressed by a….
Peri-Stimulus Time Histogram (PSTH)
How to make a PSTH?
event related neural firing
- Align spike sequences with stimulus onset (or any event) which repeated n times
- Divide the stimulus period S into N bins of size D
• Count the number of spikes (ki) within individual bin for
all trials
• Compute the histogram
what is Inter-Spike Interval (ISI)?
A temporal coding - when it fires , not by how many times.
Calculate individual firing - and the distance between
Neuronal firing pattern is considered to have a Poisson distribution: a renewal process.
traditionally, the ISI sequence is assumed to be random ! assumption behind rate coding ….. if so, the structure would be random ….