Add an Average Moving Line To Chart

Bài này minh họa các bước để thêm một line thể hiện biến thiên trung bình hay trend của một biến số theo thời gian.

Dữ liệu được sử dụng là Covid19 ở tỉnh Đồng Nai từ ngày 17/8/2021 đến 12/10/2021.

Package TTR cần thiết để tạo nên data biến thiên trung của một variable.

library(tidyverse)
library(ggthemes)
library(plyr)
library("dplyr")                        
library("readxl")
library("data.table") 
library(TTR)

df <- read_excel("E:/CovidDongNai/Donacovid21.xlsx")
df <- df %>% mutate_at(vars(date2), as.Date, format="%Y-%m-%d") #character to Date

Để vẽ biểu đồ biến thiên của số ca nhiễm covid mới (newcases) và sau đó thêm vào đường biến thiên trung bình trong 7 ngày, tôi dùng function MAS của TTR để tạo ra biến số biến thiên trung bình MAS trước.

Lưu ý là từ ngày thứ n = 7 thì MAS mới có dữ liệu.

df$MAS <- TTR::SMA(df$newcases, n = 7)
head(df$MAS, 15)

 [1]       NA       NA       NA       NA       NA       NA 686.0000 682.4286
 [9] 707.5714 745.7143 771.4286 694.1429 687.1429 640.8571 689.8571

Vẽ biểu đồ số ca nhiễm covid mới và đường biến thiên trung bình n = 7 ngày.

pl <- ggplot(df, aes(x=date, y=newcases)) +
  geom_line(color="#69b3a2", size=1) +
  ylim(0,1500) +
  geom_hline(yintercept=1000, color="orange", size=.5) +
  geom_label( 
    data=df %>% filter(newcases > 1000), # Filter data first
    aes(label=newcases)) +
    ggtitle("So Ca Nhiem Moi Covid19 o Dong Nai 
            tu 17/8/2021 den 12/10/2021") +
    ylab("So Ca Benh") + 
    xlab("Ngay")

pl + geom_line(aes(x=date, y=MAS, color = "steal blue"))+
   scale_colour_discrete(labels = "Trung bình")

Warning: Removed 6 row(s) containing missing values (geom_path).

Đường màu đỏ thể hiện sự biến thiên trung bình 7 ngày của số ca nhiễm mới covid ở Đồng Nai.