WinBUGS

  • あちこち、つまみ食いをしているうちに、こんなサイトへと流れ着きました。
    • その名も『RとWinBUGS』
  • RとWinBUGSで、MCMC-ベイズをしましょう、と。
  • 確率密度分布。
  • MCMCで無記憶(1回前までの情報のみ)で確率密度を更新。
  • そのためには、更新するときに、分布が変化しないことが必要
  • 共役事前分布…一度、自分でメモしている…(これも無記憶性か)
    • 二項分布での事象にβ分布
    • ポワソン分布での事象にガウス分布(ポワソン分布とガウス分布を畳み込むと負の二項分布になる)
      • Rで glm.nb (負の二項分布によるglm)を使う・・・