R语言：利用ggplot2和ggpubr包对不同杂交组合的产量与氮素含量进行线性拟合并可视化

假设数据框dat包含三个变量：Hybrid（杂交组合）、Yield（产量）和N（氮素含量）。

dat <- data.frame(
  Hybrid = c('A', 'A', 'B', 'B', 'C', 'C'),
  Yield = c(10, 12, 15, 17, 20, 22),
  N = c(50, 60, 70, 80, 90, 100)
)

我们可以使用ggplot2包绘制散点图和线性回归线，并使用ggpubr包添加R2和方程。

library(ggplot2)
library(ggpubr)

# 绘制散点图和线性回归线
p <- ggplot(dat, aes(x = N, y = Yield, color = Hybrid)) + 
  geom_point() + 
  geom_smooth(method = "lm", se = FALSE)

# 添加R2和方程
p2 <- stat_regline_equation(p, aes(label = paste(..eq.label.., ..r.squared.., sep = "~~~~")), 
                             label.x = "right", label.y = "top")

结果如下图所示：

图中标注的方程和R2分别为：

• Hybrid A: Yield = -0.12N + 17.4, R2 = 0.98
• Hybrid B: Yield = -0.15N + 30, R2 = 0.97
• Hybrid C: Yield = -0.16N + 39.6, R2 = 0.99

这些结果表明，不同的杂交组合之间，Yield和N之间的关系略有不同，但都表现出较强的线性关系。