Ширина полосы переменной ширины

Я хотел бы создать диаграмму области/столбца в R, похожую на это: Plot from David MacKay's book "Sustainable Energy" (сюжет из книги Дэвида Маккей (отлично) "Устойчивая энергия" )

Я, честно говоря, даже не могу найти подходящее имя для такого сюжета. Кажется, это гистограмма с переменными ширинами. Определите мощный коммуникационный инструмент.

ОТВЕТЫ

Ответ 1

Вы можете сделать это с базовой графикой. Сначала мы укажем некоторые ширины и высоты:

widths = c(0.5, 0.5, 1/3,1/4,1/5, 3.5, 0.5)
heights = c(25, 10, 5,4.5,4,2,0.5)

Затем мы используем стандартную команду barplot, но указываем, что пространство между блоками равно нулю:

##Also specify colours
barplot(heights, widths, space=0, 
        col = colours()[1:6])

Поскольку мы указали ширину, нам нужно указать метки оси:

axis(1, 0:6)

Чтобы добавить линии сетки, используйте функцию grid:

##Look at ?grid to for more control over the grid lines
grid()

и вы можете добавить стрелки и текст вручную:

arrows(1, 10, 1.2, 12, code=1)
text(1.2, 13, "A country")

Чтобы добавить квадрат в верхнем правом углу, используйте функцию polygon:

polygon(c(4,4,5,5), c(20, 25, 25, 20), col="antiquewhite1")
text(4.3, 22.5, "Hi there", cex=0.6)

Все это дает:

enter image description here

Кроме того: в показанном графике я использовал команду par для настройки нескольких аспектов:

par(mar=c(3,3,2,1), 
    mgp=c(2,0.4,0), tck=-.01,
    cex.axis=0.9, las=1)

Ответ 2

Вдохновленный кодом сообщение в блоге, о котором я упоминал выше,

df <- data.frame(x = c("Alpha", "Beta", "Gamma", "Delta"), width = c(25, 50, 75, 100), height = c(100, 75, 50, 25))
df$w <- cumsum(df$width)
df$wm <- df$w - df$width
df$wt <- with(df, wm + (w - wm)/2)

library(ggplot2)
p  <- ggplot(df, aes(ymin = 0))
p1 <- p + geom_rect(aes(xmin = wm, xmax = w, ymax = height, fill = x))
library(grid) # needed for arrow function
p1 + geom_text(aes(x = wt, y = height * 0.8, label = x)) + 
     theme_bw() + labs(x = NULL, y = NULL) + 
     theme(axis.ticks = element_blank(),axis.text.x = element_blank(), 
     axis.text.y = element_blank(), legend.position = "none") + 
     annotate("text", x = 120, y = 83, label = "a Beta block") + 
     geom_segment(aes(x = 100, y = 80, xend = 80, yend = 75), 
     arrow = arrow(length = unit(0.5, "cm")))

beta blocker