sheephead

sheephead

zero-inflatedなシカ密度

ベイジアンクリギングネタのまたつづき。 あの後、いろいろプロットして遊んでみると、シカ密度ってば0頭/km2なところが相当多いことに気づいてしまいました。 ちなみに、plot(density(deer))ってやって実行するとこんな感じ。
 http://blog-imgs-40.fc2.com/m/y/u/myuhe/deer.png
 う〜ん、0ばっかり。ま、逆に0がない方が恐ろしいわけで良いことではあるんだけど、これだと単純にポアソンでぽーんってわけにはいかないんでしょうな。

 まずは、kuboさんのところで探してみると、こういう0がやたら多いポアソン(zero-inflated poisson)というのも、winbugsでどうにかできるようです。
 お次は世界中のバグサー(bugs-er、つまりwinbugs使ってるひとのことを言うんだろう)の試みを調べるべくwinbugsのMLアーカイブの中から探してみると、世界中には同じような悩みを持つ人も少なからずいるようで、何件か見つかります。 
BUGS Archives -- 2005 (#86)
BUGS Archives -- March 2008 (#6)
 アプローチは結構違うような気がするんですけど、とりあえず隠れマルコフ使ってやっているのかな、と。

 ただ、未だに隠れマルコフがなんなのかが理解できていないとゆー致命的な問題は残すわけですが。 kuboさんところに、Bayesian Models for Categorical Dataなる本にzipモデルのことが書いてあるということが書いてあったので、まずはそれ読んでみますか。ただ、値段がね・・・ちとお高い。あまり、買いまくるのもあれなので、webcatで調べると、都合の良いことに九大図書館に蔵書があるよう。

 来週はどうせ九大に行くんだし、コピーしてくか。
  Data Analysis Using Regression and Multilevel/hierarchical Modelsって本にも書いてありますな。 こっちは、なんと熊本県大にあるそうな。む〜いい加減カード作りにいかねばね。 にしても、まだ相当時間かかりそうですな。

 って右往左往してたら、ODさんはすでに手持ちのデータでマッピングまで終わっていらっしゃいました。
 む〜これが噂に聞くマルチゲイツ級AIなのか


diigo it このエントリをはてなブックマークに登録 このエントリを del.icio.us に登録 このエントリをlivedoorクリップに登録 この記事をPOOKMARKに登録する このエントリをSaafブックマークへ追加

2008年06月30日 | 日記 | トラックバック:0 | コメント:0

コメント


コメントの投稿







管理者にだけ公開する


トラックバック

トラックバックURL
http://myuhe.blog68.fc2.com/tb.php/861-a5f87b12


follow myuhe at http://twitter.com

CalendArchive

≪ ≫
- - - 1 2 3 4
5 6 7 8 9 10 11
12 13 14 15 16 17 18
19 20 21 22 23 24 25
26 27 28 29 30 31 -

最近はこんなの聞いてます

今週のベスト5

Tree-Comment

  • 10/15  にょき
  • 09/21  さおりん
  • 09/10  myuhe
  • 09/10  なか
  • 08/22  myuhe
  • 08/20  さおりん
  • 08/19  myuhe

プロフィール

myuhe

Author:myuhe



物欲の神様
衝動買いの伝道師

mixiは↓
 http://mixi.jp/show_friend.pl?id=1864001
  
last.fmは
 http://www.last.fm/user/myuhe/
 

カテゴリー

  • 日記 (705)
  • 音楽 (59)
  • ニュース (28)
  • ipod (8)
  • firefox (9)
  • 写真 (15)
  • 書籍 (11)
  • 映画 (3)
  • w-zero3[es] (5)
  • 【統計】R (19)
  • X02HT (1)

タグクラウド

ブログ検索


アフィリエイト 通信販売
ホームページ アフィリエイト レンタルサーバー FC2ブログ