Вы можете настроить его под свои нужды, но вот основная концепция:
Мы определим новую двоичную переменную, которая определяет, заполнять ли столбики или нет. Мы передадим это в аргумент fill
в качестве фактора и используем scale_colour_manual
, чтобы дать им нужные нам цвета. Я также передал настройку параметра colour
на красный, чтобы вы могли видеть белые полосы на белом фоне.
#Sample data
df <- data.frame(
divergence = rnorm(10),
exons = paste("E", sample(1:20, 10, TRUE), sep = ""),
probset_id = sample(letters, 10, FALSE)
)
#Binary flag for fill
df$fill <- with(df, ifelse(divergence > .3, 1,0))
ggplot(data = df, aes(as.character(exons), divergence, fill = factor(fill))) +
geom_histogram(colour = "red", legend = FALSE) +
scale_fill_manual(values = c("1" = "red", "0" = "white"), legend = FALSE) +
facet_grid(~ probset_id, scales="free_x", space= "free") +
theme_bw() + opts(title="GRP78/HSPA5 (HRN vs LR)")
Будет генерировать что-то вроде этого. Я забыл set.seed()
перед сохранением изображения, чтобы ваши мили могли отличаться, но посчитайте это доказательством концепции: