相関関数を求める

2種類の細胞のスパイクがどれだけ同期しているのか、相関関数を求めて類似性を見たいと思っています。

添付された２つのファイルには、2種類の細胞（SPKC01aとSPKC01ｂ）のスパイクの数を1ミリ秒ごとに記録したものです。行がミリ秒（計４００ミリ秒）列は試行回数（計13回の試行）を示しています。単純にその2つのデータの相関を計算したいのですが、どの様に行ったら良いでしょう？

論文で以下の構文で相関関数を求めるらしいのですが、この公式でMATLABで行うと、異なった細胞の種類でも全て同じような結果になってしまいます。おそらくxcorrを用いて計算すれば良いと思うのですが、単純にSPKC01aとSPKC01ｂをｘ、ｙとしてｘcorr(x,y)として計算しても期待通りの結果が得られません。

Answer 1

[Saito]

今試してみましたのは、

Xs=sum(SPKC01a);

Ys=sum(SPKC01b);

D=finddelay(Xs, Ys);

[C9, lag9]=xcorr(Xs, Ys);

t9=lag9(I9);

plot(lag9, C9, [t9, t9], [-0.5, 1];

これで正しいのでしょうか？