R绘制COVID-19新增病例趋势图

1
2
3
4
5

rm(list = ls())
options(digits = 4)

library(readxl)
covid <- read.csv("owid-covid-data.csv")

1
2

country <- read_excel("国家和地区代码.xlsx",
                     col_names=TRUE)

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19

covid1 <- subset(covid, 
                 subset = (iso_code=="BRA"))

covid2 <- transform(covid1, 
                    peo_vac=people_vaccinated/10000, 
                    peo_fvac=people_fully_vaccinated/10000,
                    ana_new=new_cases,
                    ana_dea=new_deaths,
                    low=0
                    )

country1 <- subset(country, 
                   subset = (X6!="NA"), 
                   select = c(X2,X3,X6))

library(dplyr)
anadata1 <- left_join(covid2,
                      country1,
                      by=c("iso_code"="X6"))

1
2
3
4
5
6

mav <- function(a,n=3){
stats::filter(a,rep(1/n,n),sides = 1)
}
anadata2 <- transform(anadata1,
                      mean_new=mav(ana_new,7),
                      mean_dea=mav(ana_dea,7))

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25

library(ggplot2)
gtitle <- paste(anadata2$X2,"_",anadata2$X3,sep='')[1]
p <- ggplot(anadata2,aes(x=as.Date(date))) + 
  geom_col(aes(y=ana_new,fill="g_col")) + 
  geom_line(aes(y=mean_new,color="g_line"),size=1) +
  ggtitle(gtitle) +
  labs(x=NULL,y=NULL) +
  scale_x_date(date_label="%y/%m/%d",
               date_breaks = "3 month",
               minor_breaks = "1 month") +
  scale_fill_manual(breaks = c("g_col"), 
                    values = c("#cad5e5"), 
                    label = c("New Case")) + 
  scale_color_manual(breaks = c("g_line"),
                     values = c("blue"), 
                     label = c("Monving Average")) +
  theme(plot.title =element_text(hjust = 0.5, vjust = 0.5), 
      legend.position = "bottom", 
      legend.title = element_blank(), 
      legend.background = element_blank()) 

p + theme(panel.background=element_rect(fill='transparent', 
                                      color='gray'),
        legend.key=element_rect(fill='transparent', 
                                color='transparent'))