「ぎょーむ日誌」目次に戻る
|
KuboWeb top に戻る
|
twilog
|
atom
ぎょーむ日誌 2011-07-01
苦情・お叱りは, たいへんお手数かけて恐縮ですが, 久保 (
kubo@ees.hokudai.ac.jp
) までお知らせください.
本日
(
kubolog20110701
) |
次の日
|
1 日前
|
7 日前
|
31 日前
|
365 日前
|
top
2011 年 07 月 01 日 (金)
0620 起床. ラジオ体操. 朝飯の準備. 朝飯. コーヒー. 0820 自宅発. 曇. 0840 研究室着.
昨日でひとまづ 「統計数理」原稿かきぎょーむが終了したわけだが, その間に放置されてたあれこれにとりくむ. Inbox にいろいろ堆積 ……
おや, 深谷君から逆 Wishart 分布についてのメモが …… うーむ,
Wishart 分布
はめんどくさい.
なンでこんなのを検討する必要があるかというと,
多変量正規分布
共役事前分布なんだよね. 以前の
ぎょーむ日誌
参照.
ちょっと文献などいろいろ見てみる.
昼飯.
事前分布たる (逆) Wishart 分布についての, 現時点での意見:
階層事前分布である多変量正規分布の事前分布 (無情報なやつ) は Wishart じゃなくてもいい んじゃないの,と考えているのでこいつの parameter 指 定問題をそんなにいっしょうけんめい考えなくてもいい かもという気もする
たとえば階層事前分布の無情報な超事前分布として, この分布を使ってよいのだろうかという点も気になっている. Gelman (2006) Prior distributions for variance parameters in hierarchical models によると,事前分布が 正規分布だとした場合,逆ガンマ分布 (逆 Wishart の一次元版) を使うのはまずいといった指摘があります.ちょっと実験してみると (きちんと実験していないのですが), たしかに標準偏差の事前分布を一様分布にしたほうがマシな結果になっているような気も ……
多変量正規分布の分散・共分散行列の無情報事前分布としては, 一様分布ってのはどうなんだろうか? 計算速度がどうなるのか,予想つかないんだけど …… 対角要素は
dunif(0, 1.0E+4)
でその他は
dunif(-1.0E+4, 1.0E+4)
で生成される対称行列 (WinBUGS の多変量正規分布関数に渡す場合はこれを逆行列にする).
いろいろかたづけて,
岩波統計モデル本
修正に復帰しようとする努力 ……
1730 研究室発. なぜかすすきので晩飯. 2040 帰宅.
[今日の運動]
朝のラジオ体操.
[今日の食卓]
朝 (0720): 米 0.5 合. タマネギ・ブナシメジ・油揚の味噌汁. タマネギさつまあげ. コマツナごまあえ.
昼 (1240): 研究室. ぱん吉のパン. 今日は弁当休養日.
晩 (1820): 南 3 西 4 の「魚菜」.
[魚菜]
つぶ貝,ゴーヤおひたし,シマアジ,カキの山椒煮,ホタテの梅煮などなど, おいしくいただきました.
本日
(
kubolog20110701
) |
次の日
|
1 日前
|
7 日前
|
31 日前
|
365 日前
|
top
KuboLog
|
KuboWeb