ぎょーむ日誌 2006-05-08
2006 年 05 月 08 日 (月)
-
0815 起床.
朝飯.
コーヒー.
0945 自宅発.
晴.
よーやく冬は終わりつつある,
という気温.
1000 研究室着.
-
論文原稿よみ.
隣の講義室から東さんの独特な語り口調が聞こえてくる.
-
年輪解析論文についてちょっと調べてみる
……
と今までまったく気づかなかったんだけど
カルマンフィルター
(Kalman filter)
応用した研究は意外と多い.
私はいままで年輪解析まわりはアヤしいヒトたちばかり,
と思っていたんだけど,
内在する欠点を自ら改善していくという点において
しごくまっとうな学問分野かも,
と認識をあらためぎみ
(しかし例によって本邦のひとびとは
世の中から取り残されぎみという気もするけど).
-
たとえば
van den Brakel & Visser (1996)
ってもう 10 年前にでた論文だけど,
従来方式について私がまさにダメだと思ってた点なんかも
きちんと指摘されていて,
The first disadvantage of this approach is that response functions
are time-independent.
Therefore, it is not possible to detect slowly changing reactons to
climate variables.
For instance,
due to a slowly changing groundwater level or industrial air pollution
it is possible for trees to react with a slowly increasing water need,
which could be detected by a slowly increasing response function for
monthly precipitation variables.
A second disadvantage is that trend and climate variables are not
modeled simultaneously with the tree-ring chronology.
Therefore, the model cannot ``find'' the optimal partition of the
systematic variation over trend and climate variables.
これを解決するために状態空間モデル
(state space model)
を使い,
観測データからカルマンフィルターのアルゴリズム
つかってパラメーター推定しました,
というしごくまっとうなもの.
-
私は
「年輪解析のほうではだれもこういうモデリングやっていまい」
と憶測して自分でそういう方法とか定式化せんといかんのかねえ
とか考えてたんだけど,
じつはとっくの昔にこのあたりの問題は解決していたのである.
-
このあたりの問題,
なぜかしら 1980 年代末-1990 年代初頭の
Biometrics
にいくつか取っかかりの論文が発表されている
……
が,
JSTOR 内の古い Biometrics 論文はここからは見れないので,
そのうち機会をみて
学内の図書館
にでも取りにいきますかね.
-
そしていつものことながら
R
おそるべしで CRAN には
sspir
package なるカルマンフィルターを glm() 的に使える package
なんかも登録されている.
-
プリンター問題.
現在このへんで使えるカラー印刷できるレイザープリンターは
7F の EPSON LP-S6500 だけになってしまった.
A802 の CANON LBP2810 はモノクロ印刷なら問題ないけど,
カラーはだめ.
近日中にカートリッジ交換で修繕する予定.
-
遅い昼飯くおうとすると,
低温研での授業おわった院生たちがどやどやと
……
テイブルのすみで食事しつつ皆さんの活発なる会話を拝聴しているうちに
思いついた conjecture:
お茶部屋雑談のハナシの無秩序さ・乱暴さはその場で形成される話者グラフ
のエッジ数 (いや exp(エッジ数)?) に比例する.
-
北大ネットの spam filtering,
そろそろ能力的な限界に到達しつつあるようで,
近ごろまた spam が増えつつある
(私のところに来たのは
Gmail
その他で始末されるわけだが).
いまのところ,
北大を通りぬける spam 数は毎日 10-20 通ぐらい,
というところか.
-
いかん,
お勉強だけで一日が終わってしまった
……
1940 研究室発.
2010 帰宅.
晩飯.
-
[今日の運動]
-
[今日の食卓]
- 朝 (0900):
米麦 0.6 合.
ハクサイ・エリンギ・豆腐・サケのシチュー.
- 昼 (1420):
研究室お茶部屋.
米麦 0.5 合.
ハクサイ・エリンギ・豆腐・サケのシチュー.
- 晩 (2040):
米麦 0.7 合.
ハクサイ・エリンギ・豆腐・サケのシチュー.