使用ggproto自定义Geom函数时映射异常及多序列比对绘图报错排查

1. 先明确数据映射与stat层的输入逻辑 #

自定义stat的compute_group方法接收的是ggplot2经过分组、映射处理后的分组数据框，而非原始传入的df。你需要确保：

• 传入的aes(position, character)中，position是连续型数值（比对位点的坐标），character是离散型序列标识（如序列ID）
• 若无需额外分组，需在StatMsawind中设置required_aes，同时调用时显式指定group=1，避免自动分组打乱数据结构

2. 修正StatMsawind的compute_group核心逻辑 #

「wrong sign in 'by' argument」通常是生成区间时步长方向与起止值不匹配导致的，结合数据结构问题，以下是修正后的stat定义：

StatMsawind <- ggproto("StatMsawind", Stat,
  required_aes = c("x", "y"), # 对应映射中的position(x)和character(y)
  default_aes = aes(color = "black", size = 0.5),
  
  compute_group = function(data, scales) {
    # 先按位点坐标排序，保证数据有序
    data <- data[order(data$x), ]
    
    # 提取当前分组的唯一位点和序列ID
    unique_pos <- unique(data$x)
    unique_seqs <- unique(data$y)
    
    # 生成每个位点的左右边界（确保start < end，避免步长符号错误）
    segments <- expand.grid(
      y = unique_seqs,
      xstart = unique_pos - 0.5,
      xend = unique_pos + 0.5
    )
    
    # 生成每个序列的上下边界（将字符型序列ID转为整数计算位置）
    segments$ystart <- as.integer(factor(segments$y)) - 0.5
    segments$yend <- as.integer(factor(segments$y)) + 0.5
    
    segments
  }
)

关键修正点 #

• 明确required_aes约束映射关系，避免数据列混乱
• 强制按x（位点）排序，保证区间生成的顺序正确性
• 固定位点区间为pos±0.5，确保xstart < xend，从根源避免步长符号错误
• 将字符型序列ID转为整数型，保证y轴位置计算的连续性

3. 修正geom_msawind的几何层定义 #

确保几何层正确关联自定义stat，并选择合适的底层几何对象（如geom_segment或geom_rect）：

geom_msawind <- function(mapping = NULL, data = NULL, stat = "msawind",
                         position = "identity", na.rm = FALSE, show.legend = NA,
                         inherit.aes = TRUE, ...) {
  layer(
    geom = GeomSegment, # 若需绘制矩形块可替换为GeomRect
    stat = StatMsawind,
    data = data,
    mapping = mapping,
    position = position,
    show.legend = show.legend,
    inherit.aes = inherit.aes,
    params = list(na.rm = na.rm, ...)
  )
}

4. 调用时的验证与调整 #

运行时需指定正确映射，若需保留原始序列顺序，可手动设置y轴离散刻度：

ggplot(df, aes(x = position, y = character, group = 1)) +
  geom_msawind() +
  scale_y_discrete(limits = unique(df$character))

排查验证步骤 #

• 在compute_group开头添加print(str(data))，查看接收的数据结构，确认x、y列的存在与类型
• 单独提取compute_group内的代码，用小批量测试数据运行，验证xstart/xend/ystart/yend的生成逻辑是否正确

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