How to do spectral estimation of a point process

5 次查看（过去 30 天）
Bence Laczó
Bence Laczó 2023-4-12
评论： Star Strider 2023-4-12
I would like to analyse neural firing in the time-frequency domain. I have found that several programs (Chronux, nSTAT) exist for this problem, but they are used for analysis of neural data over multiple trials. I would like to analyse data in which I have only one recording from one particular position. I have found that the spectrum of point process can be expressed as the Fourier transform of its autocovariance function. Can somebody tell me which command should be used for it in Matlab?

请先登录，再回答此问题。

回答（2 个）

Star Strider
Star Strider 2023-4-12
For time-frequency analyses, I usually use the pspectrum function with the 'spectrogram' option.
  4 个评论
显示 2更早的评论
Bence Laczó
Bence Laczó 2023-4-12
My data are spike trains. It contain zeros and ones (where there is a spike) or I also have the interspike interval values in milisecs.
Star Strider
Star Strider 2023-4-12
Ran in:
If I remember correctly, that is a Poisson process, so perhaps something like this —
ISI = poissrnd(rand(1,250)*10)*1E-3;
Fs = 1E3;
figure
pspectrum(ISI, Fs, 'spectrogram')
colormap(turbo)
I have never analysed spike trains, so I don’t know the accepted procedures for analysing them.
.

请先登录，再进行评论。

Image Analyst
Image Analyst 2023-4-12
So you have a signal emanating from a point? So you have time series data. I'd just use fft or pwelch. It doesn't matter if your signal came from only one location.
  1 个评论
Bence Laczó
Bence Laczó 2023-4-12
编辑：Bence Laczó 2023-4-12
I have interspike interval data or spike locations which are marked by ones in an array of zeros.
Now I see what was misleading. I have neuronal firing data from different positions in the subthalamic nucleus. There was no specific task during the recording only resting state data.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Spectral Measurements 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by