[FieldTrip] Using PPC method
awutz at mit.edu
Mon Aug 28 14:49:54 CEST 2017
the spike timestamps in the fields ".timestamp" and ."time" should correspond to all trials that occur within the respective trial periods. In the spike-field coherence tutorial, if you compare the data structure for the entire recording (there called "spike") with the trial-cut data structure (there called "spikeTrials"), you can see that there are less spikes in the latter. This is because those trials, which are not within the specified trial periods, are left out.
The tutorial worked great for me. Unfortunately, it stops at the point where inferential stats have to be computed. I understand that I have to break up the relationship between spike timestamp per trial with lfp phase per trial. This will require random resampling of the lfp-phases. Intuitively, I would shuffle the trial vector for the lfp and then calculate the ppc for each random partition (using the same number of spikes). Let me add my own questions here.
1) It is unclear to me what happens when I have more than one spike per trial (e.g. 3 spikes in trial 1). Should I construct the resampling distribution with as many draws from each trial as there where in the original data (e.g. 3 draws from trial 1 but at different time points)?
2) I wonder whether only the trial vector should be shuffled (preserving the time points per trial compared to the original effect but now from a different, random trial). In other words, if I had in the original data 3 spikes which happened in the trials [1 2 3] and at the times [-.5 0 1], should the random partition look like this:
partition: trial = [3 1 2], time = [-.5 0 1]
or should the time points be taken completely at random? I hope, I was able to be clear enough with my explanation.
From: fieldtrip-bounces at science.ru.nl [fieldtrip-bounces at science.ru.nl] on behalf of 马天阳 [tianyangma2013 at gmail.com]
Sent: Monday, August 28, 2017 5:05 AM
To: fieldtrip at science.ru.nl
Subject: Re: [FieldTrip] Using PPC method
I am also a new learner using PPC method. In your reply to Julia, I don’t quite understand why there are same number of values in timestamp and time, since timestamps is the spike times from the beginning of the recording to the end and time is relative to trigger in each trial.
In my case, I have timestamp of spikes and LFP data 0.5 s before and after stimulus in each trial . Do you think it is possible I could use the code by somehow manipulating the format of input?
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the fieldtrip