ぎょーむ日誌 2004-11-(21-30)

2004 年 11 月 21 日 (日)

• 0930 起床． 朝飯． コーヒー．
• 怠業 …… してる場合ではないんだが． すでに 1320． 自宅発の北大構内走 1330-1430． 晴． 走ってるあいだに昨日議論したダケカンバ葉群分布の推定を， MCMC 法に帰着させる方法がまとまってくる． 体重 73.2kg． 重い．
• 1600 自宅発． 曇． 1610 研究室着．
• また苫小牧な問題にとりくむ． 苫小牧ボスは土日ナシ状態のよーで ……
• 昨年の M2 で今は相模国で働いてる大澤君がきたので雑談．
• 2020 研究室発． 大澤君その他合計 9 名で晩飯を食いにでる． 北大東方に集中分布している 「スープカレー」屋たちにはことごとく却下される (スープ切れ & 満員御礼)． 北 14 の新しそうな中華料理屋に入る． 2200 帰宅．
• 明日は朝から苫小牧へ．
• [今日の運動]
  • 北大構内走 1330-1430． ストレッチング．
• [今日の食卓]
  • 朝 (1030): 米麦 0.7 合． ネギ・豆腐の電磁波蒸し． ダイコンのキムチの素づけ．
  • 昼 (1450): 米麦 0.7 合． ネギ・豆腐の電磁波蒸し． 丸天． ダイコンのキムチの素づけ．
  • 晩 (2050): 北 14 西 4 に新しくできた (と思われる) 学生むけふう中華料理屋「瀋陽飯店」． これまた学生向けコース料理 (6 品) と思われる 800 円コース． かなりおてごろか． 味つけは塩がききすぎなんだけど．

2004 年 11 月 22 日 (月)

• 0600 起床． 朝飯． コーヒー． 0635 自宅発． 曇． 0645 JR 札幌駅着． 0649 快速エアポートに乗って同発． 0724 南千歳で下車． 0730 東室蘭いき普通に乗って同発． 0752 苫小牧着．
• さて， ここから研究林までどう行くか， だが． 出頭時間まで一時間ほどあるので， 歩いてみることに． 地図がないのでやや不安なのだが …… 0852 苫小牧研究林着． 日浦さんによると， 駅からここまで歩く偉業を達成したのは， 竹中さんについで二人目の快挙だそうである．
• 0900 年輪な勉強会 …… といっても鍋嶋さん・日浦さん， 信州大の安江さん， そして私の 4 名である．
• まず 1 時間ほど苫小牧ボスこと鍋嶋さんによる， 苫小牧モデルの現状説明． いつものごとく要点を押さえたものである． 過去一週間ほどのあいだに「劇変」した 推定計算方式とそこから出てくる結果のいくつかについては， 私が少しばかり補足説明． 安江さんのご指摘によると， やはり年輪年代学方面でも似たようなモデリングは 検討されているよーで ……
• つぎに 安江さん による年代気象学の紹介． 解説の題材としては， 道北のヤチダモとシベリアのカラマツ． たいへんわかりやすいものだった．
• 環孔材なヤチダモについては 孔圏 (毎年の年輪の内側，道管がとおってるところ) と孔圏外 (外側) のどちらの挙動が直径成長を決めているか， というもの． なお， 環孔材といっても開葉期に新しい道管が機能してるわけではなく， しばらく時間がかかって使えるようになるらしい． つまり， 最初は道管以外のところに水を通しているらしい． なお， 道管は毎年つかいすて．
• クラスノヤルスク北方のツラで行われているカラマツの 成長解析では斜面の南北で応答してる気象要因が異なるのでは， といったあたりを調べるもの． デンドロバンドをつけてみると， 直径成長は 6 月に局限されているらしい． 成長速度おそく， 100 年たっても DBH 10-15cm ぐらい． 100 年に一度ぐらいの火事でで， 地表のコケ類消失し， 夏場の永久凍土層が深くなり， 更新可能になるとのこと．
• 年輪な世界でよく使われてる手法．
  1. cross dating: 年輪の年を特定する
  2. 「標準化」 : 年輪の時系列データの ノンパラメトリック回帰 (spline?) を行う
  3. 標準化曲線からのずれと， 気象要因の相関を調べる
  問題は 2. と 3. だな． 二重推定． また 3. は検定をつかっているんだけど， 多重性とか考慮していなくて， とにかく相関あるところを列挙する． 推定の計算の前提はどこも線形かつ等分散正規分布． ちかごろは， このあたりを改良する手法がいろいろと提案されつつあるらしい．
• 1200 ごろ終了． 昼飯にでる． 1330 研究林もどる．
• 鍋嶋さんが安江さんを研究林内に案内されるので， 私も金魚のふんのごとくついていく． 苫小牧ボスぢきぢきの御案内で， 畏れおおいことです …… 研究林南東から反時計回りにまわるようなコース． ミズナラの伐採・施肥実験区， ミミズ enclosure 実験区 (TOEF 的人海戦術で莫大な量のミミズが除去されていたり)， 広大なるエゾシカ包囲柵 (ただし内部には「少なくとも一頭」ぐらいしかいないらしい)， 台風によって倒壊した帯状植林区画 (シラベなんかも植えられていたり ……そして急激にまばらになったちょー高密度区の 生き残り個体は ぱいぷ樹木 的な枯れあがりでよかった)， 倒れ flux tower， そしてクレーンサイトで 現在のモデリング対象個体の生き残りを見物したり．

[ミミズ除去実験区]

苫小牧研究林は なぜかしら日本国内でもかなりミミズ密度の高い森林らしい． このミミズの働きを調べるため， ミミズのぞいた実験区を防御壁で囲ってる様子． かなりたいへんな作業だったそうで．

• 日没がすっかりはやくなってしまった …… 1650 研究林にもどる． 院生部屋で年輪雑談． Ultimate Tree-Ring Web Pages なるデータベイスがある， と教えていただく． 追記． Fritts や Shashkin が関与したという， 非線形な (?) 年輪解析モデルについて， あとでこのあたりで調べておくこと．
• 年輪勉強した皆さんに加えて， TOEF PD の豊田さんとその師匠で こちらに調査に来ている横国大の 金子さん とともに， 苫小牧駅南側二条通の小さな居酒屋「釧路」 (なんて名前だ)． 1900 ごろ到達． 土壌動物ばなしなど楽しく．
• そろそろ帰りの汽車が無くなってくるような気がしたんで， 中座して 2030 居酒屋発． 2040 JR 苫小牧駅． 2048 普通に乗って同発． 札幌まで 1410 円． この時間帯だと札幌方面いき普通は 1 時間に 1 本ぐらい． 2116 南千歳で快速エアポートに乗り換え． エアポート混んでる． 2151 札幌着． 2210 帰宅．
• [今日の運動]
  • しいて言えば， 朝の苫小牧市縦断
• [今日の食卓]
  • 朝 (0610): 米麦 0.7 合． ネギ・豆腐・卵の電磁波蒸し．
  • 昼 (1250): 幌内川ぞいの居酒屋， 五合 (はんじょう) の昼定食． シマホッケ．
  • 晩 (1910): 苫小牧駅南の居酒屋「釧路」． 刺身・鍋．

2004 年 11 月 23 日 (火)

• 0840 起床． 朝飯． コーヒー． 怠業 …… いかん怠業してる場合ではない．
• 1110 自宅発． 曇． 1120 研究室着．
• 自由集会 関連のメイル書き． たしかに面白いものになりそうだけど， なかなかたいへん．
• 時刻はすでに 1350か． 昼飯．
• さっさと統計学授業の準備を …… と思っているのにメイル書きぎょーむつづく． 母子里ダケカンバ関係の 拡散型 GaAsP フォトダイオード G1118 の方向特性の図 (in 腐れわーど) を 浜松ホトニクス→加藤京子さん経由で送っていただいたのだけど …… ヘンてこな図なんだけど， よくよく見るとこれって単なる cos(θ) にすぎないんでは? という気がしてきた．
• ぢりぢりたるメイル書きばかりで， 進捗感ナシ． お茶部屋に逃げると， 種子の処理に疲れきった女性大学院生が 「私は 1g 未満の種子がこの世に存在することが許せません」 と立腹しておられるし (しかし， 伺うところかなりたいへんそうな作業のようで)．
• 1700 ごろから井田君の売りこみ面接の内容を拝聴することに …… 結局， こういうのって 「どう改良すればよいか」 を検討するため 「そのあたりのもっと詳しい内容を」 ということになり 「それって他にもいろいろ場合わけできるぢゃん」 とかハナシが限りなく発散していくのが常である．
• で， 発散を押さえこんで， なんとかコメントをまとめるとすでに時刻は 20 時すぎか．
• またメイルぎょーむ． 日本動物行動学会大会 の準備でおいそがしそうな粕谷さんをさらにいろいろとわずらわせる．
• 間瀬さんにもとっとと返信しなければならないんだけど， 返信文案考えると， なぜかしら統計学勉強が始まってしまい， 文章はまったく進捗しないという悪しき状況である．
• 明日の授業の準備がぜんぜんススんでいない． 昨年の使いまわしは避けられないけれど， せめてもうちょっとマシにしたい， という気分なんですが． ばててきたのでいったん撤退． 2200 研究室発． 2210 帰宅． またいろいろとメイル書き．
• 2400 すぎに晩飯． 授業の準備はどこにいったのやら．
• [今日の運動]
  • 運動してない ……
• [今日の食卓]
  • 朝 (0930): リンゴ． ヨーグルト．
  • 昼 (1350): 研究室お茶部屋． 「北欧」ベイグルサンドイッチ．
  • 晩 (2410): 米麦 1.0 合． ダイコンのキムチの素づけ． ネギ・ニンジン・豆腐のミートソース． これは残り野菜をきざんでレトルトパウチドミートソース で煮たもの． いままで全力運転しかできなかった壊れガスレンジで， 不完全ながら火力調整やる方法を発明した． さてさて．

2004 年 11 月 24 日 (水)

• 0730 起床． 朝飯． コーヒー． 0845 自宅発． 晴． 0855 研究室着．
• さっそく授業の準備を．
• まーた北大ネットが不調． 本日は学内 → 学外が途絶． DDoS 攻撃でもくらって腐れ proxy server が死んだんだろうな． 毎回毎回おなじよーなしくじりばかりで進歩のない連中， というか．
• 上の不調は昼ごろ復旧． 昼飯．
• ひたすら授業の準備 …… と思ったら， 三浦さんの ThinkPad が virus に感染? 北大の mail gateway server とやらもアテにならんな． それから北大の「強制通過型」proxy server が腐れていて， ときどき情報をとおさなくなる (denied 表示をだす)． これは動作がなぜか確率論的で， しばらく待ってから再アクセスすれば OK だったりする． めちゃくちゃだ．
• 1550 授業の準備完了． ファイルのアップロードとか． 1610 教室である C104 に移動． 液晶プロジェクターの準備．
• 1630 すぎ， M1 もっぱらあいての 統計学授業， 開始． 統計学的データ解析で何ができるか， 乱数と推定が重要， 割算すんな， というのがメインか．
• 1755 授業終了． ばてた． M1 あいての最初の統計学授業の内容としては， 昨年よりはマシだったかも …… 内容を削減したんで． 次は来週の月曜日． 配布資料は 35 部ぐらい必要だな (次回はこのリレー式講義の最終回なんで)．
• なんか今日は方々からやたらと統計こんさるメイルがくる． ひととーり返信を試みる．
• しかし途中でばててきたので 2000 研究室発． 2030 帰宅． 体重 73.8kg． 重い． 晩飯．
• なんかメイル書きがいつまでたっても終了しない． 自由集会とかその他とか．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0750): 米麦 0.7 合． ネギ・ニンジン・豆腐のミートソース．
  • 昼 (1255): 研究室お茶部屋． 米麦 0.7 合． ネギ・ニンジン・豆腐のミートソース．
  • 晩 (2210): 米麦 0.9 合． キャベツ・セロリ・エノキダケのカレー．

2004 年 11 月 25 日 (木)

• 0730 起床． 朝飯． コーヒー． 0830 自宅発． 雪． 0840 研究室着．
• 輪読会の予習． この第 3 章というのは， 話題限定ひきかえに丁寧な説明を実現していて， この本の中ではまっとうなほうだ．
• 1000 より輪読会 …… と思ったら， 小菅せんせーが出てきて黒板に何やら書いている． どうやら先週に担当した章で重要な部分を隠しすぎたことを 糾弾されたらしく， その「補足」とのことである． で， これまた一見すると隠匿内容を提示してるかにみせかけて， じつはホントに重要なポイントは隠している， ということがわかった． 私は小菅せんせーのこういう高等隠匿戦術には慣れているので， そのあたり指摘してみた． 状況不利とみるやただちに欺瞞説明を放りだし， いつものごとく見事なる遅滞・後退運動にて， まんまと逃げられてしまった． うーむ．
• で， Life-history evolution in plants 第 3 章， 担当は平林さん． 最適な対食植防御なるものが， ゲイム理論的状況では成立しない， というのがハナシの中心． しかし， 自分の損得勘定ばかり熱心な ひどく利己的な個体たちの足のひっぱりあい， という凶々しい状況について計算している， ということが信じられなかった善良なる大学院生 はそのあたりを理解しそこなっていたので， 私が横からあれこれと補足: 「だから， みんなが真面目に防御してるんだったら自分はさぼっても良くて ……」
• 1230 終了． 生態学会 大阪大会 の申し込み． まだ 500 人か． これから明日の夜までに 1000 人ぐらい申しこみそうだな．
• 1250-1330 北大構内走． 晴． 昼飯． 昼飯後に郵便局に大阪大会の金を払いにいく． 6500 円 + 70 円． 今回はなんとなく懇親会から逃げる， ということで． 6000 円だしなぁ．
• うう． 29 日 (月) は予備審 2 本立か． すけぢゅーりんぐが …… というか， その日の夕方も統計学授業なんですよね． 前日までに準備しろってことか． うーむ．
• さて， 自分の参加申し込みはやったけど， 自由集会の要旨作文はまだだ． ということで， 粕谷さん文章を改竄して要旨をでっちあげてみる． メイルで関係者に送信．
• 統計こんさる業 …… のはずが， こういう図を描いて遊んでいるだけだったり．

• といったことだけで終ってしまい， 1950 研究室発． 2010 帰宅． 体重 74.0kg． あいかわらず重い． 晩飯． ばててたのか， すぐに寝てしまった …… が夜中に目がさめてしまった．
• [今日の運動]
  • 北大構内走 1250-1330． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0750): 米麦 0.7 合． キャベツ・セロリ・エノキダケのカレー．
  • 昼 (1340): 研究室お茶部屋． 米麦 0.8 合． キャベツ・セロリ・エノキダケのカレー．
  • 晩 (2120): 蕎麦． ニラ・シイタケ・卵の炒めもの．

2004 年 11 月 26 日 (金)

• 0740 起床． 寝不足ぎみ． 朝飯． コーヒー． 0840 自宅発． 晴． 0850 研究室着．
• 苫小牧直径成長モデル， 推定計算プログラム in R の改造． 光合成気象値・代謝気象値を別 list にする推定プログラムで， これらの組み合わせを同時に選抜するルーチンを追加した． 試験運転もかねて， 「ねんのため」 計算開始．

  ……
  #     AIC = 2210.97
  # min AIC = 2071.37
  # --- wv.mb set = 7/27 ---
  # --- wv.ps set = 39/54 ---
  temp.ps: best=16.0/width=6.0
  ppfd: factor=0.02
  vpd: factor=1.00
      1998     1999     2000     2001     2002     2003     2004 
  26.65381 23.59560 22.93142 24.22508 25.75016 24.64961 25.81366 
  # AIC = 2218.24, np = 4+8 (scale.log const ps.prev mb6)
    scale.log       const     ps.prev         mb6 
  -0.53700289  1.22804230  0.01813444  0.15193820 
  #     AIC = 2218.24
  # min AIC = 2071.37
  # --- wv.mb set = 7/27 ---
  # --- wv.ps set = 40/54 ---
  temp.ps: best=17.0/width=6.0
  ppfd: factor=0.02
  vpd: factor=1.00
      1998     1999     2000     2001     2002     2003     2004 
  26.64290 23.88847 23.50196 25.19602 25.54211 23.64804 24.48687 
  ……
  

  てなかんぢで進行中． 27 × 54 = 1458 気象値セットしか調べてないのに， たぶん 7 時間ぐらい必要． 次はこの「同時選抜の高速化」 をやらんといかんね ……
• 苫小牧直径成長モデルの wiki page (すでに十数MB になってる) のりすとら少々． ときどきりすとらせんと見通しが悪くなる．
• いろいろと助言いただき， 今日中に自由集会要旨は提出できそう． うっかり忘れそうでコワい． お， 今回は講座の液晶プロジェクター (EPSON ELP-710, 1000 ANSI Lm) を運搬しなくてもいいかも．
• 苫小牧直径成長モデルの推定計算， やはりバグがでたか． 修正に 20 分ほど費す． 再計算開始 1147． 問題は R プログラミング作法にありそうな気がする．
• 本日の午前中は， かとーオフィスとなりの講義室から， 研究科長の何だかヘンなとーくが聞こえっぱなし． うーむ．
• というふうに昼飯前の時間はすぎていってしまった． 1230 北大構内走発． 曇天． 1310 もどる． 昼飯．
• で， いろいろやるべきことがあるにもかかわらず， それらを放りだして統計学勉強に没頭してしまう．
• 甲山さんのもちこんできたシュート構造三次元化こんさる． こういうのって， 発生過程を考慮した特徴づけをやらんと無意味だと思うんだが …… ぢつは現物をあつかってる人たちのほうがそういうのは無視して， どうでもいい平面投影角とかばかり見てしまうんだよね．
• Vine Linux 3.0 が最近更新されんなぁ， と思っていたら， 時代はすでに 3.1 に移行していた． /etc/apt/sources.list を変更して sudo apt-get update && sudo apt-get dist-upgrade すると， えらくたくさんの rpm が入れかわってしまった． X.Org まわりとか．
• 1700 から研究科のセミナーにでる． 本日は露崎さんでくだんの西オーストラリア話． あとから気づいたんだけど， 途中で示された 「DCA ではうまく天然林を特徴づけられない」 なる図は分類アルゴリズムがでたらめであるからではなく， 逆にその「特徴づけ」計算とやらの方法論的限界をあらわしていたわけだ． つまり， 植生学分野でよく使われている加重線形和では 植生なヒトたちの考えるタイプわけルールを記述できてない． 理由を簡単に言えば， そのルールが非線形であるからだ．
• 1847 自由集会 の要旨を送る． 自由集会ペイジも更新．
• 苫小牧計算， 6.5 時間ぐらいで終了か． 高速化せんといかんね．
• 2020 研究室発． 2030 帰宅． 晩飯．
• [今日の運動]
  • 北大構内走 1230-1310． ストレッチング．
• [今日の食卓]
  • 朝 (0810): 米麦 0.7 合． キャベツ・セロリ・エノキダケのカレー．
  • 昼 (1330): 研究室お茶部屋． 米麦 0.7 合． キャベツ・セロリ・エノキダケのカレー．
  • 晩 (2100): 米麦 1.0 合． キャベツ・セロリ・エノキダケのカレー．

2004 年 11 月 27 日 (土)

• 0900 起床． 朝飯． コーヒー．
• 怠業もーど …… にひたってるわけにもいかんので， 1420 自宅発の北大構内走． 雪． イチョウがいまごろ落果して， 腐敗臭ただよう． 1505 帰宅． 体重 74.0kg． すごく， 重い． 昼飯．
• 1610 自宅発． 雪． 1620 研究室着．
• いかん． 何も進捗しない． 何しにきたんだ．
• ホントに仕事がススまんので， 書籍の発注リストなど作ってみる． これはかなりの へっぽこモードでもススめることのできる作業である．
• …… だけで時間が経過してしまい 2220 研究室発． 2240 帰宅． 冷蔵庫の製氷室の氷とり． 以前はこれはピッケルで氷塊を粉砕していたわけだが， さいきんは手口を変えた． 氷のしたに洗面器をおき， 電気を切って数時間放置しとけば， とくに水難状態にもならずに氷がとれる． 今日は室温が低かったんで最後の氷塊はとれてなかった． これは手で外す．
• 晩飯． もう真夜中ぢゃん．
• [今日の運動]
  • 北大構内走 1420-1505． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0930): うどん．
  • 昼 (1530): うどん． 今日はうどんの日です (うそです)．
  • 晩 (2400): 米麦 0.8 合． モヤシ・セロリ・ネギ・シイタケ・卵の炒めもの． 壊れガスレンジは数日前に直った． たぶん配管系がいかれてたみたいなんだけど， 無理矢理つかっているうちに， どうにかなったようだ． また， 明日になったら壊れていても不思議ではないんだが．

2004 年 11 月 28 日 (日)

• 0900 起床． 朝飯． コーヒー．
• 怠業 …… ばかりしてるわけにもいかんので， 他のことはさておき， とどこおっていたメイル返信． まずは間瀬さんにいろいろと教えていただいた御礼 …… で， とても trivial な内容なのに 1 時間以上を費す． うーむ， われながら遅い．
• これまたとどこおっていた， glmmPQL() に関する質問 (のひとつ) の回答を検討する …… 結果に Log-likelihood とか出てまぎらわしいな． しかも summary() とると， ご丁寧にも AIC & BIC 計算している． ふーむ． ともあれ返信．
• glmmPQL() は一般化線形モデルを quasi-likelihood (QL) ふうに 定式化して， lme() (in nlme) に渡してるだけなんで計算とかは確実っぽく見える (う? PQL の pennalized はどこにいった?)． そうか， なんでパラメーター推定値に関して t 検定やるのか， といった疑問については nlme 方面を調べないとわからないわけだな． いや， 最大 QL 推定値の分布ってわかってるのかしらん?
• ともあれ， 滞納なメイルは …… あ， アカマツ作図がまだ残ってた． とりあえずこれは明日の夜以降， ということで．
• 1205 自宅発北大構内走． 晴れたり雪ふったり． 気温はかなり低い． 1250 帰宅． 昼飯．
• 1350 自宅発． 晴． 1400 研究室着．
• 時間もなくなってきたんで， 明日の統計学授業の準備．
• 自然科学について何ごとかを修めるべく入学してきた ここ地環研大学院生は， 2 年間の修士課程の中で「データをどう数理的にとりあつかうか」 について 2 回分の講義しか受けない …… 他ならぬ私のやる これ だけ． しかも全員が受講するわけでもない．
• 受講生の大半は学部のときにも統計学とかやっていない (少なくとも「理解できる統計学授業」は受けてない) ので， 内容的には高校生に教えるつもりでやらねばならない． しかも， 高校数学を忘れてしまった高校生あいてに …… 本邦ではなぜこれほどまでに 「データ解析なんて簡単，やればできる」 という信仰が流布してしまったのか， なかなかナゾなところだ．
• で， 前回は統計モデルと推定の関係を理解しよう， と強調したんで， 今回もしつこくそのあたりのいろいろな例を挙げて， 「データの傾向とばらつきを図にする → 推定の試み」 という基本動作 を浸透させたいところなんだけど …… まあ， 「検定」の説明します， と予告したのでその準備をせんといかん．
• 「検定」などという難しいことは， (「理解できる統計学授業」を一回しか聞いたことない) 彼らにはまだまだ時期早尚というべきなんだけど …… その指導教官たちは 「データを統計学的に解析しろ」 とは言わずに 「検定にかけろ」 などと乱暴な 「指導」しやがることが多いので， 私としては先まわりにしてテをうたずにはいられない．
• で， その「検定」の説明であるが …… 昨年は連続確率分布から得られた標本に対する検定なるものを説明した． ありていに言って， 連続確率分布は説明するのもされるのも難しい． そこから得られる標本統計量の分布もまた連続分布なので， さらに難しい． このふたつについて院生が理解できるか， もしくは理解できたふりができなければ， 「検定」なんぞというしち面倒なしろものの説明は さらに難しくなる．
• ということで， 今年は離散な事例での説明を試みてみよう． 離散確率分布の密度関数は「度数グラフ」として簡単に説明できるんで， 確率なるけったいな概念と観測値たるカウント数の対応が断然つけやすくて， 説明するのもされるのもラクなのである． また R でプログラム作り …… 今ごろやってて間に合うのか?
• 1620． なんか進捗しないと思ったら， かとーオフィス内， けっこう寒くなっていた． ガスストーヴつける．
• すでに夜． そっかー， 離散確率分布の場合， 危険水準◯% の棄却域， というのを設定するときに 5% とかそういうのを 厳密には図示はできんわけだな． おもしろい．

[コイン投げ 20 回 × 一万回]

むろんこれはシミュレイションなんぞやらずとも， 単純な二項分布になるってのはわかるわけだが (R だと dbinom(0:20, 20, prob = 0.5)) …… 試行回数が 1000 ぐらいだと， ぜんぜん二項分布っぽくなくて面白い． で， 10000 とかにするんだったらむしろ barplot() すべきだったかも， と後になってから思いついたり．

• どうも EPS 形式で出力させるのはマズそうなので， jpeg() 出力と jpeg2ps 連携わざでいくことに． jpeg() では日本語が出力されんか． じゃ， 敵国語で， と調べてみる …… コインのオモテは obverse (俗語では heads)， ウラは reverse (tails) …… いや軸ラベル問題は LaTeX で片づけるか．
• ぢりぢりと進捗させるが前途遼遠． はらへったので撤退． 2040 研究室発． 寒い． 2050 帰宅． 晩飯．
• 晩飯後も準備はつづく． 明日はなんだかいろいろありそうなんで， あまり時間がない． ということで， 本日おいつめられてるわけ．
• 関係ないことにふと気づく． 先日， 提出した自由集会の要旨， 自由集会 web page の URL 書くの忘れてた …… ま， いいか．
• 準備が終わらん．
• 時刻 2732 (0332)， けっこうでかい地震． 震度 3 ぐらいかな． ラジオによると釧路あたりは震度 5 強だとか ……
• ときどき余震． お， 震度地図 でたか (tenki.jp で更新されてる)．
• さーて， 授業の準備のほうは …… 7 割ぐらいできたか? もう寝たいなぁ． 時刻すでに 2830 すぎ．
• [今日の運動]
  • 北大構内走 1205-1250． ストレッチング．
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0930): リンゴ．
  • 昼 (1320): 米麦 0.7 合． モヤシ・セロリ・ネギ・シイタケ・卵の炒めもの．
  • 晩 (2150): 米麦 1.0 合． セロリ・ニラ・エノキダケ・豆腐・卵の炒めもの．

2004 年 11 月 29 日 (月)

• 0830 起床． ねむい． 朝飯． コーヒー． 0920 自宅発． 晴． 0930 研究室着．
• また授業の準備のつづき …… しかしアタマが動かん．
• 1020 より平尾君予備審． 1140 ごろ終了．
• また授業準備のつづき．
• 1330 より浦口さん予備審． 1510 ごろ終了 …… と思ったらさらに 20 分ほど拘束された．
• また授業準備のつづき． 1610 ごろどうにかこうにか終了． 配布資料の印刷を依頼 (この目的のために M1 石橋君を雇用していた)． PDF ファイルアップロード．
• 1630 統計学授業， 全 2 回のうち 2 回目． 検定なるものは何なのかを解明する内容． 1750 終了． おわった． ばてた．

[いかにもニセっぽいコイン]

本日の検定説明にもちいた． このニセな雰囲気が説明のうえでたいへん重要なのである． 統計学で計算量や考える量を減少させるのは， こういった主観的な部分だったり， という気もする．

• お茶部屋で予備審終了おいわいケイキのおしょうばんにあずかる …… とつぜん牧さんあらわれる． 今月いっぱいで閉店になる地環研前焼き鳥屋 「おばちゃん」 に顔をだすため． 明朝すぐに東京にもどるとのこと． あいかわらず律義だ．
• 2030 研究室発． 雪． いきなり 積雪 20 cm ぐらいか? 2040 帰宅． 晩飯． ばたばたした一日だった．
• [今日の運動]
  • うんどう休養日．
• [今日の食卓]
  • 朝 (0840): 米麦 0.7 合． セロリ・ニラ・エノキダケ・豆腐・卵の炒めもの．
  • 昼 : 食ってない．
  • 晩 (2100): 米麦 0.7 合． セロリ・ニラ・エノキダケ・豆腐・卵の炒めもの． 昼に食わなかった弁当．

2004 年 11 月 30 日 (火)

• 0620 起床． 朝飯． コーヒー． 0750 自宅発． 晴． 0800 研究室着．

[いきなり積雪]

昨日の朝はまったく雪なんか積もっていなかったんだけど．

• 全国的に修論・D 論しーづんということで， 週に何通か統計学相談メイルいただく． 今朝も一通返信．
• なーんか北大でも学術リポジトリを作ろうとしてるみたいね …… とりあえず， アンケート に答えてみる． 図書館 関係は， 北大の学術的な基盤の中ではかなりマトモに機能してる部門だと思う．
• また北大ネット不調． 学外につながったりつながらなかったり via http． 北大の学術的な基盤の中ではかなりマトモに機能してないへっぽこ部門．
• それだけじゃなく， ちかごろは接続先とはほぼ無関係にこういう ``Access denied'' しやがるし． これは北大の proxy server (強制通過型) の設定・挙動がおかしいのでこういう動作になる．

• うーむ …… さっさと次の仕事にかかるべきなんだが， 「よーやく授業おわった」 脱力からのりはびり状態で， 何もてきぱきとススまん．
• 苫小牧直径成長モデル …… 次の一手をどうするべきか判断が難しいので， お尋ねメイル送ってから， べつに急ぎでもない推定プログラム改善作業なんぞやってみたり． 名目的には， 推定プログラムの点検とでもなるのか …… とサボタージュ労働者は考えるのであった．
• そういや， アカマツ作図宿題もあったなぁ ……
• 苫小牧モデル推定プログラムの手直しだけで昼飯前の仕事終了． これは次の高速化の布石． 雪がふってるんで， しばらくは北大構内走は様子見． 昼飯． 昼飯くってるといきなり煙まみれな大学院生があらわれたり (注: 本人が燃焼していたわけではない)．
• うーむ． 苫小牧モデル作文の私が担当すべき箇所に関しては， 時間的よゆーが生じているような …… DEFCON 一段さげてよいのだろうか．
• とりあえず， 午前の続きで， 推定プログラムのアヤしげ改善のつづき． ともあれ， アヤしげな部分ごとに関数に隔離されてるんで， 困ったことにはならないんだけど．
• しかし， 気分的にはかなりラクなはずの修正作業もススまん． ということで， 先日つくった書籍リストを送って本の発注．
• 苫小牧モデル推定高速化問題， 予想があたった． この複雑なモデルのパラメーター推定には glm() のような単純なやりかたは使えない． しかしながら， 直径成長年変動の説明に使う代謝気象値・光合成気象値を 数千・数万の候補の中からしぼりこむには (特定するわけではない)， glm() みたいな計算で十分， とわかった． これでうまくいけば 7-8 時間かかっていた計算を 10-20 分ぐらいに 短縮できる．
• で， 一回試運転やってみると， ほとんどうまくいった． ただしミズナラでは AIC が 0.2 だけ大きくなってしまった． ということで， 絞りこむ候補数を 10 → 20 にしてみるか．
• 20 にしてもミズナラだけだめだ …… ちくしょう， 食植昆虫に食われまくりの変態樹木め．
• アルゴリズムを変更する． ノミネイト数固定方式はヤメ． 良さそうなやつから順番に試していく． そして， ある気象値セットが N 回にわたって「勝ちぬく」 ことができたら， その気象値を採用するという「勝ちぬき戦」型を採用する． とりあえず N = 10 ぐらいでやってみるか …… ミズナラ以外はこんなことやる必要もないんだけどね．
• ミズナラ， N = 10 でもダメか． ならば N = 20 なら?
• ダメだ． N = 30 もダメ …… そうだな， 粗い選抜であっても glm(..., family = poisson) はまずかったか? glm(..., family = quasipoisson) にしてみる． 平均 vs 分散の関係はこんなふうだし …… あれ? ミズナラでは分散がだいたい平均に比例してるぞ?

• quasipoisson 推定で使うのは擬似尤度だから aic は計算できん． そこで deviance (これも擬似っぽいな) 見て選抜かけてみる …… ダメだ． ミズナラだけはどうしてもうまくいかん． 他の 5 樹種は完ペキなのに．
• ということで， 人間に (少なくともこの私に) 悪意をもっているとしか思えんミズナラ一樹種にわずらわされるのはヤメ． AIC だって 574.5 → 574.7 と 0.2 悪化するだけだし． これ以外は以前とまったく同じ． 計算時間は 7.5 時間 → 5 分， と短縮． 90 倍の高速化 (60-70 時間を必要としてたころと比べると 800 倍以上の高速化)． すばらしいことですね．
• といった結果どもを wiki にアップロード．
• 1830 研究室発． 2000 帰宅． 体重 73.4kg． 晩飯．
• 苫小牧直径成長モデル， PDF 出力式の作図プログラム作りなおし． うう． われながら健忘症とゆーか， ちょうど 一ヶ月前 にやらかしたのと同じ間違いを連発． いやはやと思いつつバグとり続く．
• [今日の運動]
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0720): 米麦 0.7 合． ネギ卵炒飯． キャベツ・ダイコン・サバぶしの味噌汁．
  • 昼 (1240): 研究室お茶部屋． 米麦 0.7 合． キャベツ・ダイコン・サバぶしの味噌汁．
  • 晩 (2050): 米麦 1.0 合． チンゲンサイ・タマネギ・豆腐の炒めもの． キャベツ・ダイコン・サバぶしの味噌汁．