ぎょーむ日誌 2010-09-(01-10)

2010 年 09 月 01 日 (水)

• 0700 起床． 朝飯の準備． 朝飯． コーヒー． 0855 自宅発． 晴． すでに真夏日． 0910 研究室着．
• 今日は 13 時から代理授業やらんといかんので …… まあ， そのときに宣伝する予定の 統計学授業 2010 の内容を検討してみる．
• 昼飯． 授業の準備．
• 13 時から環境科学院大学院講義「植物生態学」担当者が誰もいないので代理授業， 参加者数 2 名． そもそもまったく人気のない授業で (最大参加者数 5)， しかも今は調査シーズンなので ……
• しかし二人あいてというのは， じつは， なかなかやりやすい． それぞれ北大農場と北大博物館の女性大学院生だったのだけど …… こちらも椅子に座って三人が至近距離で雑談するところから始めればよい． 「どんな研究やってます」 「植物の遺伝子の発現と形質の関係を……」 「なるほど，これこれこういうデータ解析できますよね」 「私は植生解析を」 「ああ，植生解析のデータ解析，あれこれやっかいですよね …… なになに DCA? ありゃーヒドい発明ですよね」 と， どういう分野でもデータ解析をネタにした会話が成立するので．
• というかんじで， 二項分布を使った統計モデルとその最尤推定で終了． 投影資料としては個体差だの階層ベイズモデルだの用意したんだけど， そこはカルく流してまたデータ解析雑談しつつ終了． R とか統計モデルについて質問いろいろでるようになったので， 当初のプランはどうでもよくなった．
• しかし最尤推定とかひさしぶりに丁寧に説明してみると， 今までの自分の説明とか， かなり手ぬきなような気がしてきた． 院生あいての授業ってのはかなり丁寧にやらんといかんね …… 来月アタマから始める今年度の 統計学授業， いったいどうなることやら．
• さいきん， 「自分主導で話す」「長時間話す」 といった機会がすごくすくないので， 80 分ほどしゃべったあとはアタマへろへろ状態になる ……
• へろへろ
• Ubuntu 10.04 と ThinkPad X31 の組みあわせで， Alt+Function key で suspsend/hibernation できなくなる件 …… 代替策としては， sudo pm-suspend だの sudo pm-hibernate などすればよいのだが， いちいちパスワードいれるのがめんどう． 検索 してみると， sudo visudo で /etc/sudoers の末尾に %kubo ALL=(ALL) NOPASSWD: /usr/sbin/pm-suspend,/usr/sbin/pm-hibernate などと追加すればよさそう． 依然として sudo つける必要あるけど， パスワードは不要になる．
• ちょっと涼しくなったような気がしたので， 冷房をとめて窓を全開にしてみた …… しばらくすると汗がでてきてぐったりしてしまった …… たしかにマシになったんだけど， じつは外気温はまだ 30 度ちかかったので．
• なんとか， 気をとりなおして 岩波統計モデル本 作業のつづき．
• なかなかススまん ……

• 本日の岩波本修正作業， 全 265 ペイジ (修正前のペイジ数) 中の 26 ペイジまで終了．
• 1805 研究室発． Aeon 桑園で買いもの． 1840 帰宅． 晩飯の準備． 晩飯．
• 2220 自宅発 ミニ大通公園 走． 2305 帰宅． 体重 71.0 kg．
• [今日の運動]
  • ミニ大通公園 走 45 分間．
• [今日の食卓]
  • 朝 (0740): 米 0.6 合． タマネギ・油揚の味噌汁． 鶏レバのショウガ煮． キャベツ・ネギ・ベイコンのサラダ．
  • 昼 (1200): 研究室． 食パン． 鶏レバのショウガ煮． キャベツ・ネギ・ベイコンのサラダ．
  • 晩 (1930): 米 0.8 合． ニンジン・タマネギ・ジャガイモ・鶏肉のスープ． ホッケ． ダイコン・コンニャク・さつまあげの煮物．

2010 年 09 月 02 日 (木)

• 0710 起床． 朝飯の準備． 朝飯． コーヒー． 0910 自宅発． 晴． 今日も真夏日とか． 0925 研究室着．
• なかなかぎょーむ起動しない ……
• 総務係にあれこれお願いするメイルだす． ぐったり． 総務係の担当のかたはたいへんに親切なんだけど， この方面とのやりとりするのはツラい． この逡巡はおそらくのところ， 魑魅魍魎なる北大事務体系にアタマつっこむ勇気が必要とされるためのもので， それは言ってみれば路上で劇薬かぶったので やむをえずドブ川で洗顔するがごときものなのかも．
• ようやく， 岩波統計モデル本 作業のつづき． 第 2 章の最尤推定まわり． 文章はぜんぶかきなおし， といった状況で ……
• 昼飯．
• 残念ながら認めざるをえないのだが， 冷房つけてたほうが仕事のはかどりは良い． 今日はぎりぎりで真夏日回避されたんだけど．
• 第 2 章の確率分布・最尤推定まわりの書き直し， 雪山ラッセル (ざんごう掘り) 的な進捗 ……

• 本日の岩波本修正作業， 全 265 ペイジ (修正前のペイジ数) 中の 33 ペイジまで終了． あともうちょっとで第 2 章の修正終了なんだけど．
• 1820 研究室発． 1845 帰宅． 晩飯．
• [今日の運動]
  • うんどう休養日 ……
• [今日の食卓]
  • 朝 (0820): 米 0.6 合． ニンジン・タマネギ・ジャガイモ・鶏肉のスープ． ダイコン・コンニャク・さつまあげの煮物．
  • 昼 (1230): 研究室． 食パン． ナッツ類．
  • 晩 (1920): 近所の柿八でおいしいお寿司をいただきました．

2010 年 09 月 03 日 (金)

• 0700 起床． 朝飯の準備． 朝飯． コーヒー． 0910 自宅発． 曇． 0925 研究室着． 雨になった．
• 午前は手がき書類に苦しんでオワった． 今後しばらくはもう手書きしたくない ……
• 昼飯． 雨．
• 午後もそのあたり少しだけ苦しんで ……
• 岩波統計モデル本 作業にたどりつく． 第 2 章修正をオワらせる努力 ……
• とりあえず参考文献の部分以外は， 第 2 章をひととーりなおしてみた．
• また手書き問題のつづき ……
• 第 3 章の修正作業にとりつく．

• 本日の岩波本修正作業， 全 265 ペイジ (修正前のペイジ数) 中の 37 ペイジまで終了．
• 1820 研究室発． ちょっと雨． 買いもの． 1840 帰宅． 晩飯の準備． 晩飯．
• 2210 自宅発 ミニ大通公園 走． さすがに今晩はかなり涼しい． 2250 帰宅．
• [今日の運動]
  • ミニ大通公園 走 40 分間．
• [今日の食卓]
  • 朝 (0800): 米 0.6 合． タマネギ・ブナシメジの味噌汁． ダイコン・コンニャク・さつまあげの煮物．
  • 昼 (1230): 研究室． 食パン． バナナ．
  • 晩 (1930): 米 0.8 合． ダイコン・ニンジン・ネギ・ブナシメジ・豚肉の味噌汁． コマイの干物． キュウリ・コンブのポン酢あえ． ダイコン・コンニャク・さつまあげの煮物．

2010 年 09 月 04 日 (土)

• 0710 起床． 朝飯の準備． 朝飯． コーヒー． ちょっと怠業．
• 明日からの東京出張， いろいろと情報の確認・準備とか．
• またちょっと怠業ぎみに ……
• 昼飯．
• 1250 自宅発． 晴． 今日は涼しい …… 結果として最高気温は 28 度だったけど． 1310 研究室着．
• 今日は月曜日の発表の修正作業を …… と思ったんだけど， Ubuntu 10.04 + Xfce4.6 の設定にはまりこんでしまった ……
• 私は Xfce4 をできるだけ簡単に使いたいと考えてるんだけど， Session 起動時の挙動がどうにもコントロールできなくて …… 勝手に xfdesktop だの nautilus が起動して， ですね．
• 結果からいえば， これは session の自動保存とかが問題であった． xfce4-session-setting を起動して
  1. ログアウト時の自動セッション保存をヤメる
  2. セッションタブに移動
  3. よけーなプログラム (xfdesktop とか nautilus) を削除 (開始時に起動しておいてほしいプログラムだけ残す)， 「セッションの保存」
  4. いったんログアウトして， ちゃんと意図どおりになっているかを確認， だめなら 1. にもどってやりなおし
  でなんとかなるとわかった．
• ついでに呪われ nautilus (この不死的設定はどーなってるのか?) をうっかり起動したときに， ヘンてこな「ですくとっぷ」が表示されて気分わるくならないように， gconf-editor で /apps/nautilus/preferences/show_desktop を false にしておく．
• わかってみれば単純なことばかりだが …… 時間を空費してしまった． 時刻はすでに 17 時．
• ついでにいうと， 昨日・今日と「Gnome 環境でなければマトモに起動しない (?) Ubuntu 10.04 の evince (PDF viewer) の (みかけのうえでの) 自律起動」 なるわけのわからぬ研究に時間をとられてしまった． 原理は簡単で， epdfview とか epiphany みたいな Gnome (Gtk?) なソフトウェアを起動しておけば， evince も起動可能という現象を悪用するもの． 将来のど忘れ & スクリプト消失にそなえて， その evince 起動スクリプト記録しとくと， たとえば，

  #!/bin/bash
  (epdfview &) && sleep 3 && (evince $1 2>/dev/null &) && sleep 1  && killall -KILL epdfview
  

  といった呪文でよい． 徹頭徹尾ばかばかしいものであるけど． やたらと sleep いれてるのは， こうしないとなかなかうまく動作してくれないから．

[ですくとっぷ]

とゆーことで， 左の Ubuntu 10.04 な ThinkPad X31 (5-6 年前の機械) は evince で投影資料 PDF ファイルを表示しつつ， 右側では「何もない背景」の上でターミナルひらいて， その投影資料の内容を vim で編集しているところ． X31 のスクリーンの一番下にならんでるボタン類が 「Xfce4 つかってますよ」 を示しているわけだが ……

• つかいまわし投影資料， いろいろと不備がみつかるので修正していく．
• 今日は岩波本修正作業， ススまなかった ……
• 1810 研究室発． 買いもの． 1840 帰宅． 晩飯の準備． 冷蔵庫の中はこれでほとんどカラになったかな …… 晩飯．
• 2220 自宅発 ミニ大通公園 走． 2300 帰宅．
• 明日から水曜日まで猛暑日・熱帯夜の東京か ……
• [今日の運動]
  • ミニ大通公園 走 40 分間．
• [今日の食卓]
  • 朝 (0740): モンクールのパン． ベイコンエッグ． プチトマト．
  • 昼 (1220): 自宅ちかくの 風香 でおいしい昼定食いただきました． イワシがとてもおいしかったんだけど， 白煮で 6 時間・味つけて 2 時間も煮てつくるとk ……
  • 晩 (1910): 米 0.8 合． ダイコン・ニンジン・ネギ・ブナシメジ・豚肉の味噌汁． サバの塩焼き． キュウリ・コンブのポン酢あえ．

2010 年 09 月 05 日 (日)

• 0700 起床． 今日から水曜日まで猛暑東京出張． 0720 自宅発． 晴． 涼しい． 0738 JR 札幌駅． 東口のパン屋で朝飯調達して， 0750 快速エアポートで同発． 朝飯．
• 0826 新千歳空港着． 0845 ADO14 便 16F に搭乗． 日曜日朝なのにけっこう満席． 0900 taxing 開始． 1025 羽田に着陸． わざわざ述べるまでもなく， 東京暑い． どこか熱帯の国にきた， という想像して気温については不平をいわぬことにする．
• 京急にのって， 1049 羽田空港発． 1124 日本橋着． 1128 東西線に乗りかえ て同発． 1138 早稲田着．

[蕎麦屋・三朝庵で昼飯]

地下鉄早稲田駅でてすぐ． 最初に食券を買う方式． 大ざる 800 円 …… 大もり (700 円) にするつもりだったのに， まちがえた． おそばは …… まあ， わたしはそば通ではないので， おいしくいただきました． やかんで蕎麦湯でる， ただしウスい． 中国語を話す店員さんがいたり．

• 1210 早稲田大学 7 号館 統計関連学会連合大会 会場着．

[早稲田大学]

大隈講堂． このへんには大学院生のころ， 何度かきたことある． IGBP-TEMA がらみのよくわからぬ会議に， なんの関係あるのかよくわからぬ私が出席させられて， といったあまり楽しくない ……

• 今日は チュートリアルの日， 一度ぐらいは数理統計なヒトのベイズ講義を聴いておこう， ということで駒木文保さんの 「ベイズ理論の現在」 に参加しました．
• 13 時から 3 時間． 会場には 150-200 人ぐらいの参加者が …… これぐらいの大人数になると， 私の勝手な憶測ですが， 「何がなんだか?」 なヒトが多数派だったんじゃないかな …… 駒木さんのハナシはホントに 「ベイズ理論の現在」 な数理統計もので， ここ数年にわたってご自分が証明・論文化された 内容だったので．
• 私はといえば， 何をやりたいのか， なぜそれをやりたいのか， どのような結果がえられてどのように改善されたか …… ぐらいまではわかったけど， 詳細といいますか情報幾何のハナシはほとんどわかりませんでした．
• 全体の話のくみたてはわりと簡単で，
  • 皆さん， 「無情報」な事前分布といえば Jeffreys の事前分布つかってるけど， ホントにそれでいいのか?
  • 事前分布の良さを リスク関数 (Kullback-Leibler divergence の期待値) で評価しつつ， Jeffreys な事前分布よりマシなやつを調べましょう
  といったもの．
• 前半のハナシはこういう道具だての説明， そして情報幾何を駆使して (モデル多様体上の Green 関数にもとづく事前分布がうんちゃらかんちゃらから) Stein の事前分布がつねに Jeffreys よりリスク関数が小さくなることが 証明されて， といったところまで．
• 休憩時間に質問にいって， Stein の事前分布ってのは Jefferys のとはどうちがうんでしょう， といった質問に行ってみたのだけど， あまりにも質問が素朴すぎたのか Stein 事前分布は常に…… という上述の駒木さんが証明された定理をもう一度説明してくださって ……
• ここでしつこく踏みとどまって， 途中で言及されていた正則化やペナルティーのハナシと， その Stein の事前分布の関係はどうなっているのでしょう， とお尋ねしてみたところ， たしかに関係ないわけではなく， ペナルティーの項をつけた場合と得られる結果は似ているが， ペナルティー項方式では常にリスク関数が小さくなるといった保証はなく， そのあたりが Stein とは異なる …… と答えていただいたので (言うまでもなく， このあたりは久保脳内変換後のものなので， あまりマにうけぬよう)， まあとりあえず私のあてずっぽう (これって多数の似たようなパラメーターをベイズ推定するときに， 似たようなパラメーターをことごとく無情報事前分布にするのは良くない， ってことか?) とはちょっとハナシがちがうような気がした．
• 後半はさらに新しい論文の結果の解説で， 多変数ロケイションモデル (Neyman-Scott model を一般化したもの) のベイズ推定の事前分布に関して Jefferys よりマシなものがあるか (というかあるのはわかっているんだけど， どう証明されるか)， といったハナシ． またモデル多様体とか双極空間とか …… えーと， 上記モデルは群モデルとして定義できて， 群モデルといわれても行列形式で定義されるとしかわからなかったんだけど， これの左不変事前分布と右不変事前分布がそれぞれ Jefferys の事前分布とそれより優れた事前分布に対応している …… といった内容だとおもいました …… うー
• さらに， 2009 年論文の， ローレンツ錐ガンマ分布モデルにおける Jefferys 分布に優越する事前分布のハナシは …… うー
• 最後におまけで， モデル選択規準 BIC がマズくなる条件の検討といった独立したトピックの 紹介があり …… えーと私の理解では， 標本数が同じであっても推定値の大小によってモデル選択の結果が変わる場合があり， それは BIC の高次の項をみればわかるんだけど， そうすると BIC に事前分布が含まれてしまってめんどうになり， Lindley's paradox (ここがよくわからん) とも関係するのだが ……
• まあフォローできなかった話もおおかったけど， いろいろ勉強になりました． 1602 東西線にのって地下鉄早稲田発． 1607 飯田橋着． ふらふらとさまよったり， 途中のスーパーマーケット (ポロロッカがあるはずだったのだけど， 撤退してマルエツに代わっていた) で晩飯の買いものして， 1633 本日の宿， ウィークリーマンション飯田橋 着．

[今回の宿]

簡単な自炊設備・鍋・調理器具・皿つき． タオルや歯ブラシもついてる． そうじやシーツ交換なし． 宿泊費はなぜか日変動するんだけど， 三泊で 1.8 万円ぐらい． 一人で外食するのが好きではない私としては， こういうところで晩飯たべるのが気楽なんだが …… 考えてみれば明日・明後日の晩飯は外になるなぁ．

• どーにもくつろいでしまって， ごろごろころがって本とかよんだりして， あしたの発表の修正とかする気になれない ……
• [今日の運動]
  • うんどう休養日 ……
• [今日の食卓]
  • 朝 (0800): 快速エアポート車内でパン．
  • 昼 (1150): 三朝庵のソバ．
  • 晩 (1930): 宿の自室． マカロニ． 野菜スープ． トマト． バナナ． ヨーグルト．

2010 年 09 月 06 日 (月)

• 東京出張二日目． 0630 起床． ねむい． 朝飯． コーヒー．
• 今日は 10 時までに会場につけばよいので， 自分の発表内容のみなおし．
• うおー， ぎりぎりまで作業してしまった． 0930 宿発． 地下鉄東西線にのって， 0940 飯田橋発． SUICA なのでラク． 0946 早稲田着． 0951 会場着．
• 10 時から 企画セッション：環境・生態データのモデル化と解析 (1)． 私も 11 時からいつものよーなハナシをして終了 ……
• どーにもぐったり感あるので， 撤退． 地下鉄東西線にのって， 1207 早稲田発 1212 飯田橋着．
• このふきんは， 出版社や印刷会社が多いような． TOPPAN の巨大なビルが四方を睥睨してるし．
• とにかく「暑い」ということを考えない努力． ここはどこか熱帯の国であり， 暑くてあたりまえといった空想．
• 薬局とかスーパーとかによって， 1233 宿に帰着． 昼飯．
• 午後はごろごろしてました． 雑用をちょっとだけススめた．
• 夕方， うろうろ．
• [今日の運動]
  • うんどう休養日．
• [今日の食卓]
  • 朝 (0700): 宿の自室． マカロニ． 野菜スープ． バナナ．
  • 昼 (1250): 宿の自室． 食パン． ローストビーフ． ポテトサラダ． トマト． バナナ． ヨーグルト．
  • 晩 (1830): 早稲田ちかくでいろいろおいしいものを食べました．

2010 年 09 月 07 日 (火)

• 東京出張 3 日目． 0750 起床． 朝飯． コーヒー． 0900 宿発 0909 飯田橋発 0914 早稲田
• 午前中は環境のセッションで． 「左右非対称」の角度分布は不思議なかんぢだ ……
• 今日も昼に撤退． 本屋で目的もなく一時間ほどすごす …… 本屋でこんなにのんびりするのは， すごくひさしぶりだ．
• また買いものして， 13 時すぎに宿に帰着． 昼飯．
• 今日は何か仕事しようと思ったのに， ごろごろ．
• 召集をうけてうろうろ．
• また宿にもどってごろごろ．
• また召集をうけてうろうろ．
• [今日の運動]
  • うんどう休養日， つづく ……
• [今日の食卓]
  • 朝 (0820): ウィークリーマンションの自室． 食パン． トマト． バナナ． ヨーグルト．
  • 昼 (1330): ウィークリーマンションの自室． 食パン． トマト． バナナ． ヨーグルト．
  • 晩 (2030): 江東区でおいしく晩ごはんいただきました．

2010 年 09 月 08 日 (水)

• 東京出張 4 日目． 0730 起床． ねむい． 朝飯． コーヒー． 怠業．
• だらだらと撤収の準備をして， 0935 宿発． 雨．
• 0944 地下鉄東西線に乗って飯田橋発． 0953 日本橋着． 都営浅草線に乗り換えて， 1005 同発． 1038 羽田空港着． 雨．
• 空港内うろうろ． 第二ターミナルは B1F の田辺書店以外にはめぼしい本屋がないような． また， 電力供給が三相なやつなので， 電力窃盗とかも難しい． そして使える無線 LAN がほとんどない．
• で， 保安ゲイトくぐって， いつものよーに第二ターミナルビル「岬の突端」までいくと， ISDN グレイ電話が除去されたあとブースでコンピュータに電力供給できたり …… ThinkPad X31 の部品はいろいろ交換したんだけど， 電池は交換してない 6-7 年もので， これによる稼働時間は 20 分ぐらいなんだよね．
• あいかわらず使える無線 LAN はない． スターバックスもダメだった．
• ちょっとだけ作文ぎょーむをススめる．
• 1250 搭乗． ADO 95 便 8A 席． 羽田空港豪雨． 離陸まちがツマってるので， 窓の外を観察してみると …… なんと， B 滑走路 (北東-南東のやつ) が着陸専用で， 海ぞいの C 滑走路が離陸専用に …… JAL 機が C まで来てるので A 滑走路は閉鎖してるんだろうな． 羽田空港には何回もきてるんだけど， B 滑走路が使われるところは初めて見た ……
• 14 時ごろ新千歳空港着陸． 快晴． すずしー! 空気も乾いている!
• 1419 JR 新千歳空港発． 1455 JR 札幌． 1508 JR 桑園． AEON 桑園で買いものして， 16 時前に帰宅．
• かたづけ． 晩飯の準備． 晩飯． 至近距離にある温泉， 「極楽湯」でのんびり．
• [今日の運動]
  • うんどう不足の日々， つづく．
• [今日の食卓]
  • 朝 (0820): 宿の自室． 食パン． ヨーグルト．
  • 昼 (1330): なかなか離陸しない ADO 機中で， ベイグルサンドイッチ． 空港でかったもの．
  • 晩 (1930): 米 0.6 合． ブナシメジ・豆腐のすまし汁． タマネギ・サケのホイル味噌焼． キャベツ・キュウリ・ネギ・ハムのサラダ．

2010 年 09 月 09 日 (木)

• 0740 起床． お， 今朝はかなり気温下がった …… ただし結果的には最高気温 25 度ごえなので夏日 (そして夏日がつづくという予報)． 朝飯の準備． 朝飯． コーヒー． 0910 自宅発． 晴． 0930 研究室着．

• 東京ぼけで， 家に ThinkPad X31 置き忘れてしまった …… ココロを落ちつけるために， まずはお茶をいっぱい飲む． 自宅まで取りに帰る．
• この午前中における災難， まだまだオワらない …… A 棟 7-8F DHCP 接続， 全滅．
• 原因を特定するまで時間かかってしまった． 闇ルータこと IO DATA NP-BBRS がとうとうぶっ壊れたのか …… と思ってさらに子細に調べてみたら …… 原因はそこではなく， なンと言いますか， 「犯人は探偵自身でした!」 ぱたーんだった．
• 私が東京出張のおともに連れていった， Planex のちび無線 LAN ルータ (兼 hub)， A 棟 7-8F ネットには無線 LAN hub として接続していたハズなのに ……

[勝手に IP アドレスくばってたヤツ]

Planex GW-MF54G2， なぜかしら無線 LAN として働かせていたのに， 勝手に「はーい，君はこの IP アドレス使いなさい」 とか研究室じゅうの端末に勝手に IP アドレスくばって大混乱をひきおこした． LAN 内に複数の DHCP server があるときに， なぜかしら片方だけが機能してしまうとゆー， この不思議さよ ……

• とりあえず， 原因が特定できたのでネットから外して後で調べることに． 私の研究室は院生部屋からの電波もひろえるので， 本日はそちらに接続．
• 夜に自宅でこのちび Planex を調べてみたのだが …… 原因がよくわからん． 無線 hub もーど (access point mode) でも DHCP server が動作していたのだろうか? というナゾの状態なので， とりあえずどのモードでも DHCP server が機能しないようにして， 明日また実験してみる．
• データ解析こんさる．
  • 世の中には Kenward-Roger method なるものがあり， これは linear mixed model (LMM) の (fixed effects な) 係数の Wald 統計量な検定に使われる
  • 標本数が多ければ LMM でも Wald 統計量でいーんじゃないの， といったハナシになるわけだが， 標本数が少ないときに 上記係数がゼロからずれやすいといった問題が生じるらしい
  • なぜ問題が生じるかというと --- 私の憶測だけど， ようするに random effects が大きいにもかかわらず， それを無視したモデリングをすることによって， fixed effects な係数がそのへん代替してしまうんだけど， サンプル数が少ないときでも同じような問題が生じるのかもしれない
  …… といったところかな? 蛇足ながら， SAS のなんちゃらぷろしじゃの中には上記 KR 補正法をつかった LMM (GLMM?) 推定なやつがあって， ブラックボックス SAS 使いなヒトたちもそういうのを使ってみるとゆーのが， 一部のぎょーかいでは「お作法」になってるそーで．
• 午後は雑用あれこれ．
• 今日は岩波本作業に着手できなかった …… 1810 研究室発． 買いもの． 1840 帰宅． 晩飯の準備． 晩飯．
• ちびルータのとりしらべ， とか．
• 2230 自宅発 散歩． 2315 帰宅．
• [今日の運動]
  • 散歩 45 分間．
• [今日の食卓]
  • 朝 (0820): 米 0.6 合． ブナシメジ・豆腐のすまし汁． キャベツ・キュウリ・ネギ・ハムのサラダ．
  • 昼 (1330): 研究室． 食パン． ナッツ類．
  • 晩 (1940): 米 0.8 合． ナメコ・豆腐の味噌汁． ダイコン・鶏肉の煮物 …… こういう煮物では鶏肉とかはちょっと煮たらすぐに取りだすのがいいのかも． キャベツ・キュウリ・ネギ・ハムのサラダ．

2010 年 09 月 10 日 (金)

• 0740 起床． 朝飯の準備． 朝飯． コーヒー． 0900 自宅発． 曇． 0910 研究室着．
• いきなり， 岩波統計モデル本 作業にとりくもうとするが …… あれやこれやと ……
• なんとか 11 時すぎに岩波本にとりつく．
• 昼飯．
• 昨日の Kenward-Roger 補正法のつづき， ちょっとだけ． 生態学とかまわりでの使われかたについては， Bolker et al. 2008. Generalized linear mixed models: a practical guide for ecology and evolution など， さらに LMM での type-I error について調べたのが Schaalje et al. 2002. Adequacy of approximations to distributions of test statistics in complex mixed linear models． やっぱり GLMM で使うのはヤバいんじゃないかな ……
• 岩波本作業のつづき． GLM 章 (第 3 章)， わりと最近に書きなおしたハズなのだが …… いや， 最近の書きなおしは第 4 章からで， 第 3 章の書きなおしやってたのは 今年の初め だったのか ……!!

• 本日の岩波本修正作業， 全 265 ペイジ (修正前のペイジ数) 中の 40 ペイジまで終了．
• 1820 研究室発． 家のカギわすれてうろうろ． 1900 帰宅． 晩飯．
• ぐったりぐー
• 2240 自宅発 ミニ大通公園 走． 2310 帰宅．
• [今日の運動]
  • ミニ大通公園 走 30 分間．
• [今日の食卓]
  • 朝 (0800): 米 0.6 合． ナメコ・豆腐の味噌汁． ダイコン・鶏肉の煮物． キャベツ・キュウリ・ネギ・ハムのサラダ．
  • 昼 (1230): 研究室． 食パン． ナッツ類．
  • 晩 (1940): 米 0.8 合． ネギ・ブナシメジ・油揚の味噌汁． コマツナ・油揚の煮びたし． キンメダイの一夜干し．