廿TT

譬如水怙牛過窓櫺 頭角四蹄都過了 因甚麼尾巴過不得

気温・湿度と熱中症の危険度(暑さ指数)の表

厚生労働省:職場における熱中症の予防について より暑さ指数(WBGT)の表をプロットしてみた。

f:id:abrahamcow:20160510230459p:plain

WBGT値 警戒水準
25℃未満 注意
25℃~28℃ 警戒
28℃~31℃ 厳重警戒
31℃以上 危 険

そろそろ熱中症に気をつけましょう。

library(pdftools)
library(pipeR)
library(tidyr)
library(dplyr)
library(ggplot2)
path1 <-"http://www.mhlw.go.jp/houdou/2009/06/dl/h0616-1b.pdf"
p4 <- pdf_text(path1)[4]%>>%
  strsplit("\n") %>>%
  unlist()
tab4 <-read.table(text=p4[-11],skip = 2,nrows = 20)
kion <- as.integer(rownames(tab4))
heatindex <- cbind(kion,tab4) %>>%
  gather(sitsudo,WBGT,-1)
heatindex <- mutate(heatindex,sitsudo=as.integer(sub("X","",sitsudo)))

tmp <- ifelse(heatindex$WBGT<25,1,heatindex$WBGT)
tmp <- ifelse(tmp>=25 & tmp < 28,2,tmp)
tmp <- ifelse(tmp>=28 & tmp < 31,3,tmp)
tmp <- ifelse(tmp>=31,4,tmp)

heatindex$level <- tmp

ggplot(heatindex,aes(x=kion,y=sitsudo))+
  geom_tile(aes(fill=factor(tmp)))+
  geom_text(aes(label=WBGT))+
  scale_fill_manual(values = c("1" = "cornsilk","2" = "yellow","3" = "orange","4"="tomato"))+
  xlab("気温")+ylab("湿度")+
  theme_bw(18,"HiraMaruProN-W4")+theme(legend.position = "none")

参考にしたページ:⭐️PDFの情報・文章をRでごっそり取得する - まだ厨二病