KuboLog 2011-09-08

苦情・お叱りは，たいへんお手数かけて恐縮ですが，久保 (kubo@ees.hokudai.ac.jp) までお知らせください．

2011 年 09 月 08 日 (木)

0610 起床．コーヒー．ラジオ体操．朝飯・弁当の準備．朝飯． 0750 自宅発．晴．なぜかしら，本日はまた真夏日になってしまうかも，とのことで． 0810 研究室着．

お， RStudio 紹介記事か …… イマドキの IDE 事情 ``デスクトップ & Webで動作するR言語の統合開発環境「RStudio」''．ゐんどーづ使いの院生のみなさん，もーそろそろ Tinn-R はやめたら …… いつのまにやら， (以前に紹介したときはダメだった) 日本語文字が使えるようになっているし (rstudio-0.94.102-i386.deb)．

ただし，私の Ubuntu 環境だとフォントまわりがいまいちおかしくて，エディタのカーソル位置がヘンみたいだ．この点はまたあとで調べてみよう．

データ解析こんさるメイルかき．質問は，

species ごとの個体数で優占度みたいなものを推定したい，観測データはプロット (調査区画) ごとの種ごとの個体数
N 種いるので多項分布を使った階層ベイズモデルつかいたい
しかし事前分布のディリクレ分布 (ゐんばぐすだと ddirch()) の使いかたがよくわからない

で，割算値ぎらいなわたしの回答:

WinBUGS の ddirch() はめんどうなので，おすすめし
ません．

奇妙に思われるかもしれませんが，このモデリングは
ポアソン分布でよいのです．ただしプロット差やプロ
ットごとの種差をうまくくみこんでください．

もし Agresti の和訳本「カテゴリカルデータ解析入門」
をおもちでしたら，「6.5.3 対数線形モデルとロジット
モデルとの対応関係」などを参照されてください (原著
をおもちでしたら，このへんに対応する部分を)．log
link 関数なポアソン回帰と，logit link 関数な多項分
布回帰のつながりが見えてくるのではないかと思います．

優占度なる割算値は例によって「いらない」わけですが，
世の中にはかかる悪しき割算値で結果を提示しなければ
理解できないヒトもいますので，そういう場合は，上記
のポアソン分布なモデルの事後分布を組み合わせて，優
占度なる呪われ値の事後分布を予測すればよいだけです．

岩波統計モデル本修正のつづき．第 5 章の尤度比検定の説明を終えるために，ちょっと数理統計な本を調べる．

[困ったときの Hoel 本]

Amazon には，意外にも「入門数理統計学」の新品があるようだけど，柳川・大和訳の「統計理論入門」は書誌データすらない …… 私は 1994-05-28 に福岡で買ったようだ．尤度比検定と他の検定の対応関係を調べる …… いやはや，このあたり，ぜんぜんわかってなかったので．

昼飯．正午あたりから曇天，予報されてた真夏日はめでたく回避された! 今年最後の真夏日は 9/1 と確定されたか?

第 5 章の章末をととのえて，ここはとりあえず終わったことにする．

さて，第 6 章を …… と思ったら，大学院生の学会せとぎわこんさる．

一件落着．さて，第 6 章を …… と思ったら，厚岸のとくに匿名を希望する料理人氏がとつぜんあらわれた．厚岸雑談など．

本日はそれでオワった． 1750 研究室発．買いもの． 1810 帰宅．晩飯の準備．

1825 自宅発ミニ大通公園走． 1905 帰宅．

晩飯の準備のつづき．晩飯．

[今日の運動]

朝のラジオ体操． 30 秒ぶらさがり × 3 回．ミニ大通公園走 40 分間．

[今日の食卓]

朝 (0720): 米 0.5 合．ネギ・アサリの味噌汁．キュウリ・モズク酢・カニカマボコ．キムチ．
昼 (1200): 研究室．米 0.7 合．梅干．ゆで豚肉ポン酢あえ．さつまあげ．ブロッコリー．サツマイモ．プチトマト．ミズナ．
晩 (1950): 米 0.5 合．ネギ・豆腐の味噌汁．しめさば．ナス・ピーマン・ネギ・ショウガ・豚肉の炒めもの．

ぎょーむ日誌 2011-09-08

2011 年 09 月 08 日 (木)