ぎょーむ日誌 2012-06-28

2012 年 06 月 28 日 (木)

• 0620 起床． ラジオ体操． 弁当・朝飯の準備． 朝飯． コーヒー． 0750 自宅発． 今日も晴． 朝から気温たかい …… まあ， 日陰に入れば涼しいのが札幌のいいところだけど． 0820 研究室着．
• ぐだぐだした一週間を脱出して， なんとかデータ解析したうけぎょーむにとりくむ． 「どうもうまくいかないな」 と放置していた時系列データ解析もの ……
• どうもうまくいかないな ……
• 昼飯．
• 昼飯を食いながら， あれこれ点検してたんだけど …… 食い終わった直後に， 問題点が特定できた． 送っていただいたデータ変換プログラムだ (+ 私の誤用) …… 今まで事前分布がまずいのかとか， 初期値の与えかたに工夫が必要かもとか， いろいろやっていたんだけど …… こんなところに原因があったとは! 驚きのあまりずっこけて椅子から落っこちそうになり， さらに疲労感の奇襲をうけてぐったりしてしまった．
• 気をとりなおして， 推定計算やりなおし． 一階のモデルは何も問題なし …… 二階はちょっとてこずりそうなので， とりあえずヤメ．
• そのあたりの速報じみたものを， 鍋嶋さんにおくってみる．

• 個体ごとの中央値 (オレンジ) ががたがたしているのは， 「全個体に共通する年変動 (common.y[y])」 が入っているため …… これをぬくともっとなめらかな変化になる．
• 全個体共通のトレンドは年 0.4% ぐらいの増加， か? (事後分布表, BUGS code)．

• 計算時間は 250 秒ぐらい， ThinkPad X201S， クアッドコアなので R の並列計算機能を使った …… ようするにゐんばぐすを 3 つ同時に走らせるわけで， これをやると計算時間は少なくとも半分にはなる．
• ひとくぎりついたよーな気がしたので， 「統計モデリング入門」 4 冊をかかえて， 2F の環境図書室 (北大図書館の分室) に行って， 著者寄贈．
• しかし， ここ一ヶ月ぐらいの「ぎょーむ日誌」をながめてみると， どたばたと意味もなく右往左往しているような ……
• 気をとりなおして， もうひとつの放置ぎょーむ再稼働 ……
• いろいろなことがアタマの中から脱落していて， 何が何やらわからない状態．
• しかもねむい． さいきん， ねぶそくぎみでもある．
• とりあえず， ややこしい結果をまとめる作図を試みてみる．

• ふーむ， わりと小さい genet でもヘンなことが発生している?

  > table(dlong[as.numeric(dlong$locus) == 1, "Genotype2"])
  
     g01    g02    g03    g04    g05    g06    g07    g08    g09 
      37     28     21     19     14     10      6      5      5 
     g10    g11    g12    g13    g14    g15    g16    g17    g18 
       5      5      4      4      4      4      4      4      4 
     g19    g20    g21    g22    g23    g24    g25    g26    g27 
       3      3      3      3      3      3      2      2      2 
     g28    g29    g30    g31    g32    g33    g34    g35    g36 
       2      2      2      2      2      2      2      2      2 
     g37    g38 single 
       2      2     86 
  

• まあ， これは地図と対応させんといかんな． なかなかしんどい．
• しかし今日はここまで， か …… 1740 研究室発． 買いもの． 1810 帰宅．
• 「ねむいよー」と思いながら買いものしてたら， ぼろぼろと買い忘れがあって， またコープさっぷろに行って買いもの再履修．
• すると連絡があって， 今日の晩飯は「かんたん晩飯」でよいという事態になった． つまり買いもの再履修は必要ではなかった．
• いやー， でもかんたん晩飯はかんたんだからいいですよ． かんたん晩飯の準備． かんたん晩飯．
• [今日の運動]
  • 朝のラジオ体操．
• [今日の食卓]
  • 朝 (0700): シリアル． ヨーグルト．
  • 昼 (1200): 研究室． 米 0.7 合． ダイコン・ニンジン・シイタケ・鶏肉の煮物． 野菜天． スナップエンドウ． ブロッコリー． ミニトマト．
  • 晩 (1930): 米 0.7 合． ニラ卵丼． ホウボウ の塩焼き …… 奇妙な外見とはうらはらに， くせのないたんぱくな味わいだった．