R语言绘制分组条形图：如何合并CSV中名称相似的thoracic相关分组

解决方法 #

你只需要在生成频数表之前，对bodysite_case做统一归类即可：所有包含thoracic关键字的分类统一重命名为同一个值，后续生成表格时会自动合并计数。

具体操作步骤 #

1. 在转换为小写的代码后，新增一行匹配归类逻辑：

# 原有代码
bodysite_case <- tolower(bodysite) 
# 新增：所有含thoracic的分类合并为同一个分组
bodysite_case[grepl("thoracic", bodysite_case)] <- "thoracic"

2. 后续生成table和绘图时，分类会自动合并，对应的横轴标签也可以调整为合并后的版本，修改后的完整绘图代码如下：

imagedata <- read.csv(file = "C:\\Users\\Owner\\Documents\\firstactivity\\data\\synthea_sample_data_csv_apr2020\\csv\\imaging_studies.csv")
bodysite <- imagedata$BODYSITE_DESCRIPTION
mod <- imagedata$MODALITY_CODE
bodysite_case <- tolower(bodysite) 
# 合并thoracic相关分类
bodysite_case[grepl("thoracic", bodysite_case)] <- "thoracic"
# 生成合并后的频数表
tab <- table(mod, bodysite_case)
# 绘图，注意names.arg现在只有9个分类了
barplot(tab, main = "Imaging Studies", ylab = "Frequency", 
        names.arg=c("Ankle","Arm","Chest","Clavicle","Knee","Pelvis","RUQ","Thoracic","Wrist"),
        col = c("red", "blue", "green", "yellow"), las=2, beside=TRUE)
title(xlab = "Body Site", line = 4, adj=0)
legend("top", fill=c("red", "blue", "green", "yellow"), c("CR", "CT", "DX", "US"))

原理说明 #

grepl("thoracic", bodysite_case)会返回布尔值，匹配到所有包含thoracic字符串的元素，将这些元素统一赋值为"thoracic"后，table()函数统计频数时会把所有符合条件的记录归为同一组，自然就实现了分类合并的效果。

内容的提问来源于stack exchange，提问作者Macromind101