再帰定量化解析による非線形時系列解析の手順１．アトラクタの再構成時間遅れパラメータの決定埋め込み次元の決定時系列の相空間内への埋め込み２．再帰プロットの作成半径パラメータの決定その他のパラメータの決定３．各指標の算出再帰率、決定率、最大線長、エントロピー、トレンド再帰プロット１．アトラクタの再構成時間遅れパラメータの決定埋め込み次元の決定２．再帰プロットの作成半径パラメータの決定３．各指標の算出再帰率、決定率、最大線長、エントロピー、トレンドアトラクタの再構成決定論に従う力学系の構造はｄ個の状態変数の関数として記述でき、その運動はアトラクタとして表現できる。分析対象となる系の構造が明らかでない場合には、１個の状態変数からｄ個の状態変数を復元するTakensの埋め込み定理を用いる。Takensの方法によって再構成されたアトラクタは本来のアトラクタと本質的に微分同相が保証される。ターケンスの手法により遅れ時間を用いて未知の状態変数を復元することを埋め込みという。１．アトラクタの再構成時間遅れパラメータの決定埋め込み次元の決定２．再帰プロットの作成半径パラメータの決定３．各指標の算出再帰率、決定率、最大線長、エントロピー、トレンド埋め込み遅延時間の決定埋め込みの際に用いる遅延時間 t は、時系列の平均相互情報量 I(t) が最初に極小となるt の値を用いる。平均相互情報量とは、ある二つの事象において、一方の情報を知ったとき他の情報に関して得られた情報量のこと。この指標はノイズに対して頑健であり、座標系のスムーズな変換に対して不変。以下の式により表される。 I (t )   xn , xn t P( xn , xnt ) log2 P( xn , xnt ) P( xn ) P( xnt ) 両手協応動作における時系列変化の平均相互情報量（cross average mutual information）１．アトラクタの再構成時間遅れパラメータの決定埋め込み次元の決定２．再帰プロットの作成半径パラメータの決定３．各指標の算出再帰率、決定率、最大線長、エントロピー、トレンド埋め込み次元の決定埋め込みの際に用いる次元は、時系列の幾何学的構造に着目する誤り最隣接法を用いる。時刻 n’ から始まるベクトル XnNN' (m) をm次元空間における Xn (m) の隣接ベクトルとする。このとき、Rn(m)を満たす隣接ベクトル数と xn  mt  xnNN '  mt / Rn ( m )  RT を満たす誤りベクトルの比がゼロになるときの m* を埋め込み次元の推定値とする。松葉（2000）より３次元および２次元のアトラクタに埋め込んだローレンツモデルの解１．アトラクタの再構成時間遅れパラメータの決定埋め込み次元の決定２．再帰プロットの作成半径パラメータの決定３．各指標の算出再帰率、決定率、最大線長、エントロピー、トレンド再帰とは相空間内の軌跡が、一定時間後、十分に近い位置へ再び訪れたとき、軌跡が再帰したと定義する。 DELAY = 0 DELAY = 10 X(t) X(t+t) × × RADIUS = 5 X(t) X(t+t) RADIUS = 10 X(t) X(t+t) １．アトラクタの再構成時間遅れパラメータの決定埋め込み次元の決定２．再帰プロットの作成半径パラメータの決定３．各指標の算出再帰率、決定率、最大線長、エントロピー、トレンド半径パラメータの決定半径パラメータは、半径と各指標を両対数軸プロットしたときに両者が直線的に変化する範囲で決定する半径（RADIUS）を最大距離の0.1% から30%まで0.1刻みで変化させたときの再帰率および最大線長。条件間の比較に用いた半径の値を赤の垂直線で示した。１．アトラクタの再構成時間遅れパラメータの決定埋め込み次元の決定２．再帰プロットの作成半径パラメータの決定３．各指標の算出再帰率、決定率、最大線長、エントロピー、トレンド RQAによって得られる指標指標算出法再帰率 (%RECUR) 再帰点と、可能な再帰点との比決定率 (%DETER) 対角線上の再帰点と全再系の決定論的構造帰点との比最大線長 (MAXLINE) エントロピー (ENTRONPY) トレンド (TREND) 意味系の安定性ノイズに関連連続した対角再帰点長の系の安定性最大値リアプノフ指数に関連シャノンのエントロピー系の複雑性再帰点の広がり系の定常性今後の可能性（１）：シミュレーションの利用 Ohgane, et al. (In press) Ohgane, et al. (In preparation) Hirashima, et al. (2003b) 今後の可能性（２）：心理学的な問題への学際的なアプローチ • • • • • • はやさと正確さのトレードオフ学習のスケジュールあがり動作やパフォーマンスの変動動作パターンの変化知覚と行為のカップリング Thank you おまけリアプノフ指数 l ・リアプノフ指数とは何か？相空間内の2点の指数的な時間拡散の指標誤差の指数的増大の速さ初期条件の変化に対する敏感さの度合・リアプノフ指数はどのように計算されるか？ x ( t ) を１次元空間内の時刻 t の座標とし，その写像を f(x(t)) とすると，一次元の力学系ではリアノフ指数 λ は次式のように定義される。 1 l  lim N  N N  log f ' ( x(t )) t 1 システムの相関次元が低く、かつリアプノフ指数が正であれば低次元カオスの可能性がある。・リアプノフ指数によって何が分かるか？不規則変動がノイズに由来するのか、それとも非線形性に由来するのか判別できる。このため、システムがどれくらいカオスに近いのかを判断する数値として利用できる。リアプノフ指数と分岐ダイアグラム複雑な分岐をする場所はリアプノフ指数が大きい情報量・情報量 I C.E.Shannonの情報伝達理論：情報とは確率である生起確率Pの情報量は –log2P [bit] 。・平均情報量（エントロピー）H n H   p(ei ) log2 p(ei ) i 1 あいまいさの尺度。エントロピーが大きいとき、何が起こるか予想することは難しい。・平均相互情報量 I (t) I (t )   xn , xn t P( xn , xnt ) log2 P( xn , xnt ) P( xn ) P( xnt ) ある二つの事象において、一方の情報を知ったとき他の情報に関して得られた情報量両手協応ダイナミクスにおけるアトラクタの形状は再帰定量化分析により推定できるアトラクタの位置f* ノイズの大きさQ および正の最大のリアプノフ指数 | l | f *  f  Q SDf  2|l | %RECUR  1/Q MAXLINE  l システムのモデル入力刺激同定反応選択プログラミング xBlack Box フィードバック環境冗長性・階層性をもち自律的に時間発展する複雑なシステム出力システムのモデル環境観測環境情報重力生理学的（神経-筋）活動物理学的時間発展生理物理学的相互作用外力動作パフォーマンスキーワード •ダイナミカルシステム •自己組織化 •オートポイエーシス •エコロジカルアプローチ •ベルンシュタイン問題 •観測問題 •自律性 •引き込み