Trực quan hóa biến định lượng

Phân tích dữ liệu khảo sát bằng R

Kelly McConville

Assistant Professor of Statistics

Bảng giá trị trung bình

out <- svyby(formula = ~DaysPhysHlthBad, by = ~SmokeNow, 
             design = NHANES_design, 
             FUN = svymean, na.rm = TRUE, 
             keep.names = FALSE)
out

  SmokeNow DaysPhysHlthBad        se
1       No        3.908984 0.1996290
2      Yes        4.951750 0.2346189

Biểu đồ cột

ggplot(data = out, mapping = aes(x = SmokeNow, y = DaysPhysHlthBad)) +
  geom_col() + 
  labs(y = "Số ngày sức khỏe kém\n trung bình trong tháng", 
       x = "Có hút thuốc?")

Biểu đồ cột về số ngày sức khỏe kém trung bình của người hút và không hút thuốc

Biểu đồ cột với thanh sai số

Biểu đồ cột về số ngày sức khỏe kém trung bình của người hút và không hút thuốc kèm thanh sai số

Biểu đồ cột với thanh sai số

out <- mutate(out, lower = DaysPhysHlthBad - se, 
              upper = DaysPhysHlthBad + se)
out

  SmokeNow DaysPhysHlthBad        se    lower    upper
1       No        3.908984 0.1996290 3.709355 4.108613
2      Yes        4.951750 0.2346189 4.717131 5.186369

Biểu đồ cột với thanh sai số

ggplot(data = out, mapping = aes(x = SmokeNow, y = DaysPhysHlthBad, 
                                 ymin = lower, ymax = upper)) +
  geom_col(fill = "lightblue") + geom_errorbar(width = 0.5) + 
  labs(y = "Số ngày sức khỏe kém\n trung bình trong tháng", 
       x = "Có hút thuốc?")

Biểu đồ cột về số ngày sức khỏe kém trung bình của người hút và không hút thuốc kèm thanh sai số

Biểu đồ histogram

ggplot(data = NHANESraw, mapping = aes(x = DaysPhysHlthBad, 
                                       weight = WTMEC4YR)) + 
  geom_histogram(binwidth = 1, color = "white") +
  labs(x = "Số ngày sức khỏe kém trong tháng")

Biểu đồ tần suất số ngày sức khỏe kém

Biểu đồ mật độ

NHANESraw %>%
filter(!is.na(DaysPhysHlthBad)) %>%
  mutate(WTMEC4YR_std = WTMEC4YR/sum(WTMEC4YR)) %>%
ggplot(mapping = aes(x = DaysPhysHlthBad, weight = WTMEC4YR_std)) + 
  geom_density(bw = 0.6,  fill = "lightblue") +
  labs(x = "Số ngày sức khỏe kém trong tháng")

Biểu đồ mật độ của số ngày sức khỏe kém

Biểu đồ mật độ theo ô

NHANESraw %>%
    filter(!is.na(DaysPhysHlthBad), 
           !is.na(SmokeNow)) %>%
    group_by(SmokeNow) %>%
    mutate(WTMEC4YR_std = WTMEC4YR/sum(WTMEC4YR)) %>%
    ggplot(mapping = 
           aes(x = DaysPhysHlthBad, 
           weight = WTMEC4YR_std)) + 
  geom_density(bw = 0.6, fill = "lightblue") +
  labs(x = "Số ngày sức khỏe kém trong tháng") + 
  facet_wrap(~SmokeNow, labeller = "label_both")

Biểu đồ mật độ theo ô

Biểu đồ mật độ về số ngày sức khỏe kém của người hút và không hút thuốc

Cùng luyện tập nào!

Phân tích dữ liệu khảo sát bằng R