ぎょーむ日誌 2008-09-(21-30)
2008 年 09 月 21 日 (日)
-
0910 起床．
生活周期，
ぢりぢりと復旧中．
朝飯．
コーヒー．
-
スキャナーがとどいてしまった．
[CANON LiDE 30]
6 年前発売の
古い製品．
Yahoo オークションで 1500 円だった (箱その他なし，USB ケイブルのみ)．
ヘッドが動かない
……
とちょっとアセったんだけど，
本体下にロックスイッチで固定されてたことに気づいてひと安心．
-
Linux (Vine Linux 4.2) 上でのスキャナー設定．
SANE
のたぐいの package を
apt-cache serach sane
とかで調べてから，
sudo apt-get update && sudo apt-get install あれこれ
(sane や sane-clients)
でインストール．
CANON LiDE 30 の driver もインストールされる．
-
設定は
いろいろ検索した結果
からてきとーに．
一番原始的には，
-
sudo /usr/bin/sane-find-scanner
とすると，結果出力に
found USB scanner (vendor=0x04a9 [Canon], product=0x220e [CanoScan],
chip=LM9832/3) at libusb:002:005
といった libusb 情報が表示される
(ここでは 002:005; これは usb ケイブル接続しなおすごとに変わる)
-
sudo chmod 666 /proc/bus/usb/002/005
-
sudo chmod 666 /dev/bus/usb/002/005
-
これで使えるようになる；
一番おてがるには (Photoshop もどきの free software である)
GIMP
のファイル
→
取りこみ
→
xscanimage
→
Device dialog ...
1.-3. の作業を自動化するスクリプトを書かんといかんね．
下の画像は 50 dpi でとりこんでみた例．
JPEG で保存したのでにじんでいる．
-
昼飯．
1350 自宅発．
晴．
JR 札幌駅周辺うろうろ買いもの．
駅の薬局で手あれ対策クリーム，
ついでに鼻炎のクスリ．
駅のすぐ南の Loft で指くつ下とか．
Loft の紀伊国屋書店うろうろしてみたけど，
けっきょく何も買わず．
駅のすぐ西のヨドバシで，
研究室用の紙 (1000 枚 920 円) とか CD-R (50 枚 1580 円) といった消耗品とか．
これは R 作図わざ本チェックで大量の紙を使うので，
まあ (いつかはもらえるハズの) 印税の一部で購入してみた．
1530 研究室着．
-
R の作図わざ本のやく文ちぇっく．
-
今日の R の作図わざ本
ほんやく文ファイルみなおし作業:
全 285 ペイジ中の 122 ペイジまで．
「ですまーち」第 3 章終了．
-
1900 研究室発．
1920 帰宅．
晩飯の準備．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0930):
シリアル．
- 昼 (1300):
米麦 0.8 合．
ネギ・ダイコン・ショウガ・ホッケの煮もの．
- 晩 (1940):
米麦 0.8 合．
ネギ・ダイコン・ショウガ・ホッケの煮もの．
カボチャの煮物．
2008 年 09 月 22 日 (月)
-
0700 起床．
朝飯．
コーヒー．
0835 自宅発．
晴．
0850 研究室着．
-
R 作図わざ本，
ほんやく文ファイル上のみなおし作業．
-
とりあえず 1.5 時間ほどやって，
library(lattice)
な第 4 章終了．
ひと休み．
-
今日の R の作図わざ本
ほんやく文ファイルみなおし作業:
全 285 ペイジ中の 148 ペイジまで．
-
ぜんぜんススんでいない，
car.normal() 解説
かきのつづき．
とりあえず，
生態学会誌のイヤな投稿規定に
対応させる図ファイル管理のしくみをヒネくって，
うだうだしてみる．
-
今さらながら，
\begin{figure*} ... \end{figure*}
を使えば二段組でも「ぶちぬき」な図が作れる，
とわかってしまった
……
まあ，
これでしくみが単純化されたからいいや．
-
夏休みなるものを申請しわすれていたので，
9/24 (水) - 9/26 (金)
を「なつやすみ」として申請してみる．
まあ，
どこかに遊びにいけるわけでもないんだけど．
-
昼飯．
-
三ヶ月前
に投稿したアリ論文原稿について，
まったく無反応状態がつづいているので，
いったいどうなっているんでしょうかねえ
メイルを編集部に送ってみる．
-
解説記事作文をぢりぢりと進捗させて，
1910 研究室発．
1925 帰宅．
晩飯の準備．
晩飯．
-
アリ論文に associate editor からメイルいただく．
なかなか名文．
やはり editor たるもの，
これぐらいの作文能力がないといかんのだろうなぁ
……
thank you for your inquiry. Your Ms was assigned to an Associate editor June 25. Eventually, 2 referees agreed to review (it seems, it was not easy to find agreeing referees..maybe because July/Aug this year have been full of conferences and I had the impression most of the people were traveling, traveling, traveling...). One referee submitted a review, the other has been reminded now 3 times. Please ask again this Friday.
-
[今日の運動]
-
[今日の食卓]
- 朝 (0810):
米麦 0.8 合．
ナス・ネギ・ピーマン・エノキダケ・卵の炒めもの．
- 昼 (1355):
研究室お茶部屋．
米麦 0.8 合．
ナス・ネギ・ピーマン・エノキダケ・卵の炒めもの．
リンゴ．
- 晩 (2130):
米麦 0.8 合．
サケのムニエル．
カボチャの煮物．
トマト．
2008 年 09 月 23 日 (火)
-
0800 起床．
朝飯．
コーヒー．
洗濯．
怠業．
1040 自宅発．
曇．
1055 研究室着．
-
R 作図わざ本，
ほんやく文ファイル上のみなおし作業．
-
とりあえず第 5 章終了
……
と思ったら第 6 章だった，
これは．
どうりで短いはずだ．
気をとりなおして，
第 5 章にもとりくんでみる
……
-
第 5 章もどーにかこーにか，
ひととーり読んだ．
昼飯．
人が少なくて平和だ．
-
今日の R の作図わざ本
ほんやく文ファイルみなおし作業:
全 285 ペイジ中の 220 ペイジまで．
-
car.normal() 解説
かきのつづき
……
-
なかなかススまず，
1845 研究室発．
1900 帰宅．
ちょっとうんどう．
晩飯の準備．
晩飯．
-
[今日の運動]
-
腹筋運動 30 ×
3 回．
背筋運動 30 ×
3 回．
腕立ふせ 10 ×
3 回．
スクワット 100 回．
ストレッチ．
-
[今日の食卓]
- 朝 (0910):
米麦 0.8 合．
ナス・ネギ・シイタケ・卵の炒めもの．
キュウリ．
- 昼 (1320):
研究室お茶部屋．
米麦 0.8 合．
ナス・ネギ・シイタケ・卵の炒めもの．
- 晩 (2040):
米麦 0.8 合．
ナス・ネギ・シイタケ・卵の炒めもの．
2008 年 09 月 24 日 (水)
-
0800 起床．
書類上の「夏休み」第 1 日目．
朝飯．
コーヒー．
ちょっと怠業．
1010 自宅発．
晴．
1025 研究室着．
-
R 作図わざ本，
ほんやく文ファイル上のみなおし作業．
-
なんとかひととーりみおわった
……
-
今日の R の作図わざ本
ほんやく文ファイルみなおし作業:
全 285 ペイジ中の 285 ペイジまで．
この作業には 8 日間を費やしてしまったか
……
- 2008-09-17 前書き + 5 ペイジまで
- 2008-09-18 42 ペイジまで
- 2008-09-19 76 ペイジまで
- 2008-09-20 91 ペイジまで
- 2008-09-21 122 ペイジまで
- 2008-09-22 148 ペイジまで
- 2008-09-23 220 ペイジまで
- 2008-09-24 285 ペイジまで
-
原稿チェックばいと院生募集のメイルをだしてから，
昼飯．
-
なんと，
地環研中庭の樹木が数日前の大風で倒壊していた
(しかもいまごろ気づいた)．
けっこう大きな樹木で，
てっぺんは A 棟 8F にとどくぐらいだった．
高さ 15 m ぐらい，直径 1 m ぐらい，
か?
いやー，
風がふいていた日 (土曜日?) に「めきめき，ばりばり」
とちょっと遠くから何かが壊れる音がしたので窓からまわりを
見たんだけど何が倒れたのかわからなかった．
下の写真は 8F の非常階段から見おろした構図
(9/26 にも写真あり)．
-
ちょっとネット雑用とか．
いろいろ細かい設定で時間とられる．
-
うーむ，
Advanced R Tutrial
ににぎやかしとして参加して，
講師 (これはじつにそうそうたる貌ぶれだ)
の皆さんのお世話係を，
といった依頼が
……
かなりまよったんだけど，
けっきょくこれは翌日になって引きうけることに．
しかし，
このまにあっくな内容に 1.5 万円を支はらう参加者って，
どれぐらいいるものなのかしらん?
東京圏だとそういうヒトがぞろぞろ集まるのかな?
あるいは IASC2008 にくる海外の大学院生
……
にはこの内容は入門的すぎるような気もするな．
無料 workshop ならともかく．
-
不本意ながら server 機おきばにされているこの部屋の
無停電装置が一個こわれた．
これってけっこう消耗品的側面があるな
……
いろいろと対処．
-
今日は解説文かきはぜんぜんススまなかった．
しかも〆切が一ヶ月のびたのでちょっと脱力ぎみ
……
しかし次のぎょーむもいろいろあるので，
今月中にはひととーり書いてしまって，
皆さんのご意見をうかがうのがよいだろう．
1855 研究室発．
買いもの．
1915 帰宅．
ぐったり．
晩飯の準備．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0850):
米麦 0.8 合．
豆腐・ワカメの味噌汁．
ネギ・シイタケ・卵の炒めもの．
- 昼 (1420):
研究室お茶部屋．
米麦 0.8 合．
カップ野菜スープ．
- 晩 (2150):
米麦 0.8 合．
キャベツ・タマネギ・ブナピー・豆腐・サケの味噌汁．
2008 年 09 月 25 日 (木)
-
0750 起床．
書類上の「夏休み」第 2 日目．
朝飯．
コーヒー．
0930 自宅発．
曇．
今朝の気温はかなり下がった．
札幌で 11 度ぐらい．
0940 研究室着．
-
何か夏休みらしいことを，
ということで歯医者の予約．
-
昨日のつづき．
横浜で開催される
国際計算機統計学学会 (IASC) 2008 年大会
(私は参加しない)
のオマケである
R Advanced Tutorial
について，
いろいろ説明のメイルいただく．
今回の (とある) がらみの invited session の講演者
(みなさん R 開発の重要メンバー)
の旅費補助のための tutorial うんぬんということらしいので
……
まあ，
R 開発者のみなさまへのカンパというかんぢで，
参加もうしこみ
(支出: 参加費 1.5 万円 + 自分の旅費の一部)
してみる．
ただし夜の懇親会の幹事なる任務つき
……
-
R-help でちょっと面白い表現のやりとり．
こういう「ありがちなお気楽」質問
I would like to know how to calculate the F-value and P-value for a
non-linear regression: y ~ a * X^b * W^c ?
(すごーく非線形なモデルの F 値や P 値ってどうやって計算すんの?)
に対して，
library(lme4)
だの
library(Matrix)
だのの開発者である Douglas Bates
の回答のでだし部分:
This answer may sound pedantic but it is a result of my having spent a
considerable part of my life thinking about nonlinear regression
models. You are assuming that it is meaningful to associate an
F-value and a p-value with a nonlinear regression model and it is not
quite that easy.
(久保訳) ここに書く回答は pedantic
(君が求めてるような実用的なものというより杓子定規な学術的なもの)
なんだけど，
非線形回帰モデルについてぼくの人生の少なからず時間を費やして考えた
ひとつの結果なんだよね．
そもそも
君は非線形回帰モデルと F 値だの P 値だの関連づけることに
何か意味があると思ってるみたいだけど，
それってそんなに簡単なことではないんだよ．
……
で，線形モデルとはちがって非線形モデルでは
(これで OK と言える) 帰無仮説を設定するのがしんどいので，
検定のたぐいはラクではないよ，
といったハナシ．
全文は r-help archive の
ここ
で読める．
-
たしかに，
お気楽質問にたいして pedantic な回答をせざるをえない状況はよくある．
重要なのは，
こういった回答に対して質問者がどう反応するか，
その観察ですな．
-
また解説記事かきのつづき．
-
夏休みらしさの演出にならぬものか，
と自宅に帰ってのんびりと昼飯．
雨．
-
午後も解説かき．
進捗速度がたいへんおそい．
最初の予定とはどんどんずれた内容になってる．
-
1700 から 1.5 時間ばかり，
研究院セミナーにでる．
国立環境研究所の
山形与志樹さん
で IPCC 森林ハナシと国環研の陸域メッシュ土地利用モデルのハナシ．
やはりこのあたりの分野のみなさん，
国際〆切におわれて「それいけー」てなかんじでモデリングしてるような気も．
私ていどの作文能力では幸か不幸か無理だなぁ
……
-
1950 研究室発．
雨．
2005 帰宅．
晩飯の準備．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0840):
米麦 0.8 合．
豆腐．
キャベツ・ネギ・ミズナ・ハムのサラダ．
- 昼 (1330):
研究室お茶部屋．
米麦 0.8 合．
キャベツ・ネギ・ミズナ・ハムのサラダ．
- 晩 (2130):
米麦 0.6 合．
キャベツ・ネギ・ミズナ・ハムのサラダ．
2008 年 09 月 26 日 (金)
-
0750 起床．
書類上の「夏休み」第 3 日目，最終日．
朝飯．
コーヒー．
-
歯医者．
レントゲンとられたり，
といった検査がメインで．
来週はその検査結果についての相談のため，
来院しなければならず．
-
1140 自宅発．
晴．
1155 研究室着．
(先日 も紹介した)
地環研中庭の倒木 (たぶんハルニレ)
の写真でもとってみる．
けっこうでかい．
直径 1.4 m ぐらい?
あとで聞いたら，
院生の宮田さんがさっそく「樹木てっぺんサンプリング」
にとりかかるそうで
……
-
昼飯．
-
car.normal()
解説かきのつづき
……
-
こんな図を追加してみた．
Gaussian Random Field の理解の助けにならんだろうか?
-
夕方，
また歯医者．
-
1740 帰宅．
のんびり怠業．
雨がふってきたので，
室内でうんどう．
体重 70.6 kg．
お，
ひさしぶりに 71 kg 以下になった．
晩飯．
-
明日はなぜかしら甲山さんたちと
羊蹄山 (京極コース)
を登ることに．
天気とかどうなることやら．
-
[今日の運動]
-
エアロバイク 50 分間．
400 kcal．
ストレッチ．
-
[今日の食卓]
- 朝 (0830):
米麦 0.4 合．
キャベツ・ネギ・ミズナ・ハムのサラダ．
- 昼 (1330):
研究室お茶部屋．
ぱん吉のフランスパン．
キャベツ・ネギ・ミズナ・ハムのサラダ．
- 晩 (2030):
キャベツ・ナス・ネギ・煮干の味噌汁．
ぱん吉のフランスパン．
2008 年 09 月 27 日 (土)
-
0455 起床．
今日は
羊蹄山
どうなることやら登山．
朝飯．
コーヒー．
土日に早おきすると NHK で中島裕子さんの声が聞ける．
0545 自宅発．
0620 甲山さんの車にのって北大前発．
本日のメンバーは甲山さん・庄山さん・宮田さん・私．
札幌の天気は曇ときどき晴．
[08:07 羊蹄山麓]
札幌から 77 km，
中山峠をこえて到着した京極町．
天気はあいかわらず雨がふったりやんだり．
正面は羊蹄山 (1893 m)．
けっきょく，
本日は山頂が一度も見えなかった．
9 月なのに気圧配置は冬型．
風が強く，
雲がどんどん流れる．
不順な天気のせいか，
今日はあちこちで何度も虹をみた．
[08:15 登山口着]
畑の横にある登山口．
車を止められる (しかし止めてよいのかどうかよくわからぬ)
場所があったので，
ここに車をおく．
入山届．
札幌からこれほど近い有名山なのに，
本日の入山者なし．
悪天のせいだろう．
[08:27 出発]
天気よろしくない中，
出発．
まずは畑の横の道を歩く．
山はあいかわらず雲の中．
[08:35 登山道始点]
ここからトドマツ植林地の中の道がしばらくつづく．
標高はおよそ 400 m．
本日の登る予定の標高差はおよそ 1500 m．
これはかなりのもので，
4 時間ちかくをついやすかもしれない．
[09:03 二合目]
すでに人工林はぬけて天然林 (二次林?) の中．
天気はあいかわらず晴れたり雨ふったり，
と安定しない．
[09:39 四合目]
このあたりからひたすら雨．
やや体調不調の宮田さん，
10:10 に車のカギをもって撤退．
[10:44 降雪]
五合目 (標高 1000 m ぐらい)
をこえたあたりから雨は雪になった．
9 月でこんなに雪がふるとは
……!
[10:48 本日の最高到達点]
雪はますますはげしく，
寒さ対策の装備 (厚めの手袋とか)
もってなかったわれわれは撤退を決定．
本日の最高到達点 (1200 m ぐらい)
で休憩．
後ろからかなり速いペイスで登ってきた 60 才
ぐらいのかたも「あんたたちが降りるんなら私もやめときますわ」
と撤退．
[12:27 下山]
11:00 に下山開始で 1.5 時間ほど費して下におりた．
標高の高いところは雪，
低いところは雨の中．
先におりてた宮田さんは車の中にいた．
[12:45 京極温泉]
ふきだし公園のちかくにある
京極温泉
へ．
入浴料 500 円 (ただし来月から 600 円)．
天気が良ければ露天風呂から羊蹄山が見えるはずなんだけど，
あいかわらず雲の中．
[13:40 ふきだし公園]
温泉に入ってから，
水の噴出口のある
ふきだし公園
へ移動．
名水の里，
ということで雨なのにそれなりの人数の観光客たちが水くみやってた．
[14:15 そば屋発]
新そばの季節ということで，
これまたふきだし公園の近くにある
一心庵
でおいしい蕎麦をいただく．
-
また甲山さんの車で送っていただいて，
1620 帰宅．
洗濯．
怠業．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0520):
キャベツ・ナス・ネギ・煮干の味噌汁．
ぱん吉のフランスパン．
- 昼 (1355):
京極町の蕎麦屋
一心庵
の特もり．
甲山さんにごちそうになってしまいました．
- 晩 (2200):
米麦 0.8 合．
ショウガ味のマッシュサツマイモ．
タマネギ・ブナシメジ・豚肉の炒めもの．
2008 年 09 月 28 日 (日)
-
0720 起床．
朝飯．
コーヒー．
洗濯．
怠業．
昨日の羊蹄山登山で右足小指がちょっとクツずれして，
あまりうんどうする気になれづ
……
-
とはいえ，
一日中ごろごろしてるのもナンなので，
ふと思いたって札幌市というか市の中央を東西に走る
「大通り」(公園のあるとこ) の西端にある円山にふらふらと散歩に
言ってみる．
標高 200 m ちょいなんだけど，
地下鉄だけでいける数すくない山のひとつなんで．
-
しかしぼんやりしてたので，
カメラも時計も忘れた．
地図はもっていたけど．
-
円山はゆっくりのぼっても 30 分ぐらいで山頂に到着できる．
登山口ふきんはやたらと地蔵がならんでるややぶきみなところで，
ここからやや細い尾根を南にすすむ．
尾根が細いせいか風倒木多数．
いったん少しくだってコルをすぎると，
そこからは円山の「まるい」部分への登りとなり，
山頂．
たしかに，
こりゃー展望がいいわ．
札幌市南端の札幌ドームと北端の (丘珠空港わきの) ツドームが
同時にみえる．
-
てな，
ぐだぐだした一日，
でした．
ぐだぐだ．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0800):
食パン．
モヤシ・ベイコンの炒めもの．
- 昼 (1230):
地下鉄円山公園駅すぐ近くの
「おばんざい とらや」
の昼定食．
- 晩 (2010):
米麦 0.8 合．
ショウガ味のマッシュサツマイモ．
タマネギ・ブナシメジ・豚肉の炒めもの．
レタス・ブナシメジの炒めもの．
2008 年 09 月 29 日 (月)
-
0710 起床．
朝飯．
コーヒー．
0830 自宅発．
晴．
0845 研究室着．
-
car.normal()
解説かきのつづき
……
全体のしめきりがのびたせいか，
当方の作業もおくれぎみ．
いやはや．
さっさと終わらせなくては．
-
ちょっとネット雑用．
また生物圏専攻の掲示板 (tdiary で運用)
の改造．
RSS を生成する
plugin/makerss.rb
をひねくって，
「ツッコミ」なるコメント投稿行が RSS にでないようにする．
これはコードをいじらなくてもオプション設定でなんとかなるはずなんだけど，
どうにもうまくいかなかったため．
しかし RSS にまで「ツッコミ」って，
なんとも変な仕様だなぁ
……
それとも皆さん RSS にかなり依存してる，
と考えるべきなのか?
-
昼飯．
昼飯後も解説かきぢたばたが継続する．
-
ちょっとプリンター設定まわりの試行錯誤．
2008-06-04
に CUPS による印刷にきりかえたのだが，
CANON が
ダウンロードペイジ
で親切にもくばっている LIPS LX のドライヴァー
cndrvcups-lipslx-1.70-1
に問題があるような気がする．
図の中に破線があらわれると，
そのあとも実線で描かれるべき部分が破線になってしまう，
という症状だ．
対処療法としては LIPS LX をあきらめて，
LIPS 4 の
cndrvcups-lips4-1.70-1
をインストールして使用する，
というもの．
まあ，
LIPS 4 でもぜんぜん問題ないけど．
-
解説かき，
ぢりぢりとススめているつもりなんだけど，
実際のところなかなか進捗してない．
1940 研究室発．
1955 帰宅．
ぐったり．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0740):
食パン．
ジャガイモ・ネギ・ブナシメジ・鶏肉のスープ．
- 昼 (1330):
食パン．
ジャガイモ・ネギ・ブナシメジ・鶏肉のスープ．
- 晩 (2130):
ジャガイモ・ネギ・ブナシメジ・鶏肉のスープ．
ヨーグルト．
リンゴ．
2008 年 09 月 30 日 (火)
-
0800 起床．
朝飯．
コーヒー．
0910 自宅発．
晴．
0925 研究室着．
-
また
car.normal()
解説かきのつづき
……
-
数式なんかが登場しない「とっつきやすいみかけ」
(みかけ，だけなんだが)
の解説文をめざしていたんだけど，
数式とかがないとどうにもこうにも説明がうまく書けないので，
いつもの解説文のごとく数式とかもいれてみることにしてみる．
-
そして図がどんどん増えてる．
えーと，
caption の幅は，
たとえば
\setlength{\captionwidth}{0.7\paperwidth}
といったかんぢで設定する，
と
……
-
今日は帰宅して昼飯にするつもりだったんだけど，
作業がなかなかおわらないので
-
なにやら水中いきもののデータ解析そうだんメイルが
……
-
あら，
鍋嶋さんからお子さん (もうすぐ 5ヵ月) の写真が
……
-
ということで，
解説かきはいったんほうりだして，
苫小牧樹木直径成長モデル論文の cover letter
を再チェックしてみる．
まあいいかなあ，
と思えたので鍋嶋さんにお送りしてみる．
-
共立出版から
R 作図わざ本
についての問い合わせ．
印刷出力したものの点検はまだやってないんだけど，
とりあえずここまではできてます，
ということで，
ほんやく原稿 PDF ダウンロードする方法を連絡してみた．
-
car.normal()
解説かきのつづき
……
私は空間統計学に興味あるんだけど，
けっこう基本的なところで誤解していたような気がしてきた．
-
久保の誤解?: GRF (Gaussian random field) は
多変量正規分布を事前分布とする確率場，
と考えればよい
-
本日の勉強の結果: これが正しいのかどうかは自明ではない．
上の定義があてはまる場合もある (R の
library(geoR)
とか，ゐんばぐすの spatial.exp() とか)．
少なくとも，
この解説でとりあげている CAR (conditional auto regressive)
model は「ぱっと見」の形式上はどうヒネくっても
多変量正規分布にはなりそうにはない
-
GRF はマルコフ確率場なので，
それはギブス確率場である，
というのは問題ないだろう．
-
マルコフ確率場:
有限の「近傍」の状態によってある地点の状態をあらわす
確率分布が定義される確率場．
-
ギブス確率場:
ギブス分布で定義される確率場．
Hammerseley-Clifford の定理によって，
マルコフ確率場はギブス確率場である．
またギブス確率場において「相互作用」の範囲を限定すると，
それはマルコフ確率場になる．
-
さてさて，
問題の空間おーとれぐれっしぶ (ここでは auto-norimal)
なモデルについてだが
……
Banerjee S, Carlin BP, Gelfand AE (2004)
``Hierarchical Modeling and Analysis for Spatial Data''
とか間瀬・武田 (2001) の「空間データモデリング」
の定義をみると
……
-
conditonal autoregressive model (CAR):
「ぱっと見」にはどうヒネくっても
多変量正規分布にはならない，
しかし局所状態は正規分布として与えられる
ギブス分布で定義されたギブス確率場 (== マルコフ確率場)
-
simultaneous autoregressive model (SAR):
勉強ぶそくでよくわからんけど，
CAR とだいたい同じ，
ただし確率場の定義だけで全体の平均がゼロとなっている．
これをさらにいろいろと単純化して，
各地点の分散が同じで，
spatial autoregression parameter ρ
なるもので結合の強さを表現できるとすると
……
これって多変量正規分布?!
-
うーむ，
とアタマをヒネりつつ
(後記: やはり上の CAR に関する記述のうち
「ぱっと見にはどうヒネくっても多変量正規分布にはならない」
はちょっとアヤしい)
1930 研究室発．
1945 帰宅．
うんどう．
体重 70.6 kg．
晩飯の準備．
晩飯．
-
[今日の運動]
-
[今日の食卓]
- 朝 (0830):
米麦 0.6 合．
- 昼 (1350):
研究室．
北大生協サンドイッチ．
- 晩 (2200):
米麦 0.6 合．
チンゲンサイ・ナス・ネギ・干しシイタケ・干しエビのスープ．
ヨーグルト．
リンゴ．
