ぎょーむ日誌 2013-05-17

2013 年 05 月 17 日 (金)

• 0620 起床． とら坊の朝ミルク …… 120 ml ほどふつーに飲んでもらう． 次に 20 ml ぐらいのミルクに抗生物質 (これで終了だ! ばんざーい!) を溶かしたものを …… かなり不審げながら飲んでくれた． つぎに別にとっておいた 20 ml にセキ止めクスリを …… あ， こっちは「つきあい程度」ですか， そうですか． 朝飯の準備． 朝飯． コーヒー．
• とら坊， あいかわらず「立つ練習 (つかまり立ち)」をすごく熱心にやっているんだけど， ときどきバランスをくずしてすてーんと転んでアタマをうったり． だいじょうぶかな? 鼻水とり器を使ったところ， 激怒． なだめながら， 0820 自宅発． 晴． 今日は平年なみに暖かくなるそうで． 0830 保育所着 …… あら， 今日は「ひよこ組」のみなさん， おそめの登園のようで． とらちゃんは何だか「ぽーっ」としております． 0850 研究室着．
• 22000 秒ほど費やして n.thin = 60 のサンプリング終了 …… 下の図に示してるがごとく， パラメーター (ま，何でもかんでもパラメーターなんだけど， 気分として「変数っぽくないもの」ぐらい) の収束は問題ないけれど， 「アリ密度」をあらわす z[park,year,spc] のうちいくつかはまだばらついている． これはあたりまえのことで， 事前分布から切断 I(lower,upper) を排除したので， たとえば「ある公園・ある年・ある種のアリの存在確率がほぼゼロ」 ということを表現する密度 z は -5 でも -10 でも「小さい値ならば何でもよい」状態になってしまい， とうぜんながら chain ごとにばらけてしまう．

• しょうがないので z[park,year,spc] に関しては両端切断した事前分布を復活してみるかな …… とりあえず試しに計算させてみよう． たぶん n.thin = 30 で問題ないはずだ．
• 新しくした LaTeX 環境 (TeXLive 2012) で Beamer などを使うための調査のようなことをうだうだと ……
• ちょっと試しに beamer な LaTeX ファイルをコンパイルしてみると …… LaTeX font error: I can't find file `pzdr' といったエラーがでる． ここらへんの問答 によると Zapf ding なフォントが入ってないのが問題のようで …… いちばん安直な解決は sudo apt-get install texlive-fonts-recommended かな?
• お， 統計モデリング入門 刊行一周年ということで， 岩波自然科学本アカウントで tweet してもらっている! みなさん， 今後ともよろしくお願いいたします．

• 昼飯．
• 灯油代のふりこみ， 過去二ヶ月間は 69.0 L で 7028 円． 前回の支払いは 12915 円だったので， まあ減っている …… しかしいまだに暖房をつけてるんだよね．
• 18000 秒ほどかかって切断正規な事前分布のモデルの パラメーター値のサンプリングができたけど (切断した分布だと計算が遅くなる?)， まあとくに何も改善されてないし， 結果は少し複雑化するので …… 今朝えられた結果をとりあえずのところ， 現時点での最新の結果ということにして， アリ群集論文のとりまとめにかかりますかね ……
• TeXLive 2012 と Beamer 調査のつづきあれこれ ……
• 1710 研究室発． 買いもの． 1740 帰宅． 晩飯の準備． 晩飯．
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0730): 米 0.7 合． ネギ・油揚の味噌汁． キャベツ・ネギ・ハムのサラダ． 小アジのショウガ煮．
  • 昼 (1240): 研究室． 米 0.7 合． 小アジのショウガ煮． ソーセイジ． ふかしたサトイモ． バナナ．
  • 晩 (1950): 米 0.7 合． ダイコンの味噌汁． サワラの西京味噌づけ焼き． キャベツ・ネギ・ハムのサラダ． コマツナのゴマあえ．