13

我想使用 ggplot2 绘制一些不同的数据项，使用两种不同的色标（一个连续的，一个来自两个不同 df 的离散）。我可以单独绘制我想要的任何一个，但我不能让它们一起工作。似乎您不能在同一个情节中使用两种不同的色标？我在这里和这里看到了类似的问题，这让我相信我想要实现的目标在 ggplot2 中根本不可能，但以防万一我错了，我想说明我的问题，看看是否有变通方法。

我有一些 GIS 流数据，其中附加了一些分类属性，我可以绘制（p1在下面的代码中）来获得： 在此处输入图像描述

我还有一组具有连续响应的位置，我也可以绘制（p2在下面的代码中）来获得： 在此处输入图像描述 但是我不能将两者结合起来（p3在下面的代码中）。我收到这个错误

比例错误[[prev_aes]]：尝试选择少于一个元素

注释掉该行scale_colour_hue("Strahler order") +会将错误更改为

错误：提供给连续刻度的离散值

基本上，似乎 ggplot2 对geom_path调用和geom_point调用使用相同的比例类型（连续或离散）。因此，当我将离散变量 , 传递 factor(Strahler)给scale_colour_gradientn比例时，绘图失败。

有没有解决的办法？data如果scales 函数有一个参数来告诉它应该从哪里映射或设置属性，那就太棒了。这甚至可能吗？

非常感谢和可重现的代码如下：

library(ggplot2)

### Download df's   ###
oldwd <- getwd(); tmp <- tempdir(); setwd(tmp)
url <- "http://dl.dropbox.com/u/44829974/Data.zip"
f <- paste(tmp,"\\tmp.zip",sep="")
download.file(url,f)
unzip(f)


### Read in data    ###
riv_df <- read.table("riv_df.csv", sep=",",h=T)
afr_df <- read.table("afr_df.csv", sep=",",h=T)
vil_df <- read.table("vil_df.csv", sep=",",h=T)


### Min and max for plot area   ###
xmin <- -18; xmax <- 3; ymin <- 4; ymax <- 15


### Plot river data ###
p1 <-   ggplot(riv_df, aes(long, lat)) + 
    geom_map( mapping = aes( long , lat , map_id = id ) , fill = "white" , data = afr_df , map = afr_df ) +
    geom_path( colour = "grey95" , mapping = aes( long , lat , group = group , size = 1 ) , data = afr_df ) +
    geom_path( aes( group = id , alpha = I(Strahler/6) , colour = factor(Strahler) , size = Strahler/6 ) ) +
    scale_alpha( guide = "none" ) +
    scale_colour_hue("Strahler order") +
    scale_x_continuous( limits = c( xmin , xmax ) , expand = c( 0 , 0 ) ) +
    scale_y_continuous( limits = c( ymin , ymax ) , expand = c( 0 , 0 ) ) +
    coord_map()
print(p1) # This may take a little while depending on computer speed...



### Plot response data  ###
p2 <- ggplot( NULL ) +
    geom_point( aes( X , Y , colour = Z) , size = 2 , shape = 19 , data = vil_df ) +
    scale_colour_gradientn( colours = rev(heat.colors(25)) , guide="colourbar" ) +
    coord_equal()
print(p2)



### Plot both together  ###
p3 <-   ggplot(riv_df, aes(long, lat)) + 
    geom_map( mapping = aes( long , lat , map_id = id ) , fill = "white" , data = afr_df , map = afr_df ) +
    geom_path( colour = "grey95" , mapping = aes( long , lat , group = group , size = 1 ) , data = afr_df ) +
    geom_path( aes( group = id , alpha = I(Strahler/6) , colour = factor(Strahler) , size = Strahler/6 ) ) +
    scale_colour_hue("Strahler order") +
    scale_alpha( guide = "none" ) +
    scale_x_continuous( limits = c( xmin , xmax ) , expand = c( 0 , 0 ) ) +
    scale_y_continuous( limits = c( ymin , ymax ) , expand = c( 0 , 0 ) ) +
    geom_point( aes( X , Y , colour = Z) , size = 2 , shape = 19 , data = vil_df ) +
    scale_colour_gradientn( colours = rev(heat.colors(25)) , guide="colourbar" ) +
    coord_map()
print(p3)
#Error in scales[[prev_aes]] : attempt to select less than one element

### Clear-up downloaded files   ###
unlink(tmp,recursive=T)
setwd(oldwd)

干杯，

西蒙

4

3 回答 3

13

你可以这样做。您需要告诉网格图形将一个图覆盖在另一个图之上。您必须得到完全正确的边距和间距等，并且您必须考虑顶层的透明度。总之……不值得。以及可能使情节混乱。

但是，我认为有些人可能会喜欢有关如何实现这一目标的指针。注意我使用这个 gist中的代码使顶部图中的元素透明，这样它们就不会使下面的元素不透明：

grid.newpage()
pushViewport( viewport( layout = grid.layout( 1 , 1 , widths = unit( 1 , "npc" ) ) ) ) 
print( p1 + theme(legend.position="none") , vp = viewport( layout.pos.row = 1 , layout.pos.col = 1 ) )
print( p2 + theme(legend.position="none") , vp = viewport( layout.pos.row = 1 , layout.pos.col = 1 ) )

请参阅我的答案，了解如何将图例添加到网格布局的另一个位置。

在此处输入图像描述

于 2014-06-11T15:17:04.090 回答
5

问题并不像你想象的那么复杂。通常，您只能映射一次美学。所以调用scale_colour_*两次对ggplot2没有意义。它会试图强迫一个进入另一个。

无论一个是连续的还是离散的，您都不能在同一个图表中有多个色标。包作者说他们也不打算添加这个。实现起来相当复杂，并且很容易制作令人难以置信的混乱图表。（由于类似的原因，永远不会实现多个 y 轴。）

于 2012-10-19T13:03:45.063 回答
3

我现在没有时间提供一个完整的工作示例，但还有另一种方法值得在这里提及：ggplot 中 geom_point (scale_colour_manual) 中的填充和边框颜色

基本上，结合使用 geom_point可以让您使用而不是审美shape = 21, color = NA来控制一系列点的颜色。这是我的代码的样子。我知道没有提供数据，但希望它为您提供一个起点：fillcolor

biloxi + 
  geom_point(data = filter(train, primary != 'na'), 
             aes(y = GEO_LATITUDE, x = GEO_LONGITUDE, fill = primary), 
             shape = 21, color = NA, size = 1) + 
    scale_fill_manual(values = c('dodgerblue', 'firebrick')) + 
  geom_point(data = test_map_frame, 
             aes(y = GEO_LATITUDE, x = GEO_LONGITUDE, color = var_score), 
            alpha = 1, size = 1) + 
    scale_color_gradient2(low = 'dodgerblue4', high = 'firebrick4', mid = 'white',
                    midpoint = mean(test_map_frame$var_score))

请注意对 geom_point 的每次调用如何调用不同的审美（color或fill）

于 2015-04-23T03:52:55.197 回答