如何将字符转为因子并按另一变量均值重排因子水平绘制ggplot箱线图

我不清楚自己哪里出错了，遇到了多个错误但无法排查。目标是绘制箱线图：y轴为girth，x轴为ForestDiv；需要将ForestDiv转换为因子，并按LatDec的均值对ForestDiv的因子水平重新排序。此前已绘制过基础箱线图，现在只需完成上述更新。以下是我尝试的代码及对应错误：

尝试1 #

代码：

Join_Macro_Site %>%
 ggplot(mapping = aes(FactorForestDiv ~ fct_reorder(girth, data = FactorForestDiv, ylim(0,200), .fun = mean(LatDec)))) + 
 geom_boxplot() +
 ylim(0,200)

错误：

Error in `f()`:
! Aesthetics must be valid data columns. Problematic aesthetic(s): y = FactorForestDiv ~ ....
Did you mistype the name of a data column or forget to add after_stat()?
Backtrace:
1. base (local) `<fn>`(x)
2. ggplot2:::print.ggplot(x)
3. ggplot2:::ggplot_build.ggplot(x)
4. ggplot2 (local) by_layer(function(l, d) l$compute_aesthetics(d, plot))
5. ggplot2 (local) f(l = layers[[i]], d = data[[i]])
6. l$compute_aesthetics(d, plot)
7. ggplot2 (local) f(..., self = self)

尝试2 #

代码：

Join_Macro_Site %>%
 ggplot(mapping = aes() + 
 geom_boxplot() +
 FactorForDiv <- factor(FactorForDiv, levels = data(order(mean(LatDec)), "ForestDiv", levels(ForestDiv)
 ylim(0,200)

错误：

Error: unexpected symbol in:
"FactorForDiv <- factor(FactorForDiv, levels = data(order(mean(LatDec)), \"ForestDiv\", levels(ForestDiv)
ylim"

尝试3 #

代码：

Join_Macro_Site %>%
 ggplot(mapping = aes() + 
 geom_boxplot() +
 Join_Macro_Site$ForestDiv <- factor(Join_Macro_Site$ForestDiv, levels = Join_Macro_Site[order(Join_Macro_Site$LatDec), "ForestDiv"])
 levels(Join_Macro_Site$ForestDiv)

错误：

Error: unexpected symbol in:
"Join_Macro_Site$ForestDiv <- factor(Join_Macro_Site$ForestDiv, levels = Join_Macro_Site[order(Join_Macro_Site$LatDec), \"ForestDiv\"])
levels"

尝试4 #

代码：

Join_Macro_Site %>%
 ggplot(mapping = aes() + 
 geom_boxplot() +
 Join_Macro_Site$ForestDiv <- factor(Join_Macro_Site$ForestDiv, levels = Join_Macro_Site[order(Join_Macro_Site$LatDec), "ForestDiv"])
 levels(Join_Macro_Site$ForestDiv) %>%

错误：

Error: unexpected symbol in:
"Join_Macro_Site$ForestDiv <- factor(Join_Macro_Site$ForestDiv, levels = Join_Macro_Site[order(Join_Macro_Site$LatDec), \"ForestDiv\"])
levels"

尝试5 #

代码：

Join_Macro_Site %>%
 ggplot(mapping = aes() + 
 geom_boxplot() +
 Join_Macro_Site$ForestDiv <- factor(Join_Macro_Site$ForestDiv, levels = Join_Macro_Site[order(Join_Macro_Site$LatDec), "ForestDiv"]), levels(Join_Macro_Site$ForestDiv) %>%
 ylim(0,200)

错误：

Error: Incomplete expression: Join_Macro_Site %>%
ggplot(mapping = aes() + 
geom_boxplot() +
Join_Macro_Site$ForestDiv <- factor(Join_Macro_Site$ForestDiv, levels = Join_Macro_Site[order(Join_Macro_Site$LatDec), \"ForestDiv\"]), levels(Join_Macro_Site$ForestDiv) %>%
ylim(0,200)

你的核心问题在于：

1. 未正确使用fct_reorder按LatDec均值排序ForestDiv的因子水平
2. 混淆ggplot2语法逻辑，不能在ggplot管道中直接修改原始数据框变量
3. 美学映射（aes）的参数传递错误

以下是两种可行实现方式：

方式1：ggplot内部动态排序因子 #

直接在aes中使用fct_reorder，指定排序依据为LatDec的均值：

library(ggplot2)
library(dplyr)
library(forcats)

Join_Macro_Site %>%
  ggplot(aes(x = fct_reorder(ForestDiv, LatDec, .fun = mean), y = girth)) +
  geom_boxplot() +
  ylim(0, 200) +
  labs(x = "ForestDiv")  # 保持x轴标签不变

方式2：先处理数据再绘图 #

先在预处理阶段将ForestDiv转换为排序后的因子，再传入ggplot：

library(ggplot2)
library(dplyr)
library(forcats)

# 预处理数据：排序ForestDiv因子水平
processed_data <- Join_Macro_Site %>%
  mutate(ForestDiv = fct_reorder(ForestDiv, LatDec, .fun = mean))

# 绘制箱线图
processed_data %>%
  ggplot(aes(x = ForestDiv, y = girth)) +
  geom_boxplot() +
  ylim(0, 200)

关键说明 #

• fct_reorder正确用法：第一个参数是要排序的因子变量，第二个是排序用的数值变量，.fun指定聚合函数（这里用mean计算LatDec均值）
• 不要在ggplot管道链中直接修改原始数据框（如Join_Macro_Site$ForestDiv <- ...），会破坏管道逻辑并引发语法错误
• 美学映射aes要明确指定x和y对应变量，避免错误使用公式写法（~）

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