「ぎょーむ日誌」目次に戻る | KuboWeb top に戻る | twilog | atom

ぎょーむ日誌 2013-02-(01-10)

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

- 01/31 | 02/01 | 02/02 | 02/03 | 02/04 | 02/05 | 02/06 | 02/07 | 02/08 | 02/09 | 02/10 | 02/11 -

本日 (kubolog20130201) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

2013 年 02 月 01 日 (金)

0710 起床．ラジオ体操．朝飯の準備．朝飯．コーヒー． 0820 自宅発．曇天．気温高くてすでにプラス 1 度． 0850 研究室着．
なにやらぎょーむな会議，いきなり．
それとは関係なく，先日の「検閲」のときに，ぎょーむ日誌生成システム (Perl で書いてるだけなのですが) のめんどくささ・ややこしさにイヤけがさしたので，そのあたりを少し改良 …… 20yy/kubologyyyymmdd.html を含む 20yy/mmdd.html を変更したい場合には，まず 20yy/kubologyyyymmdd.html だけを修正してから generate_log10.pl 20yy/mmdd.html とする …… すぐにど忘れする自分自身のための備忘録．
これってつまらん改訂だけど，全体かきなおしのときとかには便利かな …… いや，そういう場合には ls log/2*/[0-9]*.html > updated.txt とすれば，引数ナシで自動処理してくれるか．
みなさん，本日配布された「本学職員の退職手当の支給水準引下げについて」とやらに興味があるようなので，とりあえずモデルケイス (影響額) とかいうファイルをネット上に置いてみる …… まあ，これは以前のごとく検閲されないでしょう．公知とかいうやつ?
2/18 には高松 (香川県のほう) でアルゼンチンアリ発表をしないといけないので，アリデータ整理のつづき．
昼飯．
ちょっと雑用がらみの文献よみ， 2 時間ほど．
廿日市アリ群集データを整える作業．ここまでやってみて，調査年・調査公園・調査アリ種の組み合わせが完全にそろっているわけではなく，「ところどころぬけている」ということがわかった．ぬけている場所をうまくさがして，欠測 (NA) あつかいにしないといけない．
これはばか正直に for() とかやってたら遅くてしょうがないので， exapnd.grid() とか tapply() を駆使して高速処理．
とりあえず図示してチェック …… うまくいってるようだ．

しかしさらに詳しく調べていくうちに，調査した公園データに何か不備があるらしいとわかった …… 77 公園のはずなのに 78 水準になっている．原因がわからないので，これは「目でみて」チェックするほかない．

                    10                     13                     14                     15 
          "佐方川緑地"           "佐方上公園"             "佐方公園"         "ちびっこ広場" 
                    16                     17                     18                     19 
      "月見台第２公園"           "緑が丘公園"       "佐方４丁目公園"   "緑が丘ちびっこ広場" 
                    20                     21                     22                     23 
          "岡の下緑地"         "中央第４公園"         "三宅第１公園"         "中央第３公園" 
                    24                     25                     26                     27 
        "中央第１公園"       "海老園第１公園"        "楽々園第2公園"       "桜尾３丁目公園" 
                    28                     29                     30                     31 
      "住吉大神宮公園"         "公民館前公園"         "深江児童公園"         "対厳山集会所" 
                    32                     33                     34                     35 
  "ふじタウン第１公園"       "ふじタウン公園"   "ふじタウン第３公園"       "柿ノ浦４号公園" 
                    36                     37                     38                     39 
      "柿ノ浦第３公園"         "宮島口浦公園"       "阿品台第３公園" "地御前ハイツ第１公園" 
                    40                     41                     42                     43 
            "野坂公園"           "橋本集会所"           "上平良公園"       "陽光台第６公園" 
                    44                     44                     45                     46 
        "宮園中央公園"         "宮園中央公園"         "宮園第３公園"               "楠公園" 
                    47                     48                     49                      5 
        "四季ヶ丘公園"           "宮内天王社"           "佐原田公園"             "畑口公園" 
                    50                     52                     53                     54 
         "串戸第2公園"             "田屋公園"        "阿品台第2公園"        "阿品台第1公園" 
                    55                     56                     57                     58 
       "阿品台第7公園"        "阿品台第6公園"    "阿品ハイツ第1公園"        "阿品2丁目公園" 
                    59                      6                     60                     61 
            "鰆浜公園"           "ひかり公園"           "可愛川公園"             "可愛公園" 
                    62                     63                     64                     65 
            "郡塚公園"          "郡塚第2公園"          "郡塚第3公園"        "陽光台第1公園" 
                    66                     67                     68                     69 
       "陽光台第3公園"       "陽光台第４公園"          "宮園第1公園"        "峰高台第2公園" 
                     7                     70                     71                     72 
              "桂公園"         "御手洗川公園"          "扇園第2公園"           "榎之窪公園" 
                    73                     75                     76                     77 
        "扇園第１公園"             "後山公園"           "丸子山公園"               "堀公園" 
                    78                     80                     81                     82 
        "余田第２公園"           "六本松公園"       "玉ノ井第３公園"       "玉ノ井第１公園" 
                    83                     84                     85                     86 
      "玉ノ井第２公園"         "鑓出第３公園"         "鑓出第２公園"         "鑓出第１公園" 
                    87                      9 
            "速谷神社"           "佐方川公園"

もちろん，上のような文字列の羅列もみたってまちがいが見つけられるわけではないので，ここでもまたいろいろと工夫 …… わかった! 陽光台第６公園 は公園番号 43 であるのに 44 となっている …… はい，これは例によって ゑくせる手作業 による入力・編集 (こぴぺ) の問題ですね．
この程度の規模， 12 年間 × 77 公園 × アリ 40 種というデータであっても，こういうまちがいがあるかどうかは R で探すしかないし， R で修正するしかない．
1800 研究室発．買いもの． 1830 帰宅．晩飯の準備．晩飯．
深夜のとらちゃんねかしつけ，哺乳瓶内のミルク残量に留意しつつ，なだめすかして …… なんとか寝てくれた．哺乳瓶はほぼカラ．燃料タンクほぼカラで空港に着陸できた飛行機 …… ってかんじでしょうか．
[今日の運動]
- とらちゃんえくささいず．
[今日の食卓]
- 朝 (0740): 米 0.7 合．ニンジン・ジャガイモ・ネギ・シイタケ・油揚げの味噌汁．厚揚げソテー・ダイコンおろし．
- 昼 (1200): 研究室．米 0.7 合．キャベツ・ネギ・ハムのサラダ．リンゴ．
- 晩 (1920): ハクサイ・ダイコン・豆腐・鶏だんごの鍋．卵雑炊．

本日 (kubolog20130201) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

本日 (kubolog20130202) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

2013 年 02 月 02 日 (土)

0720 起床．朝飯の準備．朝飯．
とらちゃんにミルクあげたり，回転 (ねがえり) の練習をカントクしたり，ねかしつけたり …… 下の画像は，ねがえりダブルアクセル (アクセルって人名だったのか) 44 秒動画を ImageMagick にて convert & montage したもの．

とらちゃんかついで，コープさっぽろ植物園店で買いもの．今日も気温たかい．しかし，明日からはまたものすごく寒くなるそーで．
帰宅してかんたん昼飯．
またミルク & ねかしつけ …… この時点でかなりぐったり．
Linux の Google-Chrome で Flush plugin が「入手できない」といった状態 (たとえば YouTube とか見れなくなる) になったら …… ここでのやりとりによると $HOME/.config/google-chrome/PepperFlash の名前を変えるか削除すればよいとのこと．実際にためしたらうまくいった．
またとらちゃんミルク & ねかしつけ …… るーぷ ……
晩飯の準備．晩飯．
23 時すぎに風呂に入ってると，むこうからとら坊の泣き声が …… じつはこれは地震予知だったようで， 2317 道東で M6.4 の地震，札幌は震度 3 ぐらいの揺れでした．
[今日の運動]
- とらちゃんえくささいず．
[今日の食卓]
- 朝 (0740): 米 0.5 合．厚揚げ．
- 昼 (1330): パン．ホウレンソウ・トマト・卵の炒めもの．
- 晩 (1920): パエリア，イカ・アサリ・鶏肉・タマネギ・ニンニク・ピーマン・パプリカ．

本日 (kubolog20130202) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

本日 (kubolog20130203) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

2013 年 02 月 03 日 (日)

0750 起床 …… よくねた，というか気絶したみたいに寝てた．天気は晴．予報されてたような悪天候ではない …… 朝飯の準備．朝飯．
またとら坊の「転がる練習」につきあう．手びょうし & 方向転換担当．だんだん泣きが混入してきたので，スリングにほうりこんでねかしつけ．
とら担ひとくぎり，コーヒー．ちょっと怠業．
廿日市アリ群集時系列データの整理のつづき．巨大なデータの一部に「ゑくせる手作業によるまちがい」があるとわかったので，まずはその位置を R で特定したい．
```
> d[d$公園No. == 44 & d$公園名 == "陽光台第６公園", 1:8]
     公園No.         公園名 chip 公園面積 調査年度 調査月日 見つけ回数 count.01
4476      44 陽光台第６公園 Ito2     1170     2011  9月26日          4        0
```
まちがいは一行だけだったのか …… ということで，データ自動変換用の R スクリプトに下記の一行を追加しておく．こういう場合はデータをいじらずに，データ読み込みプログラムで修正するのが鉄則だ．
```
	d[d$公園No. == 44 & d$公園名 == "陽光台第６公園", "公園No."] <- 43
```
この修正データを使って割算値図の作図やりなおし …… とりあえず …… 「アルゼンチンアリの多い・少ないで sort した公園ごと時系列 plot」でどうだ? 何か傾向があるような気がする．

データ解析ってのは，こういう図が描ければ，まあ「半分以上は終了」と考えてよいでしょう． R 大将 Ripley さんも「あるデータ解析のやりかたを考えるのと同じぐらい，そのデータ視覚化の方法の検討は重要」 (ぎょーむ日誌 2006-12-18) と日本人 R 使いたちあいてに訓示しておられる．
昼飯の準備．昼飯．
とら担したり，怠業したり …… 天気も雪がふったり，いきなり晴れたり．

またアリ群集データを調べる．アリの種ごとの標本個体数．

> sort(tapply(ds$n, ds$spc, sum, na.rm = TRUE))

   8   22   25   30   39   10   12   18   38    5   37    4   27   19   26    2    6   14   23   35 
   1    1    1    1    1    2    2    3    3    6    6    8    9   10   12   13   15   28   75   81 
  21   28   15   36   31   32   40   13    9   16   33   34    3    7   20   24   17   29   11    1 
  92  101  113  115  261  270  288  300  311  313  443  543  574  605  958 1006 1724 1832 1946 2174

そして，「いくつの公園に出現したか?」という集計．

> s <- (ds$n > 0) & (!is.na(ds$n)); sort(tapply(ds[s, "park"], ds[s, "spc"], function(v) length(unique(v))))

 8 22 25 30 39 10 12 18 27 35 38  5 37 19  4  6  2 26 14 28 21 23 36 16 34 15 13 40 31  9 20  1  3 
 1  1  1  1  1  2  2  2  3  3  3  4  4  6  7 10 11 11 15 21 26 30 30 31 33 37 38 41 44 46 50 52 52 
24 29 11 33  7 32 17 
56 57 58 59 62 62 75

AEON 桑園に買いものにいく …… 気温はマイナス 6 度ぐらいで，すごく寒く感じる．ここしばらく，「脱」真冬日が続いていたので．
晩飯の準備．晩飯．
[今日の運動]
- とらちゃんえくささいず．
[今日の食卓]
- 朝 (0900): パンケイキ．
- 昼 (1300): 米 0.7 合，マグロ丼．エノキダケ・油揚げの味噌汁．
- 晩 (1900): 米 0.7 合．エノキダケ・油揚げの味噌汁．カレイの干物．ハクサイ・タマネギ・ブナシメジ・豚肉の炒めもの．

本日 (kubolog20130203) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

本日 (kubolog20130204) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

2013 年 02 月 04 日 (月)

0710 起床．朝飯の準備．朝飯．コーヒー． 0740 自宅発．曇天．寒い． 0810 研究室着．
アリ群集動態モデルの BUGS code 化 ……
うーむ，いきなり複雑なモデルでは無理そうか．単純化，単純化．
昼飯．
おかげさまで，統計モデリング入門の第 6 刷が 3 月にでるようです!
とりあえずアルゼンチンアリの影響だけを考慮したものを ……
試行錯誤つづく．うっかり大量のパラメーターの記録を命じてしまったので， MCMC サンプリングに時間がかかってしまった．
うーむ …… サンプルの様子をくわしくみてみると，これを単純に長くしただけではダメらしい，とわかった．何かべつの工夫が必要．

censoring and truncation わざである I(min, max) を使ってみるか …… というあたりで時間ぎれ．
1810 研究室発．買いもの． 1850 帰宅．とらちゃんだんす．晩飯．
よーやくねてくれた．
お，昨日掲載するのを忘れてた画像が …… いやー，土日はやっぱりアタマがぐるぐるしてますな．

[卵ぐるぐる]

朝のパンケイキやくときに，卵をあわだてるとよいという記述がどこかにあったので，ブレンダーでやってみた (粉とか混ぜるのもブレンダー) …… あっというまに一様化されてしまったので驚き，写真をとった次第．これは他の料理とかにも応用できそう．

[えほうぐるぐる]

なぜかしら六花亭でかくのごとき「えほうまき」が売られていたので買ってしまった …… 350 円．

[今日の運動]
- とらちゃんえくささいず．
[今日の食卓]
- 朝 (0720): 米 0.7 合．ネギ・油揚の味噌汁．
- 昼 (1210): 研究室．米 0.7 合．カタクチイワシのショウガ甘露煮．リンゴ．
- 晩 (1950): 米 0.7 合．ハクサイ・ブナシメジ・豚肉のクリーム煮．サケのムニエル．

本日 (kubolog20130204) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

本日 (kubolog20130205) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

2013 年 02 月 05 日 (火)

0640 起床．ねむい．コーヒー．朝飯の準備．朝飯 …… とらちゃん，いきなりばくはつぶつ事件．二重のかくのうようき破損．洗濯． 0840 自宅発．雪． 0900 研究室着．

[まいくろうぇいぶ加熱器具]

家で不要になったので，自分の研究室にひきとった …… 8F お茶区画には，すでにこれより良いやつが入っているので．まあ，昼めしどきは同区画は院生たちの排他的領域となることだし．ということで，今日はタクシー通学， 890 円でした．

これからイヤな GLM バナシを聞かされそうな気がしたので，ちょっと事前に検討 …… R の glm() で使う offset 項を log(t) としよう． family = poisson(link = "log") の場合，これは t に比例して平均が増大する …… といった解釈しやすいモデルとなる．
二項分布の場合は， family = binomial(link = "cloglog") つまり complementary log-log リンク関数を使えば，これまた解釈しやすいモデルになる．
それでは，もし仮に family = binomial(link = "logit") とした場合にはどうなるか? …… log(p / (1 - p)) = (線形予測子) となると考えると，これはたとえば p = A t / (1 - A t) といった式になる!! (A は exp(線形予測子)) …… といったことに初めて気づいた．
これは生態学では Holling の円盤方程式として知られる奇妙な式で …… しかし，まああまり使いみちのなさそうなハナシだなあ．
で， 1030 より研究室セミナー，本日は生駒君の修論ちょくぜん …… だいじょうぶかなああああ
昼飯．さっそくまいくろうえいぶ!
午後は廿日市アリ群集モデルの改良のつづき．姑息な BUGS censoring わざ，すなわち dnorm(m, tau)I(-3,0) といった切断正規分布モデルにすると，今までうまくいかなかった計算がうまくいった ……!!
これは観測値 (見つけた個体数) が「ゼロ」というのがやたらとある場合に，かくのごとく「下側」を censoring しとかないと「小さい値ならなんでもいいんでしょ」とばかりにてきとーな値でひっかかってしまって， MCMC 試行ごとに異なる結果が得られる …… といったほにゃららな現状があるため．
姑息な censoring わざはこのやっかいな問題をとりつくろうのに使える …… ああ，以前解決できなかった「zero-inflated な空間構造ありデータ」問題なんかも，これでうまく「ごまかし」できそうだな．
ともあれ，上位 10 種だけのアリ群集については計算がうまくできるようになった．アルゼンチンアリに対抗する「戦闘力」は種ごとにけっこうちがうみたいだ ……
Unix というか bash 上における日本語環境設定，とりあえず (冗長かもしれないが) export LC_ALL=ja_JP.UTF-8 LC=ja_JP.UTF-8 LANG=ja_JP.UTF-8 としておけば OK．
研究室のぶっこわれプリンター (EPSON LP-S6000) 問題，けっきょく (めんどくさくなった甲山さんの決定で) 新しく買い直すことになった． LP-S6000 用の消耗品があれこれたくさんあるんだが ……
また，午前の生駒問題とからんで，指導教員たる甲山さんから樹木の死亡モデルあれこれについて，いろいろと説明しなければならぬことに ……
20 種で 2600 秒ぐらい …… アルゼンチンアリぱわー，たいしたもんですね．

1730 研究室発．買いもの． 1810 帰宅．晩飯．
[今日の運動]
- とらちゃんえくささいず．
[今日の食卓]
- 朝 (0720): 米 0.7 合．ジャガイモ・タマネギ・油揚の味噌汁．
- 昼 (1210): 研究室．米 0.7 合．レトルトパウチドのトマトソース．リンゴ．
- 晩 (1920): 米 0.7 合．ジャガイモ・タマネギ・油揚の味噌汁．サンマのナンプラー煮．ジャガイモのあげむし．

本日 (kubolog20130205) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

本日 (kubolog20130206) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

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 合．ネギ・エノキダケのすまし汁．鶏肉に「からあげ粉」をつけてオーブンでグリルしてできた，からあげっぽい何か．キャベツ・タマネギ・カニのサラダ．

本日 (kubolog20130206) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

本日 (kubolog20130207) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

2013 年 02 月 07 日 (木)

0650 起床．朝飯の準備．朝飯． 0650 自宅発． 0820 研究室着．コーヒー．
廿日市アリ群集モデルの改訂．「(各アリ種の) 戦闘力の平均」の平均というパラメーターを増やすと …… やはりというべきなのか，なぜかしらというべきなのか，ともかくうまくいった．北大くらうど PC で 4-core 並列で 4400 秒ぐらい．やっぱり，アミメアリはヘンだ．

これはあやしげな結果なのだけど，どうにもよくわからないので，とりあえずデータをとった伊藤さんにメイルを書いてみる．

...
とりあえず，現時点で得られている推定結果ですが……
やはりといいますか，アルゼンチンアリの「他種をおさえる」効果はかなり明瞭み
たいです (他の 4 種と比較して)．ただし，なぜかしらサクラアリはアルゼンチンアリ
によって個体群密度が下がるといった効果がないようです．アルゼンチンアリとサク
ラアリがちゃんと同じ公園にいるかどうかは，まだチェックしていませんが，何か
サクラアリに特徴的なことがあるのでしょうか．

クロヤマトアリとサクラアリは他種にあまり影響を与えないようです．

トビイロシワアリはアルゼンチンアリを減らすようです (ホントか?)．

まったくナゾで私には理解できないのがアミメアリで (統計モデルが不適切なだけ
なのかもしれませんが)，低密度 (見つけ数が多くない状態) でも他種に大きな効果
があり，アルゼンチンアリを増やす (!?) とかクロヤマアリ・トビイロシワアリを減
らすといった推定結果が得られております．
...

あとで伊藤さんから返信いただき，サクラアリ (ant species 4) のカラダの小ささとか，アミメアリ (ant species 5) が見つかる公園の特徴などを教えていただいた ……
計算まち時間に，メタアナリシスなどの勉強，論文よみ．今日の午後のセミナーの予習．
アリ群集のほうは内積 inprod(,) を使うと計算が遅くなるかどうかの検証をやってみることに．
実験の結果，どうやら，
- inprod(,) を使うと計算がおそくなる
- inprod(,) を使うと sampling にしくじって trap 窓が出る
…… ということらしい，とわかった．

ということで，ばかばかしいが，こういう BUGS coding にならざるをえない．

mz[park, year, ant] <- (
    z[park, year - 1, ant]
    + log.r[ant]
    + r.year[year]
    - alpha[ant, 1] * expz[park, year - 1, 1]
    - alpha[ant, 2] * expz[park, year - 1, 2]
    - alpha[ant, 3] * expz[park, year - 1, 3]
    - alpha[ant, 4] * expz[park, year - 1, 4]
    - alpha[ant, 5] * expz[park, year - 1, 5]
    - alpha[ant, 6] * expz[park, year - 1, 6]
    - alpha[ant, 7] * expz[park, year - 1, 7]
)

種数を 7 に増やしている理由のひとつは， 5 種のときに何かおかしなことが生じていなかったかを点検するため．しかしこれでも alpha[,] の個数は 25 個から 49 個に増えてしまった．計算に時間かかりそう．

昼飯．
勉強のつづき．
14 時から理学部 5 号館で，ニュージーランドの中川さんのセミナー (ぎょーむ日誌 2013-01-17)， ``Tackling big questions with meta-analysis: a case study of dietary restriction and longevity''．とてもおもしろいハナシだった (らくがきメモ note20130207.pdf)．

その後，なぜかしら (環境科学院の実験棟に教授室を移動した) 木村さんの部屋で，木村さん・中川さん・野間野さんと 17 時すぎまで楽しく雑談．
廿日市アリ群集モデル，ぢりぢり発展中．
1730 研究室発．買いもの． 1800 帰宅．とら担．晩飯．
…… さてさて，帰りぎわに投げておいた計算は 6000 秒ほどで終了．うーむ，やはり第 7 アリ種「ウメマツオアリ」の挙動がおかしい ……

[今日の運動]
- とらちゃんえくささいず．
[今日の食卓]
- 朝 (0710): 米 0.7 合．ネギ・エノキダケのすまし汁．キャベツ・タマネギ・カニのサラダ．
- 昼 (1230): 研究室．米 0.7 合．鶏肉に「からあげ粉」をつけてオーブンでグリルしてできた，からあげっぽい何か．リンゴ．
- 晩 (1920): 米 0.7 合．ネギ・カブの味噌汁．ホッケ．カブシュート・油揚の煮びたし．

本日 (kubolog20130207) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

本日 (kubolog20130208) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

2013 年 02 月 08 日 (金)

0700 起床．ラジオ体操．朝飯の準備．朝飯．コーヒー． 0800 自宅発． 0830 研究室着．
アリ群集の挙動を確認 …… うーむ …… 時系列データから相互作用を抽出するのは，けっこう難しいことかも．
野村さんのデータ解析こんさる．「面積比が s1:s2 であるときに，個体数がそれぞれの領域に s1:s2 で出現する --- かもしれないし，どちらかにかたよるかもしれない」という問題を R の glm() で考えようとすると，けっこうめんどくさいとわかった ……

これはロジスティック回帰ではなく， offset 項わざ + ポアソン回帰であつかうべき問題なのかもね …… しかし本日は時間ぎれ．
まさか超幾何分布とか Polya-Eggenberger 分布をもちだすようなハナシでもないよな ……
昼飯．
なぜかしら午後は新しく購入された (予算が赤字ゾーンにとびこんでいるにもかかわらず) A810 院生部屋用のプリンター Brother HL-4570CDW の設定でかなりの時間を使ってしまった．
このプリンターは無線でつながるので，そのように設定． Ethernet ケイブルがどんどん不要になりますな．
ここまではよかったのだが，自分の Linux 用の設定ではまってしまった …… Brother も Linux 用 driver などを準備してくれていて， hl4570cdwcupswrapper-1.1.1-5.i386.deb と hl4570cdwlpr-1.1.1-5.i386.deb をダウンロードして dpkg -i …… しかしうまくいかない．
ネットで検索しても解決策がみつからない …… どうも問題は device URI に問題があるらしく， CUPS は lpd://(なんだか暗号みたいなもの)/BINARY_P1 と設定してくれたのだが，これが機能してなさそう．
CUP では修正できないようなので，あちこちひっくりかえしてみたのだが …… /etc/cups/printers.conf 内に <Printer HL-4570CDW_A810> というセクションがあり，その中に URI の設定があったので，プリンターの IP アドレスにおきかえて DeviceURI lpd://192.168.1.57/BINARY_P1 などとしてみると，これで無事に印刷出力できるようになった．
1730 研究室発．買いもの． 1800 帰宅．晩飯の準備．晩飯．
[今日の運動]
- とらちゃんえくささいず．
[今日の食卓]
- 朝 (0720): 米 0.7 合．ネギ・カブの味噌汁．
- 昼 (1240): 研究室．米 0.7 合．レトルトパウチドシチュー．リンゴ．
- 晩 (1930): 米 0.7 合．ネギ・カブの味噌汁．サバの開き．ミズナ・油揚の煮びたし．サツマイモのショウガ・バター蒸し．

本日 (kubolog20130208) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

本日 (kubolog20130209) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

2013 年 02 月 09 日 (土)

0810 起床．ラジオ体操．朝飯の準備．朝飯．コーヒー．
今日は一日中，とらちゃんと家の中にいた．
午前は，大規模なばくはつぶつ事件の対処に苦慮．
昼飯．
午後は，哺乳瓶，スリングや Beco でとらちゃんだんす，よーやく寝てくれてぐったり …… ループをぐるぐる．

[こういうかんじで]

寝られてみたり．

晩飯の準備．晩飯．
[今日の運動]
- とらちゃんえくささいず．
[今日の食卓]
- 朝 (0840): ぱん吉の角食．トマト卵やき．
- 昼 (1300): スパゲッティー，ミズナ・ニンニク・ベイコンのアーリオーリオ．サツマイモのショウガ・バター蒸し．
- 晩 (1930): 米 0.7 合．ニンジン・タマネギ・鶏肉のスープ．チンゲンサイ・マッシュルーム・干しエビの炒めもの．はんぺん．

本日 (kubolog20130209) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

本日 (kubolog20130210) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

2013 年 02 月 10 日 (日)

0820 起床．ラジオ体操．朝飯の準備．朝飯．コーヒー．
今日もとらちゃんないちにちで終了 …… だっこして AEON に行ったり．
[今日の運動]
- 朝のラジオ体操．とらちゃんえくささいず．
[今日の食卓]
- 朝 (0830): ぱん吉の角食．ニンジン・タマネギ・鶏肉のスープ．
- 昼 (1320): 極楽湯の食堂，ごま蕎麦とソースカツ丼 (小)． 500 円．
- 晩 (1920): カジカ・ネギ・エノキダケの味噌鍋．卵雑炊．

本日 (kubolog20130210) | 次の日 | 1 日前 | 7 日前 | 31 日前 | 365 日前 | top

- 01/31 | 02/01 | 02/02 | 02/03 | 02/04 | 02/05 | 02/06 | 02/07 | 02/08 | 02/09 | 02/10 | 02/11 -

KuboLog | KuboWeb