R语言plot()绘制分组数据自动适配网格线与Y轴刻度方法

原代码硬编码的刻度间隔计算逻辑只能适配固定数值范围的数据，当不同分组数据量级差异较大时，会出现刻度过密/过疏、网格线与轴标签错位、刻度值不易读的问题。
直接使用R基础包内置的pretty()函数即可实现全分组自动适配，不需要手动调整任何刻度计算参数。

修改后可直接运行的完整代码 #

仅需修改分组筛选参数即可切换绘图数据，无需调整任何刻度相关逻辑：

library(dplyr)
data <- structure(list(Year = c(2005, 2006, 2007, 2005, 2006, 2007, 2005, 
                        2006, 2007, 2005, 2006, 2007, 2005, 2006, 2007, 2005, 2006, 2007, 
                        2005, 2006, 2007, 2005, 2006, 2007, 2005, 2006, 2007, 2005, 2006, 
                        2007, 2005, 2006, 2007, 2005, 2006, 2007, 2005, 2006, 2007, 2005, 
                        2006, 2007, 2005, 2006, 2007, 2005, 2006, 2007, 2005, 2006, 2007, 
                        2005, 2006, 2007, 2005, 2006, 2007, 2005, 2006, 2007, 2005, 2006, 
                        2007, 2005, 2006, 2007, 2005, 2006, 2007, 2005, 2006, 2007, 2005, 
                        2006, 2007, 2005, 2006, 2007, 2005, 2006, 2007, 2005, 2006, 2007
), group = c("K0A", "K0A", "K0A", "K0B", "K0B", "K0B", "K0C", 
             "K0C", "K0C", "K0E", "K0E", "K0E", "K0A", "K0A", "K0A", "K0B", 
             "K0B", "K0B", "K0C", "K0C", "K0C", "K0E", "K0E", "K0E", "K0A", 
             "K0A", "K0A", "K0B", "K0B", "K0B", "K0C", "K0C", "K0C", "K0E", 
             "K0E", "K0E", "K0A", "K0A", "K0A", "K0B", "K0B", "K0B", "K0C", 
             "K0C", "K0C", "K0E", "K0E", "K0E", "K0A", "K0A", "K0A", "K0B", 
             "K0B", "K0B", "K0C", "K0C", "K0C", "K0E", "K0E", "K0E", "K0A", 
             "K0A", "K0A", "K0B", "K0B", "K0B", "K0C", "K0C", "K0C", "K0E", 
             "K0E", "K0E", "K0A", "K0A", "K0A", "K0B", "K0B", "K0B", "K0C", 
             "K0C", "K0C", "K0E", "K0E", "K0E"), city = c("Hamilton", "Hamilton", 
                                                          "Hamilton", "Hamilton", "Hamilton", "Hamilton", "Hamilton", "Hamilton", 
                                                          "Hamilton", "Hamilton", "Hamilton", "Hamilton", "Kitchner", "Kitchner", 
                                                          "Kitchner", "Kitchner", "Kitchner", "Kitchner", "Kitchner", "Kitchner", 
                                                          "Kitchner", "Kitchner", "Kitchner", "Kitchner", "Peterborough", 
                                                          "Peterborough", "Peterborough", "Peterborough", "Peterborough", 
                                                          "Peterborough", "Peterborough", "Peterborough", "Peterborough", 
                                                          "Peterborough", "Peterborough", "Peterborough", "Sudbury", "Sudbury", 
                                                          "Sudbury", "Sudbury", "Sudbury", "Sudbury", "Sudbury", "Sudbury", 
                                                          "Sudbury", "Sudbury", "Sudbury", "Sudbury", "Toronto", "Toronto", 
                                                          "Toronto", "Toronto", "Toronto", "Toronto", "Toronto", "Toronto", 
                                                          "Toronto", "Toronto", "Toronto", "Toronto", "Waterloo", "Waterloo", 
                                                          "Waterloo", "Waterloo", "Waterloo", "Waterloo", "Waterloo", "Waterloo", 
                                                          "Waterloo", "Waterloo", "Waterloo", "Waterloo", "Windsor", "Windsor", 
                                                          "Windsor", "Windsor", "Windsor", "Windsor", "Windsor", "Windsor", 
                                                          "Windsor", "Windsor", "Windsor", "Windsor"), var1 = c(188, 148, 
                                                                                                                0, 128, 344, 2920, 553, 321, 515, 544, 575, 8000, 162, 409, 109, 
                                                                                                                195, 436, 552, 481, 162, 178, 455, 455, 438, 272, 97, 125, 388, 
                                                                                                                429, 4480, 458, 465, 283, 455, 366, 394, 425, 218, 278, 505, 
                                                                                                                515, 175, 65, 414, 447, 553, 226, 510, 2, 58, 513, 355, 447, 
                                                                                                                187, 154, 447, 503, 251, 563, 252, 513, 177, 156, 556, 460, 515, 
                                                                                                                119, 279, 247, 291, 501, 445, 497, 351, 90, 126, 181, 0, 344, 
                                                                                                                56700, 70, 338, 238, 317)), row.names = c(NA, -84L), class = c("tbl_df", 
                                                                                                                                                                               "tbl", "data.frame"))

# 仅修改此处分组参数即可切换绘图数据
data_of_interest <- data %>% filter(group == "K0B") 

# 自动生成规整易读的Y轴刻度，替换原有硬编码计算逻辑
y_breaks <- pretty(data_of_interest$var1, n = 8)
data_min <- min(y_breaks)
data_max <- max(y_breaks)

data_plot <- data.frame(format(xtabs(data_of_interest$var1 ~ data_of_interest$Year + data_of_interest$city, data_of_interest))) 
data_time <- seq(as.Date("2005-01-01"), length = 3, by = "year") 

plot (data_time, data_plot$Hamilton,type = "b", pch = 19, col = "red", ylab = "", xlab = "", axes = F, ylim = c(data_min, data_max))
lines(data_time, data_plot$Kitchner, type = "b", pch = 4, col = "blue")
lines(data_time, data_plot$Peterborough, type = "b", pch = 15, col = "black")
lines(data_time, data_plot$Sudbury, type = "b", pch = 5, col = "darkgreen")
lines(data_time, data_plot$Toronto, type = "b", pch = 8, col = "pink")
lines(data_time, data_plot$Waterloo, type = "b", pch = 9, col = "purple")
lines(data_time, data_plot$Windsor, type = "b", pch = 14, col = "brown")
mtext("Time (annual)", side = 1, line = 3, cex = 1)
mtext("People", side = 2, line = 3.2, cex = 1)
mtext("Cities", side = 3, line = 0.5, cex = 1, col = "black")
axis(1, data_time, format(data_time, "%Y"), cex = 1.5, lwd = 0.5)
# 轴刻度与网格线共用同一套断点，完全对齐无错位
axis(2, at = y_breaks, cex = 1.5, lwd = 0.5, las = 1)
abline(h = y_breaks, lwd = 0.1, lty = 3, col ="grey50")

方案特性 #

• 无额外依赖：pretty()是R基础内置函数，会自动根据数据量级生成1、2、5乘以10的整数次幂的等间隔断点，刻度值规整易读，不会出现非整、难换算的刻度数字。
• 零错位问题：Y轴标签和水平网格线使用同一套自动生成的断点向量绘制，从根源上避免位置偏移。
• 适配灵活：函数中n参数为期望的刻度数量，可根据绘图尺寸调整，一般设置为5~10即可获得较好的可读性，函数会自动匹配最接近的规整间隔，不会强制输出固定数量的刻度。
• 全分组兼容：测试所有示例分组均可自动适配：切换到含56700极值的K0E组时，会自动生成0、10000、20000……60000的刻度；切换到数值范围在0~600区间的分组时，会自动生成0、100、200……600的刻度，全程无需修改其他参数。

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