errorbar

Đọc dữ liệu vào R

library(haven)

## Warning: package 'haven' was built under R version 3.6.3

men = read_dta("D:/data/bp_wc_bmi/temp/men.dta")

View(men)
library(ggplot2)

## Warning: package 'ggplot2' was built under R version 3.6.3

Mã hóa các biến phân nhóm

men$a10x1 = factor(men$a10x1, levels = c(1,2),  labels = c("WC<79","WC>=79"))
men$bmi = factor(men$bmi, levels = c(1,2,3), labels = c("<19,8","19,8-20,4","20,4<"))

tạo theme đẹp, nền trắng

 the = theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(), panel.background = element_blank(), axis.line = element_line(color="black"))

Định nghĩa p2 , gồm 3 biến.

a8x1 là huyết áp tối đa a10x1 là vòng eo bmi là BMI

p2 = ggplot(men, aes(bmi,a8x1,fill=a10x1))

Tạo biểu đồ barr

p2  + stat_summary(fun = mean, geom = "bar", position = "dodge")

Vẽ biểu đồ errorbar

p2  + stat_summary(fun= mean, geom = "bar", position = "dodge") + stat_summary(fun.data = mean_cl_boot, 
               geom = "errorbar", 
               position = position_dodge(width = 0.90),
               width = 0.2) + the + xlab("BMI tertiles(kg/m^2)") + ylab("Systolic blood pressure (mmHg)") + coord_cartesian( ylim = c(100, 180))