如何修改loess.ancova函数绘图的坐标轴与图例？

1. 修改坐标轴标签 #

loess.ancova函数启用plot=TRUE时，内部绘图逻辑不直接支持xlab/ylab参数，可通过以下两种方式修改：

方法一：绘图后用title()补充标签 #

先执行基础绘图命令，再调用title()覆盖或补充坐标轴名称：

# 生成初始绘图
loess.ancova(y=dt$response,x=dt$var1,group=dt$species,family= "gaussian",method="Speckman", plot = T)
# 自定义坐标轴标签
title(xlab = "你的X轴名称", ylab = "你的Y轴名称")

方法二：手动提取拟合结果自主绘图（灵活性更高） #

关闭plot=T，提取拟合数据后用基础绘图或ggplot2完全自定义：

# 执行分析但不生成默认绘图
fit_result <- loess.ancova(y=dt$response,x=dt$var1,group=dt$species,family= "gaussian",method="Speckman", plot = F)

# 用基础绘图实现自定义
# 绘制原始散点
plot(dt$var1, dt$response, col = dt$species, pch=16, xlab="自定义X轴名称", ylab="自定义Y轴名称")
# 添加分组拟合曲线
for(i in seq_along(fit_result$group)){
  lines(fit_result$x, fit_result$fit[,i], col=i, lwd=2)
}

2. 修改图例 #

针对默认绘图的图例修改 #

如果保留函数默认绘图，可直接用legend()函数覆盖原图例：

loess.ancova(y=dt$response,x=dt$var1,group=dt$species,family= "gaussian",method="Speckman", plot = T)
# 自定义图例：位置、标签、样式均可调整
legend("topright", legend=c("物种1", "物种2", "物种3"), col=1:3, lwd=2, title="分组名称")

针对自主绘图的图例修改 #

在方法二的自主绘图代码末尾，添加自定义图例即可，示例：

# 接方法二的绘图代码
legend("bottomleft", legend=levels(dt$species), col=seq_along(levels(dt$species)), lwd=2, pch=16, box.lty=0)

内容的提问来源于stack exchange，提问作者Gonzalo de Quesada