• ベストアンサー

マルコフ連鎖モンテカルロ

マルコフ連鎖モンテカルロ法で、既約で非周期性を満たすマルコフ連鎖の持つ不変分布を、推定パラメータの目標分布となるように推移核を構成する際に、ギブス・サンプラーアルゴリズムを用いることを考えた場合、大きな流れとして何故、推移核をギブスサンプラーで与えられる式を導入することによって、推定したいパラメータの目標分布が不変分布に収束していくのかその原理がよくわかりません。 参考書として、「東洋経済」から出版されている「ベイズ計量経済分析」を用いています。 誰か教えていただけないでしょうか?

質問者が選んだベストアンサー

  • ベストアンサー
回答No.2

http://www.ism.ac.jp/~iba/iwaPDF/altproof.pdf ここのサイトにも書いてありますね。。

malakunn
質問者

お礼

ありがとうございます。 ネットは自分も結構調べたのですが、このサイトは見つけることができませんでした。 熟読してみます。

その他の回答 (3)

回答No.4

。。。遷移とうってるつもりが変移になってました。。。。 orz

malakunn
質問者

お礼

いえ、大丈夫ですw 意味は理解できますから。 親切にいろいろ教えていただいて本当にありがとうございます。

回答No.3

リンク先の条件B*は 変移確率πによる変移を1回繰り返して、xからx'に有限の確率で到達出来る というものです。

回答No.1

任意の分布P_0(x)から出発して変移を一回した後の分布が P_1(x) = \sum_{x'}P_0(x')\pi(x' → x) は P_1(x)=_cP(x)+(1-c)R_1(x) と書けます P(x)は定常分布です。0<c<1はP_0によらない定数です。 R_1(x)は確立分布とみなせるので再度計算します。 計算するとまた同じように確立分布が出ます。 この変形を繰り返すとm回変移した分布P_m(x)は P_m(x)=(1-(1-c)^m)P(x)+(1-c)^mR_m(x) 0≦R_m(x)≦1, \sum_{x}R_m(x)=1 |R_m(x)|≦1だからm→∞のときP_m(x)→P(x)

関連するQ&A

  • 定常なマルコフ連鎖の意味

    定常なマルコフ連鎖の意味 定常なマルコフ連鎖とは定常分布を持つ、非周期的・再帰的・既約なマルコフ連鎖という意味でしょうか? つまり、初期状態に依存しない極限があると判断しても良いのでしょうか?

  • マルコフ連鎖の質問です

    至急解答をお願いいたします!!! マルコフ連鎖に関する問題です 下のマルコフ連鎖について 1、定常分布 2、既約になるための必要十分条件 3、p1=p2=1/2、q1=0、q2=1のとき(1)、(3)の再帰性を調べよ という3問がわからず困ってます。 <マルコフ連鎖> r:(1)から(1)に戻る確率 (r=1-p1-q1) p1:(1)から(2)へ移る確率 p2:(2)から(3)へ移る確率 q1:(1)から(3)へ移る確率 q2:(3)から(2) 1-p2:(2)から(1) 1ーq2:(3)から(1) 説明がわかりにくくて申し訳ありません。 かなり困っているので、よろしくお願い致します。

  • マルコフ連鎖(確率過程)に関する問題です!

    マルコフ連鎖に関する質問です。 下のマルコフ連鎖について 1、定常分布 2、既約になるための必要十分条件 3、p1=p2=1/2、q1=0、q2=1のとき(1)、(3)の再帰性を調べよ という3問がわからず困ってます。 <マルコフ連鎖> r:(1)から(1)に戻る確率 (r=1-p1-q1) p1:(1)から(2)へ移る確率 p2:(2)から(3)へ移る確率 q1:(1)から(3)へ移る確率 q2:(3)から(2) 1-p2:(2)から(1) 1ーq2:(3)から(1) 説明がわかりにくくて申し訳ありません。 かなり困っているので、よろしくお願い致します

  • マルコフ過程(確率過程)に関する質問です!

    マルコフ連鎖に関する質問です。 下のマルコフ連鎖について 1、定常分布 2、既約になるための必要十分条件 3、p1=p2=1/2、q1=0、q2=1のとき(1)、(3)の再帰性を調べよ という3問がわからず困ってます。 <マルコフ連鎖> r:(1)から(1)に戻る確率 (r=1-p1-q1) p1:(1)から(2)へ移る確率 p2:(2)から(3)へ移る確率 q1:(1)から(3)へ移る確率 q2:(3)から(2) 1-p2:(2)から(1) 1ーq2:(3)から(1) 説明がわかりにくくて申し訳ありません。 かなり困っているので、よろしくお願い致します。

  • MCMCによるサンプリング

    メトロポリス法やギブスサンプリングについて勉強しているのですが、 任意の確率分布を発生させるのにマルコフ連鎖を利用していますが、 わざわざこういう方法を取るメリットは何なのでしょうか? マルコフ連鎖を使うことのメリットでなくてもギブスサンプリングなどのメリットが知りたいです。 確率分布にしたがってサンプルさせたいなら 乱数を発生させて、とる求める各状態の確率の大きさに 0~RAND_MAXを分割して a < rand() < bのときはこの状態を取る b < rand() < cのときはこの状態を取る・・・・ というようにやればいい気がしてしまいます。 よろしくおねがいします。

  • マルコフ連鎖モンテカルロとは?

    よく,マルコフ連鎖モンテカルロ法(MCMC法)というのを聞くのですが,色々読んでも具体的にイメージできません。 マルコフ連鎖,モンテカルロと単独なら理解できるのですが。。。 掻い摘んでいうと,どういう手法なんでしょうか?

  • ガウス分布に対するベイズ推定

    計算・演算が苦手で嫌いな理系大学生です. 質題にもありますが,ガウス分布に対するベイズ推定の計算がほとんど分らない状態で困っているので,この場をお借りして質問させていただきました. 問題は以下の通りです. p(x|u)=N(x|u,b^-1) D={x_1,x_2,・・・,x_N} :観測データ集合 p(u)=N(u|m,a^-1) :uの事前分布 uは未知パラメータ,b,aは精度,m,uは平均,Nはガウス分布とする. このときのベイズ推定でパラメータ推定分布p(u|D)と推定分布p(x|D)を求めよ というものですが,計算が複雑で計算力もないため答えが求まりません. ちなみに答えは p(u|D)=N(u|Um,g^-1) Um=1/(N+λ)*(Σ^N_n=1*x_n+λ*m) g=a+N*b p(x|D)=N(x|Um,a^-1+b^-1) となるらしいです. この計算ができる方にご回答いただけたら幸いです. 詳細な説明もできるだけ含みで,ご回答よろしくお願い致します.

  • マルコフ連鎖

    ある投手が直球、カーブ、その他変化球をランダムに選択して投球している。直球のつぎに直球、カーブ、その他変化球を投げる確率は、それぞれ1/3,1/2,1/3、カーブのつぎに直球、カーブ、その他変化球を投げる確率は 、それぞれ1/2,1/4,1/4、その他変化球のつぎに直球、カーブ、その他変化球を投げる確率は、それぞれ1/2,1/3,1/6、この投手の投球行動を有限マルコフ連鎖でモデル化するとして、 十分時間が経ったとき、この投手が直球、カーブ、その他変化球を投げる割合を求めるにはどうしたらよいのでしょうか?

  • マルコフ連鎖について

    大学入試問題の「確率と漸化式」の問題を解いているとそのほとんどがマルコフ連鎖になっているような気がします.理論的に扱いやすいからだと思われます.Googleのページランクの技術もマルコフ連鎖が使われているということを聞いたことがありますので,実用的にも広く使われている理論のようです. 大学の先生方にマルコフ連鎖を高校数学で教えようとか,そういう意図があるのでしょうか.少ないとは思いますが,誰かご存知の方お願いします.

  • マルコフ連鎖について

    状態空間S={0,1,2,3}上のマルコフ連鎖の推移確率行列P=(0 1/2 0 1/2) (0 2/3 1/3 0) (0 0 0 1 ) (0 0 1/2 1/2)で与えられている (1)状態の再帰性と一時性について調べよ (2)状態の周期について調べよ この問題がわかる方いらっしゃいますか?