ぎょーむ日誌 2003-03-24
2003 年 03 月 24 日 (月)
- 0730 起床.
朝飯.
コーヒー.
0830 自宅発.
晴.
0840 研究室着.
- だらだらと
……
ぎょーむ日誌書きなど.
何をやってんだか.
まあ,
りはびりということで.
- 冨沢さんが霧多布から来ておられるので,
アカエゾマツ - ハンノキ議論など.
- 昼飯食ってからも
自由集会ペイジ
の補強作業などのんびりやってたら,
「2 本目の修論」
を書いてる M2 井田君からマルハナバチ統計相談.
- で,
またまた分散分析 (under 等分散・正規分布) かよぉぉお
……
とはいえ,
内容は面白そうなので Perl と R を駆使して取り組んでみる.
- 割算・クラスわけという「生態学的情報隠蔽」を重ねた結果から遡及して,
原データに近づいていくにつれ,
徐々にいろいろなことが明らかになっていく.
- まずはハチのとびかたによる花粉移動評価のシミュレイションだが
……
ちょっとでも複雑な計算をやりたい場合は,
ゑくせる使うな
ということ.
一個体一個体の数字をよくよくアタマのなかで計算してみると,
どうもつじつまのあわぬところが出てきた.
呪われスプレッドシートにたちかえって調べてみると,
何列にもわたって面倒な式を指定してるからだよ
……
ということで,
けっきょくは Perl でシミュレイションを作りなおしてみた.
こちらのほうがよほど簡単で見とおしがよく応用がきく.
- その間にも三中さんからメイルいただき,
ばてによって欠落した私の記憶を補完していただいた.
ありがとうございました.
- Perl ハチを Mac OS X 上でも動作するようにする
……
井田君の G4 PowerBook にインストールされてる
Apple 社謹製の XDarwin とやら,
ちと腐れていないか?
xterm の挙動とかがヘンだ.
なぜか readline ライブラリとの連携がうまくいってない.
死ぬほど使いにくい.
- R によるシミュレイション結果の解析と作図
……
ほーら,
割算で消されてしまった構造があぶりだされてきたでしょうが.
- あいかわらず R 初心者な私はしょっちゅう行き詰まったりして,
ISwR
とか
Crawley 本
とか調べまわる
……
それでもわからないと,
最終兵器かとー先生にお願いするわけで.
なるほど,
dataframe から特定の名義変数をふくむ行を選択するのは
data[data.type %in% "A",]
(data
内の type
という列が
"A"
である行だけ抽出する)
と書けばよろしいわけですか.
%in%
なんて二項演算子はどこにも書いて
……
あ,
R で help(match)
すると解説と使用例がでてきた.
- で,
最終的な作図も当然ながら R で片づけたいところだが
……
今週末には就職のため札幌を離れる井田君にこれ以上の
試行錯誤に時間をとらせるわけにもいかないようで
(しかもおひすと postscript 系の図は相性が良くない)
……
計算結果をテキスト出力してあとは好きなように,
ということで.
ええ,
ゑくせる作図でも何でも気にせず片づけちゃってくださいよ.
- りはびりのつもりがいきなり修論締切まぎわ的なちょー高速解析で終わった
1 日であった.
ばて.
2120 研究室発.
2130 帰宅.
体重 71.4kg.
つくばやせ + マルハナバチやせ,
というところか.
1kg ぐらい落ちてるよ.
- [今日の運動]
-
一種のやみあがり状態,
ということで運動せづ
……
しかしぼちぼちとまた体力の涵養に
つとめなくては
- [今日の食卓]
- 朝 (0800):
米麦 0.7 合.
チンゲンサイ・タマネギ・マイタケの炒めもの.
- 昼 (1240):
弁当.
研究室お茶部屋.
米麦 0.7 合.
チンゲンサイ・タマネギ・マイタケの炒めもの.
- 晩 (2240):
米麦 0.7 合.
ジャガイモ・タマネギ・ブナシメジのシチュー.