|
rmeddis@0
|
1 function PSTH =UTIL_makePSTH(inputData, dt, PSTHbinWidth)
|
|
rmeddis@0
|
2 % UTIL_PSTHmaker accumulates values into bins.
|
|
rmeddis@0
|
3 % No corrections are applied
|
|
rmeddis@0
|
4 % usage:
|
|
rmeddis@0
|
5 % PSTH=UTIL_makePSTH(inputData, dt, PSTHbinWidth)
|
|
rmeddis@0
|
6 %
|
|
rmeddis@0
|
7 % arguments
|
|
rmeddis@0
|
8 % inputData is a channel x time matrix
|
|
rmeddis@0
|
9 % PSTH is the reduced matrix, the sum of all elements in the bin
|
|
rmeddis@0
|
10 %
|
|
rmeddis@0
|
11
|
|
rmeddis@0
|
12 [numChannels numDataPoints]= size(inputData);
|
|
rmeddis@0
|
13
|
|
rmeddis@0
|
14 % Multiple fibers is the same as repeat trials
|
|
rmeddis@0
|
15 % Consolidate data into a histogram
|
|
rmeddis@0
|
16 dataPointsPerBin=round(PSTHbinWidth/dt);
|
|
rmeddis@0
|
17 if dataPointsPerBin<=1;
|
|
rmeddis@0
|
18 % Too few data points
|
|
rmeddis@0
|
19 PSTH=inputData;
|
|
rmeddis@0
|
20 return
|
|
rmeddis@0
|
21 end
|
|
rmeddis@0
|
22
|
|
rmeddis@0
|
23 numBins=floor(numDataPoints/dataPointsPerBin);
|
|
rmeddis@0
|
24 PSTH=zeros(numChannels,numBins);
|
|
rmeddis@0
|
25
|
|
rmeddis@0
|
26 % take care that signal length is an integer multiple of bin size
|
|
rmeddis@0
|
27 % by dropping the last unuseable values
|
|
rmeddis@0
|
28 useableDataLength=numBins*dataPointsPerBin;
|
|
rmeddis@0
|
29 inputData=inputData(:,1:useableDataLength);
|
|
rmeddis@0
|
30
|
|
rmeddis@0
|
31 for ch=1:numChannels
|
|
rmeddis@0
|
32 % Convert each channel into a matrix where each column represents
|
|
rmeddis@0
|
33 % the content of a single PSTH bin
|
|
rmeddis@0
|
34 PSTH2D=reshape (inputData(ch,:), dataPointsPerBin, numBins );
|
|
rmeddis@0
|
35 % and sum within each bin (across the rows
|
|
rmeddis@0
|
36 PSTH(ch,:)=sum (PSTH2D,1);
|
|
rmeddis@0
|
37 end
|
|
rmeddis@0
|
38
|
