KuboLog 2006-05-11

苦情・お叱りは，たいへんお手数かけて恐縮ですが，久保 (kubo@ees.hokudai.ac.jp) までお知らせください．

2006 年 05 月 11 日 (木)

0730 起床．朝飯．コーヒー． 0915 自宅発．雨．よーやく地環研前のエゾヤマザクラ開花．しかし，あとひと月ちょいで夏至なんですが． 0930 研究室着．

あらら，またしても web server の HDD に関する警告が (なぜかしら今朝だけはでていないけど) ……

May 10 04:08:45 hosho kernel: hda: dma_intr: status=0x51 { DriveReady SeekComplete Error }
May 10 04:08:45 hosho kernel: hda: dma_intr: error=0x40 { UncorrectableError }, LBAsect=71358237, high=4, low=4249373, sector=27532912
May 10 04:08:45 hosho kernel: end_request: I/O error, dev 03:02 (hda), sector 27532912

/dev/hda そろそろ買いかえどきかなあ．ひっこしが面倒だ．

苫小牧直径成長論文原稿，いつまでも放置するのはまずい，と取りくみはじめる．

あまり進捗ないまま昼飯．日浦さんがきておられたので高所作業保険の分担金 (この研究室ぶん) をおわたししておく．

私はときどき「自分で考えたみょーなコトバ」でネット検索して遊んでみることあるんだが …… まさか「重力波キック」なるコトバが見つかるとは．そしてどういう原理で双子ブラックホール (が一体になったもの?) が銀河中核部から蹴りだされるのかわからん．

さらにこの「重力波キック」は日本どくとくのコトバか? というのも原文みてみても gravitational wave kicks なる語はない．そもそもそういう語は検索しても見つからない …… が， gravitational kicks というコトバはあるみたいだな．「重力きっく」? 別モノか? 双曲線軌道 (hyperbolic orbit) とかとは関係なさそうだし．謎．

昨日「なぜかしら教科書はいつまでたっても出版されないけど」と書いてしまったせいなのか，ずっと前に Amazon に注文してた James S. Clark 教科書のうち最初の一冊，現在もっともらでぃかるな生態学データ解析本，

Clark JS and Gelfand A (eds). 2006. Hierarchical Modelling for the Environmental Sciences Statistical methods and applications. Oxford University Press.

が届いた．

[環境科学のための階層モデリング]

詳しい内容は出版社のペイジに．ペイジ数は 205 と多くはないが (厚さも 1cm 弱)，内容はひじょうに「アツい」ものがある．これまでにない本・時代の変わり目，ということで「画期的」なのかも．

とにかく最初から最後までひたすら Bayesian 推定 & MCMC 計算な生態学データ解析．最初のほうに入門的な内容． Gibbs sampler の導出例だの， WinBUGS のコード例だの．なんと MCMC++ なる MCMC 計算用の C++ クラスライブラリもあるのか．この本の contributor の一人， Kent Holsinger 作．

あとは Bayesian データ解析の研究紹介．この本の contributor 数は 30 人ぐらい，なんだけど．全体四部構成で

Introduction to hierarchical modelling (階層モデリング入門)

ベイズ推定，グラフィカルモデリング， Gibbs sampler，計算ソフトウェアに関する入門的な内容

Hierarchical models in experimental settings (実験状況下での階層モデルたち)

FACE (Free-Air CO₂ Enrichment -- 野外二酸化炭素高濃度化) 実験状況下での植物の成長データ解析

Spatial modelling (空間モデリング)

空間相関を考慮した生物の「いる・いない」問題 (分布推定の問題)

Spatio-tempo modelling (時空間モデリング)

前の部に続いて，今度は生物分布の時間変動の問題を「隠れた状態」つかった階層ベイズモデルであつかう

いやはや，スゴい本ですね．

1500 ぐらいから夜にかけて，塩寺さんと熱帯林葉っぱの生存時間解析． library(survival) にある Cox 比例ハザードモデルがうまく適用できてるみたいだな …… 推定されたパラメーターのもとで，任意の生残曲線を出力するには survfit() にあたえる data を工夫する．

そして統計モデリングに関するメイル議論．空間モデリングで使う Gaussian Random Field (GRF) について．このあたり，アタマの中で GRF が動いてくれてるような気がする，つまり明確にわかっているようでいながら，ぢつは GRF に関する重要な性質がよくわかっていない．とくに GRF つかったベイズ推定における事後分布はいったい何なんだろうか，というのはうまく説明できないような．いや，私のアタマの中での図示のしかたが悪いだけなのかな．

2115 研究室発． 2130 帰宅．晩飯．

[今日の運動]

意外にも，というか今週はうんどう不足がつづく …… なンかだらだらした日常なのかしらん?

[今日の食卓]

朝 (0820): 米麦 0.7 合．サバ塩ゆで．ニラ・海藻・豆腐のスープ．
昼 (1310): 研究室お茶部屋．米麦 0.5 合．サバ塩ゆで．ニラ・海藻・豆腐のスープ．
晩 (2230): 米麦 0.7 合．卵やき．ニラ・海藻・豆腐のスープ．

ぎょーむ日誌 2006-05-11

2006 年 05 月 11 日 (木)