Подтвердить что ты не робот

Различные легенды и цвета заливки для граненых ggplot?

Извините, что не включены никакие примеры данных для моей проблемы. Я не мог найти способ легко создать пример файла формы. Надеемся, опытные пользователи ggplot могут видеть, что Id нравится делать из приведенного ниже описания.

Я получил:

  • A data.frame X с информацией о выборках (plotid, var1, var2, var3, var4,...)

  • Файл формы полигона Y с пространственной информацией для графиков выборки

Импорт шейп файла Ymaptools) и fortify ing как data.frame Z (ggplot2) отлично работает. melt ing X до X_melted работает одинаково хорошо. merge -ing Z и X_melted до mapdf также работает.

Это означает, что теперь мы имеем data.frame в длинной форме с пространственной информацией и var1, var2, var3,...

Теперь я хочу построить этот фрейм данных следующим образом:

pl1 <- ggplot(mapdf,aes(long,lat),group=group)
pl1 <- pl1 + geom_polygon(aes(group=group,fill=value),colour="black")
pl1 <- pl1 + facet_grid(variable ~ .)
pl1 <- pl1 + coord_equal(ratio = 1)
pl1

Результат - хороший сюжет с одной панелью для каждой переменной. Карты панелей идентичны, но цвет заливки изменяется со значениями переменных. До сих пор все работает как шарм... с одной проблемой:

Переменные имеют разные значения min и max. Например, var1 идет от 0 до 5, var2 от 0 до 400, var3 от 5 до 10 и т.д. В этом примере легенда для заполнения цвет идет от 0 до 400. var2 красиво рисуется, но var1 и var3 находятся в основном в одном цвете.

Можно ли использовать другую легенду для каждой панели фасета? Или это просто не возможно (пока) возможно с facet_wrap или facet_grid в ggplot?

Я мог бы сделать отдельные графики для каждой переменной и объединить их с видовыми экранами, но там много переменных, и это будет большой работой.

Или может быть другой пакет или метод, который я мог бы использовать для выполнения того, что Id нравится делать?

И помощь будет очень оценена.:)

Изменить: С помощью описания ggplot2 -package я построил пример, который иллюстрирует мою проблему:

ids <- factor(c("1.1", "2.1", "1.2", "2.2", "1.3", "2.3"))
values <- data.frame(
id = ids,
val1 = cumsum(runif(6, max = 0.5)),
val2 = cumsum(runif(6, max = 50))
)
positions <- data.frame(
id = rep(ids, each = 4),
x = c(2, 1, 1.1, 2.2, 1, 0, 0.3, 1.1, 2.2, 1.1, 1.2, 2.5, 1.1, 0.3,
0.5, 1.2, 2.5, 1.2, 1.3, 2.7, 1.2, 0.5, 0.6, 1.3),
y = c(-0.5, 0, 1, 0.5, 0, 0.5, 1.5, 1, 0.5, 1, 2.1, 1.7, 1, 1.5,
2.2, 2.1, 1.7, 2.1, 3.2, 2.8, 2.1, 2.2, 3.3, 3.2)
)

values <- melt(values)
datapoly <- merge(values, positions, by=c("id"))

p <- ggplot(datapoly, aes(x=x, y=y)) + geom_polygon(aes(fill=value, group=id),colour="black")
p <- p + facet_wrap(~ variable)
p

Панель справа показывает различные значения для var2 на карте. На панели слева, однако, все полигоны имеют один и тот же цвет. Это логично, потому что для всех панелей используется только один цветовой градиент. Могу ли я использовать другой цветовой градиент для каждой панели?

4b9b3361

Ответ 1

В настоящее время для каждого графика может быть только одна шкала (для всего, кроме x и y).

Ответ 2

С поддержкой сетки

align.plots <- function(..., vertical=TRUE){
#http://ggextra.googlecode.com/svn/trunk/R/align.r
  dots <- list(...)
  dots <- lapply(dots, ggplotGrob)
  ytitles <- lapply(dots, function(.g) editGrob(getGrob(.g,"axis.title.y.text",grep=TRUE), vp=NULL))
  ylabels <- lapply(dots, function(.g) editGrob(getGrob(.g,"axis.text.y.text",grep=TRUE), vp=NULL))
  legends <- lapply(dots, function(.g) if(!is.null(.g$children$legends))
                    editGrob(.g$children$legends, vp=NULL) else ggplot2:::.zeroGrob)

  gl <- grid.layout(nrow=length(dots))
  vp <- viewport(layout=gl)
  pushViewport(vp)
  widths.left <- mapply(`+`, e1=lapply(ytitles, grobWidth),
                        e2= lapply(ylabels, grobWidth), SIMPLIFY=F)
  widths.right <- lapply(legends, function(g) grobWidth(g) + if(is.zero(g)) unit(0, "lines") else unit(0.5, "lines")) # safe margin recently added to ggplot2
  widths.left.max <- max(do.call(unit.c, widths.left))
  widths.right.max <- max(do.call(unit.c, widths.right))

  for(ii in seq_along(dots)){
    pushViewport(viewport(layout.pos.row=ii))
    pushViewport(viewport(x=unit(0, "npc") + widths.left.max - widths.left[[ii]],
                          width=unit(1, "npc") - widths.left.max + widths.left[[ii]] -
                                                 widths.right.max + widths.right[[ii]],
                          just="left"))
    grid.draw(dots[[ii]])
  upViewport(2)
  }
}



p <- ggplot(datapoly[datapoly$variable=="val1",], aes(x=x, y=y)) + geom_polygon(aes(fill=value, group=id),colour="black")
p1 <- ggplot(datapoly[datapoly$variable=="val2",], aes(x=x, y=y)) + geom_polygon(aes(fill=value, group=id),colour="black")
align.plots( p,p1)

Ответ 3

Возможно, немного неортодоксально, но вы можете попробовать факторировать свою "ценность". Например:

p <- ggplot(datapoly, aes(x=x, y=y)) + geom_polygon(aes(fill=factor(value), group=id),colour="black")
p <- p + facet_wrap(~ variable)
p

ggplot2 использует факторы для создания легенд. Поэтому, если вы можете добавить столбец, который принимает "значение" и разбивает его на факторизованные диапазоны, вы можете заменить "значение" на диапазоны.

Создайте столбец, например "f":

    id variable        value   x    y f
1  1.1     val1   0.09838607 2.0 -0.5 0.09-0.13
2  1.1     val1   0.09838607 1.0  0.0 0.09-0.13
3  1.1     val1   0.09838607 1.1  1.0 0.09-0.13
4  1.1     val1   0.09838607 2.2  0.5 0.09-0.13
25 2.1     val1   0.13121347 1.0  0.0 0.13-0.20

...

Затем используйте:

p <- ggplot(datapoly, aes(x=x, y=y)) + geom_polygon(aes(fill=f, group=id),colour="black")
p <- p + facet_wrap(~ variable)
p

Вам нужно будет указать категории, которые вы хотите, что может потребовать много времени. Но, по крайней мере, график выйдет, как вы этого хотите. В основном, вы перекодируете данные в другой столбец. Вот несколько примеров:

http://www.statmethods.net/management/variables.html

Ответ 4

Рискуя указать очевидное, кажется, что вы должны раскрашивать проценты вместо исходных значений. Затем ваши измененные значения и ваша легенда будут идти от 0 до 1.