ぎょーむ日誌 2007-05-25

2007 年 05 月 25 日 (金)

• 0830 起床． コーヒー． 朝飯． 0930 自宅発． 曇． ちょっと寒い． 0945 研究室着．
• うだうだとしてから， またぢりぢりとトドマツ原稿にとりくむ． 今日は (discussion の内容も限定されてきたので) モデル説明まわりの作文改善， とか． しかし作文しようとすると， いろいろトドマツ質問メイルをいただいたり．
• うーむ， 散布カーネルな確率分布の平均値って， ``tail'' 部分に敏感だろう， と． そして調査 plot を広げれば広げるほどますます 「遠方に散布されてしまった子供」が見つかったりするので， それにつれて 「平均散布距離」はいよいよ延伸されてしまうかも， と．
• そして「random effects とは何か?」 の定義に関して「一般化した」説明をまたひとつ思いついた． 平均値にもとづく説明はヤメて， 「分布の中央値は変えないけれど， ばらつきを変えるのが random effects」 というのはどうだろう? これに対して fixed effects は中央値を変える， と． おそらくこれで consistent なハズなんだけど， わかりやすくないので (つまりふつーは確率分布を中央値で考えないので) あまり多くの人から支持されないような気もする．
• 作文なおしは着手できぬまま昼飯．
• R-News の最新号 (Volume 7/1, April 2007) に ``A New Package for Fitting Random Effect Models'' として library(npmlreg) (Nonparametric maximum likelihood estimation for random effect models; library(help = npmlreg)) が紹介されてるな (引用太字は久保による):

  Random effects have become a standard concept in statistical modelling over the last decades. They enter a wide range of applications by providing a simple tool to account for such problems as model misspecification, unobserved (latent) variables, unobserved heterogeneity, and the like. One of the most important model classes for the use of random effects is the generalized linear model. Aitkin (1999) noted that ``the literature on random effects in generalized linear models is now extensive,'' and this is certainly even more true today.

  で， ここでいう nonparametric maximum likelihood は random effect parameters の分布を決めずに (たとえば glmmML() だとこれを正規分布と仮定している)，

  The random effects distribution can be considered as an unknown mixing distribution and the NPML estimate of this is a finite discrete distribution. This can be determined by fitting finite mixture distributions with varying numbers of support points, where each model is conveniently fitted using a straightforward EM algorithm.

  …… ということで， EM algorithm つかった GLMM の推定 package なのである． ふーむ， またこんど使ってみるか ……
• 8 月のアリ研究会， 事務とりまとめの三浦さんからのキビしい追及 (その背後に executive order?) から逃れられなくなり， 「アリデータの統計モデリング」 なるいいかげんな講演タイトルをだしてしまう (しめきりは 5/31 だったはずなのに……)． さっさとトドマツと地衣類繁殖の作文ぎょーむをかたづけてしまわなくては ……
• と思ったら原稿チェックのすばやい練さんから 「ここが変ですよ」 ご指摘一覧が． うう， 私が昨日追加した作文も含まれている ……
• 修正してアップロード． 時刻はすでに 1700 ちかい．
• トドマツ原稿修正のほうはあまり進捗せず …… 1800 研究室発． 雨． 1825 帰宅． 晩飯の準備． 晩飯．
• Wired Vision のりにゅーあるが進捗している (Hotwired からのお知らせ)．
• [今日の運動]
  • うんどう休養日
• [今日の食卓]
  • 朝 (0850): パンケイキ． ニラ・豆腐のトマトスープ．
  • 昼 (1330): 研究室お茶部屋． 食パン． チーズ． リンゴ．
  • 晩 (2020): 米麦 0.7 合． ネギ・ジャガイモ・ニンジン・ブロッコリー・ブナピー・ 干しエビ・ホタテのクリームシチュー． 豚タン・タマネギ・プチトマト・レタスのサラダ． キュウリ．