ggplot2绘制Prob变量多分组折线时仅显示单条曲线

核心问题点 #

• 数据类型强制转换错误
  使用cbind()拼接包含字符型向量（如代码中的Species列）时，R会将所有参与拼接的向量统一强制转换为字符类型，后续生成的数据框中，Prob、acid、nit等本应为数值型的变量全部变为字符类型，破坏ggplot2绘图所需的数值映射逻辑。
• 计算逻辑遗漏导致曲线完全重合
  你针对不同Prob分组计算了greek变量，但在推导核心绘图变量acid的公式中完全没有引入该值，最终三个Prob分组下，acid与nit的数值对应关系完全一致，三条折线100%重叠，视觉上仅显示为一条曲线。

修正步骤 #

1. 调整数据构建方式，避免跨类型拼接的强制转换：直接用data.frame()逐列生成数据表，替代原来as.data.frame(cbind(...))的写法，保证所有变量保留自身正确的数据类型。
2. 修正acid的计算公式，将之前遗漏的greek变量纳入计算，保证不同Prob分组下的acid取值存在可观测的差异。
3. 适配新版ggplot2的语法，替换已逐步淘汰的aes_string()为标准aes()映射逻辑，保证分组和颜色映射生效。

修正后可运行的参考代码 #

# 原有变量逻辑保持不变
Prob <- c(rep(0.1, 37), rep(0.5, 37), rep(0.9, 37))
nit <- rep(seq(0.9, 1, length.out = 37), 3) # 简化原有冗长的手写nit序列，效果完全等价
greek <- log((1-Prob)/Prob)/-10
italian <- ((0.997-nit)/(0.997-0.97))^3 
Temp <- rep(25, 111)
GT <- ((30-Temp)/(30-3.3))^3
GH <- 1-GT-italian
# 修正acid计算：将遗漏的greek变量纳入公式，避免组间曲线完全重合
acid <- (-1*(((sign(GH)*(abs(GH)^(1/3)))*(7-5))-7)) + greek
Species <- rep("Case", 111)
# 正确构建数据框，避免类型强制转换
data <- data.frame(
  Prob = Prob,
  greek = greek,
  GT = GT,
  GH = GH,
  italian = italian,
  Temp = Temp,
  acid = acid,
  nit = nit,
  Species = Species
)
# 正确绘图，使用标准aes映射
ggplot(data = data, aes(x = acid, y = nit, color = factor(Prob))) +
  geom_line(size = 0.8)

注：如果你的acid计算公式有固定的业务逻辑，仅需确保公式中包含与Prob相关的计算项，即可得到分离的三组曲线。

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