R图:按组标记

qaxu7uf2  于 2023-11-14  发布在  其他
关注(0)|答案(1)|浏览(94)

我正在处理的数据是一个聚类数据,在一个组中有多个观察结果,我生成了一个毛毛虫图,并希望为每个组(zipid)标记,而不是每一行,我当前的图形和代码看起来像这样:

text = hosp_new[,c("zipid")]
  ggplot(hosp_new, aes(x = id, y = oe, colour = zipid, shape = group)) +
  # theme(panel.grid.major = element_blank()) +
  geom_point(size=1) +
  scale_shape_manual(values = c(1, 2, 4)) +
  geom_errorbar(aes(ymin = low_ci, ymax = high_ci)) +
  geom_smooth(method = lm, se = FALSE) +
  scale_linetype_manual(values = linetype) +
  geom_segment(aes(x = start_id, xend = end_id, y = region_oe, yend = region_oe, linetype = "4", size = 1.2)) +
  geom_ribbon(aes(ymin = region_low_ci, ymax = region_high_ci), alpha=0.2, linetype = "blank") +
  geom_hline(aes(yintercept = 1, alpha = 0.2, colour = "red", size = 1), show.legend = "FALSE") +
  scale_size_identity() +
  scale_x_continuous(name = "hospital id", breaks = seq(0,210, by = 10)) +
  scale_y_continuous(name = "O:E ratio", breaks = seq(0,7, by = 1)) +
  geom_text(aes(label = text), position = position_stack(vjust = 10.0), size = 2)

字符串
卡特彼勒图:
x1c 0d1x的数据
每种颜色代表一个区域,我只想要一个标签/每个区域,但不知道如何删除重复的标签在这个图。任何想法?

juud5qan

juud5qan1#

关键是让geom_text只为每个zipid返回一个值,而不是多个值。如果我们希望每个zipid标签位于其组的中间,那么我们可以使用id的平均值作为每个标签的x坐标。在下面的代码中,我们使用stat_summaryh(来自ggstance包)来计算标签x坐标的平均值id,并为每个zipid返回一个标签。

library(ggplot2)
theme_set(theme_bw())
library(ggstance)

# Fake data
set.seed(300)
dat = data.frame(id=1:100, y=cumsum(rnorm(100)), 
                 zipid=rep(LETTERS[1:10], c(10, 5, 20, 8, 7, 12, 7, 10, 13,8)))

ggplot(dat, aes(id, y, colour=zipid)) +
  geom_segment(aes(xend=id, yend=0)) +
  stat_summaryh(fun.x=mean, aes(label=zipid, y=1.02*max(y)), geom="text") +
  guides(colour=FALSE)

字符串
x1c 0d1x的数据
您也可以使用faceting,正如@user20650所提到的。在下面的代码中,panel.spacing.x=unit(0,'pt')删除了facet面板之间的空间,而expand=c(0,0.5)在每个面板的侧面添加了0.5个单位的填充。这些共同确保了刻度线之间的恒定间距,即使是跨facet。

ggplot(dat, aes(id, y, colour=zipid)) +
  geom_segment(aes(xend=id, yend=0)) +
  facet_grid(. ~ zipid, scales="free_x", space="free_x") +
  guides(colour=FALSE) +
  theme_classic() +
  scale_x_continuous(breaks=0:nrow(dat), 
                     labels=c(rbind(seq(0,100,5),'','','',''))[1:(nrow(dat)+1)], 
                     expand=c(0,0.5)) +
  theme(panel.spacing.x = unit(0,"pt"))

相关问题