R Markdown

Load library

library("readxl")
library("dplyr")
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library("tidyr")
library(forcats)
library(ggplot2)

read data

da.c <- read_excel("chung192021.xlsx")

data transformation

dat_longer <- da.c %>% pivot_longer(cols =  c(2:4),  names_to = "Nam",  values_to = "So.luong") %>% glimpse()
## Rows: 69
## Columns: 3
## $ Syndrome <chr> "Yêu thống", "Yêu thống", "Yêu thống", "Tý", "Tý", "Tý", "Hạc…
## $ Nam      <chr> "year.2019", "year.2020", "year.2021", "year.2019", "year.202…
## $ So.luong <dbl> 145, 130, 66, 90, 120, 30, 60, 51, 15, 53, 32, 45, 50, 36, 18…
l <- c("Yêu thống", "Tý","Lạc chẩm","Hạc tất phong","Tọa cốt phong","Thất miên","Kiên tý","Hư lao","Đầu thống","Huyễn vựng","Khái thấu","Vị quản thống", "Ma mộc","Vị nghịch","Hiếp thống","Phong chẩn","Kiên thống","Khẩu nhãn oa tà","Bán thân bất toại","Nhĩ chứng","Tỵ uyên","Kinh loạn", "Chứng lâm")



levels(dat_longer$Syndrome)
## NULL
dat_longer$Syndrome = factor(dat_longer$Syndrome, levels = l)
dat_longer$Syndrome
##  [1] Yêu thống         Yêu thống         Yêu thống         Tý               
##  [5] Tý                Tý                Hạc tất phong     Hạc tất phong    
##  [9] Hạc tất phong     Lạc chẩm          Lạc chẩm          Lạc chẩm         
## [13] Tọa cốt phong     Tọa cốt phong     Tọa cốt phong     Kiên tý          
## [17] Kiên tý           Kiên tý           Thất miên         Thất miên        
## [21] Thất miên         Hư lao            Hư lao            Hư lao           
## [25] Vị quản thống     Vị quản thống     Vị quản thống     Vị nghịch        
## [29] Vị nghịch         Vị nghịch         Đầu thống         Đầu thống        
## [33] Đầu thống         Khái thấu         Khái thấu         Khái thấu        
## [37] Ma mộc            Ma mộc            Ma mộc            Tỵ uyên          
## [41] Tỵ uyên           Tỵ uyên           Nhĩ chứng         Nhĩ chứng        
## [45] Nhĩ chứng         Khẩu nhãn oa tà   Khẩu nhãn oa tà   Khẩu nhãn oa tà  
## [49] Phong chẩn        Phong chẩn        Phong chẩn        Kiên thống       
## [53] Kiên thống        Kiên thống        Chứng lâm         Chứng lâm        
## [57] Chứng lâm         Huyễn vựng        Huyễn vựng        Huyễn vựng       
## [61] Bán thân bất toại Bán thân bất toại Bán thân bất toại Hiếp thống       
## [65] Hiếp thống        Hiếp thống        Kinh loạn         Kinh loạn        
## [69] Kinh loạn        
## 23 Levels: Yêu thống Tý Lạc chẩm Hạc tất phong Tọa cốt phong ... Chứng lâm
levels(dat_longer$Syndrome)
##  [1] "Yêu thống"         "Tý"                "Lạc chẩm"         
##  [4] "Hạc tất phong"     "Tọa cốt phong"     "Thất miên"        
##  [7] "Kiên tý"           "Hư lao"            "Đầu thống"        
## [10] "Huyễn vựng"        "Khái thấu"         "Vị quản thống"    
## [13] "Ma mộc"            "Vị nghịch"         "Hiếp thống"       
## [16] "Phong chẩn"        "Kiên thống"        "Khẩu nhãn oa tà"  
## [19] "Bán thân bất toại" "Nhĩ chứng"         "Tỵ uyên"          
## [22] "Kinh loạn"         "Chứng lâm"
dat_longer$Syndrome <- fct_rev(dat_longer$Syndrome)

Vẽ biểu đồ

# Plotting the Data in ggplot2
ggplot(dat_longer,  aes(x = Syndrome,  y = So.luong,  
                          fill = Nam,  label = So.luong)) +
  geom_bar(stat = "identity") + geom_text(
    size = 2,  position = position_stack(vjust = 0.5), colour = "black")  +
  coord_flip()+
  theme_bw()+
  ylab("Số lượng")+
  xlab("Các chứng YHCT")+
  guides(fill=guide_legend(title="Nam"))+

  
  scale_fill_discrete(labels = c("Nam 2019", "Nam 2020", "Nam 2021"))