如何在R绘图中添加其他研究的对数方程对比曲线？

步骤1：转换另一研究的方程形式 #

你拿到的对比方程是log₁₀(X) = -0.4254 + 3.2857log₁₀(Y)，需要将其转换为Y（卵数）关于X（体长）的函数，才能在图上绘制曲线：

1. 移项整理：log₁₀(Y) = (log₁₀(X) + 0.4254) / 3.2857
2. 两边取10的幂：Y = 10^[(log₁₀(X) + 0.4254)/3.2857]
   也可简化为指数形式：Y = 10^(0.4254/3.2857) * X^(1/3.2857)，两种形式计算结果一致。

步骤2：在Base R中添加对比曲线 #

直接在现有绘图代码后添加curve()函数，用不同颜色/线型区分两条曲线：

x <- mydata$length
y <- mydata$EstEggs
log_model <- lm(y ~ log(x))
coefficients <- coef(log_model)
cat("Equation: y =", coefficients[1], "+", coefficients[2], "*log(x)\n")

# 绘制原始散点和本研究拟合曲线
plot(x, y, main = "Fecundity", xlab = "Length (cm)", ylab = "Number of Eggs", pch = 16)
curve(coefficients[1] + coefficients[2] * log(x), add = TRUE, col = "black", lwd = 2)

# 添加另一研究的对比曲线
curve(10^((log10(x) + 0.4254)/3.2857), add = TRUE, col = "red", lwd = 2, lty = 2)
# 或用简化指数形式，结果相同
# curve(10^(0.4254/3.2857) * x^(1/3.2857), add = TRUE, col = "red", lwd = 2, lty = 2)

# 可选：添加图例区分曲线
legend("topleft", legend = c("本研究拟合曲线", "其他研究对比曲线"), 
       col = c("black", "red"), lwd = 2, lty = c(1, 2))

可选：用ggplot2实现（更美观的可视化） #

如果偏好ggplot2，代码如下：

library(ggplot2)

# 定义对比方程的函数
compare_func <- function(x) {
  10^((log10(x) + 0.4254)/3.2857)
}

ggplot(mydata, aes(x = length, y = EstEggs)) +
  geom_point(shape = 16) +
  # 添加本研究拟合曲线
  geom_smooth(method = "lm", formula = y ~ log(x), se = FALSE, color = "black", linewidth = 1) +
  # 添加对比曲线
  stat_function(fun = compare_func, color = "red", linewidth = 1, linetype = "dashed") +
  labs(title = "Fecundity", x = "Length (cm)", y = "Number of Eggs") +
  theme_bw() +
  # 配置图例
  scale_color_manual(values = c("black", "red"), 
                     labels = c("本研究拟合曲线", "其他研究对比曲线")) +
  guides(color = guide_legend(title = "曲线类型"))

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