Mittelpunkt diskreter divergierender Skala in ggplot2

Question

Ich arbeite mit einigen Gitterdaten und habe Probleme mit diskreten divergierenden Skalen zu arbeiten. Insbesondere, wie Sie den Mittelpunkt so einstellen, dass er nicht in der Mitte des Bereichs liegt. Dies ist ein reproduzierbares Beispiel zu bekommen, was ich meine:

library(ggplot2) 
grid <- expand.grid(lon = seq(0, 360, by = 2), lat = seq(-90, 0, by = 2)) 
grid$z <- with(grid, cos(lat*pi/180) - .7) 

ggplot(grid, aes(lon, lat)) + 
    geom_raster(aes(fill = cut_width(z, .1))) + 
    scale_fill_brewer(palette = "RdBu") 

Hier ist die Mitte der Skala ist nicht die Kluft zwischen positiven und negativen Werten. Ich weiß, dass ich eine kontinuierliche Skala verwenden könnte, aber ich finde, dass weniger Farben helfen, was ich zu zeigen versuche.

Gibt es eine Möglichkeit, den Mittelpunkt in einer diskreten Skala zu verschieben? Andere Alternativen, die das gleiche Ergebnis erzielen, sind ebenfalls willkommen.

Answer 1

Das Problem ist, dass Ihre Schnittpunkte nicht symmetrisch um 0 fallen und sich direkt auf Ihre Farben beziehen. Ein Ansatz ist die manuelle Ihre Schnittpunkte so eingestellt, dass sie drehen sich um 0. Dann, nur um sicher zu machen, um nicht ungenutzt Ebenen in der Legende zu löschen:

zCuts <- 
    seq(-.7, 0.7, length.out = 10) 

ggplot(grid, aes(lon, lat)) + 
    geom_raster(aes(fill = cut(z, zCuts))) + 
    scale_fill_brewer(palette = "RdBu" 
        , drop = FALSE) 

Wenn Sie bereit sind, mit einem gehen stattdessen Gradienten solcher diskreten Farben, können Sie scale_fill_gradient2 die standardmäßig Zentren bei 0 und liegt zwischen zwei Farben verwenden:

ggplot(grid, aes(lon, lat)) + 
    geom_raster(aes(fill = z)) + 
    scale_fill_gradient2() 

Oder, wenn Sie wirklich die Interpolation von Color Brewer wollen, können Sie das limits Argument in scale_fill_distiller setzen und stattdessen einen Gradienten bekommen. Hier setze ich sie auf + und - der Bereich um 0 (max(abs(grid$z)) erhält die größte Abweichung von 0, ob es die min oder die max ist, um sicherzustellen, dass der Bereich symmetrisch ist). Wenn Sie mehr als die 11 verfügbaren Werte verwenden, ist, dass wahrscheinlich der beste Weg zu gehen:

ggplot(grid, aes(lon, lat)) + 
    geom_raster(aes(fill = z)) + 
    scale_fill_distiller(palette = "RdBu" 
         , limits = c(-1,1)*max(abs(grid$z)) 
         ) 

Wenn Sie mehr Farben wünschen, ohne einen Gradienten zu tun, sind Sie wahrscheinlich zu müssen gehen konstruiere deine eigene Palette manuell mit mehr Farben. Je mehr Sie hinzufügen, desto weniger unterscheiden Sie zwischen den Farben. Hier sehen Sie ein Beispiel, bei dem zwei Paletten zusammengeführt werden, um sicherzustellen, dass Sie mit unterschiedlichen Farben arbeiten.

zCuts <- 
    seq(-.7, 0.7, length.out = 20) 

myPallette <- 
    c(rev(brewer.pal(9, "YlOrRd")) 
    , "white" 
    , brewer.pal(9, "Blues")) 


ggplot(grid, aes(lon, lat)) + 
    geom_raster(aes(fill = cut(z, zCuts))) + 
    scale_fill_manual(values = myPallette 
        , drop = FALSE)