У меня есть график линии GGPlot2
с нечетным числом граней. Я хотел бы добавить boxplot маргинального распределения значений x в пустом пространстве. BoxPlot должен быть горизонтальным и делиться общей осью X с другими гранями. Поскольку BoxPlot по умолчанию GEOM является вертикальным, требуется CORD_FLIP ()
. Из-за этого я не верю, что можно включить данные BoxPlot в том же DF, что и другие аспекты, использующие вариабельную переменную фиктивную фактору для FaceTing.
Использование GRID
, я могу идентифицировать пустой ViewPort и вставьте boxplot, но я хочу иметь линию осей X. Ответы на подобные вопросы (см. здесь или здесь здесь ) Предлагаем использование align_plots
в пакете GGEXTRA
, но я не верю, что будет Работайте с границей. Я включил простой воспроизводимый пример ниже. Если я получу эту работу, мне также придется редактировать пустую панель Grob для создания новой этикетки, соответствующего другим граням. Мы ценим любые предложения.
library(ggplot2)
#generate df for faceted line graphs
df <- data.frame(x = rep(1:100, times=7),
facet_var = as.factor(rep( 1:7, each=100)),
y = runif(7*100)
)
#create faceted line graphs
p <- ggplot(data = df, aes(x, y)) +
geom_line() + facet_wrap( ~ facet_var, ncol=2)
#generate df for boxplot
xdata <- runif(1000, min = 0, max = 100)
boxdf <- data.frame(x=xdata, group=rep(1,length(xdata)))
#create boxplot removing axes and margins
q <- ggplot(data = boxdf, aes(as.factor(group),x)) + geom_boxplot() +
coord_flip() + labs(x=NULL) +
opts(axis.text.x = theme_blank(), axis.title.x=theme_blank(),
axis.text.y = theme_blank(), axis.title.y=theme_blank(),
axis.ticks = theme_segment(colour = "white"),
panel.margin = 0, plot.margin = unit(rep(0,4), "lines")
)
print(p)
downViewport("panel-14-5")
print(q, newpage=F)
Редактировать: После полезного ответа Kohske я попытался адаптировать код для разных х пределов и разрывов. Вот код с лимитами X и разрывы изменены для диапазона (0,80). Возможно, я не хватаю что-то в коде, который должен быть изменен с ограничениями.
library(ggplot2)
df <- data.frame(x = rep(1:80, times=7),
facet_var = as.factor(rep( 1:7, each=80)),
y = runif(7*80)
)
# label for marginal plot
df <- rbind(df, data.frame(x = NA, y = NA, facet_var = "Boxplot wow"))
p <- ggplot(data = df, aes(x, y)) +
geom_line() + facet_wrap( ~ facet_var, ncol=2) +
# set limits for adjustment
coord_cartesian(xlim = c(0, 80)) +
#scale_x_continuous(breaks = 1:4*20)
opts()
xdata <- runif(1000, min = 0, max = 80)
boxdf <- data.frame(x=xdata, group=rep(1,length(xdata)))
q <- ggplot(data = boxdf, aes(as.factor(group),x)) + geom_boxplot() +
# set breaks and limits for adjustment
coord_flip(ylim = c(0, 80)) + labs(x=NULL) +
scale_y_continuous(breaks = 1:4*20) +
# opts for full region drawing:
# see https://kohske.wordpress.com/2010/12/25/drawing-on-full-region-in-ggplot2/
opts(
legend.position = "none",
panel.margin = unit(0,"null"),
plot.margin = rep(unit(0,"null"),4),
axis.ticks = theme_blank(),
axis.text.x = theme_blank(),
axis.text.y = theme_blank(),
axis.title.x = theme_blank(),
axis.title.y = theme_blank(),
axis.ticks.length = unit(0,"null"),
axis.ticks.margin = unit(0,"null")
)
print(p)
# remove unused panel
grid.remove("panel-14-5")
downViewport("panel-14-5")
print(q, newpage=F)