ぎょーむ日誌 2013-02-06

2013 年 02 月 06 日 (水)

• 0650 起床． 朝飯の準備． 朝飯． コーヒー． 0800 自宅発． 0830 研究室着．
• 廿日市アリ群集動態の試行錯誤のつづき …… とりあえず公園差・年差をいれてみたが， 公園差は余計だということがわかった． 考えてみれば各集団の動態は初期値を「自由」に選べるようにしているので， ここに公園差がすでにくみこまれているわけで．

• アルゼンンチンアリの戦闘能力 (alpha) はあいかわらずのようで．
• さて， 他の dominant なアリたちの密度効果みたいなものを組みこみたいのだが ……
• いろいろ試行錯誤してるうちに， 腸内細菌叢 (gut flora, bowel flora) のデータが送られてきた …… とりあえずデータを点検して作図． やっぱり， rat 個体間で差があるし， 1 rat からとられたサンプル間にも大きなばらつきがある． 二重の overdispersion．

• とりあえず図を見れるようにして， 「しばらくまってください」 というメイル．
• 昼飯．
• アリ群集動態の試行錯誤つづく． 40 種のうち dominant な 10 種を選んで， こいつらの戦闘力行列 (10 × 10) を推定させようとしたのだが …… ゐんばぐすでは「時間がかかりすぎる」と判明した．
• ということで， どんどん簡単にしていく試行錯誤．
• とりあえず「攻撃側」が主要 2 種の場合は， なんとか MCMC サンプリングできる， とわかった． ただし， これは「一方的に攻撃」といった単純な状況．
• この単純な状況で， 攻撃がわの種を 5 ぐらいに増やせないかな ……?
• 攻撃側 5 種， 全 10 種の場合． モデルの再検討が必要． 北大くらうど PC で 4-core 並列でやって 3350 秒ぐらい．

• 結果をみていくと， dominant 5 種の中にはアルゼンチンアリに「対抗」 していそうなやつも ……
• 本来ならモデルをもう少し改善すべきだが， 今日はもう時間がないので， ただたんに種数を減らし (10 種 → 5 種)， サンプリング間隔をちょっと長くして様子をみる …… といった姑息な計算を命じておいて撤退準備． ああ， 無能．
• 1730 研究室発． 買いもの． 1800 帰宅． 晩飯の準備． 晩飯．
• 今晩もとらちゃんにぶんぶんふりまわされて (ぐったり) …… よーやく寝てくれた．
• さてさて， 北大くらうどに投げておいた計算はどうなったかな …… 4600 秒ほどで終了していた．
• あいかわらずモデルに何か問題ありそうなかんじなのだが …… 「5 番目の species」 の挙動が何やらあやしい， というのは言えるかもしれない．

• 「5 番目のアリ species」 ってなんだよ， と調べてみると ……

  > v.antname[as.numeric(levels(dss$spc)[1:10])]
   [1] "アルゼンチンアリ"     "クロヤマアリ"         "トビイロシワアリ"     "サクラアリ"          
   [5] "アミメアリ"           "ハリブトシリアゲアリ" "ウメマツオオアリ"     "ルリアリ"            
   [9] "オオズアリ"           "ハリナガムネボソアリ"
  

  アミメアリ …… ああ， あの (永住的な) 巣を作らずに集団移動している連中か． いかにもあやしそうだな． くわしいことは伊藤さんにお尋ねしてみますか． ちょっと検索 してみると， アルゼンチンアリ vs アミメアリとか， アミメアリの凶暴さとか， いろいろでてくるな ……
• [今日の運動]
  • とらちゃんえくささいず．
• [今日の食卓]
  • 朝 (0730): 米 0.7 合． タマネギ・油揚の味噌汁． イワシの丸干し．
  • 昼 (1220): 研究室． 米 0.7 合． レトルトパウチドカレー． リンゴ．
  • 晩 (1920): 米 0.7 合． ネギ・エノキダケのすまし汁． 鶏肉に「からあげ粉」をつけてオーブンでグリルしてできた， からあげっぽい何か． キャベツ・タマネギ・カニのサラダ．