ggplot position_dodge с несколькими барами ошибок
У меня есть пять точек данных для каждой перекрестной классификации двух категориальных переменных. Я пытаюсь добавить некоторый равномерный интервал между панелями ошибок, чтобы они не перекрывались в ограненном графике ggplot2, но не дали результатов. Данные очень похожи на это...
library(ggplot2)
library(dplyr)
df0 <- iris %>%
group_by(Species) %>%
mutate(long_sepal = ifelse(Sepal.Length > mean(Sepal.Length),
yes = "long", no = "short")) %>%
group_by(Species, long_sepal) %>%
mutate(petal_rank = order(Petal.Width)) %>%
filter(petal_rank <= 5)
> df0
# Source: local data frame [30 x 7]
# Groups: Species, long_sepal [6]
#
# Sepal.Length Sepal.Width Petal.Length Petal.Width Species long_sepal petal_rank
# (dbl) (dbl) (dbl) (dbl) (fctr) (chr) (int)
# 1 5.4 3.9 1.7 0.4 setosa long 1
# 2 4.6 3.4 1.4 0.3 setosa short 1
# 3 5.0 3.4 1.5 0.2 setosa short 2
# 4 4.4 2.9 1.4 0.2 setosa short 3
# 5 4.9 3.1 1.5 0.1 setosa short 4
# 6 5.4 3.7 1.5 0.2 setosa long 3
# 7 5.8 4.0 1.2 0.2 setosa long 4
# 8 5.2 4.1 1.5 0.1 setosa long 2
# 9 5.5 4.2 1.4 0.2 setosa long 5
# 10 4.5 2.3 1.3 0.3 setosa short 5
# .. ... ... ... ... ... ... ...
Код участка до сих пор.
ggplot(data = df0,
aes(x = long_sepal, y = Petal.Width,
ymin = Petal.Width-0.05,
ymax = Petal.Width+0.05)) +
geom_pointrange(position = position_dodge(width = 0.4)) +
facet_wrap(~ Species, scales = "free")
Я попытался возиться с аргументом ширины, но я получаю один и тот же график независимо от значения.
2 ответа
Решение
Я считаю, что вам нужно добавить фактор-переменную в df0
это позволит вам увернуться от наблюдений внутри long_sepal
,
# Create new grouping column by which we can dodge
df0$fac <- factor(unlist(sapply(as.vector(table(df0$long_sepal)), seq_len)))
# Assign fac to "group = " inside the aes()
... aes(x = long_sepal, y = Petal.Width, group = fac, ...
Я использовал position = position_jitter()
прежде чем избежать наложения точек, я думаю, что это может быть решением для вас.
Пытаться:
ggplot(data = df0,
aes(x = long_sepal, y = Petal.Width,
ymin = Petal.Width-0.05,
ymax = Petal.Width+0.05)) +
geom_pointrange(position = position_jitter()) +
facet_wrap(~ Species, scales = "free")