How to do spectral estimation of a point process

Question

Bence Laczó 2023-4-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1945804-how-to-do-spectral-estimation-of-a-point-process

评论： 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?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Star Strider 2023-4-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1945804-how-to-do-spectral-estimation-of-a-point-process#answer_1214374

For time-frequency analyses, I usually use the pspectrum function with the 'spectrogram' option.

4 个评论
显示 2更早的评论隐藏 2更早的评论

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 2023-4-12

在 MATLAB Online 中打开

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.

.

请先登录，再进行评论。

Answer 2

Image Analyst 2023-4-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1945804-how-to-do-spectral-estimation-of-a-point-process#answer_1214459

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 个评论
显示 -1更早的评论隐藏 -1更早的评论

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.

请先登录，再进行评论。