ぎょーむ日誌 2011-07-21

2011 年 07 月 21 日 (木)

• 0640 起床． 朝飯・弁当の準備． 朝飯． コーヒー． 0810 自宅発． 晴． まあ今日も夏日を回避できそーな涼しさで (実際には，最高気温が 25.2 度だったので夏日)． 0830 研究室着．
• ぢりぢりと， 岩波統計モデル本 第 7 章 GLMM 章の修正．
• 今まであっちでちょこちょこ， こっちでがりがりと文章・図を削っていって， 合計 20 ペイジ以上の短縮をしたのだけど …… この章は図をつけくわえたほうが良さそう．
• ちょっとめんどうなので， 時間かかる．
• 昼飯．
• その後もぢりぢり．
• こんなかんじの「分布を混ぜるとは?」 を説明する図． R と LaTeX マクロの複合作図わざ．

• ようするに GLMM によって生成される混合な確率分布の概念図にすぎないのだが …… 「えーっ? なんでこんなのわざわざ図示するの?」 と思われるヒトもいるかもしれないけれど， GLMM ゆーざーたちに対する， さりげなくかつ執拗なる調査観察の結果， じつは少なからぬヒトたち， このへんわかってないと判明したため．
• ありがちな状況:
  「論文投稿したら GLMM 使えっていわれたから， R でやって再投稿したら accept されたんだけど ……いまいちハッキリしないんだよね」
  「どのあたりがハッキリしませんかね」
  「えーと，ほら GLMM って『個体のランダム効果』なんかをモデル式に加えて， なんか誤差項が正規分布みたいな? …… だけど GLMM でも GLM と同じようにポアソン分布のままだから， それでいいんだよね?」
  「いえいえ，個体差の値ごとに異なるいろいろなポアソン分布を， 個体差の重みづけて混ぜあわせて，ポアソンとはちがう分布にしてるんです」
  「えっ，それって正規分布になるんだっけ?」
  …… といったよーな， まるでかみあわない対話をこれまで何度くりかえしてきたことか．
• まー， GLM は直線回帰のアナロジーで (やや) ブラックボックス的に「理解」して， その「外側」に何かパーツをくっつけると GLMM になるとゆー把握 …… なんでしょうかね．
• 数式がよくわからなかったら， 図を描いてみればいいと思うんだけど …… 自分で試行錯誤して作図するヒトもあまりいないので， 過保護主義な岩波本に図を追加しますか， とゆーかんぢで．
• この図でもわからない読者はけっこういるだろうけど， 何人かは「あ，そうか」となるかもしれないので， 「分布を混ぜる」図に 2 ペイジを消費．
• ばててきたし， 撤退時刻だからそろそろ …… とうだうだしていたら， Eli さんにつかまり今日もいきなり PC 修理ぎょーむ．
• うー， これは尋常いちようならざるブチこわれかたですよ …… ファイルシステムがやられてる?
• うまくいかないので， また明日以降ということで終了． 1840 研究室発． 1900 帰宅． おそくなってしまった． 晩飯の準備． 晩飯．
• あとかたづけなどやって， 少し休んでから， また図の改良をすこしばかり．

• 明日はどたばたとした一日になるのだろうか ……
• [今日の運動]
  • うんどう休養日になってしまった ……
• [今日の食卓]
  • 朝 (0720): 米 0.5 合． ネギ・タマネギ・ホタテの味噌汁． ラッキョウ． キュウリのピクルス．
  • 昼 (1230): 研究室． 米 0.7 合． カブ・ニンジン・タケノコ・コンニャク・鶏肉の煮物． イカ一夜干し． ゴボ天． プチトマト． キュウリ．
  • 晩 (2000): 米 0.5 合． ネギ・タマネギ・ホタテの味噌汁． シマホッケ開きのグリル． ダイコンおろし． ニンジン・ズッキーニ・牛肉のオイスターソース炒め．