ぎょーむ日誌 2013-02-03

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 合． エノキダケ・油揚げの味噌汁． カレイの干物． ハクサイ・タマネギ・ブナシメジ・豚肉の炒めもの．