[論文レビュー] How (not) to assess the importance of correlations for the matching of spontaneous and evoked activity

József Fiser, Máté Lengyel+3 more|arXiv (Cornell University)|Jan 28, 2013

Neural dynamics and brain function参考文献 5被引用数 7

ひとこと要約

この論文は、オクンらが神経活動の相関が付随的であると主張するのを挑戦し、そのサーローに基づく分析が概念的および統計的に誤りであることを示している。最大エントロピー・モデルと発散度指標を用いて、元のデータにおける静止状態と刺激誘発状態の活動パターンの一致が、低次の統計量だけでは説明できないことを示し、高次の相関が観察された神経パターン一致に不可欠であることを確認している。

ABSTRACT

A comment on `Population rate dynamics and multineuron firing patterns in sensory cortex' by Okun et al. Journal of Neuroscience 32(48):17108-17119, 2012 and our response to the corresponding reply by Okun et al's (arXiv, 2013).

研究の動機と目的

• オクンらの主張、すなわち静止状態と刺激誘発状態の神経活動パターンの一致が、集団の発火レートの揺らぎといった単純な統計的性質によるものであるという主張を反証すること。
• オクンらのサーローに基づく分析手法における概念的および統計的欠陥を特定し是正すること。
• 単一ニューロンの発火レートや集団スパiking数を越えた高次の相関が、観察された神経活動パターンの一致を説明するために不可欠であることを示すこと。
• オクンらが提案した単純なサーローでは説明できないことから、ベルケスら（2011年）の元の発見を検証すること。
• 低次の統計構造と高次の統計構造の両者による一致を区別できる、より厳密な発散度ベースの比較法を提唱すること。

提案手法

• 単一ニューロンの発火レートと、同時に活性化するニューロン数の分布を保持するが、機能的結合性は仮定しない最大エントロピーのサーロー・モデルを構築すること。
• Kullback-Leibler 発散度を用いて、真の神経活動分布（P, Q）とそれらのサーロー近似（P̃, Q̃）およびそれら同士を比較すること。
• 発散度指標 D[P||Q] を用いて、静止状態（S）と刺激誘発状態（M）の一致を評価し、D[P̃||Q̃] と比較することで、その一致が低次の統計量に起因しているかどうかをテストすること。
• オクンらのサーロー手法を用いて、元のベルケスらのデータセットを再分析し、その主張の妥当性を検証すること。
• 合成のガウス混合例を用いて、サーローに基づく比較の概念的限界を説明する。
• 真の分布 P とサーロー分布 Q̃ の間の発散度 D[P||Q̃] を新たに導入し、P と Q の一致が低次の統計構造か高次の統計構造かを区別できるようにすること。

実験結果

リサーチクエスチョン

• RQ1 静止状態と刺激誘発状態の神経活動パターンの一致は、単一ニューロンの発火レートや集団スパiking数といった低次の統計的性質によって完全に説明可能か？
• RQ2 オクンらが提案したサーロー・モデル（低次の統計量のみを保持）は、観察された神経活動パターン一致を再現するのに十分か？
• RQ3 オクンらの分析における、サーローの静止状態と刺激誘発状態の分布の類似性は、特に V1 において、実際のデータの統計的性質を的確に反映していると言えるか？
• RQ4 発散度ベースの比較は、静止状態と刺激誘発状態の一致が高次の相関に起因するか、それとも単に低次の統計量に起因するかを区別できるか？
• RQ5 相関の真の役割は何か、特に感覚皮膚における統計的最適性の文脈で。

主な発見

• オクンらのサーロー・モデルは、元の V1 データの統計的性質を正確に再現できず、発展に伴い真の分布とサーロー分布との間の発散が顕著に増加する（静止状態活動では Spearman の ρ = 0.68、p = 0.005）。
• 正しい発散度指標を用いて再分析した結果、真の分布とサーロー分布の間の発散が顕著であるため、静止状態と刺激誘発状態の一致は低次の統計量だけでは説明できないことが示された。
• オクンらが用いたサーローに基づく比較法は本質的に誤りであり、サーロー分布間の一致が低次の構造か高次の構造かを区別できないため、相関の重要性を評価するのに情報が得られない。
• 著者らは、非ガウス分布を用いた合成例を通じて、真の分布 P と Q の間の発散度 D[P||Q] が小さく、かつ同様にサーローの発散度 D[P̃||Q̃] もほぼ同じであっても、その一致が低次の統計量に起因するとは限らないことを示した。
• 成熟個体では、真の静止状態と刺激誘発状態の一致（D[M||S] ≈ 100）が保持されているが、真の分布とサーロー分布の間のクロス距離（例：D[M||S̃Okun] ≈ 200）は顕著に大きいことから、高次の構造が不可欠であることが確認された。
• 提案された発散度比較 D[P||Q̃] は、高次の統計構造による一致と低次の統計構造による一致を区別できるが、オクンらが提唱した手法とは異なり、その区別ができない。