厚労省はCOVID-19関係のデータをオープンデータとして公開するようになった。これをRでプロットしてみる:
df = read.csv("https://www.mhlw.go.jp/content/pcr_positive_daily.csv", fileEncoding="UTF-8-BOM")
plot(as.POSIXct(df[[1]]), df[[2]], type="h", xlab="月", ylab="感染者数", lwd=3, col="orange")
さらに死亡数を重ねてみる:
df2 = read.csv("https://www.mhlw.go.jp/content/death_total.csv", fileEncoding="UTF-8-BOM")
points(as.POSIXct(df2[[1]])[-1], diff(df2[[2]]), type="h", lwd=3, col="red")