This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
analysis:course:week10 [2013/11/24 11:52] mvdm [Visualizing the results] |
analysis:course:week10 [2018/07/07 10:19] (current) |
||
---|---|---|---|
Line 95: | Line 95: | ||
</code> | </code> | ||
- | Next, we plot the spikes of a single cell (cell 5, which we met in the previous module) where the rat was when each spike was emitted: | + | Next, we plot the spikes of a single cell (cell 7, which we met in the previous module) where the rat was when each spike was emitted: |
<code matlab> | <code matlab> | ||
% get x and y coordinate for times of spike | % get x and y coordinate for times of spike | ||
- | iC = 5; | + | iC = 7; |
spk_x = interp1(Range(sd.x),Data(sd.x),Data(sd.S{iC}),'linear'); | spk_x = interp1(Range(sd.x),Data(sd.x),Data(sd.S{iC}),'linear'); | ||
spk_y = interp1(Range(sd.y),Data(sd.y),Data(sd.S{iC}),'linear'); | spk_y = interp1(Range(sd.y),Data(sd.y),Data(sd.S{iC}),'linear'); | ||
Line 124: | Line 124: | ||
spk_binned = ndhist(cat(1,spk_x',spk_y'),[SET_nxBins; SET_nyBins],[SET_xmin; SET_ymin],[SET_xmax; SET_ymax]); | spk_binned = ndhist(cat(1,spk_x',spk_y'),[SET_nxBins; SET_nyBins],[SET_xmin; SET_ymin],[SET_xmax; SET_ymax]); | ||
- | imagesc(spk_binned') | + | imagesc(spk_binned'); |
- | axis xy; colorbar | + | axis xy; colorbar; |
</code> | </code> | ||
Line 143: | Line 143: | ||
tc = spk_binned./(occ_binned .* (1./VT_Fs)); % firing rate is spike count divided by time | tc = spk_binned./(occ_binned .* (1./VT_Fs)); % firing rate is spike count divided by time | ||
- | pcolor(tc'); shading flat | + | pcolor(tc'); shading flat; |
- | axis xy; colorbar; axis off | + | axis xy; colorbar; axis off; |
</code> | </code> | ||
Line 176: | Line 176: | ||
figure; | figure; | ||
pcolor(tc'); shading flat; axis off | pcolor(tc'); shading flat; axis off | ||
- | axis xy; colorbar | + | axis xy; colorbar; |
</code> | </code> | ||
Line 253: | Line 253: | ||
tc(isinf(tc)) = NaN; | tc(isinf(tc)) = NaN; | ||
- | pcolor(tc'); shading flat | + | pcolor(tc'); shading flat; |
- | axis xy; colorbar; axis off | + | axis xy; colorbar; axis off; |
</code> | </code> | ||