ぎょーむ日誌 2013-11-14

2013 年 11 月 14 日 (木)

• 0550 起床 …… ここ 2 日ほど「とら黎明泣き」がとまっているような?? …… コーヒー． 朝飯・弁当の準備． 朝飯． とらめし． 0810 自宅発． 路上の雪はほとんど消えた． 0820 保育所着． 0850 研究室着．
• また bootstrap べんきょう ……
• ま， bootstrap なんていろいろ使い道があるんだけど， たとえば信頼区間 (confidence interval, CI) を構成するなんてのは「いかにも bootstrap 法むけ」 な問題かもしれない． ここしばらく勉強していて， この問題は (あたりまえのことなのかもしれないが) 二つの部分にわけられるのでは …… と (よーやく) 気づいた．
• bootstrap で CI を構成するときには， 少なくとも次のふたつの問題を考えなければならない:
  • 何を，どのようにリサンプリングするのか?
  • リサンプルされた数値から，どのように CI を構成するのか?
  これらは，いちおう別々の問題としてあつかえそうだ．
• 後者の CI の構成に関しては， bootstrap 法の教科書・解説に列挙されているように，
  • パーセンタイル法 (変換不変)
  • bootstrap-b 法 (バイアス補正)
  • BCa 法 (変換不変 + バイアス補正)
  といった方法がある …… カッコ内に記したようにそれぞれ特徴がある． 手持ちの日本語本の中では， 「計算統計 I ― 確率計算の新しい手法」 (統計科学のフロンティア 11 -- たぶんすでに絶版) の第 I 部「ブートストラップ法入門」 がうまく解説していると思った (第 I 部著者の一人，Wang Jinfang さん の新著「ブートストラップ入門」 も入手するべきかなあ……)．
• このあたりを R でやる方法については， hoxo_m さんの解説 とかにあります． 信頼区間の「評価」ではこういう「被覆確率」なんかをよく使うようだ …… 蛇足ながら， 信頼区間というものがよくわかっていないヒトは， ここらへんでつまずいている．
• さてさて， 自分では何も確認していないんだけど (先行研究の文献調査とかやってませーん & Efron and Tibshirani 1993 も持ってませーん状態)， 上記の CI 構成法はいずれも 「もとになる材料である数値をどのようにリサンプリングするべきか」 を明示的には指定していない …… はず．
• ということで， そのあたりを検討・改善しましょう， 何かいいことあるかもね …… というのがいま読んでいる論文原稿．
• いやー， 私の読解では， ですが， こういうことをあまりきちんと書いていないように思えたので， bootstrap の勉強をしつつ構想の意図みたいなものもおしはからねばならぬ状態でして ……
• ともあれ， 「どのように，なにをリサンプリングするのか」 は重要な問題です …… とはいえ， bootstrap なんてのは非常に集中的に研究された問題のようだから， すでに数理統計な専門家がだれかが， この論文原稿であつかっている問題について証明つきで何ごとかを結論しているのではなかろーか?
• 上記のようなことを， 論文原稿よみのとりあえずの報告としてメイル書き ……
• ちょっとひと休みして， 昼飯．
• メイルかきのつづき ……
• やっと終わった …… ありゃまあ， 時刻はすでに 1450． いやはや．
• いったん bootstrap 論文原稿からは離脱しよう， まだ完全には読みおわっていないんだけど．
• さて， 次に何をやるかだが …… 先週あげてみたやることりすと， はこんなものだが …… だいたいのしめきり順:
  • bootstrap な論文原稿よみ (さっさとやらなくては)
  • 根の論文原稿よみ (さっさとやらなくては)
  • 「統計モデリング入門」のまちがい一覧の更新 (さっさとやりたいね)
  • 論文原稿の査読報告かき (10/19 まで)
  • 「統計数理」の企画書かき (月末でもいいのかな……その前に執筆者の皆さんに連絡しなくては)
  • 生態学会広島大会の口頭発表編成のおてつだい (例年どーり 11 月末か?)
  • 生態学会さーばー管理情報の共有 …… これはあまり難しく考える必要はなく， 何かてきとうなファイル共有 (編集 permission ができるやつ) にすれば問題なかろう (年内になんとかすればいいのかな?)
  • 自分の論文原稿の改定とアルゼンチンアリ論文かき (このあたりを一番さきにやりたい!!! …… んだけど)
  次は根の論文原稿よみをやるべきなんだけど， ぐずぐずしているうちに査読報告かきの締め切りがせまってきたので …… この査読ぎょーむはちょっとイヤなのだが …… まあ， revision なので先に着手して「てきとー」に終わらせてしまいますか． いやはや ……
• などとうだうだしていたら， すかさず bootstrap 論文原稿の著者から返信がきた． むちゃくちゃに仕事が速いヒトだなぁ ……
• 返信どんどんきていろいろ教えてもらっていますが ……
• 撤退時刻． 1700 研究室発． 保育所でとら回収． 1740 帰宅． とら風呂． とらめし． 晩飯．
• とりあえず， わりとあっさり寝てくれた． 20 時すぎ．
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0630): 米 0.7 合． ネギ・油揚げの味噌汁．
  • 昼 (1200): 研究室． 米 0.7 合． ソーセイジ． ブロッコリー． ミニトマト．
  • 晩 (1830): 米 0.7 合． ネギ・マッシュルームの味噌汁． ニラ・エリンギ・豚肉の炒めもの．