KuboLog 2006-10-20

2006 年 10 月 20 日 (金)

0820 起床．う． 9 時間以上ねてしまった．朝飯．コーヒー． 0920 自宅発．雨． 0935 研究室着．

なぜかしら朝はアタマが三次元計算むきになっているらしく，昨日の志水さんの地衣類問題に関して，三次元ぎみなモデリングを試みてみる．とりあえず， R で作図．

[大雪山主稜線上の地衣類]

…… が好む岩面? 志水さんデータを使った作図．北緯 43° において日当たりが悪いところはダメそう．しかしながら，日当たり良すぎるところもやはり ……

本日も江川さんとサロベツデータ攻防 (いつの間に「攻防」?!) …… とはいえ，ゑくせるファイルという自由のない専制国家から生態学的データたちを救出する，というデータ解析の第一段階は終了しつつある．

あれこれ作業にひとくぎりついたので，北大構内昼飯調達の旅にでる．研究室お茶部屋にもどって昼飯．やはりおそめの昼飯くってた亀山さんいわく「ツガザクラ作文『すらんぷ』だなんてトンでもない．絶好調ですよ，絶好調 …… ってぎょーむ日誌に書いといてくださいよ」うんぬん，さらに「ぎょーむ日誌」利用による情報操作の試みあったんだけど，のこりの部分は読者の想像力で補完できるだろうから省略．

モデルのあたりをつけるか，ということで R の library(Matrix) の lmer() に志水さん大雪山地衣類データをほうりこんで，いろいろと試行錯誤してみたんだが …… ダメだ，二重・三重の random factor をくみこむと計算にかなり時間かかり (10 分ぐらい)，しかも途中で「難しくて計算できません」と投げ出してしまう．この程度の問題でナゲるなよなぁ …… まあ一般化線形混合モデル (GLMM) 計算は「すぐに難しくなってしまう」と考えるべきなんでしょうなあ．ということで，こいつも Bayesian おくり，と (なんか Bayesian が悪いことみたいだな)．

そしてサロベツデータ解析第二段階． Perl から R へ．

[まずは作図]

江川さんデータから作図 (ただし左の図はちょっと作図プログラムミスが含まれてる，と気づいた)．これの前段階 (R によるまた別のデータ変換) でちょっと苦闘してしまった．それを横で見てた江川さんが「R ってかなり難しいんじゃないか」と不審の目で見るので，それを払拭するデモンストレイション，「いやー，ここから先は簡単ですよ，簡単」などと冷汗ながしつつ宣伝文句を口ばしりつつこの作図プログラムを作った次第．もともとサロベツデータは (指導教官・露崎さんの趣味で) ややこしい実験設定のややこしい構造になっており， pdf() で全データを出力させると左のような図が 4 ペイジ続くことになる．

…… 「R って難しいんでは」との疑念いだきつつある院生を説得しつつ，ややこしい構造のデータを操作するプログラミング，というのは脳のあちこちの機能をまるちたすく的に使うらしく，かなりばててしまった． 2000 研究室発． 2015 帰宅．体重 72.4kg．まあ，これは一時的な脱水的体重減少だろう …… 冷汗ながしすぎたか? 晩飯．

いかん．夕方からの数時間で乏しい言語能力がスリきれてしまったかんぢだ．アタマ麻痺してメイルが書けん．

[今日の運動]

カラダは動かすべきだったか

[今日の食卓]

朝 (0840): 米麦 0.5 合．レトルトパウチドカレー．
昼 (1420): 研究室お茶部屋．食パン．
晩 (2035): スパゲッティー．豆腐トマトソース．ワカメスープ．

ぎょーむ日誌 2006-10-20

2006 年 10 月 20 日 (金)