如何用ID列排序ggplot中Variables分面？现有代码失效求解决

问题背景 #

有一份长格式tibble数据（Data_Penelitian2），共540行，包含Zonation、Variables、Value、ID列，涉及3个因子水平与20个变量。使用ggplot绘制了Value与Zonation的小提琴图+箱线图+抖动点图，按Variables分面时默认按字母排序，希望借助ID列对Variables分面进行排序（避免手动排序），但尝试的facet_wrap代码未生效，原代码如下：

calc_stat <- function(x) {
  coef <- 1.5
  n <- sum(!is.na(x))
  # calculate quantiles
  stats <- quantile(x, probs = c(0.1, 0.25, 0.5, 0.75, 0.9))
  names(stats) <- c("ymin", "lower", "middle", "upper", "ymax")
  return(stats)
}

windowsFonts(Times=windowsFont("Times New Roman"))

p <- ggplot(Data_Penelitian2, aes(x=Zonation, y=Value)) + 
    geom_violin(trim=FALSE) + 
    stat_summary(fun.data = calc_stat, geom="boxplot",alpha = 0.2, width = 0.2 ) +  
    geom_jitter(width = 0.3, alpha = 0.5, aes(color = Zonation))+
    facet_wrap(factor(Data_Penelitian2$Variables~., levels = unique(Data_Penelitian2$ID)),scales="free", ncol=4) +
  # scale_fill_viridis(discrete = TRUE,  option = "turbo", alpha=0.7)+
   theme_minimal()+ theme(text=element_text(family="Times")) +
    theme(legend.position = "none")

解决方案 #

原代码的问题在于错误地在facet_wrap内部处理因子排序，ggplot的分面顺序由分面变量的因子水平决定，正确做法是先将Variables转换为按ID排序的因子，再进行分面：

步骤1：构建按ID排序的Variables因子水平 #

假设每个Variables对应唯一的ID，先提取两者的唯一映射并按ID排序：

# 加载dplyr（如果未加载）
library(dplyr)

# 获取按ID排序的Variables顺序
var_order <- Data_Penelitian2 %>%
  distinct(Variables, ID) %>%  # 保留每个Variables对应的唯一ID
  arrange(ID) %>%              # 按ID升序排列
  pull(Variables)              # 提取排序后的Variables列表

# 将Variables转换为因子，水平为排序后的顺序
Data_Penelitian2$Variables <- factor(Data_Penelitian2$Variables, levels = var_order)

步骤2：修改ggplot的facet_wrap代码 #

直接使用处理后的Variables变量分面即可，无需在facet_wrap内额外处理：

windowsFonts(Times=windowsFont("Times New Roman"))

p <- ggplot(Data_Penelitian2, aes(x=Zonation, y=Value)) + 
    geom_violin(trim=FALSE) + 
    stat_summary(fun.data = calc_stat, geom="boxplot", alpha = 0.2, width = 0.2 ) +  
    geom_jitter(width = 0.3, alpha = 0.5, aes(color = Zonation))+
    facet_wrap(~Variables, scales="free", ncol=4) +  # 直接使用排序后的因子变量
    # scale_fill_viridis(discrete = TRUE,  option = "turbo", alpha=0.7)+
    theme_minimal() + 
    theme(text=element_text(family="Times"), legend.position = "none")

说明 #

• 核心逻辑是通过因子水平控制分面顺序：ggplot分面时会遵循因子的levels顺序排列，因此提前将Variables的因子水平设置为按ID排序的结果，就能实现预期的分面顺序。
• 原代码中facet_wrap(factor(Data_Penelitian2$Variables~., levels = ...))的写法不符合语法，facet_wrap的第一个参数应为分面公式（如~Variables），而非直接创建因子。

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