Having trouble with the spike triggered average

5 views (last 30 days)
I keep getting empty matrices when trying to determine a spike triggered average. Below is my code. Can anyone help??
spkr=[NaN,spikes(1:end-1)]; trig_spk=find(spikes>0 & spkr<=0); x1=1:length(spikes);
y2=zeros(1,length(trig_spk)); figure;plot (x1,spikes,'g',trig_spk,y2,'.') dog=[]; cat=[]; for i=1:length(trig_spk) try dog(i,:)=binarized( trig_spk(i) : trig_spk(i)+4000-1); cat(i,:)=spikes(trig_spk(i) : trig_spk +4000-1); end; end end;

Answers (4)

Star Strider
Star Strider on 16 Mar 2015
I keep getting empty matrices too, because I don’t have your data.
Care to let us know what your ‘spikes’ data are?

Holly Hunsberger
Holly Hunsberger on 16 Mar 2015
Attached is the data file. Thanks!

Holly Hunsberger
Holly Hunsberger on 16 Mar 2015
Also, the window is 200ms. (200ms before a spike and after)
  2 Comments
Star Strider
Star Strider on 16 Mar 2015
When you post it, will your data also have a time variable?
Image Analyst
Image Analyst on 16 Mar 2015
After you click "Choose file..." remember to click the "Attach file" button or else it won't attach.

Sign in to comment.


Holly Hunsberger
Holly Hunsberger on 17 Mar 2015
Time wasn't in the original data file, but I figured it out to be t=1:length(stimulus); t=t./sr;

Categories

Find more on Electrophysiology in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!