%
%  Psth.m 
%  コメントは% を使います、各行の% 以降 はコメントとして無視されます。
%    
%  またスクリプトや関数の最初の行からのコメントは
%　>> help 関数名
%  を実行したときに表示されます。
%  ここに関数の使い方などを書き込んでおけば、他人や未来の自分（しばしば作
%  った本人が関数の使い方を忘れます）がスクリプトを直接読まずに使い方がわ
%  かって便利です。
%   
%   Psth.m : dataからラスター、PSTHを描画する。
%
%    Input :  data ; column1- event code, column2 - Time(msec) 
%

align = data(data(:,1)==64,2);
allspikes = data(data(:,1) ==1|data(:,1)==2 ,2)

buf = [];
for i = 1:length(align)
   spikes = allspikes-align(i);
    spikes = spikes(spikes > -1000 & spikes < 1000);
    buf=[buf ; ones(size(spikes))*i spikes];
end

[spcount, bins]  = hist(buf(:,2),[-995:10:995]);
sprate = spcount/(10*10^-3)/length(align);

subplot(2,1,1) 
plot(buf(:,2), buf(:,1), '.');
subplot(2,1,2)
bar(bins, sprate,1);