在R中使用group_by调用hist()

# 首先，确保你已经安装并加载了dplyr包
if (!requireNamespace("dplyr", quietly = TRUE)) {
  install.packages("dplyr")
}
library(dplyr)

# 假设我们有一个名为df的数据框，其中包含变量group和value
df <- data.frame(
  group = rep(c("A", "B"), each = 50),
  value = c(rnorm(50, mean = 5), rnorm(50, mean = 10))
)

# 使用group_by和do来为每个组绘制直方图
library(ggplot2) # 需要加载ggplot2包来使用autoplot
df %>%
  group_by(group) %>%
  do(autoplot(~value, data = ., bins = 10))

# 或者，如果你想保存每个组的直方图为独立的图形文件
library(gridExtra)
plots <- df %>%
  group_by(group) %>%
  do(plot = autoplot(~value, data = ., bins = 10))

# 将所有图形保存到pdf文件中
pdf("histograms.pdf")
grid.arrange(grobs = plots$plot)
dev.off()

library(ggplot2)
ggplot(df, aes(x = value, fill = group)) +
  geom_histogram(bins = 10) +
  facet_wrap(~group)