Я знаю, что опоздал с этим, но я думаю, что нашел довольно простое решение.
Если вы посмотрите на исходный код floating.pie()
(например, позвонив getAnywhere(floating.pie)
), вы заметите, что он использует очень простой, но эффективный подход: рисование круговых сегментов в виде многоугольников. Если все, что вам нужно из ваших гистограмм, это столбцы (без меток, осей и т. Д.), Вы можете следовать тому же подходу и написать свою собственную функцию. Вот быстрая и грязная версия:
# the function
mapbars <- function (x, xllc = 0, yllc = 0, barwidth=1, maxheight=10){
# calculate how long each bar needs to be
bars <- (x/max(x)) * maxheight
# get some quick colors
col <- rainbow(length(x))
for(i in 1:length(x)){
# figure out x- and y coordinates for the corners
leftx <- xllc + ((i-1) * barwidth)
rightx <- leftx + barwidth
bottomy <- yllc
topy <- yllc + bars[i]
# draw the bar
polygon(x=c(leftx, rightx, rightx, leftx, leftx),
y=c(bottomy, bottomy, topy, topy, bottomy),
col=col[i])
}
}
x
для значений, представляемых столбцами
xllc
и yllc
укажите положение нижнего левого угла левой панели в любой системе координат, которую вы используете
barwidth
и maxheight
используются для масштабирования размера баров
Вот демоверсия с базовым sp
сюжетом. Я не думаю, что я работал plotrix
раньше, но исходя из того floating.pie
, как работает, я бы предположил, что это также должно работать plotrix
.
library(sp)
library(maptools) # just for easy access to a background map
# load some country borders as a background
data("wrld_simpl")
plot(wrld_simpl)
# zoom on a bit …
mexico <- subset(wrld_simpl, NAME=="Mexico")
plot(mexico, axes=TRUE)
# data for the bars
x1 <- c(4, 7, 1, 2)
# plot
plot(mexico, axes=TRUE)
mapbars(x=x1, xllc=-110, yllc=20, barwidth=.5, maxheight=5)
legend(x="topright", pch=22, col="black", pt.bg=rainbow(x1), legend=c("foo", "bar", "baz", "foobar"))
# add another one:
x2 <- c(9, 21, 64, 45, 33, 43, 12, 7)
mapbars(x=x2, xllc=-100, yllc=25, barwidth=.2, maxheight=2)
Результат выглядит так:
ggsubplot
пакетом, но теперь он устарел и не будет работать (как вы упомянули). Возможно, этот пост может послужить отправной точкой: stackoverflow.com/questions/36063043/…