ぎょーむ日誌 2011-10-02

2011 年 10 月 02 日 (日)

• 0720 起床． 朝飯． コーヒー．
• 怠業．
• また Cardamine な統計モデリングのつづき ……
• あまりにもうまくいかないので， GeoBUGS の spatial.***() な関数にもテをだしてみる ……

• 予想された結末ではあるが， spatial.***() な統計モデリングは， サンプリングが 100 倍ぐらい遅くてつかいものにならない ……
• 昼飯の準備． 昼飯．
• 時間ぎれなので， Cardamine したうけはいったん中止． 無念．
• 次のしたうけぎょーむに着手するまえに， 岩波統計モデル本 おくれのおわびメイル (複数) かき， その他メイルかき．
• 放置してた事務雑用．
• しんどい．
• ラジオ体操． ひとやすみ．
• そして次の「危険な状態になっているしたうけぎょーむ」 マラリア蚊の統計モデリングの準備． 104 行 (これはふつー) × 211 列 (これは異常!) な CSV ファイル．
• これは空間構造あるんで， まずは地図づくりだ． 調査地は このへん． どうでもいいが， この Google Maps の Victoria Lake ケニア側の地図はちょっとデータおかしいな． 湖面に意味不明なかくかくたちが．

• 自分の作図にはふつーに shapefile とか探したほうがいいか …… 検索すると， Kenya の Administrative areas - Latitude and longitude coordinates - ESRI Shape file format KEN_admin_SHP.zip …… ダウンロードして， ほどいて R で作図．

  ken <- readShapeSpatial("KEN")
  

  これはケニヤの行政区全体が入っている．

  plot(ken)
  

  ふーむ， 場所によっては細かく分割されてるのね． 部族ごととか?

• で， Victoria 湖畔の調査地はどのへんかな， と sp オブジェクトをまさぐってみると，

  < ken@data[541:541,]
       ID     LBL FIP MMT_ID SHORT__FRM          LONG_FRM  ADM0   ADM1     ADM2  ADM3  ADM4
  540 541 KEN-541  KE    KEN      Kenya Republic of Kenya Kenya Nyanza Homa Bay MBITA GEMBE
         ADM5 STL.0 STL.1 STL.2 STL.3 STL.4 STL.5
  540 MIITINE   121     6     3     3     -     -
  

  このへんだとわかったので， 作図．

  plot(ken[541:553,], xlim = c(34.1, 34.3), ylim = c(-0.463, -0.46), col = "#aaffaa", bg = "#aaaaff", axes = T)
  points(d$Lon, d$Lat, pch = 20, col = "#ff400080")
  

  赤丸が調査した家．

• すでに日没時刻． ちょっとうんどう． 1850 自宅発 ミニ大通公園 走． かなり寒い． 雨ぱらぱら． 1920 帰宅．
• 外で晩飯．
• [今日の運動]
  • ラジオ体操． 腕立ふせ 10 × 3 回． 腹筋運動 30 × 3 回． ミニ大通公園 走 30 分間．
• [今日の食卓]
  • 朝 (0740): ぱん吉の角食． チーズ． ヨーグルト．
  • 昼 (1220): スパゲッティー． レトルトパウチドのカニトマトソース． キュウリ・タマネギ・トマトのサラダ．
  • 晩 (1950): ミニ大通公園の西端ちかくの Agog で玄米プレート， ラムのグリル．