相互相関解析

塩崎 博史
独立行政法人理化学研究所 脳科学総合研究センター
DOI：10.14931/bsd.1087　原稿受付日：2012年5月14日　原稿完成日：2012年8月16日
担当編集委員：藤田 一郎（大阪大学 大学院生命機能研究科）

英語名：cross-correlation analysis

　相互相関解析とは、二つの時系列信号の類似度を評価することで、信号間の関係（同期している、片方が遅れている、無関係である等）を検討する方法である。類似度の評価に相互相関関数を用いることから、その名がついた。神経科学の分野においては、主に細胞間の機能的結合を推定する目的で、同時計測した二つの神経細胞の活動に対して相互相関解析が行われる。

解析方法

　相互相関関数の定義は学問分野、研究者によって異なる。ここでは神経科学の分野でしばしば用いられる定義について述べる。

図１．相互相関関数の例
二つの仮想的な細胞のスパイク活動から計算した相互相関関数。平均10ヘルツで活動電位を発射する二つの細胞から5分間にわたり同時に活動を計測する実験をシミュレートし、相互相関関数を求めた。あるビンにおいて細胞が発火している状態を $X(t)=1$ で、発火していない状態を $X(t)=0$ で表し、計算を行った。A、二つの細胞が同期して活動する傾向がある場合、相互相関関数は時間差0にピークを持つ。B、2つの細胞の活動の間に何の関係性もない場合、相互相関関数は平坦となる

　ある二つの神経細胞の活動（例えば活動電位の発生タイミング）を同時に計測したとする。計測期間を $T$ 個のビンに区切り、 $t$ 番目のビンにおけるある細胞の活動を $X(t)$ で、もう一つの細胞の活動を $Y(t)$ で表す。このとき $X(t)$ と $Y(t)$ の相互相関関数（相互相関ヒストグラム） $C_{XY}(\tau )$ は、次のように定義される。

$C_{XY}(\tau )=\sum _{t=1}^{T}X(t)\ Y(t+\tau ),$

ここで $\tau$ は $X$ と $Y$ の間の時間差（time-lag）を表す。相互相関関数は、細胞 $X$ の活動と細胞 $Y$ の活動との時間的関係性を反映する（図1）。

図2．相互共分散関数の例
二つの仮想的な細胞のスパイク活動から計算した相互共分散関数（青線）。相互相関関数（黒線）には、細胞活動間の共分散（相関）に由来する成分と、平均発火率の変化に由来する成分が含まれる。細胞活動が独立である場合に期待される相互相関関数（赤線）を差し引くことで、共分散に由来する成分を抽出することができる。この例では、二つの細胞の活動の間に正の相関がある。

　神経活動はしばしば確率過程としてモデル化される。この場合、相互相関関数の値は複数の統計量（細胞活動の平均や共分散）を反映する。例えば、二つの細胞の活動が独立、つまり共分散が0であっても、両細胞の活動の平均が時間的に同じように変化すると、相互相関関数は時間差0で最大となる場合がある。実験データから計算した相互相関関数と、二つの細胞の活動が独立である場合に期待される相互相関関数（帰無仮説）の差を取ることで、（いくつかの仮定のもとに）二つの細胞の活動が独立かどうかを統計的に検討することができる。この差は相互共分散関数 $Cov_{XY}(\tau )$ と呼ばれる。

$Cov_{XY}(\tau )=\sum _{t=1}^{T}{\bigl \{}X(t)Y(t+\tau )-{\mu }_{X}(t){\mu }_{Y}(t+\tau ){\bigr \}},$

ここで ${\mu }_{X}(t)$ と ${\mu }_{Y}(t)$ は $t$ 番目のビンにおける細胞 $X$ と細胞 $Y$ の活動の平均を表す。相互共分散関数は、細胞 $X$ の活動と細胞 $Y$ の活動の間の相関（共分散）の度合いを表す（図２）。関数 $Cov_{XY}(\tau )$ のことを相互相関関数と呼ぶ場合もあるので、注意が必要である。細胞活動の平均の変化が感覚入力によって引き起こされている場合、細胞活動の平均の積（帰無仮説）を信号相関、相互共分散関数をノイズ相関と呼ぶ。なお、上のように定義した相互共分散関数の値は計測時間や平均活動度の違いによって変化する。異なる実験間で結果を比較するために、相互共分散関数を総スパイク数、細胞活動の平均、分散等で割ることで正規化する場合がある。

　実際の実験においては、細胞活動の平均 ${\mu }_{X}(t)$ 、 ${\mu }_{Y}(t)$ は未知であるので、帰無仮説は計測した細胞活動をもとに設定しなければならない。帰無仮説の設定は、同じ感覚刺激を繰り返し与えることで得た細胞活動データを用いて行う。通常、片方の細胞活動データの試行番号をランダムに並べ替えたり、試行番号を１つずらすことで、細胞活動が独立であった場合に得られるであろう相互相関関数を計算する。前者の帰無仮説設定法により得た相互共分散関数をシャッフル補正相互相関ヒストグラム（shuffle-corrected cross-correlogram）、後者の方法により得た相互共分散関数をシフト補正相互相関ヒストグラム（shift-corrected cross-correlogram）と呼ぶ。なお、スパイク活動の相関の有意性検定については、上記以外にも数多くの帰無仮説設定方法が考案されている。

　相互相関解析はスパイク活動の他、膜電位、膜電流、カルシウム応答などに適用されている。

解釈

　相互共分散関数の形状から、神経回路の機能的結合関係を推定することができると考えられている。例えば、ある二つの細胞の活動から計算した相互共分散関数が時間差0に幅の狭い大きなピークを持つ場合（図2）、二つの細胞は共通の興奮性入力を受け取っていると考えられる。相互共分散関数のピークの位置、幅を分析することにより、細胞間の興奮性結合や抑制性結合を推定することも可能である。細胞 $X$ から細胞 $Y$ への興奮性結合の強度を定量化するために、細胞 $X$ のスパイクの後どのくらいの割合で細胞 $Y$ がスパイクを発射したか（efficacy）や、細胞 $Y$ が発射したスパイクのうちどのくらいの割合が細胞 $X$ のスパイクの後に発生したか（contribution）といった指標が用いられる。

　相互相関解析は機能的結合を間接的に推定する方法であるため、結果の解釈には曖昧性が残る可能性が指摘されている。例えばシミュレーションによって、異なるメカニズムで働く神経回路から同じようなピーク位置、幅を持つ相互相関関数が得られる場合があることが示されている。