KuboLog 2013-05-09

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

2013 年 05 月 09 日 (木)

0600 起床 …… なぜかしら，本日とら坊は朝ねぼうでこの一時間後に起床．朝飯の準備．とら坊の朝ミルク & おむつ替え．朝飯．コーヒー． 0820 自宅発．晴． 0830 保育所着．今朝はあひる組 (一歳クラス) の子たちがひよこ組にあそびにきていたんだけど，体格はとらちゃんよりひとまわり小さい子でも，きわめて安定した二足歩行をしているというだけで「ぜんぜん別種のいきもの」に見えるなあ …… 0900 研究室着．

unix まめちしき: tar には --newer DATE というオプションがあり DATE より新しい日付のものだけを選択的に tar に追加する．圧縮版である tgz でも OK．

unix まめちしき: ファイルの変更日付 (でいいのかな?) をあらわす文字列をとりだす方法 …… おそらく他にもたくさんある中のひとつ．

  $ echo `ls -l --time-style=+%Y-%m-%d ファイル名` | cut -d' ' -f6
  2013-05-09 # ← 出力例

いきなり大学院生 (カトリンさん) のぶっこわれ PC の HDD さるべいじ作業にかりだされる ……

[かなり難しい PC 解体作業]

HP ProBook 4420s とかいうやつで，外観からは HDD 取り出し不可能に思えたのだが …… 不幸中のさいわいなのか， HP の web site にわりとくわしい解説があったり， YouTube に解体作業動画があって助かった． HDD とりだすだけなら上記動画ほどあちこち外す必要はなく， HP の解説ペイジどおりにやればよい．それでもけっこう難しいけどね．

HDD とりだしに成功したところで，いったん釈放されたので，昼飯．

午後はファイルコピーの作業手順を確立，その作業を自分でやってもらうようにしたところで，また一時的に釈放された．

アリ群集モデルの改善作業に復帰．とりあえず， density (連続値) みたいなものを occurrence ({0, 1}) に置き換えたモデルを試してみるか．

北大くらうどにもっていって計算させてみる．

待ち時間がけっこう長い．

二項分布とかの link 関数に complementary log-log というやつがあるが (そしてそれを今のアリ群集モデルで使っている)， complementary ではない log-log リンク関数もあるようだ (アイルランドのリムリック大学のサイト内の解説)．

一時間まっても終わらんな．やはり潜在変数を離散値にしたモデルってのは，ゐんばぐすは苦手そうだ ……

というか，そろそろゐんばぐすとかヤメたい …… 以前は R 内では bugs オブジェクトを使いたい (mcmc クラスのオブジェクトとかではなく) などといった理由づけをしていたんだけど，まあ，これはいまや相互に変換可能になりつつあるしな …… 完全に変換できるわけではないが．

4600 秒ほどで (3 プロセス並列化) 終了．結果は，やはり良くないな．何か単純なまちがいが残存しているように見える．

うーむ …… アリ種間相互作用をあらわすパラメーター alpha[defense, offence] の事前分布を「もと」にもどしてみるか …… すなわちあのヘンタイ的な切断正規分布 dnorm(mean, tau)I(-0.01,) に ……

推定計算やりなおし．また 5000 秒ぐらいかかるんだろう ……

3 月に発売された丸山さんたち (あ，あの小松さんも共著者の一人か!) のアリの巣の生きもの図鑑 (The Guests of Japanese Ants) の Amazon れびゅー，とてもアツいですなあ …… みなさんホントにムシ好きだ．新聞の書評とかでもいろいろとりあげられている …… ネットで掲載写真の一部がみられるが，たしかに写真ながめるだけでもすごい本だという思えてくる (私はとくにムシ好きではないけど)．

やっぱり c-log-log ではなく log-log をつかうべきだったかな?

計算はなかなか終わらないので撤退することに． 1720 研究室発． 1740 帰宅．晩飯の準備．晩飯．

こそだてすとれすか何なのか …… 自分ではまったく気づいてないんだけど，私は寝ているときにウナされていて，とつぜん「気をつけてー!!」とか叫んだり，ばりばりと歯ぎしりとかしているらしい …… なんかばててるみたいですね ……

[今日の運動]

とらちゃんえくささいず．

[今日の食卓]

朝 (0700): 米 0.7 合．ダイコン・豚肉のスープ．ゴボウ・ニンジンのきんぴら．
昼 (1200): 研究室．米 0.7 合．紅茶豚．ゴボウ・ニンジンのきんぴら．ミニトマト．リンゴ …… さすがに季節はずれになって (しかも安いやつなので) 内部が変質していてあまりおいしくない．
晩 (1920): 米 0.7 合．ダイコン・豚肉のスープ．カマスの干物．カブシュート・ネギ・ショウガ・ブナシメジ・卵の炒めもの．オカヒジキのカラシマヨネイズあえ．

ぎょーむ日誌 2013-05-09

2013 年 05 月 09 日 (木)