ぎょーむ日誌 2005-05-(21-31)

2005 年 05 月 21 日 (土)

• 0700 起床． 朝飯． コーヒー． 怠業． 1200 自宅発． 晴． 1220 研究室着． 昼飯． フェアバンクス 近郊の調査地設定からもどってきた女性大学院生 (北海道民) アラスカ話などいろいろとうかがう． なんでも
  1. 北海道と内地の「類似度」
  2. 北海道とアラスカの「類似度」
  を比べると， 2 のほうが「ぜんぜん近い」そうだ． いわく， アラスカ人が北海道弁をしゃべれば， もはや両者に何のちがいもない， と．
• アカマツ原稿の数式 (10 個弱) をわーど98 for 旧 MacOS にて入力する作業． 甲山さんから PowerBookG3 かりて作業する． なんか， 数式おぶぢぇくとって扱いがむちゃくちゃにヘンだな ……

[PowerBookG3]

OS は MacOS 8.6 が動いておりました． いまさらながら， Equation Editor にはいらいらさせられる． たぶん使ったのは 6-7 年ぶりだ．

• 数式入力， ファイルアップロード， メイル連絡 …… で 1500 ごろ終了． ふう．
• 1540 北大構内走． 晴． 気温高い． 樹木の開葉はじまっている． 1625 もどる． ばてた．
• アカマツ原稿， 投稿用の図の準備 …… これは以前に作った自動作図プログラム (in R) の修正などなど， で簡単に片づいた． 作業はカンタンだけど空腹なので時間かかる． LaTeX に \includegraphics でよみこむ． PDF ファイル化． アップロード． メイルで連絡． 1850 終了． 下の図は R の pch で指定する作図シンボルの一部 (example(points) で表示される)．

• はらへったので撤退． 1920 研究室発． 1940 帰宅． 体重 74.4kg． 晩飯．
• [今日の運動]
  • 北大構内走 45 分間．
• [今日の食卓]
  • 朝 (0730): 米麦 0.5 合． キュウリ． ハクサイ・ニラ・トウフの炒めもの．
  • 昼 (1220): 研究室お茶部屋． 「北欧」そうざいパン (カレー & 明太子)．
  • 晩 (2030): マカロニ． ニンニク豆板醤ソース． サバ酒蒸し．

2005 年 05 月 22 日 (日)

• 1020 起床． 朝飯． コーヒー． 怠業．
• Yahoo! オークションで落札したジョギングシューズとどく． 新品 (というか未使用というか)． とあるスポーツ用品店が出品． どうもこういうのは「2004 年版」とかあるらしく， 「カタおち」になるとオークションに安くでまわったりするんだよね ……

[asics GEL-1090-SW (TJG856)]

トレイニング用， つまりクッションが厚いタイプ． これで走ってみると， たしかに足への負担がへったような気にさせられる． 送料込 5800 円．

• 1420 自宅発北大構内走． 晴． 1530 帰宅． 体重 73.8kg． 昼飯．
• 1630 自宅発． 曇． 1645 研究室着．
• 「補償」のほうのアカマツ原稿はひとくぎりついたんで， その残骸たるモデリング部分のほうもなんとか …… と思うわけだが， それとは関係ない論文あれこれを読んでみたり．
• 2100 研究室発． 2130 帰宅． 晩飯．
• [今日の運動]
  • 北大構内走 70 分間．
• [今日の食卓]
  • 朝 (1040): コッペパン．
  • 昼 (1610): マカロニ． レトルトパウチドカレー．
  • 晩 (2250): 米麦 0.8 合． 豆腐． タマネギ・ニンニク・パセリ・豚肉のトマト煮こみ．

2005 年 05 月 23 日 (月)

• 0700 起床． 朝飯． コーヒー． 0815 自宅発． 曇ちょい雨． 0830 研究室着．
• 道新 (北海道新聞) 一面はやはり新函館までとどく新幹線着工 (完成は 10 年後とか --- 札幌まで届くかどうかは未定)． 採算性アヤしい辺境新幹線ごときに かくのごとく一喜一憂する背景には， 「内地」に関する北海道の不思議な二律背反があるように見える．
• アカマツ「補償」アロケイションな原稿が 小林さんのところでまとまりつつある一方， その「残骸」とでも言うべきアヤしげ成長モデル …… さるべいじ作業にとりくんでみる． まさにいまやらねば永久に放置してしまいそうなんで． はい， コレのことです．

• あ． 更新すんの忘れてた．
• アカマツモデリング， 説明の順番いれかえたりいろいろ削除したり …… かなり乱暴にりすとらしてみる．
• あいかわらず雨． 昼飯．
• 北大ネット， DNS もしくは proxy server まわりがおかしい． 動作が不安定．
• ふーん， 世の中には Environmental and Ecological Statistics などというぢゃーなるもあるのか …… いまヒネってるアカマツモデルは投稿できそうにないけど． というのも， これはそれほど statistical な見かけではないため．
• BibTeX おそるべし …… 文献データベイスの書式まちがえると， かなりスサまじい無限ループにおちこむ． Ctrl+c は通用せず Ctrl+z で make を suspend して killall でやっつけた． しかも defunct すなわちゾンビになったりするし …… 親プロセスの make を fg にして終了すると死んでくれたけど．
• で， ぢりぢりと修正作業つづける …… うちに， ばててしまったので， 現実逃避のため Wikipedia 上をうろうろするうちに， なぜかしら George Dantzig の項にいきあたる． つい十日前に亡くなっていたとわかった． ふーむ．
• 1900 研究室発． 雨あがってる． 1915 帰宅． 晩飯．
• [今日の運動]
  • うんどう休養日になってしまった
• [今日の食卓]
  • 朝 (0720): 米麦 0.5 合． タマネギ・ニンニク・パセリ・豚肉のトマト煮こみ．
  • 昼 (1330): 研究室お茶部屋． 米麦 0.5 合． タマネギ・ニンニク・パセリ・豚肉のトマト煮こみ．
  • 晩 (2030): 米麦 0.8 合． コンブあえもの． チンゲンサイ・ピーマン・ネギ・卵の炒めもの．

2005 年 05 月 24 日 (火)

• 0850 起床． どうも， あいかわらずヘンだ． 朝飯． コーヒー． 0930 自宅発． 曇． 0945 研究室着．
• 朝から北大ネットは絶不調． まあ， よくあることだが． そして例によって何の役にもたたない掲示が 「外からは見えない」 場所に何やら書いてる:

  	2005/05/23午後〜24日にかけて、ネットワーク遅延が頻発しております。
  	ウィルス、ワーム、スパイウェア等に感染・活動している端末が学内に数台存在していると
  	思われますが、原因の特定が困難な状態です。
  	各ユーザ様におかれましては、管理端末のウィルスチェックなどの徹底をお願い致します。
  

  毎度毎度のことなんだから， さっさと不審 request 発信源ぐらい特定しろよ． すげー中央集権的に TCP/IP の流れを制御してるだから， それぐらい簡単にできるでしょ．
• M1 の大学院生からいきなり (調査地間移動を最小にする) 「巡回セイルスマン問題」 をとけ， と言われてあわてる． 検討してみると， 必ずしも「最短経路」でなくてもよいことがわかった． ひと安心． 蛇足ながら R の最適化関数 optim() の example の最後に， この問題を SANN (simulated annealing) で近似的に解くデモンストレイションがある (example(optim) で下のような結果が表示される)．

• 午前中は本日の「(文献紹介の) セミナーを聞く準備」に費される． つまり論文よみ． (先週とおなじく) かかる準備が必要となり理由は …… 紹介者が内容を正確に要約してくれるかどうか， かなりアヤしいところがあるから． 時間はとられるけれど， ダマされるよりはマシ， ということで．
• またメイルサーヴァー不調か …… くだんの「北大ネット隠しざんげペイジ」 (ざんげではなく単なるぐちかもしれん) では

  	なお、この遅延に関連してメール配信処理にも支障を来す場合があります。
  	メールアプリケーションで送信が完了している場合でも、北大メールサーバ内で
  	timeoutを起こしている可能性がありますので、ネットワーク遅延解決後に
  	配信確認をお願い致します。
  

  といつものごとく趣旨不明なことを． 送受信どちらもダメ， ということで Mail Gateway Server とやらが腐れてるんでは (かとー説)． いやはや． 昼飯．
• 1300 文献紹介の セミナー， 本日は M2 石橋君， 樹木のカタチと機能の関係について． やはり予想どーり当方の十分なる予習は必須だった (全論文に関して)． それはともかくとして， このあたりの理解 (観測データとモデルの対応と解釈) は現在に至ってもいぜんとして難しいまま， という気がする．
• その後， 工藤さん主催による大学内 (実験室などでの) 安全講習会． 独法化以降， こういうのが義務づけられている． 1530 ごろ終了．
• 今年度修論情勢イラク化の予感に耐えきれず対策を講じてみる． 牛原さんの屋久島の葉っぱデータはとりあえずゑくせるという牢獄から救出し， 本人にも R の使いかたはひととーり教えた． ということで次は， 現時点で「唯一」(!) 解析可能なデータセットをもつ実吉さんからデータファイルを借り受ける．
• いや， 最終的には本人に R でデータ解析してもらうわけだけど …… そもそも， のべ反復数 1000 超しかも時間構造ある発芽実験データなんで ゑくせるではぜんぜん無理． で R にもっていくためには， この呪われファイル形式からデータを「救出」せんといかんわけで， これがけっこう手間がかかる …… で， 今回のも， やはり Perl のような強力なケガれ言語なしでは， とうてい整形不可能なほどに ……
• R もなかなかすぐれた文字列処理能力をもつけれど， 生態学研究者個々人のきわめて独創的な表現能力に対応するためには， Perl のほうがむいているように思われてならない． 端的に言えば， 毒をもって毒を制する． R ではお上品すぎるのである． そして院生には R はススめるけれど Perl は (もはや) おススめではない (「とりあえずまずは R を」政策)． ということで， このあたりのケガれ作業は私の独壇場となる．
• しかしながら， こういう救出作業はなんべんやってもなかなか 「うまく問題解決しつつある」 といった気になれない． 個々のワザは改善してるんだろうけど …… なんというべきか， 万引き多発が止まらぬスーパーマーケットの警備員が 「アヤしい客を見つけるのはうまくなったんだがなぁ」 とタメイキをついてるよーな状況か?
• さいきん， 昼飯の炭水化物量が少ないせいか太陽高度とともに急速に血糖値さがる． はらへったので撤退． 1910 研究室発． 1940 帰宅． とりあえずリンゴなどかじる．
• 本日もうんどうしてないんで自宅発北大構内走にでるか …… と検討するも， なかなかふみきれず逡巡したまま晩飯の準備．
• 2210 自宅発北大構内走． 夜の北大構内はぢつはランナー密度が高い． みなさん， 人目をしのんでこっそり走るのである． また， こんな時間に集団舞踊する北大学部生男女も． これは北辺の若者たちが「ようやく春が来た」とはしゃいでるわけではなく， YOSAKOI ソーラン祭り の準備． 2255 帰宅． 体重 74.6kg． うーむ． 晩飯．
• [今日の運動]
  • 北大構内走 45 分間．
• [今日の食卓]
  • 朝 (0910): 米麦 0.5 合． コンブあえもの．
  • 昼 (1220): 研究室お茶部屋． 米麦 0.5 合． チンゲンサイ・ピーマン・ネギ・卵の炒めもの．
  • 晩 (2320): 米麦 0.6 合． ニラ卵炒飯． ダイコン・コンブ・サンマの味噌汁．

2005 年 05 月 25 日 (水)

• 0750 起床． 夜中， なんかノドがいたくて， 何度もうがいなどしていた． 朝飯． コーヒー． 0855 自宅発． 曇． 0910 研究室着．
• 北大ネット， DNS まわりがいぜんとしておかしいような気がする． しかしメイルやりとりはできてるみたいだ．
• お， Crawley 先生の ``Statistics: An Introduction using R'' (Wiley, Amazon) がすでに先月末に出版されてるな．
• ノドはかなりマシになってきたけど， 今度はアタマが少しふらふらするよーな．
• 北大ネット， やはりヘンだ．
• アカマツ， モデリングのほうの原稿なおしにとりくむ． いつもごとく， じつに匍匐前進な進捗で．
• やはり風邪ぎみなのか， 寒気がするんでジョギング用ジャージ上を着用している． 北大構内走やめ． 昼飯． すでに調査しーづんに突入してるんで， みなさんおいそがしそう．
• 昼飯のあともアカマツ原稿． ぢり．

[Gnuplot な図]

いやはや， このあたり Gnuplot 使ってたころに作ったもんで． しかも gnuplot 4.x は 3.x とは一部互換性がなくなってる， とわかった (あの劣悪な with まわり， おそらく …… 後記: 翌日のぎょーむ日誌参照)．

• やはり体調よくないので早めに撤退． 1850 研究室発． 1910 帰宅． 体温 37.0°C で予想してたほど悪くない． いきなり寝てしまう． どうやら， 昨晩「うがい」のため何度もおきたせいで寝不足状態， だったということにすぎないよーで …… われながらまぬけだな． 2200 ごろおきて晩飯． これでまた生活周期がずれるんではないか?
• いま手もとでりすとらしてるほうの アカマツ論文はかなりモデリングよりなんだが …… だったら森林衰退がどうこうとかより (これもむろん書きますけど)， むしろモデリングとしてはどういった趣旨なんだ というあたりを強調してみせるのが 「役にたたない科学」 を標榜している私には似つかわしいやうに思われる． そこでこの苦しいモデリングのすろーがんを打ちだすなら， としばらく考えて …… ``data oriented modeling'' とでもフれてまわれば， 少しはぽぢちぶなかんぢだと錯覚させることができるんでは， とゆー気になってきた． 少なくとも， 冗談としてはおもしろい …… この点はとても大切である．
• まあ， 検索すればないわけではないけど， 学術なところではあまり使われてないよーで (modeling だけでなく modelling まで含めるともう少し増える)．
  • "data oriented modeling" : 85 件 (www.google.co.jp)
  • "data oriented modeling" statistics : 16 件 (www.google.co.jp)
  • "data oriented modeling" : 4 件 (scholar.google.com)
• やはり眠くならない． そしてノドはやられぎみのまま．
• [今日の運動]
  • まぬけ休養日
• [今日の食卓]
  • 朝 (0820): 米麦 0.5 合． ダイコン・コンブ・サンマの味噌汁．
  • 昼 (1350): 研究室お茶部屋． 米麦 0.5 合． ニラ卵炒飯． ダイコン・コンブ・サンマの味噌汁．
  • 晩 (2240): 米麦 0.8 合． チンゲンサイ・ピーマン・ネギ・ショウガ・エノキダケ・豆腐・ホタテの炒めもの． ホタテは貝柱「以外」のパーツが一パック 100 円で売られていたもので …… まあ， 亜鉛の補給にはなるでしょう． ネギ・ショウガが無意識のうちに混入してるのは， 昼ごろ M1 女性大学院生から 「ネギ・ショウガを食べてないから 風邪ぎみなんかになったりするんですよ!」 と反論の余地なく断定されてしまったため， と今ごろになって気づいた．

2005 年 05 月 26 日 (木)

• 0730 起床． ノドはしだいに回復しつつある． しかしよく眠れなかったんで， アタマがふらふらする． 朝飯． コーヒー． 0900 自宅発． 晴． 0915 研究室着．
• 北大ネット， また腐れている． この「朝になったらコワれる」 って何か北大内の人間の挙動と関係ありそうな． しかし， 管理者の連中は状況を何も把握できてないようだ． あいかわらず．
• いっぽう私自身は …… アタマがややふらふらするけど， 体温は平熱ということでこれは寝不足だろう． そして今朝 0900 の札幌の気温は 12.5°C で特に寒いわけでもないはず． しかし， かとーオフィス内はサムく感じるので， 家からもってきたセーターを着こむ．
• みなさん， 調査にでかけてしまったので留守番状態．
• またアカマツ原稿の修正作業にとりくむ．
• 道北の 豊富町 (人口 5000 人弱) の警察 (駐在所: 0162-82-1703) からいきなり電話． サロベツ湿原の横の同じ場所に何日も前から車が置いてあって， こちらへの連絡先とかが書いてあるんだけど， だいじょうぶなんですか， との問いあわせ …… それは調査に行ってる院生たちのもので， と説明して， とりあえず調査隊長の西村さんに連絡してみる． とくに何か変わったことあるわけでなく， 毎日たまたま同じ場所に止めてるだけとのこと． また駐在所に連絡して問題ありませんと伝える． やっぱ， ヨソから来たのは目だつんでしょうなぁ．
• しかし警察から電話あったときは， また誰かが調査中に交通事故とかやったのか， と思ってしまった． ここらへんの院生たちは車での移動量がなかなかすごくて， ですね． たとえば， ここ札幌からサロベツまでおよそ 300Km の道のりである． これは東京-名古屋間 (東名高速) の距離とほぼ同じ．
• いかん， またアタマがへろへろしてきた． 体温上昇してるわけでもない (と思う) のにどういうわけかしらん． まあ， 単なる睡眠ぶそくなんだろう． 何度も起きてしまっただけでなく， 眠りも浅かった， と．
• 昨日の Gnuplot 問題， 今までの version では

  	plot ... with points ps 2 7
  

  とゆーふーに書いても OK だったんだけど， 4.x からは ';' expected とエラーになる． この場合は，

  	plot ... with points ps 2 pt 7
  

  としなければならぬ， とわかった． まあ， もう Gnuplot なんかめったに使わんと思うけど (R で作図するんで)．
• あまり進捗せず昼飯． あいかわらずふらふら． いまこの近辺ではやってる 結核 では， と長谷川宣教師殿からおどされる． 結核の症状 とやらをみると， 「体重減」 とあり， これは満たしてないよーな． まあ， 長びくようなら医者にいく， ということで．
• 甲山家の奥様が来てプリンター購入こんさる依頼． 当方， ススんでいなかったのですぐに仕事をほうりだして， ネット上を検索してみる． ふーむ， こんにちではスキャナつき複合機は 三万円台 なのか ……
• アカマツつづき． しかし 1620 ごろばてきったので， 中断する． 残りの時間はくだんの発芽データ救済， という Perl ケガれ仕事に着手してみるか．
• しばらくやってると， これはかなり難しい変換作業だとわかってきた． 反復の中に時間経過「○日目に○個発芽」とゆー nested な構造が含まれるため．
• さらにいろいろな自動チェックしながら変換していくフィルター (in Perl) を作って元データ (in CSV) を流しこむと …… ばんばんエラー警告をだしてくる． これはあちこちに不正規なるメモ書きがあったりするため． やっぱ， ゑくせるってぜんぜんダメぢゃん (ちょっと複雑な構造のデータを扱うには， すなわち生態学の研究に用いるには)． さらに， 入力した本人に確認せんことにはススめないんだけど …… しかしながら， 本人なぜかいないんだよね．

[カエデ]

地環研前． 開葉し花もつけてる． シュートの構造とか気になる (よくわかってないんだけど)．

• ということで， 本日もはやめに撤退． 1830 研究室発． 1900 帰宅． 体温 37.2°C よくないかんぢだ． 晩飯の準備． 晩飯．
• どうも， 体調いよいよ悪化しつつあるんでは …… とりあえず， 風邪薬など服用してみる．
• 深夜， 眠れんので体温測定． 時刻は 2440． 体温 38.0°C …… なんでクスリのんで横になってるのに， こうなるのか? たすけてー
• [今日の運動]
  • へろへろ休養日
• [今日の食卓]
  • 朝 (0750): 米麦 0.4 合． チンゲンサイ・ピーマン・ネギ・ショウガ・エノキダケ・豆腐・ホタテの炒めもの．
  • 昼 (1250): 研究室お茶部屋． 米麦 0.5 合． チンゲンサイ・ピーマン・ネギ・ショウガ・エノキダケ・豆腐・ホタテの炒めもの．
  • 晩 (2030): 米麦 0.6 合． ニラ卵． タマネギ・ニンジン・ショウガ・豆腐のカレー．

2005 年 05 月 27 日 (金)

• 0900 起床． 昨晩から本日午前ぐらいまでが最悪だったのかも． 大げさにいうならばあたかもノド全体が炎症おこしたみたいな． 一日中ねていた．
• 1000 体温測定 37.9°C
• 1315 体温測定 37.3°C …… 風邪薬のせいか， のどはすこしマシになる．
• 1850 体温測定 37.3°C
• 2240 体温測定 36.9°C
• [今日の食卓]
  • 朝 (1000): ヨーグルト．
  • 昼 (1330): シリアル．
  • 晩 (2240): 米麦 0.4 合． タマネギ・ニンジン・ショウガ・豆腐のカレー．

2005 年 05 月 28 日 (土)

• 0800 起床． 朝はなぜかノドいたい． しかし朝飯くって風邪グスリのんだらなおった． おそるべし， 風邪グスリ．
• 0810 体温測定 37.1°C
• 1100 体温測定 36.7°C
• 1330 体温測定 36.7°C …… ま， 体温に関してはほぼ平熱にもどったかな． ノドについてもかなりマシ．
• 1530 風邪グスリもつきたので， 買い出しにでかけてみる． 嗚呼， クスリ依存症． ついでに本でも， ということで JR 札幌駅へ． 10 分ほど歩くだけなんだが， 過去 40 時間以上ほとんどねていたんで， 直立歩行に違和感ある．
• 駅のクスリ屋で風邪薬 (新ルル A ゴールド, 30 錠, 924 円, 使用期限 2007.12) 買ったとこころ， テキはこちらが弱ってるのをすばやく見ぬいて， 「これもいっしょにのまれると回復が速いのですが」 とよくわからぬ 液体 小ビン三本もつけようとする． 一秒間にわたるニラみあいのすえ， ここは気あいで負けてしまって 1024 円よぶんな出費が生じてしまった． うーむ， この社会で病人とはカネを支払わされる立場なのですな …… いや， 私のよーな「不慣れな」病人が狙われるのか?
• 帰宅してまた寝る．
• 2130 ごろ起きて晩飯． かなり回復してきた． おそらく， クスリばかりのんでるせいなのか， アタマがぼーっとしている． これってきっと脳によくないんだよ． しかし晩飯後もなんとなく服用してしまう． にわか病人はヨワい．
• まあ， ともかく近ごろ札幌あたりではやってる 結核 ではなかったよーで．
• [今日の食卓]
  • 朝 (0830): タマネギ・ニンジン・ショウガ・豆腐のカレー．
  • 昼 (1340): タマネギ・ニンジン・ショウガ・豆腐のカレー．
  • 晩 (2210): スパゲッティー． バジルソース． ネギ・ニンニクのオムレツ．

2005 年 05 月 29 日 (日)

• 0810 起床． ノド・体温問題ナシ …… とゆーことでおかげさまで風邪はほぼ治ったよーだ． ただし， 三日以上も風邪グスリをのみつづけてしまったので， 体調いささか悪いよーな気がする． カゼには効くけど， こりゃ一種の毒ですよ． いまさらながら． あと， 睡眠が不規則だった影響もありそう． 朝飯． コーヒー． ヤク物摂取はやめ．
• 静養と称する怠業． 1110 自宅発． 晴． 1125 研究室着． やっぱクスリぼけしてるよーな． アタマぼーっとする．
• りはびりと称して， とりあえず木曜日にほうりなげた発芽実験データの整形作業にとりくんでみる． まあ， こういう単純ぱずるみたいなやつのほうがよさそうなんで． こないだ苦闘してた「よけーなコメント」 はそこだけ記録しとけばいいわけで ……
• しかしながら， 今度は「あるはずのデータがない」 問題につまってしまって投げ出す． なかなか難物ですよ．
• かとーオフィスのソファにひっくりかえって， 統計学本などながめる …… Stein 統計量とかって面白そうなんだけど． これのよしあしってさー， そういう variance の期待値みたいなんで評価してよいのだろうか． 不偏性とかはすでに備えている， というつもりなのか， それはまた別の問題とゆーことなのか ……
• なんかアタマの中の統計学な部分が励起されたんで， そのあたりネット上でうろうろと調べてみる． やはり， というべきか Hidden Markov Random Field (隠れマルコフ確率場) という考えかたはすでに 10 年以上まえから使われていたか …… ただし， 画像解析の分野で． そしてこれも当方の予測どーり， Markov Chain EM algorithm (MCEM 法) でパラメーター推定とかやっとるね． まさか， モデル選択のやりかたまで検討されてる， とは思わんかったけど． 画像な分野ってヘンな連中ばかりですよ．
• このへんの分野では Markov Random Field のことを (Bayesian inference の) prior (事前分布) と呼んだりすることもあるのか． たしかに， prior みたいな位置づけになるなぁ． 言われてみれば． 私は mixed model からの連想でランダム変量の分布などと呼んでた．
• そこからの連想でさらに調べると， 世の中にはすでに On Estimation and Prediction for Spatial Generalized Linear Mixed Models (Zhang, H. 2002. Biometrics 58:129-136) みたいな論文もあるわけで． これは空間構造ありの一般化線形混合モデルの推定とかやりました， というもの． 生態学って， ひょっとして世の中から取り残されてるんぢゃ ……
• この論文， ``a Monte Carlo version of the EM gradient algorithm'' を開発しましたとのことで， 私が思いつくようなことは皆さんとっくに使ってるんですなぁ． ここでいう the EM gradient algorithm ってのはうだうだと 今月はじめ に「一歩だけススめばいいぢゃん EM algorithm」 として検討してたアレと同じ …… しかも， 決定論的なほうはどうも十年ほど前から使われてたよーで (Lange 1995 とか)．
• というふうに， りはびりと称して検索したり論文ダウンロードばかりでオワりつつある． しかし上のごとくよけいな知識をアタマにいれたら， 風邪薬の毒はだいぶぬけていったかんぢだ． 2000 研究室発． 2030 帰宅． 体重 74.8kg． 病床太り． 晩飯．
• [今日の運動]
  • 北大構内走 50 分間．
  • 腹筋運動 30 × 3 回． 腕立ふせ 5 × 3 回．
• [今日の食卓]
  • 朝 (0830): ヨーグルト．
  • 昼 (1230): 研究室お茶部屋． 「北欧」サンドイッチ．
  • 晩 (2120): 米麦 0.7 合． 乾燥野菜・ナメコ・ワカメ・豆腐・シジミの味噌汁．

2005 年 05 月 30 日 (月)

• 0800 起床． 朝飯． コーヒー． 0945 自宅発． 晴． 1000 研究室着． まだ風邪グスリが体内に残存してるかんぢでアタマがぼーっとする． ノド・セキも多少の問題あり． ぢつに， 生体防御・免疫系とゆーのはびみょーなもんですなぁ ……
• アカマツ (モデリングのほう) 原稿の修正にとりくむ． 見なおしてみると， まだまだ修正が必要ですなぁ． ``data oriented modelling'' なるすろーがんは冗談としては面白いが， 前面に出すとちょっと持ちこたえられない， というのがわかってきた． うーむ．
• 昼飯．
• 東北大の河田さんみえる． 1300-1900 工藤さん部屋で モデリング相談． エゾノツガザクラ vs アオノツガザクラ， (大雪山系) 頂上決戦 ……
• 詳細を書くにはこのペイジの余白は少なすぎるので， また後日にでも． ナガレとしては
  1. 工藤さん (ツガザクラぷろぢぇくと代表) がツガザクラづの懇切丁寧なる説明
  2. 河田さん (モデル設計主任) がすごくたくさんの遺伝子を仮定した， ややこしー進化モデルを構築していく
  3. かかるモデルのプログラミングは死にそうとの 危機感つのってきた久保 (下請け業者) が「そんなのやだやだ」を正当化する口実 を真剣に検討しはじめる …… すると， 観測された現象を説明可能 (ぢゃないかもしれないけど知ったことかよ的その場しのぎ) な最単純モデルを発見
  4. とりあえず 3. の安直なやつを「ぷろとたいぷモデル」 として久保が作ってみて， それを動かしてみながら河田さんがその次の改良を考える
  で収束した． 最単純モデルとは， 工藤さんの説明にでてきた「環境傾度に沿った浸透交雑」 (あるかもよ，とゆーのは本日初めて聞いたんだが) なる概念を悪用して， エゾノツガザクラがひどく「だらしない」 & アオノツガザクラは「身持ちがかたい」 生きものだと決めつけると， うまく ``hybrid'' (とゆーことになってる連中) 形成の非対称性 (とそれが環境傾度で変化する不思議さ) が生じるでありましょう， というもの．
• 1910 研究室発． 甲山さん・高田さんも合流して 5 人でひだまり庭で晩飯． その後， ASPEN ホテルロビーで河田さんの統計こんさる …… 2250 帰宅．
• [今日の運動]
  • やみあがり休養日
• [今日の食卓]
  • 朝 (0830): 米麦 0.5 合． 乾燥野菜・ナメコ・ワカメ・豆腐・シジミの味噌汁．
  • 昼 (1210): 研究室お茶部屋． 米麦 0.5 合． 乾燥野菜・ナメコ・ワカメ・豆腐・シジミの味噌汁．
  • 晩 (1930): 北 11 西 3 の「ひだまり庭」．

2005 年 05 月 31 日 (火)

• 0730 起床． 朝飯． コーヒー． 0840 自宅発． 晴． 0855 研究室着．
• セミナーの予習． たくさん論文があってたいへんだ．
• 1300 セミナー …… のはずなんだが， 発表者たる北村君がいない． 集まったひとびとは口々に「何かまた不幸なことが」と懸念している． というのも， これまで彼が何かひと前で発表しようとするまさにその日になると， 交通事故・病院送りその他もろもろの事件が発生して中止になるからだ． とりあえあずセミナー中止の連絡をだす．
• さて本日はいかなる凶報が， と待ち受けていると …… 30 分ほどして連絡あり． 予想どーり 事故だった (もはや誰も驚かない)． 今回は人身ではなく物損． 現在は某所に抑留中． いやー， それは不幸中の幸いだったねー， と安心させるところが北村君のスゴさと言うべきなのかもしれない．
• とゆーことで脱力して昼飯． お茶部屋で亀山さん他に昨日のツガザクラ安直モデルの説明など． なんか， さっさと計算せんといかんよーな．
• で， さっさとアカマツ作業に復帰すべきなんだが …… 上のツガザクラ問題が気になって集中できない． モデルの中での抽象化された遺伝子流動はよくわかってるんだけど， その「実装」みたいなのが思いつかないため．
• しょうがないので， この問題にいったん専念してみる．
• 自殖だの花粉制限だの訪花昆虫だのめんどうな考えがアタマの中 ぐるぐるとまわる ……
• えーと， 隣花受粉の一番安直なモデルは幾何分布， か． しかしここで隣花受粉は現象にどういうカタチで顕現しているのか …… などなど面倒あれこれあって， ですね．
• これってつまり 「ある花の結実数は何で決まるか」 という単純な問いを出発点にすればいいだけぢゃん …… とゆー簡単なことに気づいたときにはすでにばてきっていた． 撤退． 2020 研究室発． 2050 帰宅． 晩飯．
• [今日の運動]
  • やみあがり休養日
• [今日の食卓]
  • 朝 (0800): 米麦 0.6 合． チンゲンサイ卵炒飯．
  • 昼 (1330): 研究室お茶部屋． 「北欧」エビチリドッグ．
  • 晩 (2140): 米麦 0.7 合． チンゲンサイ・タマネギ・ニンジン・ブナシメジのカレー．