Добавьте метки подгрупп к графику дрожания в ggplot2
Вопрос
У меня есть почти коробочный график, похожий на jitter-plot:
dt <- rbind(se,cb,cb.se)
qplot(ds, size, data=dt, geom="jitter", colour=root, facets = test ~ .)
Я бы хотел поместить итоговую метку для каждой группы в середине графика - например, итоговый размер здесь:
aggregate(list(size=dt$size), list(dt$ds, dt$test), sum)
Group.1 Group.2 size
1 b217 se 9847
2 c10 se 97296
3 c613 se 21633
4 c7 se 207540
...
Я пробовал использовать + geom_text(aes(x=ds, y=128, label=sum(size)), size=2)
чтобы добавить метки, но я получаю одну и ту же метку в каждой позиции - как я могу получить сумму только этого раздела данных?
Редактировать: Вот где я сейчас нахожусь - может быть, я просто иду в неправильном направлении
data <- rbind(se,cb,cb.se)
labels <-ddply(data, c("ds", "test"), function(df) sum(df$size))
ggplot(data=data, aes(x=ds)) +
geom_jitter(aes(y=size, colour=root)) +
geom_text(data=labels, aes(x=ds, y=600, label=V1), size=3) +
facet_wrap(test ~ .)
Этот код не работает - я получаю undefined columns selected
ошибка...где-нибудь.Может быть, это из-за множества data=
разделы?
Решение
Поскольку вы не предоставляете выборочные данные, я продемонстрирую решение с использованием случайных данных.
set.seed(1)
n <- 100
dat <- data.frame(
ds = sample(paste("x", 1:8, sep=""), n, replace=TRUE),
size = runif(n, 0, 250),
root = sample(c(TRUE, FALSE), n, replace=TRUE),
test = sample(c("se", "cb", "cb.se"), n, replace=TRUE)
)
head(dat)
ds size root test
1 x3 163.68098 TRUE cb.se
2 x3 88.29932 TRUE se
3 x5 67.56504 FALSE cb
4 x8 248.17102 TRUE cb
5 x2 158.37332 TRUE cb
6 x8 53.30203 FALSE cb.se
p <- ggplot(dat, aes(x=ds, y=size)) +
geom_jitter(aes(colour=root)) +
facet_grid(test~.)
Создайте фрейм данных, содержащий данные метки.Обратите внимание на использование summarize
.Это говорит ddply
чтобы создать новый столбец в data.frame
labels <- ddply(dat, .(ds, test), summarize, size=round(sum(size), 0))
head(labels)
ds test size
1 x1 cb 193
2 x1 cb.se 615
3 x1 se 274
4 x2 cb 272
5 x2 cb.se 341
6 x2 se 1012
p + geom_text(aes(x=ds, label=size, y=128), data=labels, size=2)
Другие советы
Взгляните сюда.Это может быть полезно Добавление прямых меток к ggplot2 и решетчатым графикам