ぎょーむ日誌 2004-01-(21-31)
2004 年 01 月 21 日 (水)
-
0740 起床．
朝飯．
コーヒー．
0850 自宅発．
雪．
0900 研究室着．
-
本日の午前のたいはんは
Annals of Botany
近年の目次全部読み & 内容ひろい読みに費されてしまった．
じつは樹木個体にまつわるモデリングにおいて，
私の考えつきそーなことはふらんすの
INRA
な人たち (どういうくくりかただ)
のそれらとぶつかることが多いのかも，
というような．
-
あいまに雑用．
釧路大会公式アカウントあてにきたメイルを何人かに転送する設定．
じつはもっと早くに設定する予定だったんだけど，
昨年中は例によって北大ネットの不調で「転送設定不可能」
状態が続いていたんだよねえ
……
-
午後は甲山さんが講読してる
Functional Ecology
を 2 年分ほど．
やはり「ぱらぱらと全部ながめる」
のは紙媒体のほうがラクですなぁ．
お，
QAIC などというマニアックな手法を使ってたりするぞ．
-
しかし御聖断くだって，
Ecological Modelling とゆーことに．
うーむ．
-
原稿みなおし，
とか．
-
粕谷さんが面倒感あふるる Poisson 分布なハナシを
biometry
に投稿されたので，
現実とーひのため回答を試みる
……
うーむ，
条件とかが明示されてなくて難しいぞ．
で，
あとから粕谷さんがフォローするところによると，
私のやりかた (無理やり一般化線形モデルに帰着させる)
は方向性が違ったようで．
それはうすうすわかってたんだけど，
ほかにやりようないぢゃん，
という気がするンだが．
-
1950 研究室発．
2010 帰宅．
体重 74.2kg
……
これはどう考えてもつじつまがあわん．
運動不足のせいで別の平衡点に移ろうとしているのか．
はたまた塩分とりすぎか?
-
二次元離散分布問題をうだうだと考え続けてしまう．
2400 消灯．
-
[今日の素読]
-
[今日の運動]
-
[今日の食卓]
- 朝 (0810):
米麦 0.7 合．
卵焼．
- 昼 (1330):
弁当．
研究室お茶部屋．
米麦 0.7 合．
チンゲンサイ・タマネギ・ショウガ・豆腐の炒めもの．
- 晩 (2130):
米麦 0.7 合．
ハクサイ・ブナシメジ・豆腐・煮干・昆布の味噌汁．
2004 年 01 月 22 日 (木)
-
0800 起床．
朝飯．
コーヒー．
0900 自宅発．
曇．
気温高い．
0910 研究室着．
-
いきなり Matt の群集生態学問題にまきこまれる．
niche 間距離
……
niche location は無視して niche overlapping だけ考慮した
distance を定義せよ
……
って．
うだうだぎろんしながら，
Solensen 類似度の連続分布版をでっちあげてみる．
うーむ．
-
南米大陸からの依頼にしたがって講座ペイジの一部つくりなおし．
いやはや．
-
さて，
まだあれこれやりのこしたままなんだけど，
次の一手も検討せねばならない．
順番つけるなら
-
苫小牧直径成長問題
-
アカマツ
-
パイプ樹木,
M1 対策,
などなどなどなどなどなど
てところになるのかなぁ．
明示的な締め切りはないのに順番はありそう，
という状況．
どーなることやら．
-
上のよーに書くと，
難問山積感があるので，
作らんといかんプログラムの種類でも列挙してみるか．
てきとーに
……
……
とこころみに書いてみたが，
そもそも
stepAIC()
は汎用である．
つまりこいつに適合するような
(model fitting 用の)
計算関数を定義できるかどうか，
がポイントなのか?
-
うーん，
それは
dropterm() & addterm()
に対応してないといけないわけで
……
うをー，
step() だの
drop1() だの
add1() だの
もあるぞ．
何がどうなってるのやら．
-
MASS 本
つらつらと眺める．
そりゃあ，
まあ線型でなければならんなぁ．
えーと，
非線型版をつくらんといかんのだっけ?
どうなんだろ．
-
おお，
世の中には
cpanflute
なるものがあるのか，
と横道にそれてみる．
で，
なかまさんは
cranflute
作られたわけだ．
-
ひさしぶりに北大構内雪上走．
昼飯．
-
苫小牧直径問題かんがえてたら，
ちょーど鍋嶋さんからメイルいただく．
うーむ，
いろいろススんでるよーで．
そして degree of freedom の限界にいどむモデリングですなぁ
……
などと，
のんびり論評してるばあいでもないのか．
えーい，
どうしたものやら．
-
と思ったら，
また Matt 群集ハナシにトラップされる
……
こんどは Pasoh かよ．
なんといいますか，
もひとり大学院生が増えたような．
かなりばてる．
-
1920 研究室発．
1940 帰宅．
-
仕事ススまん．
うう．
苫小牧直径問題，
どうにかせんといかんのだが，
いまいち霧中模索状態だ．
2300 消灯．
-
[今日の素読]
-
[今日の運動]
-
北大構内走 1310-1340．
雪上を走るのは，
やはりあまりよろしくない．
ストレッチング．
-
腹筋運動 30 ×
3 回．
-
[今日の食卓]
- 朝 (0810):
米麦 0.7 合．
ハクサイ・ブナシメジ・豆腐・煮干・昆布の味噌汁．
- 昼 (1400):
弁当．
研究室お茶部屋．
米麦 0.7 合．
ハクサイ・ブナシメジ・豆腐・煮干・昆布の味噌汁．
- 晩 (2030):
米麦 1.0 合．
ニラ・卵の炒飯．
2004 年 01 月 23 日 (金)
-
0730 起床．
朝飯．
コーヒー．
0820 自宅発．
晴．
0830 研究室着．
-
1000 より
論文セミナー．
本日は長谷川宣教師のモジュール神学のおしえの布教．
ところで，
紹介されてる論文なんだけど，
よく考えられた観測やってるのとは対照に，
データ解析はそのデータ構造 & やりたいことに合致してない．
「まず成長がきまって，余り資源で繁殖では」
と本文で言ってるんだから，
花のカウントデータをシュート成長で説明するような
一般化線型混合モデルがばっちり使える状況なのに
……
たいへんもったいない．
図を見ての憶測だけど結果も変わるだろう．
実際にやっておられることは，
花数→伸長量で「なんでもかんでも等分散正規分布」わーるど，
しかも何故かしら枝次数がランダム変量にされた ANCOVA
……
なぜ生態学者は対象を観察するときと同じぐらいに
注意ぶかくデータのカタチの観察しようとはしないのだろうか．
-
終了後，
モジュール神学の聖典のひとつっつーことになってるらしい
White (1979)
The plant as a metapopulation
なる文献をダウンロードして読んでみる
……
なーんだ，
この題名にある ``metapopulation'' なる用語は
階層化複数集団のことではなく，
「共有された性質 (遺伝子など) をもつ要素からなる集団」
という意味だと明記されているな
(meta- には shared の意味アリとのこと)．
しかも樹木というより，
草本なクローナル植物を念頭においてるよーな．
-
北大構内雪上走．
体重うんぬんはともかく，
カラダを動かさないとアタマが不調になるんで．
院生に聞いたところによると，
冬の早朝には北大構内はクロスカントリースキーの
練習やるヒトたちでいっぱい，
だそーで．
やっぱ積雪期はそっちのほうがいいのかな．
-
もどって昼飯．
ユニさんから東
カリマンタン
北部のダヤック族の集落で
お酒すすめられたハナシなどうかがう．
-
夜まで
……
苫小牧直径データのモデリングはどうなるんだよぉ
(これを放置することは背筋の冷たくなる恐怖感を伴う)．
しかし，
群集パラメーター推定のアヤしさはさておき，
Pasoh 樹木たち 26 万本の直径→明るさ関係の推定結果
(log-normal 分布，平均・分散それぞれは exp(三次関数)
とかでいーんぢゃないの
……
と言ったらマにうけられてしまった)
見てると勉強にはなりますなぁ．
おもに
R
の挙動に関する面でだけど．
-
ということで，
これもまた他の下請け仕事にも役に立つんですよ，
きっと
……
と誰に弁明しているのだろうか?
-
今日もかなりばてた．
1820 から講座内パーティー．
ユニさん歓迎 & 岳明さん昇進祝い．
セミナー室で鍋．
-
2410 研究室発．
2420 帰宅．
-
2500 消灯．
-
[今日の運動]
-
北大構内走 1300-1330．
ストレッチング．
-
[今日の食卓]
- 朝 (0740):
米麦 0.7 合．
- 昼 (1350):
弁当．
研究室お茶部屋．
米麦 0.7 合．
- 晩 (1830):
講座内パーティー．
2004 年 01 月 24 日 (土)
-
0820 起床．
朝飯．
コーヒー．
-
北海道ハワイ化計画
の進捗の遅れに関して，
平八郎提督よりおしかりのメイルいただく．
また新しい神算鬼謀というべき秘策が
……
ところで先日、長野県の田中知事が、県名を信州にするなどと申し立てて世間の失笑を買ったという風説はシベリアまでは届いておりますまいが、そんなことがあったのである。
貴殿もお察しのとおり、これは単に県名がどうのこうのという問題ではないのだ。長野県は一気にアメリカの州になるべく運動を始めているのだ。そのうち北アルプス、中央アルプス、南アルプスが西ロッキー、西アパラチア、西シエラネバダ山脈などと改称される事は火を見るより明らかである。表銀座・裏銀座も、表ウオール街・裏ウンターデンリンデンになってしまうことも容易に予想されることである。さらに言えば、松本も信濃ツェルマットと改められることも考えられよう。
それに比べ北海道はどうなっているのかとおもへば、全く進展が見られんのである。しかもせっかくの兵力が中東に裂かれつつある。全然いかんのである。案ずるに、知事にやる気がないのであろう。リコールすべきである。
ここで秘策がある。次の知事にクリントンをあてるのだ。名付けて、「なんとなく、クリントン」作戦だ。これで田中康夫がずっこける間に一気に話を進めることが出来そうである。
今がチャンスだ、久保君。
-
「道北侵攻食い止め専門部隊」
だった第二師団もハワイ化のために中東に抽出されちゃうことだし
(ただしハワイ化構想抵抗勢力に対する欺瞞工作のため，
あたかもオランダの手下であるかのように偽装されている)，
いまは「雌伏」のときでは，
と愚考を申し上げたところ，
[揮毫]
いただく．
午海斎なる雅号のよーで．
-
1640 自宅発．
曇．
途中，
ICI スポーツで山スキーのシールの留具
「ワイドストレッチャー」
(ゴム製，幅 80mm，2 個 1 組)
税込み 1701 円購入．
1700 研究室着．
-
さっそくシール修繕．
ついでにスキーを 7F 実験室にもっていき，
板の末端にノコギリで「コ」の字型のシール固定溝を作る．
なかなかに悲惨だった
奥手稲
山行の戦訓を反映している．
-
ぢりぢりと苫小牧直径問題の手さぐり
……
-
2230 研究室発．
2240 帰宅．
体重 74.0kg．
ふーむ
……
-
晩飯．
2430 消灯．
明日はなぜか
春香山
にいく予定．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0830):
米麦 0.6 合．
卵炒飯．
- 昼 (1340):
米麦 0.6 合．
レトルトパウチドもの麻婆春雨．
- 晩 (2310):
米麦 0.7 合．
チンゲンサイ・モヤシの炒めもの．
2004 年 01 月 25 日 (日)
-
0700 起床．
朝飯．
コーヒー．
0810 自宅発．
快晴かつ寒い．
0820 研究室着．
-
春香山
に山スキーに行くんで装備の準備．
あ，
ヘッドランプ忘れてた
……
2000 年に春香山に
行ったとき
は日没ぎりぎりまで行動してたんだが．
-
0900 に集合のはずなンだが
……
0908 隊長たるかとー先生よりお電話
「昨晩のみすぎたんで，
いま起きたばかりです」．
0930 研究室発．
平尾君の車に荷物を積んで出発．
ちょっと我が家によってもらう．
ヘッドランプをとってもどる．
0948 北 36 の麻生自動車学校でかとーさんと合流．
ちょーど出てきたところだった．
小樽方面へ．
-
桂岡のセヴンイレヴンで買い物．
住宅地を通りすぎてすぐのところにある駐車場着．
時刻は 1040．
うーむ，
先行登山者，
とゆーか先行スノーモービル乗りたち多数のようで．
1055 出発．
スキーにシールをつけて林道を登る．
快晴．
[始点 (標高 124m)]
なんかしばらくは除雪されてた．
そっから先は先にいったヒトたちがふみかためていて，
ですね．
ラッセルとかぜんぜんなかった．
[林道 (標高 124-570m)]
二日酔いに苦しみながら登高してくるかとーさん．
それはともかく
……
このルート，
以前に一度は往復したことあるはずなのに
何も思い出せん．
[土場 (標高 570m)]
林道終点．
石狩湾がよく見える．
湾の向うは暑寒別岳．
遠くに夕張山地・大雪・日高らしきものまで．
じつに穏やかな冬の日．
ここに来たのは覚えてるなぁ．
さすがに．
クリックすると写真拡大．
[銀嶺荘 (標高 700m)]
休日のみかもしれんけど冬もやってる山小屋．
春香山は札幌から近くて，
スキーには手軽な雪山なんで，
この時期でも登山者が多い．
-
1220 土場 (林道終点) 着．
ちょっと昼飯くってから，
1240 同発．
1320 銀嶺荘着．
1330 同発で，
1400 春香山山頂着．
[春香山山頂 (標高 906m)]
3 時間ほどのゆったりとした登高で山頂着．
石狩湾だけでなく石狩平野にひろがる札幌市もよく見える．
上の写真右はおつかれかとーさん．
左の写真手前は平尾君．
-
1430 山頂発．
雪の状態よく私以外の二人はくるくるとターンしながら，
頂上直下の斜面を滑べり降りていく．
私は例によって斜滑降 + キックターンでぢりぢりと．
1450 銀嶺荘．
スノーモービル通過あとにはまると面倒．
1520 土場．
おお，
前回の半分の時間でここまで戻ったか．
-
ここから林道なんだけど，
ここもまた今回はラク．
前回はがちがちに凍結してたんだけど，
今回はそんなに雪面が硬くないので，
ボーゲンなる私が使用できる唯一の非破壊的制動技術が使えるため．
かとーさんによると，
いいとこばかりのテレマークスキーもこういう広くない林道は苦手で
(ボーゲンのごとき制動ができないため)，
くたびれるとのこと．
たしかに，
こういう場面だけは，
歩行に不便なほどに硬い山スキー靴は便利かも．
-
あと，
スキー下手な人にとって背負ってる荷物が軽い，
というのはたいへんに重要．
コケる回数がぜんぜん違ってくるし，
コケたときの起き上がりやすさもぜんぜん違う．
-
1550 駐車場．
あいかわらず晴．
年末の奥手稲のあれに比較すると，
実に平和な山行でした．
天気も穏やかだったし．
1600 同発．
1645 研究室着．
-
1830 研究室発．
1840 帰宅．
体重 73.8kg．
晩飯．
-
ばててるのに
Linux カーネル構築試行錯誤にはまってしまう．
2500 消灯．
-
[今日の食卓]
- 朝 (0715):
米麦 0.7 合．
チンゲンサイ・モヤシの炒めもの．
- 昼:
山に登りながら氷砂糖なめたり．
- 晩 (2000):
米麦 1.2 合．
チンゲンサイ・ハクサイ・ショウガ・マイタケの炒めもの．
2004 年 01 月 26 日 (月)
-
0800 起床．
朝飯．
コーヒー．
0910 自宅発．
晴．
0920 研究室着．
-
さーて，
累積量がなぜか減ってくれない下請け仕事を片づけていきましょうかね
……
と思ったら，
また Linux カーネルまわりの問題
(不要不急)
にトラップされてしまった．
ひとくぎりついたのが 1100．
ひたすらにまぬけだ．
-
苫小牧直径成長問題に取り組む．
とりあえず R の関数でも作るかなぁ
……
-
北大構内走．
曇天．
走り終ったとたん雪がふってきた．
昼飯．
-
R
で推定プログラムを作ろうとしてんだけど，
ススまない．
計算式，
とゆーか確率分布のカタチは頭の中ではっきりしてんだよね．
べつに難しいもんでもないし．
プログラミングに必要と思われるていどのメモ書き
(といっても数式が列挙されてるだけなんだが)
もやった
……
しかしススまない．
思うにその数式を構成してる記号たちの定義，
notation があいまいなままなんだよねえ．
それが気になってススまないのである．
不思議なことに．
-
とゆーことで，
というか
……
まあ，
いそがば回れっつーことで，
メモを直すところから着手．
-
ふーむ．
実は世の中にはふたとーりの幾何分布があることを発見．
-
確率変数
x
が非負の値となる
P(x) = p (1 - p)x
-
確率変数
x
が正の値のみとなる
P(x) = p (1 - p)x - 1
R
やたいてい本では 1. なんだけど，
(統計学本いろいろ出してる)
蓑谷さんは 2. 派だな．
というのも，
幾何分布の定義が「成功確率 p のベルヌーイ試行において，
初めて成功するまでの回数」
となってるんで．
そりゃそれが本来の定義かもしれないけど，
ゼロから始まるやつでないとあれこれ不便な場合も多々あるわけで．
-
とうだうだと．
-
1810 研究室発．
駅の南の銀行 ATM へ．
帰りは本屋をふらふらと．
2000 帰宅．
体重 73.8kg．
-
昨日，
PDA 工房から新しい Clie T600C が届いた．
前の中古品はなぜかしら輸送途中に完全に壊れてたらしく，
今回はまた別の (前よりも状態よさそうな)
中古品を送ってくれたんで
……
で，
これと Linux の HotSync でコケてるのである．
カーネルに
visor モジュールを
自動的に読ませるところまではできたんだけど，
そこから先がうまくいかない．
-
しかし今晩はそこで時間とられてたわけではなく，
なぜか当家の炊事用流しのまわりの掃除で夜がふけていった．
すくなくとも一年以上は放置してあったので，
なかなかに楽しくない状況になってたもんで
……
-
2500 消灯．
-
[今日の運動]
-
北大構内走 1300-1330．
ストレッチング．
-
[今日の食卓]
- 朝 (0820):
米麦 0.7 合．
チンゲンサイ・ハクサイ・マイタケの炒めもの．
- 昼 (1350):
弁当．
研究室お茶部屋．
米麦 0.7 合．
チンゲンサイ・ハクサイ・マイタケの炒めもの．
- 晩 (2200):
米麦 1.2 合．
ハクサイ・タマネギ・ホタテの炒めもの．
2004 年 01 月 27 日 (火)
-
0800 起床．
朝飯．
0900 自宅発．
晴．
0910 研究室着．
-
うーむ，
幾何分布の簡単な計算ができん．
ばかだ．
-
などとやってたら，
釧路大会シンポジウム申し込み．
申し込みペイジ
→
jes51@ees.hokudai.ac.jp
→ あちこち，
という経路になってる．
それはいいんだが，
釧路大会実行委員会メンバー
とやらの誰が返信すべきなのかわからず，
とりあえず関連しそうなヒトたちにつっこみメイル．
-
新型 virus 発生，
とゆーことで，
本日は感染メイルが多数とどく．
-
北大構内雪上走．
昼飯．
-
夕方ごろになって，
自分の阿呆さかげんに気づく．
幾何分布に関するごくごく簡単な計算やってるだけなのに，
記憶の中にある統計量と一致しない，
という問題だったんだが
……
なんべんやっても同じ答えになる計算が正しく，
記憶のほうが間違ってた．
正確にいうと，
昨日のぎょーむ日誌に書いてた「二種類の幾何分布」
の区別がついてなかった．
これはたんにゼロから始まるか否か，
だけではなく，
平均・分散といった統計量もちょーっと違って来るんだよね．
そんなことにも気づかなかったとは
……
-
疑問は氷解したけど，
脱力したんでお茶部屋で意外性追求の雑談．
「マイルド」 (非マニアック)
パズルファンには意外と女性が多いんではないか，
北大ではない某大学の生態学研究室のせんせーがたは
意外と計算機ゲイム好きなんでは，
というような．
-
計算のツメ．
両側幾何分布 (two tailed geometric distribution か?)
は分散がふつーのやつの 2 倍になります，
という結果．
これって当たりまえっぽいんだけど，
証明に至る過程
(場合わけして個別に計算していき最後に組み合わせる)
はあまり直感的でないものになってる
……
計算の見とおしがよろしくなかったんだろうか?
-
1940 研究室発．
2000 帰宅．
晩飯．
-
こちらのペイジ
など参考にしながら，
とりあえず赤外線ポート経由で，
PalmOS 機と Linux Thinkpad の間で通信できるようにした．
要点は
-
BIOS でちゃんと IrDA 使えるように設定しとく
-
irda-utils
は必須
(version
0.9.14-9.src.rpm
かこれより新しいものを探す)
--
なンで vine にはないんだろうねえ?
-
/etc/sysconfig/irda
で
DEVICE=/dev/ttyS1
と設定
--
これはナゾ
-
起動は
/etc/init.d/irda start
でよい
--
起動したら
/sbin/ifconfig
で irda0
が出るかどうか確認
-
port は
/dev/ircomm0
--
あらかじめ mknod & chmod
しとく
しかし，
まあ赤外線通信はおそいのう．
-
USB 経路を探索してるうちに夜がふけてしまった．
2500 消灯．
-
[今日の運動]
-
北大構内走 1300-1330．
ストレッチング．
-
[今日の食卓]
- 朝 (0820):
パンケイキ．
- 昼 (1350):
弁当．
研究室お茶部屋．
米麦 0.7 合．
ハクサイ・タマネギ・ホタテの炒めもの．
- 晩 (2130):
米麦 1.0 合．
ニラ・ニンニク・卵の炒めもの．
2004 年 01 月 28 日 (水)
-
0800 起床．
朝飯．
コーヒー．
0900 自宅発．
曇．
0910 研究室着．
-
えーい，
朝から
……
Linux Thinkpad のファイルシステムを，
とりあえず
ext2 から
ext3
に移行してみる．
ext3 は遅くて好きではないんだけど，
このちょー高速 Thinkpad ならなんとかなるだろう，
ということで．
このペイジ
など参考にしながら
tune2fs & mkinitrd
使って
……
と思ったら，
久保設定カーネルでは起動 ramdisk 使えないようになってたんで，
朝からカーネル再構築．
cat /proc/mounts
で ext3 でマウントされてると確認．
-
ぢゃーなりんぐファイルシステム，
ってなんとなく server 業者なヒトたちのために作られたもの，
と思ってたんだけど，
ここ数日の体験からカーネルやその周辺を
やたらとぶち壊してしてまう
(もの好きなというか奇特なというか)
開発者連中のために作られたんでは，
という気がしてきた．
-
そして今日も多数の virus メイル来信．
-
ほぼ一日，
複合分布問題のメモ書き
(というほどには長くない)
に費されてしまった．
これよーやく関数コード書けるか?
-
査読依頼．
呪われわーぷろソフトウェアなフォーマットのファイルを
飛ばしまくるのやめてくだされ．
OpenOffice
の解毒能力が向上してるとはいえ
(数式オブジェクトがうまくほどけてるのには感心させられた)，
図とかはまだまだなんだよねえ．
-
1630
Trendy セミナー
，
本日はこの建物 1-2F (地学系)
地球生態
講座の Gautam さん，
ネパール低地林 Shorea
な木の人間による利用と個体群動態的な関連について．
まあ，
生態学っつーより林学 & 文化人類学だな．
ethnosilviculture
(これ新造語っぽいな
……
起源はethnobotany
あたりだろう)
とか言ってるし．
興味ぶかくはあるけど，
樹木成長の解析としても文化人類学としても
まだまだきっちり仕上げる余地はありそう．
-
1840 研究室発．
1900 帰宅．
体重 74.2kg．
うーん，
なぜか重い．
晩飯．
-
2300 消灯．
-
[今日の運動]
-
北大構内走 1300-1330．
ストレッチング．
-
腹筋運動 30 ×
3 回．
-
[今日の食卓]
- 朝 (0810):
米麦 0.7 合．
卵焼．
- 昼 (1350):
弁当．
研究室お茶部屋．
米麦 0.7 合．
コンブ．
- 晩 (2000):
米麦 1.0 合．
ハクサイ・タマネギ・コンブ・煮干・油揚の味噌汁．
鶏レバ総菜．
2004 年 01 月 29 日 (木)
-
0800 起床．
朝飯．
コーヒー．
0900 自宅発．
晴れてるのに雪．
0910 研究室着．
-
Poisson + 両側幾何分布の
R
関数書き．
-
昼すぎになってよーやく両側幾何分布の密度分布関数
dttgeom()
と両側幾何分布乱数
rttgeom()
ができた．
[両側幾何分布]
r^abs(x) (1 - r) / (1 + r)
てなところで．
-
北大構内走．
走ってるあいだにも晴れたり雪がふったり．
昼飯．
-
北米大陸東側を制圧して無事に戻られた松木親分から，
北大山スキー生態学者どもに
豊羽鉱山
周辺スキー登山のお誘い
(目標，無意根山?)．
緊張しつつ「参加させてください」請願をだすと，
「今度はリベンジ（？）ということで、万全の体制で望めたらいいですね」
とのお返事．
眼光紙背に徹すれば，
「今度しくじったら，どうなるかわかってるんでしょうね……」
という意味に違いない
……
そうに決まってる．
-
いつものことながら，
R の作画問題でてこずる．
これって便利なのか不便なのか．
-
とりあえず複合分布のパラメーターを
R の
optim()
で推定できるような example を作る．
推定結果は悪くない．
-
さらに，
樹木の成長量を説明するモデルに
適用できるように拡張していく．
[離散成長量モデル]
のパラメーター推定．
これは乱数で生成した偽造データを使っている．
気のせいか，
連続成長量モデルよりパラメーター推定が難しいような気がする．
また，
真のパラメーター値に依存して，
系統的な偏りが発生しているような気がする．
そのあたりはデータの偏りにも
とうぜんながら，
推定値の間には相関が生じてしまっている
……
やはりパラメーター推定は難しい．
-
一年ちょい前
にも同じよーなことやってましたなぁ
……
このときは連続分布 (ガンマ分布 + 正規分布)，
数値最尤推定計算は Perl でやって，
作図は gnuplot してた．
たしかに何もかも R でできるんなら，
大学院生などにとってはラクなことなのかもしれない．
-
しかしながら，
ですね．
今回も計算手順をよく考えて工夫してるつもりなんだけど
……
R 上での数値計算はやや遅いんでは，
と思うわけですよん．
いや，
ぜったいに遅い．
Perl と比較してもなお．
まあ，
私が工夫してるつもりでも，
それがちっとも R に最適化された計算手順になってない可能性は
かなり
あるんだけどね．
-
う?
いつのまにか 8F に煙が．
いや 7F もか?
……
どうも 6F の某講座実験室でぼやかなンかがあったよーで．
とりあえず拡大してはいない．
1800 ごろだったかな．
大学院生たちはエレベータで避難
(煙が有毒っぽいかんぢだったんで)．
私はめんどーだったんで，
かとーオフィスの窓をぜんぶあけて居座る．
-
推定計算がうまくいく場合とうまくいかない場合，
というのをいろいろ生成してみる．
計算ぢりぢりぢり．
-
2120 研究室発．
2130 帰宅．
体重 74.2kg．
うーむ．
晩飯．
-
2400 消灯．
-
[今日の素読]
-
仕事のちょーし良くても悪くても復旧できぬこの不思議さよ
……
やはり朝にやるべきなんだろうな
-
[今日の運動]
-
北大構内走 1330-1400．
ストレッチング．
-
[今日の食卓]
- 朝 (0810):
米麦 0.7 合．
ハクサイ・タマネギ・コンブ・煮干・油揚の味噌汁．
- 昼 (1410):
弁当．
研究室お茶部屋．
米麦 0.7 合．
ハクサイ・タマネギ・コンブ・煮干・油揚の味噌汁．
- 晩 (2230):
米麦 1.0 合．
ハクサイ・タマネギ・コンブ・煮干・油揚の味噌汁．
2004 年 01 月 30 日 (金)
-
0800 起床．
朝飯．
コーヒー．
0850 自宅発．
曇．
0900 研究室着．
-
1000 論文セミナー ，
本日は奥田君が bivariate な
(と論文中では書いてるが，そういってよいのか?)
空間統計学論文で，
個体のランダム配置とランダム名札づけを的確に区別せよ，
って内容．
そりゃ，
区別は重要ですけどね
……
動物集団についてはわからんのだけど
(しかし昆虫などでは距離法はあまり有効ではあるまい)，
植物の集団に関してはランダム配置生成のモデルって，
あまり使いどころないんぢゃないの?
「一様」に見える場所でもじつは住適不適があるから，
標本数が増えるとたちまち Poisson 分布からずれると思うんだよねえ
(むしろ空間解像度のあらい区画法むきなのか?)．
ほとんどの用事はランダム名札づけ (random labeling)
ですむんじゃないでしょーか．
-
北大構内雪上走．
北 12 「北欧」パンで昼飯調達．
ついでに ICI スポーツで山スキーシール糊 (グルー)
のチューブ 1 本 1260 円購入．
もどって昼飯．
-
本日午後は PlotNet なる
ゑくせるファイル交換同人の集まりをそのへんでやってるよーなんだけど，
うかつに出席すると死ぬような雑用ばかり押しつけられることになるので逃げる
(夕方すぎ，
気の毒なるかとー先生はまさにそういう境遇におかれていた)．
苫小牧直径成長モデリング問題に関して，
おそるべき勢いで解析していく鍋嶋さんメイルにお返事かきなど．
-
大澤君とシウリザクラ雑談．
苫小牧シウリザクラの現時点において残存してる
未解決かつ重要なナゾとして，
永光さん解明された遺伝子型を投影してみると，
genet なるものが排他的なかたまりを形成してるよーに見えるんだけど，
これってホントにそうなのか?
そうだとしたら何故なのか?
種子由来の個体はホントに巨大 genet 内部に定着できないのか?
……
というあたり．
じつはこれは大澤調査において，
各 ramet に対する
「所属 genet」
なる名札づけがまっとうだったかどうか，
を左右する現象なんだよねえ．
-
7-8F でヒマそーに見えるのは私だけであったために，
昨日の「A 棟けむり蔓延事件」に関する説明うんぬんに出席せよ，
との指令を雪野さんからいただく．
これは抗命不可能ってやつですね．
1600 事故現場たる A606 実験室へ．
説明された事故の経緯としては
-
北大のあちこちの実験室で draft chamber 新型化の工事
がいっせいにおこなわれている
(あとで日浦さんに聞いたところによると，
全体で数億円を費しているとか)
-
これは地環研事務ではなく大学本部施設課の担当だった
-
工事の日程がおしていた
-
古いダクトを交換する作業と新 chamber 導入準備のための
天井鉄骨切断工事を同時並行してやっていた
(ふつーはやってはいかんらしい)
-
鉄骨切断の火花が旧ダクト内に飛びこみ，そこに残留して
いた「何か」(不明) に引火して小火災が発生した
(ふつーはそうならないように防災シートなるものをつける
んだけどやってなかった)
-
工事現場には消火器がなかったので 6F 備え付けの消火器
2 本を消費して火災消火した -- われわれの研究室まで
ながれてきた「けむり」らしきものの正体はこの消化剤
のものらしい (?)
-
事件発生は 17:30 で地環研事務の人間がおらず，本部のヒト
は全館アナウンスのやりかたがわからなかったので，連絡
がまったくないまま事件が終った
対策，のつもりらしいこととしては，
火花がとびちっても火災にならんよーに，
今後は防災シートなどをつけるようにする，
とか．
最後の質疑応答で誰も指摘しなかったんで，
「今回の事故発生時に棟内にいる人たちに一切の連絡なかったので，
時間外 (夜間 or 休日)
にも全館アナウンスができるように手配してもらえませんか」
と私が言ってみた．
すると，
今まで大学本部事務の責任の糾弾に熱心だった地環研事務官がやや狼狽して，
その件に関しては打ち合わせましょうと本部事務官に協調を求めたのであった．
ゲイム理論的，
というべき状況なんだろか?
-
あと，
講座が注意すべきこととしては，
実験室などで工事がおこなわれるときに，「火花が飛び散る」
「消化剤ばらまかれる」事態が発生しうるので，工事前には
可然物 危険物のたぐいは撤去し，消化剤をかぶったら困る
機械はビニールシートなどで保護すること，
ということのよーで．
-
小菅せんせーがまた隠匿してくれましたよ
……
8F 闇ネットの IP アドレスに使えないものがある，
と思ったら小菅せんせーの隠匿計算機がこっそり使用していたのであった．
-
biometry におけるやりとり．
粕谷さんメイルから引用させていただくと:
久保さん：
三中さんはなぜ「ぶーとすとらっぷ」ではなく「ぶーつすとらっぷ」
と表記されるのでしょうか? 何か謂れ因縁故事来歴物語がありそう
な気がします．
三中さん：
単なる日本語の感覚的〈慣れ〉の問題でしょうが，「ブート」と言われたよりも
「ブーツ」の方が〈靴っぽい〉語感がありません？（主観的なのは重々承知）
私も以前から気になっていました。他の分野ですが、bootlegもブ−ツレッグとブ−
トレッグと両方見られるようです。統計にあまり関係なくなってすみません。
私は外来語であっても日本語としての語感を重視する，
という三中さんの立場に意表をつかれた思いであった．
蛇足ながら検索してみると，
今日ではブートレッグのほうが使用例が多く見つかるよーな
……
-
1800 すぎ，
PlotNet 会議も終ったよーなんで，
遠方から参加された皆さんと JR 札幌駅北西地下の居酒屋へ．
-
2020 研究室にもどる．
山スキーシールの糊 (グルー) はりかえの試み．
うまくいくかどうか．
こちらのペイジ
など参考にしてアイロンの準備なんかはしてみたんだけど．
どうも古いグルーを取り除くのが難しい．
技法の問題かアイロン温度の問題か道具の問題か．
とりあえず新しいグルーをのせてアイロンで圧延するだけにしよう
……
これまた難物で，
「糊がついてもくっつかないシート」
が必要である．
手近にあるいろんな素材をためしてみた．
最良だったのは，
ガムテイプを紙に並列に張りつけていってシート状にしたもの．
これって糊面につけてもキレイにはがせる．
-
PlotNet 二次会後の清野さん・西村さんがお茶部屋にみえる．
とくに清野さんは絶好調モードであり，
2400 まぢかに妙技「終電の魔術師」をひさびさに披露してくださった．
-
2530 研究室発．
2540 帰宅．
-
[今日の運動]
-
北大構内走 1250-1320．
ストレッチング．
-
[今日の食卓]
- 朝 (0810):
米麦 0.7 合．
ハクサイ・タマネギ・コンブ・煮干・油揚の味噌汁．
- 昼 (1410):
研究室お茶部屋．
「北欧」ベイグルサンドイッチ．
- 晩 (1820):
PlotNet 宴会につきあう．
2004 年 01 月 31 日 (土)
-
0830 起床．
朝飯．
コーヒー．
-
慢性週末無気力症候群
……
時間がおそろしく浪費されてるよーな．
いや，
これも休養のうちなのかしらん?
-
spamassassin
なるすぱむフィルター使ってんだけど
……
これって安直な仕掛けのわりにはそこそこには見分けてくれる．
しかし非すぱむなメイルも削除されたりする．
HTML メイルなんかは無視できない確率で誤排除されてる
(排除されても復旧の可能性はあるんだが，
「油断」してると復旧不可能状態になる)．
まあ，
相手かまわず HTML なメイル送るってのもなぁ．
-
この問題って，
経路上の最終防衛線たる自製メイルフィルターを少し改善すれば
うっかりミスによる誤消去の可能性を減らせるな，
と気づいたんで，
コードをなおしてみる．
すこしは生産的なことができたか?
-
怠業といえば
……
長谷川宣教師殿いわく「ついに免責」
とのことで，
http://www.freeml.com/message/morphy-members@freeml.com/0000036
だの
http://ecos-h8.sourceforge.jp/sogwipo/pukiwiki.php?MorphyWiki
だの
……
ぢつに教訓的な事件でありましたなあ．
-
1620 自宅発．
雪．
古本屋だの Homac だのふらふらと．
1830 研究室着．
装備を整える．
地図は
地形図閲覧システム
からダウンロード．
画像処理．
-
メイル書きとか．
計算メモ書き，
ほんのすこしだけススめる．
-
2030 研究室発．
北 12 生協で買い物．
2100 帰宅．
-
ということで，
明日は
無意根山
「方面」に豊羽鉱山から入山だ．
メンバーは綿貫夫妻，松木さん，原拓史さん，
というじつにコワい顔ぶれだ．
明朝 0800 完全装備にて集合
……
と思ってたら 0730 に変更との連絡あり．
うー，
山頂を狙ってるな，
このヒトたち．
さっさと寝るべきなんだが，
日ごろの悪しき習慣でなかなか眠れん．
すでに 2420 すぎ．
-
[今日の食卓]
- 朝 (0850):
パンケイキ．
- 昼 (1330):
サツマイモ．
リンゴ．
- 晩 (2130):
米麦 1.0 合．
ハクサイ・卵の炒飯．
