【第776回】 メールとか打ち合わせ会議とか（2014年9月4日）
- 6:00起床。2番のバスで出勤。
- 今日も英文を含むメールへの返事と打ち合わせ会議で時間が過ぎていく。合間にレポート採点の続き。
- それにしても全米オープンで錦織選手がベスト4に入ったのは素晴らしい快挙だと思う。
- 殺虫剤を撒いた日より後にサンプル（ドライアイスベイトトラップ）した蚊からデングウイルスが検出されたという報道からすると，あの程度の殺虫剤散布では不十分だったのだろう（このトラップによっても減っているはずとはいえ）。となれば，患者から吸血した蚊の寿命が尽きる9月末くらいまでは，代々木公園周辺で新規感染が続く可能性はゼロではないので，当面蚊が好みそうな場所は封鎖するという判断は合理的だと思う。あと，取材に行く人は，日本で普通に売られているスプレーの虫除けではなく，DEETをたっぷり含んだtropical strengthの昆虫忌避剤（液体の）を塗って行くべきだろう。
- 16:30からの某重要会議は無事終了。
- ふとpyramidパッケージでオーバーレイ描画させるためのアイディアを思いついて関数pyramidf()を作ってみた。par()からframeのデフォルトを受け取る部分がxaxpとyaxpではいけないのか，ずれてしまう問題が解決しないので，とりあえずオーバーレイには下記のように空plot()しておいて，同じframeを指定する必要がある。この問題が解決したらpyramidパッケージに組み込んで1.4版にするつもり。実行例としては，上でリンクしてあるpyramidf.Rを実行した後で下記のコードを実行（fmsbパッケージは日本人口を使うために呼び出すだけで，関数を使うわけではない）すると，その下の画像を描くことができる。
library(fmsb)
plot(c(0,100), c(0,100), type="n", frame=FALSE, axes=FALSE, xlab="", ylab="")
# overlay
pyramidf(data.frame(Jpopl$M1960, Jpopl$F1960, Jpopl$Age), frame=c(10, 90, 60, 100),
Clab="", Cstep=10, Laxis=0:3*5e+5, AxisFM="d")
pyramidf(data.frame(Jpopl$M2010, Jpopl$F2010, Jpopl$Age), frame=c(10, 90, 60, 100),
Clab="", Lcol=NA, Rcol=NA, Cstep=10, Laxis=0:3*5e+5, AxisFM="d")
# pararell
pyramidf(data.frame(Jpopl$M1960/10000, Jpopl$F1960/10000, Jpopl$Age), frame=c(0, 45, 0, 50), Clab="",
Cstep=10, Laxis=0:3*50, AxisFM="d", Llab="1960年男性（万人）", Rlab="1960年女性（万人）")
pyramidf(data.frame(Jpopl$M2010/10000, Jpopl$F2010/10000, Jpopl$Age), frame=c(55, 100, 0, 50), Clab="",
Cstep=10, Laxis=0:3*50, AxisFM="d", Llab="2010年男性（万人）", Rlab="2010年女性（万人）")
- この関数とmaptoolsを使えば，地図上に自治体単位で人口ピラミッドを書き込むことも可能になるはず。年齢区分をある程度大雑把にしないと見えないかもしれないが。
- 18:30から30分ほど，博士論文題目提出届について等，Alimさんから相談を受けた。レポートを見るのはもう少しで終わりそうだが，20:00までだと少し無理か。
- 直通終バスで帰宅。挽肉とピーマンとナスの炒め物を作り，レトルトご飯を電子レンジで温めて晩飯にした。軽井沢ビールが美味だったのだけれども，少しく酔ってしまったので，今夜はもう仕事は無理だな。
△Read/Write COMMENTS
▼前【775】（今日は名谷で会議いくつか（2014年9月3日） ） ▲次【777】（院生統計解析相談と成績処理（2014年9月5日） ） ●Top