7
$\begingroup$

I want to characterize a binary time-series signal x (derived from neuron action potential data) in the frequency domain.

Should I use the FFT of the original signal x, or the power spectrum (FFT of the autocorrelation of x, which is defined as $k(n)=E[x(t)\cdot x(t-n)]$ )?

Also, does anyone know any literature which explicitly talks about this, or at least took sides?

$\endgroup$
9
  • 4
    $\begingroup$ What is your ultimate objective with "characterizing" the data? Definitely read Spikes when you have a chance. The short answer from there is that it depends on what you are looking to get out of the signal. If you're looking to drive a device or something of that nature, the FFT is going to be computationally burdensome, and may not give you sufficient information. $\endgroup$ May 28, 2014 at 21:43
  • 3
    $\begingroup$ As an aside, for papers to look at, George Gertstein is the guru of the analysis of spike trains. $\endgroup$ May 28, 2014 at 21:43
  • 3
    $\begingroup$ I'm using "information" as a loaded term. I was referring to both measures that you are talking about, sorry if that wasn't clear. Neither is going to tell you anything meaningful. If you're locking to a specific event (sensory input or specific motor output - e.g., an EMG of a muscle) the frequency variation will tell you something, but a "free form" recording really isn't. You'd be better off examining local field potentials (most extracellular recording systems can separate these out in the spike sorting processing) if you are trying to tease out frequency changes. $\endgroup$ May 29, 2014 at 20:02
  • 2
    $\begingroup$ Otherwise, this is really just an exercise in signal processing, to be frank. $\endgroup$ May 29, 2014 at 20:03
  • 2
    $\begingroup$ (not trying to be an overbearing curmudgeon here, just offering some constructive criticism. I'm happy to see a methods question on here!) $\endgroup$ May 29, 2014 at 21:00

2 Answers 2

2
$\begingroup$

For binary data (spike-trains), people use a spike-train autocovariance function (same as autocorrelation but subtracting the mean firing rate). It's commonly used to study brain oscillations (e.g. Engel & Konig, 1991). There is a good description of the math in the Dayan & Abbott (page 27).

Dayan, P., & Abbott, L. F. (2001). Theoretical neuroscience (Vol. 806). Cambridge, MA: MIT Press. Engel, A. K., & Konig, P. (1991). Interhemispheric synchronization of oscillatory neuronal responses in cat visual cortex. Science, 252(5009), 1177.

$\endgroup$
1
$\begingroup$

As this answer explains, the Fourier transform of the autocorrelation of $x$ is the same as the magnitude squared of the Fourier transform of $x$. The autocorrelation of x can be expressed as $x(t)*x(-t)$, where * means convolution and the Fourier transform of the autocorrelation is then $X(w)X^*(w)=|X(w)^2|$, which of course is the Fourier transform of the original signal $x(t)$.

As for references, see this answer.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.