ぎょーむ日誌 2005-12-21

2005 年 12 月 21 日 (水)

• 0840 起床． 朝飯． コーヒー． 0945 自宅発． 晴． 1000 研究室着．
• 昨日の北大ネット途絶に関して …… 管理者連中はまた「なかったこと」にしてるよーで． あるいはホントに途絶したことに気づいてなかったとか?
• 苫小牧直径成長問題， 樹木個体 random effects の積分のやりかたを変える， とゆー昨日の改造プログラムの試運転． 最尤推定はうまくできるけど， 事後分布の作図がうまくいかない． 試行錯誤しつつばぐとり． 正午前にまっとうそうな結果が得られるようになった．
• こんなかんぢで． 二ヶ月前の 10/27 ごろに計算してたのとほぼ同じなんだけど， 事前分布の「下の端」まで計算するようになった， というのが相違点である．

• 1130 北大ネットまたしても原因不明の通信途絶． 少なくともこの建物 (地環研 A 棟) は完全に孤立している． A 棟 hub が「ツマってる」のか， ということで 3F まで行ってまた無断強制再起動を試みたんだが …… ぜんぜん効果ナシ． ということで北大情報基盤センターとやらが管理している 地環研ルーターかそれより上流 (ただし北大内) のどこかが壊れているもよう． 1157 ごろ一瞬だけ外界につながったけどすぐに沈黙．
• 昼飯． 1300 すぎてもネットはネット復旧せず．
• かとーさんが i-mode などで調査し考察したところ， どうも A 棟だけがネットから分断されていて， その原因は (以前にも同じ事例があったんだが) A 棟内の virus 呪われゐんどーづ機暴走なんかにあるんでは， とのこと． そこでまた棟 hub を無断再起動してみたんだけど， 状況は変わらず．
• しょうがないんで内線 4494 でネットワーク掛に電話． 状況説明． やはり先方は管理者のくせにこの通信途絶に気づいてなかった． えーと， ここ 8F の hub 番号は T2208， 使用ポート番号の例としては 02 など． 調査するとのこと．
• しばらくしてから脱力返答． 通信不能なので状況わからず． ぎょーしゃさんとやらを現場に派遣する． 時間かかるけどそれでいいか? …… これって「だめだ」といったらもっと迅速に復旧してくれるわけ?
• ぼーっとしてるのも何なのでかとーさんと調査つづける． どうも A 棟 hub と 8F hub のどちらかに問題あるらしい． 8F 闇ネット内には問題なし， と実験によって確認． となると， 当方には打つテなし， か．
• 1440 きゃたつかかえた「ぎょーしゃさん」とやらが登場． 北海道ゼロックスから二名きた． ThinkPad X40 の USB から 8F hub の RS232C シリアルポートに直結 …… するまでもなく uplink ぎれ， とわかる (私はそのあたり注意してなかった)． で， 一人が 3F hub まで行ったところ， 8F につながる (12 番ポートの) ケイブルがぬけかけていた， とのこと． さしなおしたら復旧． 誰がやったんだ?
• ともあれ北大ネット途絶時の点検手順を更新．
  1. かとーオフィス内の断線を調べる
  2. A 棟 8F 内の通信状況を調べる
  3. A 棟 8F の uplink を調べる
  4. 3F hub の 12 番ポート のケイブルを点検
  5. 3F hub の電源再投入
  6. それでも動かなければ 情報基盤センター 4494 に電話する
  というところか．
• なぜかしら間瀬さんが RjpWiki 内に Gmail 宣伝 ペイジを作っておられる． 先日 このへんで書いた ``View as HTML'' 機能紹介しました， とわざわざ連絡をくださったので ……
• 苫小牧樹木直径成長問題， 本日えられた結果などアップロードしてみる． 「個体差」 (random effects) をあらわす事前分布のいれかた (積分のやりかた) を変えると推定結果がびみょーに変わった． 全体の最大化対数尤度も改善しているので， 喜ぶべきことなのかもしれないが …… まあ， ひじょーにびみょーな差位であるし ……

• イタヤカエデ (AM) なんかにおいて， 大サイズで成長が「へたれる」ようになった， といったところが見てくれに関して変わったところ． これは「私は集団内のすごくへっぽこ樹木個体です」と自称する連中の 「個性」を (雑ながら) 端まで計算できるようになったため …… と考えたいのだが， あまりにも苫小牧ボスの意向にそった推定結果なので， 改造点を中心に計算プログラムの再点検にとりくまねばなるまい． というのも， 私が「無意識」のうちにボスの意向に迎合する計算結果を出すような コード書きをやってた可能性が ……
• 変数名改善などもかねて， 一行ずつ点検． 変数名を変えたので試験運転． どちらも問題ナシ， のように見える． そもそも， すごく単純な改造だし．
• ということで， じつは大サイズではやっぱり落ちていたのか． 成長速度．
• こうなってくると， 例の hyperspecies / species / individual 三階層ベイズ推定とかやったらどうなるんだろう， と好奇心わくわけだが， これはすけぢゅーる破綻という墓穴への直通経路になっているんで， 今回はあきらめることに．
• こちらのモデリングに関してはボス指示待ちだな． うかつに動けん．
• 1910 研究室発． 買いもの． 1935 帰宅． 晩飯．
• 直径成長モデル， さっそく質問いただく． てきとーに回答しようとして …… ふと手が止まる． そういや， random effects の積分が少し手ぬきぎみだったかも． 予想としてはほとんど影響ないはずなんだが， だんだんコワくなってきた． ということで， とりあえずまずはそちらを点検．

	dx.re <- 0.05
	x.re <- seq(dx.re * 0.5, 6.0, dx.re)
	sum(dgamma(
		x.re,
		scale = scale,
		shape = 1.0 / scale # because mean == 1
	)) * dx.re

[いいかげんな積分?]

いいかげんに離散化したガンマ分布を積分した値． scale という分散を定めるパラメーター がでかいところでは 1 からのずれが大きくなる (分散 > 平均² な部分)． まあ， 不幸中のさいわいでこんなに「個体差」分散が大きいのは オオモミジ (AP) だけ， なんだけど …… まあ， きちんと補正するようにしよう．

• [今日の運動]
  • うーむ， 冬はうんどう休養日になりがち ……
• [今日の食卓]
  • 朝 (0900): 米麦 0.7 合． ダイコン・ニンジン・ネギ・ブナシメジ・タラすりみの味噌汁．
  • 昼 (1240): 研究室お茶部屋． ツナマヨパン．
  • 晩 (2030): 米麦 0.7 合． ダイコン・ネギ・チカの味噌汁． ホウレンソウ卵とじ．