function out = rate(data,range,half)
% ”­‰Ξ—¦‚̐„’θ
%
% [t-half, t+half] ‚Μ‹ζŠΤ(• 2*half+1)‚Ε”­‰Ξ—¦ r(t) ‚̐„’θ‚πs‚€
%

trial=size(range,1);		%%% ŽŽs”
len=range(1,2)-range(1,1);	%%% Ψ‚θo‚·’·‚³

out=zeros(len,1);	%%% ”­‰Ξ—¦ r(t) ‚̐„’θ’l 

for i=1:trial
  train=data(range(i,1)-half<data & data<=range(i,1)+len+half)-range(i,1);
			%%% ŽŽs‚²‚ƂɁ}half‚Ύ‚―—]•ͺ‚ɐ؂θo‚·D
  for t=1:len
    out(t)=out(t)+length(train(t-half <= train & train <= t+half));
  end
end

out = out/trial/(2*half+1);