都道府県別飲用水(水道水)環境放射能測定結果

(このページは岡山大学の垂水先生のページのミラーです)

データ

• mext-water.csv
• mext-water.xls
• 元データ
  • 文部科学省　http://www.mext.go.jp/a_menu/saigaijohou/syousai/1303956.htm　より 上水（蛇口水）の調査結果 次のファイルより作成
    • 平成23年3月22日（火曜日）採取 （PDF:90KB）
    • 平成23年3月21日（月曜日）採取 （PDF:75KB）
    • 平成23年3月20日（日曜日）採取 （PDF:186KB）
    • 平成23年3月19日（土曜日）採取 （PDF:175KB）
    • 平成23年3月18日（金曜日）採取 （PDF:171KB）

R program

library("ggplot2")

str_replace<-function(d,pat,a)  {
  sub(pat,a,d)
}

dwater <- read.csv('http://face.f7.ems.okayama-u.ac.jp/~t2/data/stat/mext-water.csv', 
    fileEncoding="sjis", skip=2, na="-", as.is=TRUE)

dwater[(!is.na(dwater[,4]))&(dwater[,4]=="不検出"),4]<-0
dwater[,4]<-as.numeric(dwater[,4])

location<-dwater$都道府県
dwater$都道府県 <- str_replace(location,"（.*）","")

dest.pref<-c(8:15,19)
dtemp<-dwater[dwater$pref%in%dest.pref,]

lct <- Sys.getlocale("LC_TIME"); Sys.setlocale("LC_TIME", "C")
ttime <- strptime(dtemp[,3], "%m月%d日")

G <- ggplot(dtemp, aes(as.POSIXct(ttime), I.131, colour=都道府県)) +
     geom_line() + geom_point() + xlab("日時(日ー月)") + ylab("I-131(Bq/kg)")

G +  facet_wrap(~都道府県)