在ggplot2中注释刻面图

Question

我正在研究下面报道的数据集（pre.sss）

pre.sss <- pre.sss <- structure(list(Pretest.num = c(63, 62, 61, 60, 59, 58, 57, 4,2, 1), stress = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L,1L), .Label = c("[0,6]", "(6,9]"), class = "factor"), time = c(1L,1L, 1L, 1L, 1L, 1L, 1L, 8L, 8L, 8L), after = structure(c(2L,2L, 2L, 2L, 2L, 2L, 1L, 1L, NA, 1L), .Label = c("no", "yes"), class = "factor"),id = c("call_fam", "call_fam", "call_fam", "call_fam", "call_fam","call_fam", "call_fam", "counselor", "counselor", "counselor")), .Names = c("Pretest.num", "stress", "time", "after","id"), reshapeLong = structure(list(varying = structure(list(after = c("after.call.fam", "after.speak", "after.send.email","after.send.card", "after.attend", "after.fam.mtg", "after.sup.grp","after.counselor")), .Names = "after", v.names = "after", times = 1:8),v.names = "after", idvar = "Pretest.num", timevar = "time"), .Names = c("varying","v.names", "idvar", "timevar")), row.names = c("63.1", "62.1","61.1", "60.1", "59.1", "58.1", "57.1", "4.8", "2.8", "1.8"), class = "data.frame")

我需要根据另一个分类变量（'压力'）的特定级别来绘制几个分类变量的计数：所以，在我的情况下，一个多面的bbble-lot会完成工作 所以我做的是以下内容：

ylabels = c('call_fam' = "call fam.member for condolences",
            'speak' = "speak to fam.member in person", 
            'send.email' = "send condolence email to fam.member",
            'send.card' = "send condolence card/letter to fam.member", 
            'attend' = "attend funeral/wake", 
            'fam.mtg' = "provide fam.meeting",
            'sup.grp' = "suggest attending support grp.",
            'counselor' = "make referral to bereavement counselor" )

p = ggplot(pre.sss, aes(x = after, y = id)) + 
    geom_count(alpha = 0.5, col = 'darkblue') + 
    scale_size(range = c(1,30)) +
    theme(legend.position = 'none') +
    xlab("Response") + 
    ylab("What did you do after learning about death?") + 
    scale_y_discrete(labels = ylabels) +
    facet_grid(.~ pre.sss$stress, labeller = as_labeller(stress.labels))

我完全按照自己的意愿获得了以下图片。

现在，我想用每个气泡标记相应数据在数据集中出现的计数。

dat = data.frame(ggplot_build(p)$data[[1]][, c('x', 'y', 'PANEL', 'n')])
dat$PANEL = ifelse(dat$PANEL==1, "[0,6]", "(6-9]")
colnames(dat) = c('x', 'y', 'stress', 'n')

p + geom_text(aes(x, y, label = n, group = NULL), data = dat)

这给了我以下错误我真的无法理解。

> p +  geom_text(aes(x, y, label=n, group=NULL), data=dat)
Error in `$<-.data.frame`(`*tmp*`, "PANEL", value = c(1L, 1L, 1L, 1L,  : 
  replacement has 504 rows, data has 46

有人可以帮我吗？

谢谢！

EM

Answer 1

此示例中仍缺少您称为贴标机功能的功能。 geom_count使用stat_sum计算参数n，即该点的观测数量。因为您可以使用此计算参数，所以您实际上不必将绘图分配给变量并提取其数据，就像使用ggplot_build一样。

这应该做你正在寻找的事情：

ggplot(pre.sss, aes(x = after, y = id)) + 
    geom_count(alpha = 0.5, col = 'darkblue') + 
# note the following line
    stat_sum(mapping = aes(label = ..n..), geom = "text") +
    scale_size(range = c(1,30)) +
    theme(legend.position = 'none') +
    xlab("Response") + 
    ylab("What did you do after learning about death?") + 
    scale_y_discrete(labels = ylabels) +
    facet_grid(.~ stress)

我添加的行计算的内容与geom_count中幕后的内容相同，但是为其提供了一个文本geom，标签映射到该计算参数n。< / p>