ぎょーむ日誌 2011-01-13

2011 年 01 月 13 日 (木)

• 0705 起床． あいかわらず雪． 朝飯の準備． 朝飯． コーヒー． 0910 自宅発． 雪． 気温零度以上になる気配のまったくない真冬日がつづく …… 歩いてる最中に視界は悪化． 50 m 先が見えない． 0940 研究室着．
• 何やら BUGS コードなメイルかきあれこれ ……
• 院生の R 往診．
• しかしこういう作業で自分のまぬけさに気づかされることもある． ifelse() 関数は vector に対応してないと思っていたのに， 対応していた． えー， これって最近の仕様変更? …… と思ったら， 2007 年の間瀬さんの 本 にわざわざ「ベクトル化条件分岐 ifelse」 とゆー小節があった …… われながらアホだった．
• ifelse() の vector 化とかこういうものです:

  > ifelse(1:9 > 5, 11:19, NA)
  [1] NA NA NA NA NA 16 17 18 19
  
  

• そして生態学会雑用はなかなか終了しない ……
• 午前中は吹雪になったり晴れたり． 昼飯．
• ぐだぐだとした大学事務雑用． このくされぎみな大学データベイス入力なシステム， ぽっぷあっぷ許可しとかないとログインできない， とゆーないすな設計で ……
• さて， 今後の予定:
  • 1/14 (金) 午後， 長時間にわたる 生態学会大会雑用会議……いーかげんに釈放してくれー
  • 1/16 (日) せんたー試験かんとく
  • 1/17 (月) M1 あいてのじゅぎょー
  • 1/19 (水) M1 あいてのじゅぎょー
  • 1/25 (火) 統数研公開講座のテキスト提出
  うーむ ……
• 来週の M1 じゅぎょーについてうだうだとした検討 ……
• セミナーの予習．
• Trendy セミナー， 本日は小野田さんで葉の力学的性質のメタ解析． オーストラリアの プロジェクト に やとわれた 小野田さんたちが， かなーり苦労されつつあちこちの葉論文の著者たちなどからデータをだしてもらって， (葉の破壊への強度) = (∝か?) (葉の厚さ) × (材密度) × (質的強度) とゆーモデルで Y = log(葉の強度) として Var(Y) = Cov(Y, log(厚さ)) + Cov(Y, log(密度)) + Cov(Y, log(質的強度)) なる「分散を共分散に分解」なることをやって， 厚さ・密度は計れるから残りは質的強度で …… とゆーハナシをきいて， ふむふむと思ったのだけど …… あとからいろいろ考えてみるとやはり少しヘンかもとゆー気もしてきた．
• まず， (葉の厚さ) × (材密度) × (質的強度) とゆーモデルが成立してなくて， たとえば一番簡単には， 葉の厚さ (mm) なんかのベキ乗といった アロメトリックなモデルだとすると， covariance 分割による「Var(Y) への貢献」みつもりは おかしいことになりそうな気がする． 厚さの a 乗だとすると a Cov(Y, log(厚さ)) となるからね ……
• そもそも分散を共分散で分割するって， だいじょうぶなのだろうか． 寄与率なる割合を Cov(Y, log(X)) / Var(Y) と計算するとのことなんだけど， これは Var(log(X)) × 相関係数(Y, log(X)) だ (後記: これはまちがい …… もっとわけのわからない量になる)． てことは， Var(log(X)) が大きくても相関係数(Y, log(X)) が大きくてもおなじ?
• まあ， 小野田さんは明日も研究室におられるようだから， 時間があればお尋ねしてみよう．
• 1830 研究室発． すごく寒い． 買いもの． 1910 帰宅． 晩飯の準備． 晩飯．
• ……
• [今日の運動]
  • 朝のラジオ体操． うんどう休養日， つづく ……
• [今日の食卓]
  • 朝 (0730): 米 0.5 合． ダイコン・ブナシメジ・豚肉の味噌汁． 生ハム．
  • 昼 (1200): 研究室． 食パン． リンゴ．
  • 晩 (2030): 米 0.5 合． キャベツ・鶏だんごのトマトスープ． マッシュサツマイモ・ショウガ・バター． チンゲンサイ・シイタケ・イカの炒めもの．