R语言base包绘制按组着色、按性别标注文本的图形问题排查

核心错误原因 #

• 因子水平映射不匹配：代码初始将df$group转为因子时，未手动指定水平与标签，默认因子水平为原始取值"0"、"1"；但后续构造颜色映射因子时，错误将levels设置为c("control", "trial group")，和现有因子水平完全不匹配，最终生成的颜色向量全为NA，无法正常着色。
• 参数类型错误：base R绘图的col参数仅接受字符型颜色值向量，直接传入因子对象无法被正确解析。
• 标签可读性问题：直接将sex的原始取值0、1作为文本标注，读图时无法直接理解含义，可按需映射为易读标签。

修正后可运行代码 #

# 提前构造group因子映射，把0/1对应到分组名
df$group <- factor(df$group,
                   levels = c(0, 1),
                   labels = c("control", "trial group"))
# 直接生成字符型颜色向量，不使用因子存储颜色值
cols3 <- ifelse(df$group == "control", "black", "blue")
# 可选：把sex的0/1映射为易读标签，不需要可跳过这步
sex_labels <- ifelse(df$sex == 0, "Female", "Male") # 可按需替换0、1对应的标注文本

# 初始化空白画布
plot(x = df$pk1,
     y = df$pk5,
     xlab = "Base Headache Severity Rating",
     ylab = "1-year Headache Severity Rating",
     main = "Headache Severity Ratings",
     type = "n")

# 添加文本标注
text(x = df$pk1,
     y = df$pk5,
     labels = sex_labels, # 若要保留sex原始0/1标注，换回df$sex即可
     col = cols3)

# 添加图例
legend(x = "bottomright",
       legend = c("control", "trial group"),
       col = c("black", "blue"),
       pch = 15,
       pt.cex = 1.5)

调整说明 #

• 所有因子转换、映射逻辑在代码开头统一处理，避免前后水平不一致的问题
• 直接用条件判断生成字符型颜色向量，符合base R绘图参数的传参要求
• 预留了sex标签的映射位置，可根据实际变量定义调整0、1对应的文本
• 若无需替换sex的原始值标注，仅需把text()函数中labels参数改回df$sex即可正常运行

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