如何在R语言中绘制帕累托前沿图（Pareto Graph）

针对你需要对比**攻击行为（有序变量）和求偶时长（连续变量）**的需求，以下是在R中绘制目标样式帕累托前沿图的具体方案：

核心思路 #

帕累托前沿的核心是识别非支配解：即不存在其他数据点，在攻击行为维度更优的同时，求偶时长维度也更优（可根据你的研究目标调整支配规则，比如攻击行为越低、求偶时长越长为更优）。我们将先筛选出这些前沿点，再用曲线连接，搭配散点展示全部数据。

具体实现步骤 #

1. 加载依赖包 #

# 首次运行时安装所需包
install.packages(c("ggplot2", "dplyr"))

# 加载包
library(ggplot2)
library(dplyr)

2. 准备数据（替换为你的真实数据） #

假设你的数据包含有序的攻击行为等级（attack_level，1-5级，1为最低，5为最高）和连续的求偶时长（courtship_time，单位：分钟）：

set.seed(123) # 保证结果可复现
df <- data.frame(
  attack_level = factor(sample(1:5, 50, replace = TRUE), ordered = TRUE),
  courtship_time = rnorm(50, mean = 15, sd = 3) %>% abs()
)

3. 筛选帕累托前沿点 #

这里假设优化目标是更高的攻击行为等级且更短的求偶时长，若目标相反，只需调整函数内的判断逻辑：

find_pareto_frontier <- function(data, ordinal_col, continuous_col, ordinal_increase = TRUE, continuous_decrease = TRUE) {
  data <- as.data.frame(data)
  ord_vals <- data[[ordinal_col]]
  cont_vals <- data[[continuous_col]]
  
  frontier_idx <- c()
  
  for (i in 1:nrow(data)) {
    current_ord <- ord_vals[i]
    current_cont <- cont_vals[i]
    dominated <- FALSE
    
    for (j in 1:nrow(data)) {
      if (i == j) next
      
      # 判断是否被其他点支配
      ord_dominant <- if (ordinal_increase) ord_vals[j] >= current_ord else ord_vals[j] <= current_ord
      cont_dominant <- if (continuous_decrease) cont_vals[j] <= current_cont else cont_vals[j] >= current_cont
      
      if (ord_dominant && cont_dominant && !(ord_vals[j] == current_ord && cont_vals[j] == current_cont)) {
        dominated <- TRUE
        break
      }
    }
    
    if (!dominated) frontier_idx <- c(frontier_idx, i)
  }
  
  # 返回按有序变量排序的前沿点
  data[frontier_idx, ] %>% arrange(!!sym(ordinal_col))
}

# 生成前沿点数据
frontier_df <- find_pareto_frontier(df, "attack_level", "courtship_time")

4. 绘制帕累托前沿图 #

ggplot() +
  # 绘制全部数据散点（灰色背景）
  geom_point(data = df, aes(x = attack_level, y = courtship_time), color = "#999999", size = 2, alpha = 0.6) +
  # 高亮帕累托前沿点
  geom_point(data = frontier_df, aes(x = attack_level, y = courtship_time), color = "#E63946", size = 3) +
  # 用平滑曲线连接前沿点
  geom_line(data = frontier_df, aes(x = as.numeric(attack_level), y = courtship_time), color = "#E63946", linewidth = 1.2) +
  # 调整坐标轴与标题
  labs(
    x = "攻击行为等级",
    y = "求偶时长（分钟）",
    title = "攻击行为与求偶时长的帕累托前沿"
  ) +
  theme_minimal() +
  theme(
    plot.title = element_text(hjust = 0.5, size = 14, face = "bold"),
    axis.title = element_text(size = 12)
  )

关键调整说明 #

• 若优化目标为更低攻击行为+更长求偶时长，调用函数时传入ordinal_increase = FALSE, continuous_decrease = FALSE即可。
• 若需要更平滑的曲线，可将geom_line替换为geom_smooth(method = "loess")（需保证前沿点数量足够支撑拟合）。

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