ぎょーむ日誌 2006-04-(11-20)

2006 年 04 月 11 日 (火)

• 0730 起床． 朝飯． コーヒー． 0850 自宅発． 晴． 本日は今年はじめて札幌市の気温が 10°C こえるそうで． 0905 研究室着．
• 昨日の新歓コンパ会場 A804， にぎやかな宴のあとというか 「あけがたまで飲んで騒いでたんだろうな」 と思わせるなかなかの惨状 …… 小山さん・工藤さん・亀山さん・井田君とかたづけ & 皿洗い． 一時間ほどかかる．
• ヤチダモ遺伝子計算， 昨日のつづき． これは明日こちらにみえる後藤さん (東大富良野演習林) データ解析のてつだい． 例によって例のごとく 「そんなの簡単にできますよ」 なる自分のセイルストークの責任をとっているのである
• めんでる計算まわりは昨日で片づけたつもりだったが …… 200 種子の「父親さがし」の関数つくりに苦闘する． しかも計算に時間かかるし． 私の ThinkPad X31 だと 4 分間ぐらい．
• しかも 200 種子のうち 99 で 「探索範囲内では父親みつからず」 だと?
• 昼飯． 活力ありあまりの _女性大学院生たちもさすがに 昨日の盛大きわまりない新歓コンパで燃えつきてしまったらしく， 本日は影もカタチもみえない …… M1 のヒトたち (の一部?) は地環研・最初の三日間せんのうセミナーとかにでてるのかな? A 棟 8F は平和で静かだ．
• ヤチダモ遺伝子計算つづき． R においてたくさんの処理を行列計算に帰着できると， すごく速く処理される (上記めんでりあんはどうやっても行列計算にならないので 時間がかかる)． しかしこれはプログラムつくってる人間のアタマがけっこうくたびれるんだよね ……
• 父親樹-母親樹の距離行列計算に疲れたので， お茶部屋でなべ洗い． 昨日の新歓コンパでつかわれたもので， 直径 80cm ぐらい．
• 計算プログラムかきのつづき …… うーむ， このヤチダモ花粉モデルの論文， すでに Molecular Ecology に accept ずみなんだが …… なんか確率計算がヘンだという気がする． まあ， このあたり気にせず書いてあるとーりにプログラムを書けばよい， のかしらん?
• しかし気になる． 母親樹があたかも花粉親であるかのように考え， このニセ花粉親がまきちらした花粉が「系外」 (plot の外) に飛びだす確率を p とする． この確率 p こそがその母親樹が系外から花粉うけとってる確率? どうしてそんな計算になるわけ?
• この件うだうだ考えてもよくわからん …… 亀山さんと相談してるうちに少しわかってきた． けっこうネが深いかも． というのも， そもそも Burczyk ふうの

  p ×(めんでる確率 | 系外) + (1 - p) ×(めんでる確率 | 系内)

  という確率 p の「いれかた」からしてぢつはスゴく ……
• と考えてるうちに M1 なヒトたちが授業からもどってきたので， あれこれとまきこまれてしまう． ThunderBird のインストール & 設定， とか．
• 違法こぴーやめましょう， ということで OpenOffice.org とか．
• えーと， MacOS X つかうかもしれない M1 に対しては …… かとーさんに教えていただいたところによると， NeoOffice のほうがいいかもよ， と ……
• 時刻はすでに 2100． いかん． しかしばててきたので， 撤退． 2210 研究室発． 2230 帰宅． 晩飯．
• ふにおちない確率計算にツマってしまってプログラミングは進捗しなかったわけで …… ここを回避して， 例題的な推定問題を構成してその計算プログラムをつくるか?
• [今日の運動]
  • 今日もうんどう休養日 ……
• [今日の食卓]
  • 朝 (0820): 米麦 0.6 合． ダイコン・ブナシメジ・豆腐・油揚・海藻・サバぶしの味噌汁．
  • 昼 (1250): 研究室お茶部屋． 新歓コンパの残りもの． 鶏なべ．
  • 晩 (2340): 米麦 0.7 合． 海藻・豆腐のスープ． コマツナあえもの．

2006 年 04 月 12 日 (水)

• 0500 起床． さくばんは， いきなり寝てしまった． 朝飯． コーヒー． やちだも Burczyk 的モデルぷろぐらみんぐのつづき． 0700 自宅発． 小雨． 0715 研究室着．
• やちだも計算プログラミング， (花粉親成功確率をあらわす) 対数尤度関数の定義は行列の転置とかにひたすら注意してれば， それほど難しくなく書ける …… しかしながら， あとから見ると

  v.likelihood <- apply(
  	t(t(phi) / phi.denominator) * matrix.mendelian,
  	2,
  	sum
  )
  

  これってたしかに意味はわかるけど， どうにかならんのかねと思わずにはいられない (しかしこういう行列計算にしてやると速い)．
• 1010 (時間ぎりぎりではあったが) ひとまず R の最適化関数 optim() によって Burczyk タイプ花粉親競争モデルの パラメーター最尤推定できるようになった． 昨日の「系外問題」のめんどうを (とりあえず) 回避するために， 「花粉親が特定されたもの」 のみを対象とする 「R 説明のための簡略化した推定問題」 になおしたためか， 得られた推定値がもと論文とはずいぶんと異なるものになってしまったけれど ……
• 1100-1630 東大富良野演習林から来てくださった後藤さんと 風媒花粉親成功モデルの検討と R での推定計算実装例に関する議論． たいへん勉強になった．
• 後述するように， 本日はばてきったので詳しくは明日にでも …… と書いて「後頭部に狙撃者の視線」的イヤな予感を感じたので 思わずうしろをふりかえってみる …… やはり 富田君 の監視下にあったか． ということで， 最小限のことだけ書いとくと
  • 昨日ひっかかっていた 「系外からの流入確率計算に関する疑問」 は必ずしもマトはずれではなく， この確率と花粉 dispersal kernel のカタチを どう関連づけるかはぜんぜん自明ではない
  • Burczyk 的な系外・系内分割方式もまた 「ホントにこれでいいのか?」 をしつこく検討してみる必要あり
  えー， つまり Burczyk 的 な多項 logistic モデルによる花粉親成功の表現， というモデルは興味ぶかくかつ有用なんでしょうが， 「もうこれで疑問の余地ナシ完ペキなる定式化」 などではなかろう， と．
• そしてど忘れせぬうちに書き記しておかねばならない …… 本日の最重要検討事項のひとつ 「このポーランド人らしき Burczyk はなんと読むのか?」
  • 後藤説: ばるざっく
  • 久保説: ぶるくつぃく
  ポーランドではよくある名前， みたいなんだけど． 折り紙研究家 として著明らしい Burczyk さんもいる．
• 後藤さん富良野に帰られると同時に M1 なヒトたちも ちかん研「これがちきゅうかんきょう科学の基本だ最初の三日間」講義 (きわめて猛烈に不評) からもどってきて， またいろいろと …… 早朝プログラミング → モデルの対話的検討 & R 議論 → 新年度年頭 M1 さぽーと (Outlook Express の設定までやってしまった …… 末世だ)， という連鎖でへろへろになる．
• そしてこういう状況で， データ解析こんさるメイル (破壊力のつよいやつ) 二通ばかりいただくと， アタマの中がホントに真っ白になった …… はい， 復旧してからお返事さしあげます．
• ふらふらと撤退． 1945 研究室発． 2000 帰宅． 晩飯． いきなり寝てしまい， そして真夜中におきてしまって， ぼーっと ……
• [今日の運動]
  • またうんどう不足すぱいらるにおちこみつつある
• [今日の食卓]
  • 朝 (0520): 米麦 0.6 合． 海藻・豆腐のスープ． コマツナあえもの．
  • 昼 (1210): 百周年記念会館のきゃら亭昼定食 680 円． 後藤さん・亀山さん・平尾君と． なんと， 亀山さん・平尾君は生態遺伝学方面における 「大御所」 だったのか!
  • 晩 (2030): 米麦 0.7 合． ネギ卵炒飯． 海藻スープ． サバ電磁波酒蒸し．

2006 年 04 月 13 日 (木)

• 0830 起床． 朝飯． コーヒー． 洗濯． 1040 自宅発． 曇． 1055 研究室着．
• 本日からかとー室長が筑波大がよいなので， A803B はさしずめ「かとー留守オフィス」とでも呼ばれるべき状態になった． かとーさん荷物はまだいろいろとあるわけで …… しかしすばやい大学院生がかとー高級椅子を「借りて」いってしまった．
• 苫小牧樹木直径モデル論文原稿にふたたび取り組んでみる． 起動状態よろしくない．
• 昼飯．
• 1328 地震． また 震度 1 ぐらいなんだけど， A 棟 8F はぐらぐらする．
• 論文原稿よみ． ぢりぢりと鍋嶋わーるどにはいりこむ． 自分で書いた部分にさしかかると， いつものことながらアタマいたくなる．
• この推定計算のわくぐみでは「樹種間比較」はかなり難しいわけだが …… ふと「全樹種まったく同じパラメーター」 という再単純モデルなら簡単に計算できることに気づく． たぶんすごくあてはまりの悪いんだろうけど …… というのも， かなり樹種差があるのに「差がない」と仮定したモデルなんで．
• A801 室「苫小牧ボス」こと鍋嶋さんに拝謁し， かかる計算の勅許をいただく． Dell desktop 機 (P4 2.4GHz) で 25 分弱で計算終了． やはり「樹種べつ」モデルの合計 AIC = 5470 にくらべて かなりよろしくないモデルだ (AIC = 5535)．

• しかしまあ， いちおうこれで「樹種間比較」最初の一歩にはふみこんでいることになる． これ以上のまぢめな計算は難しい． それでも論文原稿修正の， 多少の役にはたつかもしれぬ．
• M1 大学院生が突然かとー留守オフィスにやってくる．
  「昨日のぎょーむ日誌に関する質問あるのですが」
  「はい ……」 (いきなりなにか花粉流動まにあっくな質問だろうか……)
  「以前から疑問に思ってたのですが」
  「はい ……」 (どきどき……)
  「『今日の食卓』の『電磁波酒蒸し』って何ですか，あれは?」
  かなりずっこけつつ電子れんぢすなわち microwave oven でつくる酒蒸しであります， と回答さしあげる．
• しかしその後はあまり進捗もなく撤退． 1910 研究室発． 1930 帰宅． 晩飯．
• 昨日は敵前逃亡したこんさるメイル書き． 二通だけなんだが …… やはり時間かかる．
• さて後藤さんと昨日検討していたハナシの一部をこのへんに書いてみる． Burczyk 的モデルとはこういったかんぢで (自殖ナシの場合):

が遺伝子型 h の母親から g という種子ができる条件つき確率， それが系外花粉による場合 と父親 j の花粉による場合 ; が (あれこれとアヤしまれてる) 「系外からの花粉である確率」， そして (系内限定の) 父親間競争は

と多項 logistic 的になっている． は母親 i における父親 j の「競争能力」で樹間距離依存性 (これが dispersal kernel できまる) とか父親のサイズとか開花状況の関数である．
• で， 後藤さんたちの新モデルは が 花粉 dispersal kernel のカタチと母親樹の plot 内位置で決まるんぢゃないの， という変更をいれたもの …… 私はその確率計算についてはあまり正しいと思っていない．
• で， こういう系外・系内分割のやりかたがそもそもおかしいのでは， という意見もありえるわけで． どこからきた花粉であろうと同じように柱頭の上で競争させればいいはず， ということで …… たとえば代替案としては，

としてしまって， 競争の部分を

とでもしてみるか? は系外花粉の競争能力で， ここをどう定式化するかがとてもむずかしいんだけど (つまり自分でもどうしたらいいかわかってないんだけど)，
• 花粉 dispersal kernel のカタチ
• 母親樹の plot 内位置
• 系外樹木密度とサイズ分布
みたいなものどもに依存している， と．
• さらに花粉親樹木あるいは相性問題を random effects でくみこむ， とか， は系外遺伝子頻度に依存しているんだけどこれを EM アルゴリズムとか Bayes 推定的に構成しながら計算する， とか， このあたりのモデリングはまだまだ考えるべきことがあると思うんだよね． 嗚呼， 難問山積 ……
• ありがたいことに 富田君 が重要な難問に関する研究を一歩ススめてくれた …… Burczyk は 「ぶぁーちっく」と読むのが適当 ではなかろーか， という考究である． なるほど， 「いんちきドイツ語ふう読み」 とでもいうべき久保説「ぶるくつぃく」 よりはもっともらしいような気がする …… たしかに Wikipedia でも ``cz''は「ちゃ」音， とあるよな．
• [今日の運動]
  • またしてもうんどう休養日か ……
• [今日の食卓]
  • 朝 (0920): 米麦 0.6 合． 海藻スープ． サバ電磁波酒蒸し．
  • 昼 (1310): 研究室お茶部屋． 米麦 0.6 合． レトルトパウチドのトマトソース．
  • 晩 (2030): 米麦 0.7 合． タマネギ・豆腐・タラすりみカレー．

2006 年 04 月 14 日 (金)

• 0700 起床． ねむい． 朝飯． コーヒー． 0855 自宅発． 晴． 0910 研究室着．
• 朝からうだうだと ……
• A801 の壊れプリンター EPSON LP-2000C 修理のため担当のヒトがくる …… そしてややめんどうな経緯あり， 保守契約 3 年間というのをすることになった． 本日は「修理されたこと」になってるんだけど， たぶん来週にでもまた修理されるのであろう． ふう．
• あ， それからエラー番号をまちがって伝達していたので， 面倒が拡大したような気がする． 正しくは 6031 だった． 反省． いまさらながら．
• 新年度といえば生命保険攻撃． 毎年新手のてごわいおねーさんたちがやってくるのだけど， 手口がわんぱたーんなのでさすがに慣れてしまった．
  「何か生命保険に入っておられますか」
  「入ってます」 (うそ)
  「それはもう加入されて長いんですか」
  「はい」 (うそ)
  「えーと，それはどちらの会社に……」
  「教えません」
  「…… ここにお名前・生年月日を書いてもらうと○○のチケットを手配 ……」
  「書きません，それって『個人情報』じゃないですか」
  で終了． ぎょーむ日誌とか書いてて「個人情報」もへちまもあったもんぢゃないけど， 個人情報保護法なんかの御威光 こういうふうに私的利用してみるのは便利ですなぁ．
• も一件の雑用あるんだけど着手する気になれず， また雑用があると論文原稿なおしもなかなか …… と現実逃避的にまた別のデータ解析こんさるメイルへの返信， ついでに GLM 参照 なるペイジを作る． 前に作った GLMM 参照 と半分ぐらい重複した内容．
• 昼飯． そして雑用の午後になってしまった ……
• まずはくだんの (まず無意味と断定できる) ソフトウェアのライセンス番号管理とかいうやつ． こういう阿呆らしい作表ぎょーむは Perl プログラミングによる問題解決が有用である． 階層構造データを CSV 変換する Perl スクリプトを作り， それを Gnumeric でゑくせる表に変換する． 会計係に送信して一件落着．
• 雑用モードになったんで， 研究室雑用 マニュアルを整備する …… いろいろ追加・改訂したんだけど， まだまだってかんぢだね ……
• 本業はあまりススまぬまま 1905 研究室発． 1920 帰宅． 晩飯．
• なンか毎晩毎晩 2-3 通の統計こんさるメイルを書いてるような気がする． 本日いただいた質問は明快なものばかりなんだけど …… けっこうたいへん．
• [今日の運動]
  • しかし気づくと真夜中になっていた ……
• [今日の食卓]
  • 朝 (0820): 米麦 0.7 合． タマネギ・豆腐・タラすりみカレー．
  • 昼 (1350): 研究室お茶部屋． 米麦 0.5 合． タマネギ・豆腐・タラすりみカレー．
  • 晩 (2020): 米麦 0.7 合． ネギ・キムチ炒飯． 海藻スープ．

2006 年 04 月 15 日 (土)

• 0530 起床． 本日は大阪の祖母宅をちょっとのぞいてみる予定． 朝飯． 紅茶 …… 寝不足ぎみのときにコーヒーのむと気分わるくなるんで．
• 0630 自宅発． 0638 JR 札幌駅着． 大阪というか神戸空港往復航空券を自動チェックイン機からとりだす． 0648 同発． 0736 新千歳空港着． 1040 円． 0805 ANA の B767-300 で離陸． 0955 波間に沈みそうな神戸空港着． 1019 ポートアイランド線に乗って同発． 1037 三宮着． 320 円 阪急の駅に移動して 1050 同発． 十三で阪急宝塚線にのりかえて 1129 豊中着．
• さーてここからが難所だ． 駅の地図をぢーっとながめて経路の把握につとめる． 雨の中を 20 分弱あるいて 豊中高校 着． そこからちょっと北側の住宅地にはいりこんで 1150 祖母宅着． すでに親戚一同そろっている． こちらのおばあちゃんも 90 才こえてるのに すごく元気で， アタマもたいへんしっかりしているというか， 私などよりよほど短期記憶・長期記憶が強められている， とゆーか ……
• 昼飯とか晩飯とかで今晩は大阪泊．
• [今日の運動]
  • 明日は札幌にもどるんで ……
• [今日の食卓]
  • 朝 (0550): 米麦 0.5 合． 海藻スープ．
  • 昼 : 大阪昼飯．
  • 晩 : 大阪晩飯．

2006 年 04 月 16 日 (日)

• 0700 起床． 朝飯． 0810 中之島の宿発． 0830 阪急梅田発． 0858 三宮着． 310 円． 0910 ポートアイランド線三宮発． 0928 神戸空港着． 320 円．
• 時間あるんで空港内すこしうろうろしてみる． コンパクトな空港． 神戸市民がたくさん見物に来ている …… 昨日も展望ラウンジにヒトがぎっしりいたので， 誰か有名人でもやってくるのかと思ってたんだけど …… そうではなく「飛行場」というモノじたいが珍しくて見物にきたらしい．
• 1050 離陸． 晴． 1220 新千歳空港着． 1234 同発． 1310 札幌着． 小雨． 昼飯など調達して 1330 研究室着． 昼飯．
• うーむ， またデータ解析こんさるメイルがいろいろと ……
• あまり進捗せず 1740 研究室発． あいかわらず小雨． 1800 帰宅． いきなり寝てしまった． 夜中ちかくにおきて晩飯． また生活周期がずれそう．
• [今日の運動]
  • うーむ ……
• [今日の食卓]
  • 朝 (0750): 宿の朝飯．
  • 昼 (1350): 研究室お茶部屋． 幕の内弁当 470 円．
  • 晩 (2330): マカロニ． レトルトパウチドカレー．

2006 年 04 月 17 日 (月)

• 0750 起床． 朝飯． コーヒー． データ解析こんさるメイル書き． 1010 自宅発． 小雨． 1025 研究室着．
• 起動状態よくない． データ解析こんさるメイル書きのつづき． ネット雑用少々．
• 苫小牧樹木直径成長モデルの論文原稿修正にとりくむ．
• おっと， 今日は「何時かわからんけどプリンターとどく」日だったな (こういうオプションにしとくとよけーな金をとられない) …… いつ来るのかわからんので， 1300 前に昼飯すませておくか． 昼飯．
• R-help メイリングリスト， help(formula) よく読め， か …… たしかに I() とかいろいろなワザがあるな． そして私はじつは formula オブジェクトがよくわかっていなくて， ぜんぜん使いこなせていないんだよね． そのうち勉強しなくては．
• A704 室に新プリンター LP-S6500 が運びこまれる． 私が行ったときにはすでに院生たちが設置を終えていた．

[EPSON LP-S6500]

カラーレイザープリンターの低下価格は現在も進行中． これは 18.8 万円だった． まだ安くなるだろう． ネットワーク設定で少してこずる． IP アドレスは簡単に設定できた． しかしそれが反映されない …… 正解は「電源をいったん切る」という野蛮なものだった．

• そしてこのプリンターの保証書なんだが …… なにやらシリアルナンバーを myepson.jp というよくわからぬところに登録すると， 一週間後ぐらいに郵送されてくるらしい．
• 甲山弟子の新 M1 の宮田さん・田辺さんが (当然のことながら，と言いますか) 「(研究課題の背景に関する) 甲山さんの説明がよくわからない」 とのことなので， 三人して「甲山わーるど」に関して知っている断片をつなぎあわせて あーだこーだと議論してみる．

[甲山わーるど]

…… のごく一部． これはどちらかといえば宮田さんがらみの， 「葉量一定」とかいう法則だか信仰だかよくわからぬアレが， 成立しているとしたらどのような条件を満たさねばならぬか， といった林分の構造と「枝さき」をつなぐハナシ． ぱいぷ樹木 なんかも無関係ではない．

• ネット雑用あれこれ． 研究室雑用 Wiki． 新 M1 の hokudai メイルアカウントがでそろってきたので， メイリングリスト をぼちぼちと運用開始． いろいろと連絡など．
• 粕谷さんから連絡いただき， 粕谷さん・私でだしていた 科研費 (萌芽研究 -- といってもシウリザクラとかの sprouting とかを研究するわけではない) がどうやら通ったとのこと． 私は寄生虫のごとく粕谷さんにタカってるだけの構造なんだが …… まあともかくひと安心． 今年度の研究費はこれとあと東「アリ」科研だけなんで．
• また苫小牧樹木の世界に入りこもうとする努力 ……
• しかし稲葉さん・川合さんと繁殖生態学の進化的側面における coarse-grained と fine-grained のちがい， といった問題をいっしょに検討することになる …… えーと， いまそういうコトを書いた本の輪読会をやっておられるとのことなので． 概念は別に難しくなのだが， その本 のその章の説明がへっぽこなので苦闘させられる．
• すでに夜中ちかい． 2300 研究室発． 2315 帰宅． 晩飯．
• [今日の運動]
  • まあ調査しーづんになれば A 棟 7-8F の院生密度も下がるはずで ……
• [今日の食卓]
  • 朝 (0840): 米麦 0.6 合． ホウレンソウ・タマネギ・ショウガ・豚タンの炒めもの．
  • 昼 (1250): 研究室お茶部屋． 米麦 0.7 合． ホウレンソウ・タマネギ・ショウガ・豚タンの炒めもの．
  • 晩 (2350): 米麦 0.7 合． ネギ卵キムチ炒飯． 海藻スープ．

2006 年 04 月 18 日 (火)

• 0810 起床． 朝飯． コーヒー． 0955 自宅発． 曇． 1010 研究室着．
• 朝からプリンターどたばたとか． Canon UM-86E (LBP-2810用ETBユニット) …… なるナゾの部品だったもの (実体は紙搬送系) の交換． これで紙づまりが少しは解消されるか?
• 仕事ススまぬ …… 小雨ふるなか北大生協まで昼飯調達の旅にでる． 研究室にもどって昼飯．

[ファインマンの手紙]

北大生協書籍部によったついでに買ってしまった． なんとなくファインマン本は手をだしてしまう …… そしてこれはハズれではない． 値段は高いけど． トンデモ系・デンパ系手紙への対応がなかなか楽しい． 原著廉価版 がもうすぐでるらしいから， そちらも入手してみたくなったり．

• 1440-2000 農学部の高橋さんとお茶部屋で R データ解析相談． ここでひとつ判明したことは …… Crawley 先生の R 本， 基本的に良い本だと思うんだがときどき効率よろしくない方向に流れるね． 「効率悪い」というのはある意味「教育的」だとは思うんだけど …… p.242 あたりからはじまる負の二項分布に関する説明， 説明にまちがいあるわけでもなく丁寧なものだとは思うんだけど， 必要もない関数を定義したり glm.nb() を紹介していなかったり ……
• 一般化線形モデルなどはいわゆる多変量解析の代替になるか， というハナシ． 現在つかわれてる 「加重値きめて足し算して固有値最大化」 みたいなたぐいの (つまり線形) 多変量解析なるものどもが， 生態学方面データ解析の (まあひかえめにいって) へっぽこツールとして濫用されているわけだが …… このあたりの大半は将来的には滅亡するのはまちがいあるまい． 理由は単純で「固有値最大化」するような加重値など 生態学にとって何か興味ある数量にはなっていないんで．
• 一般化線形モデルは要因に加重値 (係数のこと) かけて (線形モデルなんで) 足し算して …… といったところが多変量の操作としては似ていて， さらにリンク関数かませて確率分布をつくって説明したいパターンに フィットするよう加重値を調節していく， というところはやみくも無意味な固有値最大化よりよほどマシ． で， モデル選択と組み合わせれば「いらない」要因は削れるし．
• しかしいまその 「加重値→……→固有値最大化」 的な (線形) 多変量解析とやらでやってることを一般化線形モデル その他あれこれといったまっとうな統計モデルで置き換えようとすると (線形多変量解析はかなり統計っぽくない)， いろいろとモメそうだな．
• たとえば ordination とやらは GLM のたぐいでは単純には置き換えられないだろうし．
• さらに， だ …… 「ひたすらむやみにあれこれ測定したあげくにそれを全部説明につかう」 みたいなイヤなお作法・根性論が美徳とされてる世界があって， ですね．
• 一般化線形モデルなんかを モデル選択でモデルからよけーな要因を削除して単純化するってのは 「予測」の役にたたんから (むしろジャマだから) いらん， みたいなハナシで …… まあ切り捨てられた要因については， それぞれ「捨てられた要因 X1 は残された要因 X2 と相関が高いんだけど， X2 よりぐしゃぐしゃしてるんで捨てられた」 みたいな説明はいちいち付けることができるわけだが．
• さらにこの生態学データ解析の多変量わーるどに Baysian なモデルが参入してくるとなると (これは「確定している未来」とでもいうべきもの) …… いったいどうなることやら．
• 2045 研究室発． 2100 帰宅． またいきなり寝てしまい， 真夜中に目がさめた．
• [今日の運動]
  • なんかヘンな生活周期にはまりこんでるな ……
• [今日の食卓]
  • 朝 (0830): ヨーグルト．
  • 昼 (1310): 研究室お茶部屋． 北大生協いなりずし 298 円．
  • 晩 (2530): クラッカー．

2006 年 04 月 19 日 (水)

• 0915 起床． うーむ， いかん …… 朝飯． コーヒー． 1030 自宅発． 晴． 1045 研究室着．
• メイル書きとか．
• XOrg-6.8.2 上の galeon-1.3.17 で slashdot.jp のペイジをみるとときどき (というか fold された投稿を展開しようとすると) X ごとオチるという豪快なとらぶるに巻きこまれる確率が高い (他のサイトではこういうことにはならない)． まあ， そんなところ見てるなよ， ということなのかもしれないが． 原因不明で興味ぶかい．
• A802 の CANON LBP-2810 カラーレイザープリンター， ETB ユニットをかえても紙づまりぎみのよーで …… 最悪の場合は定着器ユニット (UM-86F?) も交換せんといかんか?
• 苫小牧樹木直径成長モデル論文にとりくもうとする． 本日は他の用事もなくぢりぢりと入り込める．
• 昼飯．
• いまさらながら年輪解析まわりの論文あつめ． やっぱヘンだよなあ ……
• まーたまた難読人名問題． Lebourgeois， ふらんすぢん …… いろいろ調べてみた結論として， おそらくこれはほぼ確実に「ルブルジョワ」だろう． bourgeious は「ぶるぢょわ」なので．
• 1700 から一時間ほど 1F 研究「院長」室よこの会議室ですごす． 院長が院内の「助手」を全員招集したので …… その人数， 十数名． まとめてクビにするのかと思いきやそういうわけでもないらしく， 来年度から「助手」がなくなり「助教」になるんだ， という説明に 60 分も費してくれた．
• で， その不明瞭なる説明の要約を試みてみると ……
  • 名前: 助手が助教になる
  • 任期がつく
  • 新任助手は募集されず， 今後は助教が公募される
  • これまで (公式には) 「教える」仕事はなかったけれど， 助教は公式に講義・実習などうけもつ …… 単位認定権? それ持ってたら居眠りの人数が減るのかしらん?
  • 他はよくもわるくも何も変わらない …… 給料かわらず， 会議でなくてよい， 修士号・博士号まわりとは無関係， などなど
  ねんのために書いておくと， 「助教」なるみょーな名称変更は 「劇ヘン」だの 「起学」だの といったここ地環研の いつものウケねらいの奇矯なふるまいではなく， 第 162 回国会における文部科学省提出法律案 が可決されたために日本全国いっせいに変わるのである．
• 1940 研究室発． 1955 帰宅． 運動． 体重 73.2kg． うんどう不足ふとりだな． 晩飯．
• [今日の運動]
  • エアロバイク 30 分間．
• [今日の食卓]
  • 朝 (1000): 米麦 0.7 合． 塩ザケ． コマツナあえもの．
  • 昼 (1440): 研究室お茶部屋． 米麦 0.6 合． 塩ザケ． コマツナあえもの．
  • 晩 (2220): 米麦 0.7 合． ネギ卵キムチ炒飯． 海藻スープ．

2006 年 04 月 20 日 (木)

• 0845 起床． うーむ， 生活周期が補正されない． 朝飯． コーヒー． 0955 自宅発． 雨． 1010 研究室着．
• ちょっとだけネット雑用． それとは無関係なんだけど， 生物圏専攻長からこういうメイルが皆さんにとどく．

  環境生物学部門で、生物圏科学専攻を担当している
  教員の皆様：
  教育経費ですが、多分去年と同じ程度（少し減額に
  なる可能性がありますが）配分されると思います
  （去年の配分額：修士：740000円、博士：112000円）。
  今年度、すでに、教育経費で支払いが可能ですので、
  もし、教育経費で支払いをしたいとお考えの場合は、
  昨年同様、納品書等、中野さんまでお持ちください。
  また、教育経費で書籍を購入した場合は、こちらに
  納品書が来ませんので、金額を中野さんまでお知ら
  せください。
  なお、以前のもアナウンスいたしましたが、経理の
  単純化のため、学生が文献複写をする場合は、指導
  教員の研究費でお支払いください（指導教員の名前
  で依頼してください）。
  

• ふーん， これぐらい研究室にお金が来るんだと思ってたら， 翌日になって

  昨日お知らせいたしました修士学生の経費は間違いです。
  宝くじにあたった気分で買い物などしないよう、お願い
  いたします。
  

  なる連絡が． なんなんだ．
• 樹木直径成長モデル論文よみ …… 旧来的な年輪気候学的なモデリングのマズさは 「標準曲線」 なるものを
  • 根拠なく …… というか「見た目」で
  • 年変動する気象要因と無関係に
  でっちあげてしまうところだ． いつまでもこういう「ばらばら推定」 とかに無意味にコダわるのはやめて， ですね …… 年変動する気象要因が重要だと考えているなら， その標準曲線なるものを計算するときに同時に考慮してやればよいわけですよ． はい．
• このあたり私などが指摘するまでもなく， このあたりの問題点をよく考えた研究は 年輪解析方面でもぢりぢりと増えつつある． ただしたいていの連中は旧来の方法を惰性で使いつづけているように見える．
• 雨の中昼飯調達にでる． 生活周期が不規則だと炊事ローテイションに悪影響がでて， ですね …… 研究室にもどって昼飯．
• 指摘された文献まわりにひととーり目をとおしたので， それを原稿にどう反映させる・させないを考える …… われながらのろのろとした作業だな．
• 1500 すぎから農学部の 城田さん としべりあカラマツ個体ベイス林分モデルの相談． 「燃え倒木更新」 などなど， みょーな空間相関が必要とされているのかも ……
• ハナシはなかなか尽きないのだけど， 城田さんは本日は 北大構内 (!) パイプオルガン 演奏会をのぞいてみる， ということで 1800 前に終了． パイプモデラーはパイプオルガンにも興味をもつのである．
• M1 プリンタードライヴァーインストール作業 …… これでおそらくほぼ終了のはず． さて， Epson LP-9200ps2 みたいな postscript プリンターのドライヴァーはどこで入手するか? という問題． 正解は「Adobe の FTP サイト」 ではないかと思う． ゐんどーづなヒトたちは ftp://ftp.adobe.com/pub/adobe/printerdrivers/win/1.x/ winstjpn.exe が使える (2006 年 4 月現在における最新版なので)． Mac なヒトたちは …… たぶん ftp://ftp.adobe.com/pub/adobe/printerdrivers/mac/8.x/ の JPN.sea.hqx? はい， Unix なヒトたちはこんなアホな苦労は何も必要ありませんね．
• しかしあいかわらずげいつ OS は馬鹿というか …… 二回も再起動させられてばてきったので撤退． 1930 研究室発． 雨． 1945 帰宅． 運動． 体重 73.8kg． 晩飯．
• [今日の運動]
  • エアロバイク 50 分間．
• [今日の食卓]
  • 朝 (0840): ヨーグルト．
  • 昼 (1350): 研究室お茶部屋． 北大生協の韓国ふうのりまき． 338 円．
  • 晩 (2320): マカロニ． ニンジン・タマネギ・ブナシメジ・豚レバのシチュー．