モデルの複雑さをどう選ぶか

Download Report

Transcript モデルの複雑さをどう選ぶか

モデルの複雑さをどう選ぶか

銅谷賢治 [email protected]

ATR 人間情報科学研究所 科学技術振興事業団 CREST

ニューロンモデルのいろいろ

McCulloch-Pitts — 0/1 logic circuit Connectionist — sigmoid back-propagation Transfer function — linear system control/identification theory Integrate-and-fire — 1D Ornstein-Uhlenbeck process Phase equation — circles

ニューロンモデルのいろいろ

(2)

FitzHugh-Nagumo — 2D state phase space Hodgkin-Huxley — 4D state voltage clamp data more currents — ~10D AHP, bursting multiple compartments — ~100D?

morphological data, Ca 2+ imaging intracellular messengers — ???

モデルの目的のいろいろ

現象の説明 bursting のメカニズムは ?

可塑性のメカニズムは ?

現象の一般化と予測 STDP で何が起こりうるか ?

同期/非同期の条件は ?

計算方式の模索 視覚運動変換 強化学習 北野 道川 加藤 青柳 筧 鮫島

具体例

下オリーブ核ニューロンの HH 型モデル (Schweighofer, Doya, Kawato; 1999) 多くの生理データを再現 計算論的意味づけ 基底核強化学習モデル (Nakahara, Doya, Hikosaka; 2001) 行動データを再現 単純な強化学習モデルに何が足らないのか

下オリーブ核/登上線維の特性

非常に低い発火周波数 — 最高でも数 in vitro では、自発リズム活動 Hz 電気シナプス (gap junction) による相互結合 → 近隣細胞の同期発火 高速,精密な運動の学習に必要な情報を伝えら れるのか?

下オリーブ核細胞の Hodgkin-Huxley 型モデル Dendrite I K_ca Soma I h I K_dr I Ca_h [Ca++] I ld I sd gint I ls I Ca_l I Na gint gc

膜電位、イオンチャネル、 カルシウム濃度のダイナミクス Somatic potential: Low- thresho ld Ca 2+ current: H-curren t: Na + current: Delayed rectifie r K + current: Current to dend rit es: Leakage current: Dend riti c potential: High-thre shold Ca 2+ current: Ca 2+ depend ent K + current: Current to soma : Leakage current: Coup li ng current: Activation/inac tivation

k, l , q , h , n , r

: Ca 2+ depend ent K + activation: Ca 2+ dynamics:

C m

Ca _ l

h

dt

Na

K _ dr

ds ls

 

dV

g s h

 

g I Ca

Ca _ l q (V

g Na g m s

K _ dr _

l k

3  l (V

V h

) (V

n

4

s

) (g int

g l

(V /p )(V

s

V l

)

s I

3

h

s I

Na

h (V

V s

Ca

V Na

) (V

s

 

V d

)

V K

) )

I K

_

dr C m dV

Ca _ h K _ Ca

sd ld c dx

  

dt

 (g

g g

l d c

g g

int (V

I Ca

Ca _ h _

h

K _ Ca /(1 

d

f (V

d r

2  

V l

)

I

(V s (V

d K d

_ 

V Ca

) 

V K

) p)) (V

V de

) (V

d dt ds

 (

x

 (

V

)  (

s

 ([

Ca

2 

dt d

[

Ca

2  ]

dt

   3.0

x

) / ])

I

Ca

s

( _

h Ca d V

) /  

I sd

V s

) )  

s V

([

de

Ca

0.075[ ) 2

I

ld

])

Ca

 2 

I

ds

]

I c

 

I ls I app

I app

電流刺激に対する応答

モデル (Iapp=0,-0.5,-1,-1.5 m A/㎠) 実データ (Bal and McCormick 1997)

2細胞間の電気シナプス結合

逆相振動 非同期振動

情報伝達シミュレーション

眼球運動誤差信号のコーディング

gc=0 gc=0.01

gc=0.2

RL Model of 2x5 Task

(Nakahara et al., 2001) visual input DLPF PMv inverse kinematics reward prediction anterior BG SMA M1 somatosensory input posterior BG

Learning Curves

1st day 3rd day 10th day

Reversed Sequence

Lesion in Visual network

New sequence

Needs for Multiple Representations

Visual representation

Quick learning candidates are explicitly given Slow execution visuo-motor transformation on every step

Motor representation

Slow learning ill-posedness in visuo-motor mapping Quick, automatic execution

どう選べばいいのか

?

詳細なモデル 実験データをダイレクトに取り込める — 未知パラメタの方がずっと多い モデルフィットをいかに客観的にできるか 単純化されたモデル 一般的性質/可能性の探索 — 前提が間違っていると全く無意味に 具体化への道筋をいかにつけるか

Computational Neuroscience

Top-down Inference How a functional system should be?

What computation is necessary?

What algorithms/methods are available?

Bottom-up Neurobiological Constraints Behavioral data Cellular/synaptic mechanisms Molecular/genetic mechanisms