Insert zeros in a vector - Random

Hi guys,

So I have this vector with 182 elements that show me timings of a neuron's spiking activity. Each element corresponds to a time point, for e.g. 672 where 672 is in ms.

My timeline is 10s or 10000 ms.

What I need to generate is a vector with 10000 elements, where there are 1s at the times corresponding to the spike-timing matrix and 0s everywhere else.

The code I wrote shows me only one 1, that too at the last time point. Here is what I wrote -

spikes = zeros(1,length(time));
s = SpikeTimes*10;
for j = 1:length(SpikeTimes)
  for i = 1:length(spikes)
      if s(j) ~= i
          spikes(i) = 0;
      else
          spikes(i) = 1;
      end
  end
end

Could anyone please help me out? Thanks.

I need to further use the generated spikes vector to calculate firing rate using binning, a rectangular window, a gaussian window and compare and discuss the results.

Answer 1

[Matt Tearle]

Perhaps I misunderstand, but you should be able to do simply:

spikes(SpikeTimes*10) = 1;

Example:

x = [1 2 7 14 17]
y = zeros(1,20);
y(x) = 1

[Abhilash]

Nope, that didn' work. What actually worked was taking away the 'else' clause. Weird.

[Abhilash]

Oh not weird. I get why taking away the else clause worked.

[Matt Tearle]

What does SpikeTimes look like? Is this the 182-element vector of timings (672 etc)?

[Abhilash]

Yes. It has timings starting at 672 and ending at 8606

[Matt Tearle]

OK so:

spikes = zeros(...)

spikes(SpikeTimes) = 1;

[Abhilash]

Thanks

Answer 2

[Dimithra Anton]

spikes = zeros(1,length(time));

spikes = zeros(1,length(time))

[Adam Danz]

spike = zeros(size(time));

to maintain the same shape.

Also avoid using length() since you never know what dimension it will return if your data aren't a vector. Instead, use

spike = zeros(1,numel(time));

[Abhilash]

Wow this was 9 years ago!