R Markdown

loading library

library(haven); library(table1); library(dplyr); library(forcats); library(ggplot2)
## 
## Attaching package: 'table1'
## The following objects are masked from 'package:base':
## 
##     units, units<-
## 
## 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
D <- read_sav ("bsNga1807.sav")
da <- D
head(da)
## # A tibble: 6 × 84
##     STT MaLK         MaBN Hovaten  Gioitinh Namsinh Nghenghiep Noicutru Ngaykham
##   <dbl> <chr>       <dbl> <chr>    <dbl+lb>   <dbl> <chr>      <chr>    <chr>   
## 1     1 VS2635353 1338043 PHAM TH… 2 [Nu]      1961 Lao dong … HCM      44621   
## 2     2 VS2629707 1332455 TRAN TH… 1 [Nam]     1953 Lao dong … HCM      13/01/2…
## 3     3 VS2632383 1335934 NGUYEN … 1 [Nam]     1952 Lao dong … HCM      44621   
## 4     4 VS2631331 1338950 TRAN TH… 2 [Nu]      1966 Lao dong … HCM      17/01/2…
## 5     6 VS2627863 1411332 TANG TI… 2 [Nu]      1957 Lao dong … HCM      17/01/2…
## 6     7 VS2632854 1334482 HO BACH… 2 [Nu]      1962 Lao dong … HCM      44652   
## # ℹ 75 more variables: Ngayketthuc <chr>, Tongsolankhamvadieutringoaitru <dbl>,
## #   Solankham <chr>, Solandieutringoaitru <dbl>, Phanloai <chr>,
## #   Thangkhamlandau <dbl>, Chandoan <chr>, MabenhBYT <chr>, MaICD10 <chr>,
## #   Chuong <chr>, Chung <chr>, Benhdanh <chr>, Thebenh <chr>, Phongchan <dbl>,
## #   Khaithau <dbl>, Yeuthong <dbl>, Toacotphong <dbl>, Hactatphong <dbl>,
## #   Mamoc <dbl>, Thatmien <dbl>, Kienty <dbl>, Chungty <dbl>,
## #   Viquanthong <dbl>, Dauthong <dbl>, Huyenvung <dbl>, Laccham <dbl>, …
da <- da %>% mutate(
  tuoi = 2022- as.numeric(Namsinh),
  tuoi = setLabel(tuoi, "Tuoi"),
  nhom.tuoi = cut(tuoi, breaks=c(-Inf, 39, 59, Inf), labels=c("< =39","40-59",">=60")),
  nhom.tuoi = setLabel(nhom.tuoi, "Nhom tuoi"),
  Gioitinh= factor(Gioitinh, levels= 1:2, labels = c("Nam", "Nu")),
  Gioitinh = setLabel(Gioitinh, "Gioi tinh"),
  Noicutru= factor(Noicutru, levels= c("HCM", "KHAC"), labels = c("TP HCM", "Tinh khac")),
  Noicutru = setLabel(Noicutru, "Noi cu tru"), 
  Nghenghiep= factor(Nghenghiep, levels= c("Lao dong tri oc", "Lao dong chan tay"), labels = c("Lao dong tri oc", "Lao dong chan tay")),
  Nghenghiep = setLabel(Nghenghiep, "Nghe nghiep"),

  Hulao= factor(Hulao, levels= c(0,1), labels = c("Khong", "Co")),
  Hulao = setLabel(Hulao, "Hu lao"),
  Khaithau= factor(Khaithau, levels= c(0,1), labels = c("Khong", "Co")),
  Khaithau = setLabel(Khaithau, "Khai thau"),
 
  Yeuthong= factor(Yeuthong, levels= c(0,1), labels = c("Khong", "Co")),
  Yeuthong = setLabel(Yeuthong, "Yeu thong"),
  Kienty= factor(Kienty, levels= c(0,1), labels = c("Khong", "Co")),
  Kienty = setLabel(Kienty, "Kien ty"), 
  Chungty= factor(Chungty, levels= c(0,1), labels = c("Khong", "Co")),
  Chungty = setLabel(Chungty, "Chung ty"), 
  Thatmien= factor(Thatmien, levels= c(0,1), labels = c("Khong", "Co")),
  Thatmien = setLabel(Thatmien, "That mien"), 
  Toacotphong= factor(Toacotphong, levels= c(0,1), labels = c("Khong", "Co")),
  Toacotphong = setLabel(Toacotphong, "Toa cot phong"), 
  Laccham= factor(Laccham, levels= c(0,1), labels = c("Khong", "Co")),
  Laccham = setLabel(Laccham, "Lac cham"),
  Hactatphong= factor(Hactatphong, levels= c(0,1), labels = c("Khong", "Co")),
  Hactatphong = setLabel(Hactatphong, "Hac tat phong"),
  Tyuyen= factor(Tyuyen, levels= c(0,1), labels = c("Khong", "Co")),
  Tyuyen = setLabel(Tyuyen, "Ty uyen")
  )

#Phuong phap dieu tri

table1(~Khaithau, data=da)
Overall
(N=805)
Khai thau
Khong 661 (82.1%)
Co 144 (17.9%)

#Hư lao

table1(~Hulao, data=da)
Overall
(N=805)
Hu lao
Khong 687 (85.3%)
Co 118 (14.7%)

#Yêu thống

table1(~Yeuthong, data=da)
Overall
(N=805)
Yeu thong
Khong 698 (86.7%)
Co 107 (13.3%)

#Kiên tý

table1(~Kienty, data=da)
Overall
(N=805)
Kien ty
Khong 738 (91.7%)
Co 67 (8.3%)

#Chứng tý

table1(~Chungty, data=da)
Overall
(N=805)
Chung ty
Khong 747 (92.8%)
Co 58 (7.2%)

#Phuong phap dieu tri ## Điện châm

table1(~Diencham, data=da)
Overall
(N=805)
Dien cham
Có 455 (56.5%)
Không 350 (43.5%)

XBBH

table1(~Xoabopbamhuyet, data=da)
Overall
(N=805)
Xoa bop bam huyet
Có 391 (48.6%)
Không 414 (51.4%)

Chườm ngải

table1(~Chuomngai, data=da)
Overall
(N=805)
Chuom ngai
Có 295 (36.6%)
Không 510 (63.4%)

Thuốc thành phẩm

table1(~Thuocthanhpham, data=da)
Overall
(N=805)
Thuoc thanh pham
Co 196 (24.3%)
Khong 609 (75.7%)

Thuốc thang

table1(~Thuocthang, data=da)
Overall
(N=805)
Thuoc thang
Co 444 (55.2%)
Khong 361 (44.8%) 

library(compareGroups)
#Moi lien quan giua nhom tuoi va  phuong phap dung thuoc
createTable(compareGroups(Khaithau~thuocyhct,data=da))
## 
## --------Summary descriptives table by 'Khai thau'---------
## 
## _______________________________________________________________________ 
##                                          Khong        Co      p.overall 
##                                          N=661       N=144              
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯ 
## Thuoc YHCT (thuoc thang+ thanh pham):                           0.551   
##     Co                                527 (79.7%) 111 (77.1%)           
##     Khong                             134 (20.3%) 33 (22.9%)            
## ¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯¯