求助：如何在R中按性别分组计算研究领域的加权均值

嗨Gabriele，别担心，这两个问题都有很直接的解决办法，咱们一步步来：

一、计算按性别分组的加权均值 #

你遇到的核心问题是weighted.mean需要对应每个观测的权重，而aggregate在分组时会自动把每个领域列和权重的对应子集传给函数。只要把权重作为额外参数传递给FUN就行，不用手动处理子集：

# 注意把weights作为参数传给aggregate的FUN
barplot_sex_weighted <- aggregate(
  x = df_dummies[, 1:19],  # 你的研究领域变量
  by = list(sex = df$sex),  # 按性别分组
  FUN = function(col, w) weighted.mean(col, w),  # 自定义加权均值函数
  w = df$weights  # 传递权重变量
)

原理很简单：aggregate会按sex分组，每个组里的col（单个领域列）和对应的w（该组的权重子集）会被一起传入函数，这样就能正确计算每个领域在不同性别组的加权均值了。

二、保存修改后的Boxplot #

你说手动修改Boxplot后无法保存，是因为R的基础绘图系统里，手动调整（比如在RStudio里拖曳、修改）是临时的，不会被保存到对象里。正确的做法是先通过代码完成所有修改，再保存图形：

方法1：用基础绘图系统保存 #

先打开图形设备（比如png、pdf），绘制并修改图形，最后关闭设备：

# 打开png设备，设置尺寸
png("weighted_sex_boxplot.png", width = 800, height = 600, res = 100)

# 绘制boxplot并保存对象（方便后续修改）
bp <- boxplot(your_weighted_data ~ sex, data = barplot_sex_weighted, 
              main = "加权均值：研究领域 vs 性别",
              xlab = "性别", ylab = "加权均值")

# 在这里添加你的修改，比如调整颜色、添加注释
box(col = "#2c3e50")
points(...)  # 如果需要添加点

# 关闭设备，完成保存
dev.off()

方法2：用ggplot2更灵活（推荐） #

如果用ggplot2绘制，所有修改都能保存在ggplot对象里，直接用ggsave保存：

library(ggplot2)
library(tidyr)

# 先整理数据为长格式（ggplot需要）
weighted_long <- pivot_longer(barplot_sex_weighted, 
                              cols = -sex, 
                              names_to = "field", 
                              values_to = "weighted_mean")

# 绘制boxplot并修改
bp_plot <- ggplot(weighted_long, aes(x = sex, y = weighted_mean, fill = sex)) +
  geom_boxplot() +
  labs(title = "加权均值：研究领域 vs 性别", x = "性别", y = "加权均值") +
  theme_minimal() +
  scale_fill_brewer(palette = "Set2")

# 保存图形
ggsave("weighted_sex_boxplot_ggplot.png", bp_plot, width = 8, height = 6)

这样不管你怎么修改ggplot对象，ggsave都会保存最终的效果，不用再担心手动修改无法保存的问题。

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