「ドリフト拡散モデル」の版間の差分

''産業技術総合研究所人間情報インタラクション研究部門''<br>

''専修大学人間科学部心理学科''<br>

英：drift diffusion model　独：Drift-Diffusionsmodell

類義語：逐次サンプリングモデル (sequential sampling model)

{{box|text=　ドリフト拡散モデルは、刺激呈示から反応が起こるまでの意思決定プロセスを説明する数理モデルの一つである。反応選択と反応時間を説明するモデルとして、心理学や神経科学の研究において幅広く用いられている。その性質は確率過程に関する数学的理論や計算機シミュレーションにより詳細に調べられている。実際の反応データからモデルのパラメータを推定することも可能であり、パラメータの個人差や群間差も検討されている。モデルの振る舞いに対応する神経活動も観測されており、意思決定の神経基盤のモデルとしても用いられている。

==概要==

 

　ドリフト拡散モデルが適用されるのは、刺激の呈示後1秒から2秒程度の時間内に2種類の反応から一つを選択することが求められるようなタイプの意思決定である。熟考の末判断を下すような時間のかかる意思決定は適用対象ではないとされる。また、いずれかの反応が正解であり、もう片方の反応は不正解とされる場合が多いが、必ずしもその限りではない。適用される意思決定のタイプとしては、呈示された刺激が以前にも呈示されたものであるか新規なものであるかを判断する[[記憶]]課題<ref name=Ratclif1978 />、明るさや動きなどの刺激の物理的特徴について判断する[[知覚]]的判断<ref><pubmed> 12812276</pubmed></ref><ref><pubmed>15066392</pubmed></ref>、その他、[[語彙]]判断<ref><pubmed>14756592</pubmed></ref>、数量的判断<ref><pubmed>24446719</pubmed></ref>、文章処理<ref><pubmed>20431690</pubmed></ref>、選好判断<ref><pubmed>20835253</pubmed></ref>など多岐にわたる。

　ドリフト拡散モデルは、[[逐次サンプリングモデル]] (sequential sampling model) の一種である。逐次サンプリングモデルでは、刺激が呈示されると生体は時間経過とともに感覚入力やこれまでの経験の記憶など、判断に必要な情報 (エビデンス) を蓄積していき、その蓄積が境界 ([[閾値]]) を越えた時に反応が出力されると仮定する。蓄積されるエビデンスはノイジーであり、確率的な揺らぎを伴うと仮定される (刺激自体がノイズを含む場合もあるが、その入力情報を伝える[[神経細胞]]の活動にもノイズが含まれると考えられる)。

　行動課題を実施した際に反応までにかかる時間は、'''図1'''に示すように、(1)刺激の読み込み、(2)エビデンスの蓄積、(3)反応 (ボタン押しなどの運動) に分解することができる。(1)刺激の読み込みと(3)反応は、判断に関わる過程ではないので、それにかかる時間は[[非決定時間]] (non decision time) と呼ばれる。(2)エビデンスの蓄積にかかる時間は、[[決定時間]] (decision time) と呼ばれる。ドリフト拡散モデルをはじめとする逐次サンプリングモデルを用いることで、非決定時間の推定と決定時間の生成に関わるパラメータの推定を行うことができる。

　ドリフト拡散モデルにおけるエビデンスの蓄積過程は、開始点 <math>z</math> から始まり、一定の[[ドリフト率]]に従ってエビデンスが蓄積される ('''図1''')。そして、上下いずれかの境界までエビデンスが蓄積されると反応が出力される。'''図1'''の場合、上側の境界<math>a</math>に到達すると反応Aが出力され、下側の境界0に到達すると反応Bが出力される。開始点から境界<math>a</math>または0までの距離が長いほど、エビデンスの蓄積にかかる時間が長くなる。また、開始点から境界までの蓄積過程における速度は、ドリフト率に依存する。反応Aが正解である場合は、ドリフト率は正の値をとると仮定される。ドリフト率が大きいほど、境界まで到達する時間は短くなる。開始点、ドリフト率、境界、そして非決定時間がドリフト拡散モデルの振る舞いを決定する主なパラメータである。

　刺激が呈示されてから、刺激情報の読み込みや反応の準備に必要な時間が経過してからエビデンスの蓄積が開始され、<math>x</math>が変化する。エビデンスの蓄積過程は以下の式のように連続時間上で定義される確率過程である、[[ウィーナー過程]] ([[ブラウン運動]]) に従うとする。

:<math>dx = v dt + \sigma dW</math>

ここで、<math>dx</math> は微小な時間間隔 <math>dt</math>の間の<math>x</math>の変化を表す。<math>v</math>はドリフト率である。反応Aが正解である場合は<math>v > 0</math>、反応Bが正解である場合は<math>v < 0</math>とする。<math>\sigma dW</math>は平均が0で分散が<math>\sigma^2 dt</math>となる正規分布に従うホワイトノイズを表す。

　ウィーナー過程は<math>dt \to 0</math>という極限における連続時間上で定義されるが、計算機上でシミュレーションする場合は、離散時間で近似する必要がある。例えばシンプルな近似法としては以下のようなものがある。微小な時間幅<math>\Delta t</math>を考え (例えば<math>\Delta t = 0.001</math> (秒)とする)、平均0, 分散1の[[wj:標準正規分布|標準正規分布]]に従う[[wj:正規乱数|正規乱数]]<math>\epsilon</math>を用いて、この時間幅<math>\Delta t</math>あたりの変数<math>x</math>の変化量を、以下の式で決定する。

:<math>\Delta x = v \Delta t + \sigma \epsilon \sqrt{\Delta t} </math>

この式で<math>x</math>を更新していくことによりエビデンスの蓄積過程をシミュレートできる。'''図2'''の各軌道はこの方法により得られたものである。

　標準的なドリフト拡散モデルでは、開始点とドリフト率、および非決定時間は、試行間で変動すると仮定される<ref name=Ratclif1978 />。ドリフト率の試行間変動は、刺激に対する注意の変動などに対応すると考えられ、標準偏差<math>\eta</math>の正規分布に従って変動すると仮定される。この変動を仮定することで、正反応より誤反応の方が反応時間が長くなるという現象が説明可能となる。これは、ドリフト率が小さくなる試行において、ノイズの影響が強くなるため誤反応が起こりやすくなり、かつ反応時間が長くなるためである。

　開始点の試行間変動は区間<math>[z-s_{z}, z+s_{z}]</math>の一様分布に従うと仮定される。これは、ある特定の刺激がどの程度呈示されやすいかについての期待が試行間で変動することを表現する。この変動により、誤反応が起こる試行で反応時間が短くなることが説明できる。これは、開始点が誤反応側の境界に寄っているときに、早い時間帯で誤反応が起きやすくなるためである。一度正反応の方に<math>x</math>が近づけば、開始点の影響はなくなるため、遅い時間帯では開始点の影響による誤反応は生じにくい。また、非決定時間も区間<math>[T_{er}-s_{t}, T_{er}+s_{t}]</math>の一様分布に従うと仮定される。

　以上のように、標準的なドリフト拡散モデルのパラメータは、開始点の平均 (<math>z</math>)、開始点の試行間変動 (<math>s_{z}</math>)、ドリフト率の平均 (<math>v</math>)、ドリフト率の[[wj:標準偏差|標準偏差]] (<math>\eta</math>)、境界 (<math>a</math>)、非決定時間の平均 (<math>T_{er}</math>)、非決定時間の試行間変動 (<math>s_{t}</math>) の7つとなる。

==選択確率および反応時間分布の理論解==

　上述のドリフト拡散モデルにおいて、各パラメータを固定した場合 (試行間変動は仮定しない場合) には、それぞれの選択肢を選ぶ確率、およびその反応時間の分布は以下のように求められる <ref name=Ratclif1978 />。下側の境界 (0) に到達し、反応Bが起こる確率は、ドリフト率<math>v</math>、境界<math>a</math>、開始点<math>z</math>の関数として、

:<math>P(v, a, z) = \frac{e^{-2va/\sigma^2}-e^{-2vz/\sigma^2}}{e^{-2va/\sigma^2}-1}</math>

で与えられる。ただしドリフト率<math>v</math>が0だった場合はこの確率は

:<math>P(v, a, z) = 1 - z/a</math>

となる。反応Bが起こり、かつ非決定時間を除いた反応時間が<math>t</math>となる条件付き確率密度は

 

で与えられる。境界<math>a</math>に到達し反応Aが起こり、かつその反応時間が<math>t</math>となる確率密度は、上の式において<math>v</math>を<math>-v</math>で, <math>z</math> を<math>a -z</math>で置き換えることで得られる。'''図2'''の上下の曲線はこれらの式により得られた条件付きの確率密度関数である。シミュレーションにより得た反応時間のヒストグラムも試行数が増えるにつれてこの分布に近づいていくことがわかる。

 

 

* Fast-dm<ref><pubmed>18183889</pubmed></ref>: Windowsで動作するソフトウェア。[[最尤推定]]、[[Kolmogorov-Smirnov]]、&chi;²最小化などの方法が可能である。

* DMAT<ref><pubmed>18411528</pubmed></ref>:ドリフト拡散モデル用MATLABツールボックス。

* EZ2<ref><pubmed>17546727</pubmed></ref>: ドリフト拡散モデルから試行間変動の推定を除いたシンプルなEZ拡散モデルを用いてモデルフィッティングするRパッケージ。

* HDDM<ref><pubmed>23935581</pubmed></ref>: ドリフト拡散モデルを用いて[[階層ベイズ推定]]を行うPythonパッケージ。

* hBayesDM<ref><pubmed>29601060</pubmed></ref>: 意思決定課題に対して階層ベイズ推定を行うRパッケージ。ドリフト拡散モデル以外にも様々な意思決定のモデルを用いることができる。

 

　それぞれのモデルや推定方法には仮定がおかれていることがあり、モデルフィッティングに用いるデータがその仮定に合っているかどうかは事前に確認する必要がある。各種推定法に関する専門家による推奨については、Boehmらの論文<ref name=Boehm2018><b>Boehm, U., Annis, J., Frank, M. J., Hawkins, G. E., Heathcote, A., Kellen, D., Krypotos, A.-M., Lerche, V., Logan, G. D., Palmeri, T. J., van Ravenzwaaij, D., Servant, M., Singmann, H., Starns, J. J., Voss, A., Wiecki, T. V., Matzke, D., & Wagenmakers, E.-J.(2018). </b><br>Estimating across-trial variability parameters of the Diffusion Decision Model: Expert advice and recommendations. <i>Journal of Mathematical Psychology</i>, 87, 46–75</ref>にまとめられている。

　ドリフト拡散モデルでのモデルフィッティングにあたっては、十分なデータ数が必要になる。特に反応時間の分布の情報を用いてパラメータ推定をする方法では、試行数が100程度の場合はドリフト拡散モデルの試行間変動性にかかわるパラメータの推定が真値からずれることが示されている<ref><pubmed>18229471</pubmed></ref>。試行間変動性のパラメータの推定にはできるだけ多くの試行数が必要になるが、試行数を増やすと参加者の動機づけが低下する、疲れの影響が出る、などの問題も生じる。また、そもそも試行数を増やすことが難しい実験状況も多い<ref name=Boehm2018 />。そこで、パラメータ推定にあたり、参加者集団のパタメータの分布も仮定した階層ベイズ推定を行うことで、各参加者の試行数は少なくとも安定した推定をする方法も提案されている<ref><pubmed>23935581</pubmed></ref><ref><pubmed>29601060</pubmed></ref>。

[[Image:DDM_z_vs_v.png|thumb|520px|<b>図3.反応時間分布に及ぼすドリフト率 (左) 開始点パラメータ (右) の影響</b><br>破線は比較のための基準となるモデル (<math>v = 1.0</math>, <math>z = 0.5</math>, <math>a = 1.0</math>) を表す。縦軸のスケールが異なるが、左右のパネルで同じ分布である。左のパネルではドリフト率を大きくした場合 (<math>v = 2.0</math>) 、右のパネルでは開始点を高くした場合 (<math>z = 0.7</math>) の結果をそれぞれ実線で表している。]]

 

 

　例えば、開始点パラメータ<math>z</math>を増加させることと、ドリフト率<math>v</math>を増加させることはいずれも反応Aの選択確率を増加させ、その平均的な反応時間を短くする効果があるが、その反応時間分布の形状に与える影響が異なる。'''図3'''の左のパネルでは、開始点パラメータ<math>z</math>を固定し、ドリフト率を増加させている (実線が増加後)。この場合、反応Aの確率が高くなり、速い反応時間の密度が増加するため平均反応時間は短くなるが、反応時間の分布のピーク (最も密度が高くなる位置) はほとんど変化しない。一方、開始点パラメータ<math>z</math>を<math>a</math>に近づけた場合 ('''図3右''') は、反応Aの反応時間分布のピークが早い時間帯にシフトし、分布の歪みが大きくなっている。

 

　ドリフト拡散モデルのパラメータの推定値を利用して、選択に関するプロセスの個人差に影響する要因も検討されている。代表的な事例は加齢の効果について調べたものである。例えば高齢者は一般に多くの認知課題において若年者に比べて反応が遅くなることが示されているが、ドリフト拡散モデルを適用して検討した研究では、高齢者の反応には長い非決定時間、そして境界の間隔が大きいという特徴はあるものの、ドリフト率には若年者との違いはないということが報告されている<ref><pubmed>19962693</pubmed></ref>。

　一方で、幼児では境界パラメータが大きく、非決定時間が長いことに加え、ドリフト率も比較的小さいことが示されている<ref><pubmed>22188547</pubmed></ref>。また、[[注意欠如・多動症]] ([[ADHD]]) や[[読字障害]] ([[dyslexia]]) を有する若年者はそうでない統制群に比べ、ドリフト率が低い傾向があることを示した研究もある<ref><pubmed>20926067</pubmed></ref><ref><pubmed>22010894</pubmed></ref>。[[知能指数]] ([[IQ]]) に関しては、高IQ群は低IQ群よりドリフト率が2倍程度高いという結果が報告されている <ref><pubmed>19962693</pubmed></ref> <ref><pubmed>21707207</pubmed></ref> 。

==神経活動との関連==

　主に[[サル]]を被験体とした単一細胞レベルでの神経活動記録により、エビデンスの蓄積過程に対応する神経活動が検討されてきた。例えば視線でターゲットを選択することで反応する意思決定課題においては、ターゲットの方向へのサッケード時に選択的に活動する[[LIP野]] ([[lateral intraparietal cortex]]) の細胞は刺激の呈示とともに徐々に活動が増加し、その発火率 (単位時間あたりの活動電位の数) がある閾値に到達したときに[[サッケード]]反応が起こるということが観測されている。その振る舞いはドリフト拡散モデルをはじめとする逐次サンプリングモデルと対応付けられて議論されている<ref><pubmed>8570606</pubmed></ref>。

 

[[Image:逐次サンプリングモデルの図.png|thumb|520px|<b>図4.逐次サンプリングモデルの種類</b><br>Ratcliff et al.,2016<ref name=Ratclif2016 />を元に一部改変。]]

 

 

[[Image:LBAの概要.png|thumb|420px|<b>図5.線形弾道蓄積モデルにおける反応と反応時間の生成過程</b>]]

 

　ドリフト拡散モデル以外の代表的な逐次サンプリングモデルとして、[[線形弾道蓄積モデル]]<ref><pubmed>18243170</pubmed></ref>がある。'''図5'''にあるように、線形弾道蓄積モデルは、ドリフト拡散モデルと類似しているが、エビデンスの蓄積の基準が絶対的なことと確率的ではない点が異なる。

　ドリフト拡散モデルでは、反応はエビデンス蓄積が上の境界と下の境界のどちらに到達するかで決まる相対的なものであった。一方、線形弾道蓄積モデルでは、それぞれの反応は独立してエビデンスの蓄積を行って、最終的に先に閾値(<math>b</math>)に到達した反応が出力される（'''図5'''の場合、先に<math>b</math>に到達した反応Aが出力される）。

　エビデンスの蓄積が始まる点を開始点(<math>a</math>)と呼び、選択肢で同一のこともあるが、異なることもある。開始点の位置の違いは、エビデンスの蓄積の前に存在する選択肢に対するバイアスとして解釈される。ドリフト拡散モデルと同様にエビデンスの蓄積の速さはドリフト率(<math>d</math>)が決めるが、蓄積過程は線形かつ非確率的である。各試行のドリフト率(<math>d</math>)は、平均<math>v</math>、標準偏差<math>s</math>の正規分布に従い、各試行の開始点(<math>a</math>)は、0から<math>A</math>(開始点の上限)の一様分布に従う。決定時間は、<math>(b-a)/d</math>で求めることができ、非決定時間 (<math>\tau</math>)は、全試行で一定とする。

　<math>a</math>と<math>d</math>は、推定するパラメータではなく、<math>v, b, A, s, \tau</math>が推定するパラメータになる。線形弾道蓄積モデルは、ドリフト拡散モデルよりも推定するパラメータが少なく、選択肢が３つ以上ある状況にも適用できるという利点がある。

==強化学習モデルとの統合==

　ドリフト拡散モデルは試行内の刺激呈示から反応出力 (選択) までのプロセスを表現するモデルであるが、試行間の選択傾向の変化を表す他の数理モデルと組み合わせることもできる。例えば、報酬に基づく学習のプロセスを表現する代表的なモデルである[[強化学習モデル]]と組み合わせた枠組みが提案されている<ref><pubmed>27966103</pubmed></ref><ref><pubmed>25589744</pubmed></ref>。