死亡者数

2021年7月27日まで

データの出所: https://www.mhlw.go.jp/content/death_total.csv

`

2020年5月9日以降

データの出所: https://covid19.mhlw.go.jp/public/opendata/deaths_cumulative_daily.csv

library(ggplot2)
library(stringr)
library(scales)
d1 <- ggplot(dat,aes(x=Date,y=死亡者数))+
  geom_line(color="grey20")+
ylim(0,70000)+
labs(title = "新型コロナウイルス死亡者数",subtitle = str_c(as_of,"までの結果"),caption = "厚生労働省オープンデータ\nhttps://covid19.mhlw.go.jp/public/opendata/deaths_cumulative_daily.csv",x="日付") +
scale_x_date(breaks = datebreaks,labels=date_format("%Y %b"))+
theme_minimal(base_family="HiraKakuProN-W3")+
theme(axis.text.x=element_text(angle = 30,hjust = 1))
d1

1日の増加を計算

# すべてにゼロを入れる。
dat$前日 <- numeric(nrow(dat)) 

# 2日目から最後の日までのところに前日のデータを入れる。
dat$前日[2:nrow(dat)] <- dat$死亡者数[1:nrow(dat)-1]

dat$増分<- dat$死亡者数 - dat$前日
head(dat)
tail(dat)
# データを2日目からのものにする。行数1減。
dat <- dat[2:nrow(dat),]
library(ggplot2)
library(stringr)
library(scales)
d1 <- ggplot(dat,aes(x=Date,y=増分))+
  geom_line(color="grey20")+
ylim(0,500)+
labs(title = "新型コロナウイルス死亡者数",subtitle = str_c(as_of,"までの結果"),caption = "厚生労働省オープンデータ\nhttps://covid19.mhlw.go.jp/public/opendata/deaths_cumulative_daily.csv",x="日付", y="死亡者数") +
scale_x_date(breaks = datebreaks,labels=date_format("%Y %b"))+
theme_minimal(base_family="HiraKakuProN-W3")+
theme(axis.text.x=element_text(angle = 30,hjust = 1))
d1

移動平均を計算

平均2 <- vector("numeric",nrow(dat))


for (i in 7:nrow(dat)) {
  平均2[i] <- (dat[i,4] + dat[i-1,4] + dat[i-2,4] + dat[i-3,4] + dat[i-4,4] + dat[i-5,4] + dat[i-6,4])/7
}
平均2 <- as.numeric(平均2)
dat$平均2　<-  平均2

dat$Date <- as.Date(dat$Date)
as_of <- dat$Date[nrow(dat)]
datebreaks <- seq(as.Date("2020-05-01"),as.Date("2023-02-01"),by ="4 month")

head(dat,10)
tail(dat,10)
library(ggplot2)
library(stringr)
library(scales)
d2 <- ggplot(dat,aes(x=Date,y=平均2))+
  geom_line(color="grey20")+
ylim(0,500)+
labs(title = "新型コロナウイルス死亡者数",subtitle = str_c(as_of,"までの結果"),caption = "厚生労働省オープンデータ\nhttps://covid19.mhlw.go.jp/public/opendata/deaths_cumulative_daily.csv",x="日付", y="死亡者数（移動平均）") +
scale_x_date(breaks = datebreaks,labels=date_format("%Y %b"))+
theme_minimal(base_family="HiraKakuProN-W3")+
theme(axis.text.x=element_text(angle = 30,hjust = 1))
d2

新規陽性者数（移動平均）

library(readr)
pcr_positive_daily<- read_csv("https://covid19.mhlw.go.jp/public/opendata/newly_confirmed_cases_daily.csv")

# saveRDS(pcr_positive_daily,file="./data/pcr_positeve_daily.RDS")
# pcr_positive_daily <- readRDS("./data/pcr_positeve_daily.RDS")

dat <- pcr_positive_daily[,c(1,2)]
colnames(dat)[2] <- "新規陽性者数"
dat$Date <- as.Date(dat$Date)
as_of <- dat$Date[nrow(dat)]
平均 <- vector("numeric",nrow(dat))


for (i in 7:nrow(dat)) {
  平均[i] <- (dat[i,2] + dat[i-1,2] + dat[i-2,2] + dat[i-3,2] + dat[i-4,2] + dat[i-5,2] + dat[i-6,2])/7
}
平均 <- as.numeric(平均)
dat$平均　<-  平均

dat$Date <- as.Date(dat$Date)
as_of <- dat$Date[nrow(dat)]
datebreaks <- seq(as.Date("2020-01-01"),as.Date("2023-02-01"),by ="4 month")

head(dat)
tail(dat)
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.