-
Notifications
You must be signed in to change notification settings - Fork 0
/
PeriEventRaster_zzg.m
118 lines (97 loc) · 4.96 KB
/
PeriEventRaster_zzg.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
function [raster rastertime ts_origin]=PeriEventRaster_zzg(fileName,DataParam,RefTS,timerange,range,bin_width)
% range=[-2 4]; range是考虑event前后几秒的数据
% fileName='D:\zzgprocessingdata\MUA\PZ Neuron activity\1109\7dat041415001\7dat041415001.nex'; %dirname
% DataParam=['TETSPK05a';'TETSPK09a';'TETSPK17a';'TETSPK21a';'TETSPK21b';'TETSPK25a';'TETSPK29a';'TETSPK45a';'TETSPK49a';'TETSPK53a';'TETSPK57a';'TETSPK61a'];%
% DataParam=DataParam(2,:);
% timerange=[1;2000];
% RefPara='EVT20'; %Nex文件中代表event的时间点，比如agression的时间等等
% bin_width=0.01;
fileName='D:\zzgprocessingdata\MUA\PZ Neuron activity\1109\7dat041415001\7dat041415001.nex';
timerange=[1;2000];
range=[-2 4];
bin_width=0.01;
RefTS='EVt';
[EVT20 DataType DataTime]=GetDataNex(fileName,RefTS,timerange);
DataParam=['TETSPK05a';'TETSPK09a';'TETSPK17a';'TETSPK21a';'TETSPK21b';'TETSPK25a';'TETSPK29a';'TETSPK45a';'TETSPK49a';'TETSPK53a';'TETSPK57a';'TETSPK61a'];%
DataParam=DataParam(2,:);
% [EVT20 DataType DataTime]=GetDataNex(fileName,DataParam,timerange);
[raster rastertime ts_origin]=PeriEventRaster_zzg(fileName,DataParam,RefTS,timerange,range,bin_width);
[Data DataType DataTime]=GetDataNex(fileName,DataParam,[0;100000]);
if isempty(RefTS)
raster=[];time=[];
ts_origin=[];
return
end
ts_origin_temp=[];
bin_num=round(diff(range)/bin_width);
raster=zeros(length(RefTS),bin_num);
CorrectIn=1:length(RefTS);
switch DataType
case 0
case -1
for i=1:length(RefTS)
temp_s=RefTS(i)+range(1);
temp_o=temp_s+bin_num*bin_width;
temp_ts=Data(find(Data>=temp_s&Data<temp_o));
% temp_ts=[temp_ts;temp_s;temp_o];
ts_origin(i).data=temp_ts-RefTS(i);
% 20151208 新添code 用来将firing frequency 低于0.5hz的trial删除掉
% for i=1:length(ts_origin)
% if numel(ts_origin(i).data)/(range(2)-range(1))>0.5
% ts_origin(i).data=ts_origin_temp(i).data;
% end
% end
% %
if isempty(temp_ts)
else
[raster(i,:),temp_n]=hist(temp_ts,bin_num);
raster(i,1)=raster(i,1)-1;
raster(i,length(raster(i,:)))=raster(i,length(raster(i,:)))-1;
% end
clear temp_index temp_ts temp_n
end
end
rastertime=(range(1)+bin_width/2):bin_width:(bin_width*(bin_num-1)+range(1)+bin_width/2);
otherwise
adfreq=DataType;
if bin_width==1/adfreq
invalidIn=[];
for i=1:length(RefTS)
s=1+floor((RefTS(i)+range(1)-DataTime(1))/bin_width);
o=s+bin_num-1;
if s<=0||o>length(Data)
invalidIn=[invalidIn;i];
continue
else
raster(i,1:bin_num)=Data(s:o);
end
end
CorrectIn(invalidIn)=[];
raster(invalidIn,:)=nan;
elseif bin_width>1/adfreq
for i=1:length(RefTS)
invalidIn=[];
bin_width_num=round(bin_width*adfreq);
raster_temp_index=zeros(bin_width_num,bin_num);
raster_temp=zeros(bin_width_num,bin_num);
temp_s=floor((RefTS(i)+range(1)-DataTime(1))*adfreq)+1;
temp_o=temp_s+(bin_num-1)*bin_width_num;
raster_temp_index(1,:)=temp_s:bin_width_num:temp_o;
if temp_s<=0||temp_o>=length(Data)
invalidIn=[invalidIn;i]
continue
end
raster_temp(1,:)=Data(raster_temp_index(1,:));
for j=2:bin_width_num
raster_temp_index(j,:)=raster_temp_index(j-1,:)+1;
raster_temp(j,:)=Data(raster_temp_index(j,:));
end
raster(i,:)=mean(raster_temp);
end
raster(invalidIn,:)=nan;
CorrectIn(invalidIn)=[];
else
'bin width should not be smaller than sampling period'
end
rastertime=(range(1)+bin_width/2):bin_width:(bin_width*(bin_num-1)+range(1)+bin_width/2);
end