厚生労働省:職場における熱中症の予防についてより暑さ指数(WBGT)の表をプロットしてみた。
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でごっそり取得する - まだ厨二病