2016-02-09 68 views
3

我想使用ggplot2創建散點圖矩陣。部分地,我應付了我的問題。我設法使用ggplot2創建三部分圖:下三角形 - 散點圖,對角變量名稱和上三角形相關係數(R^2)。下面我給出了一段我的數據+代碼生成我的圖表。 我的數據:R:使用ggplot2的散點圖矩陣,主題因方面面板而異

f = 
C xval V yval corr V1 
1 1622 1 1622 1  2ng 
1 1622 2 1639 0.997 2ng 
1 1622 3 1584 0.992 2ng 
1 1622 4 1549 0.99 2ng 
1 1622 5 1541 0.993 2ng 
1 1622 6 1543 0.994 2ng 
1 1622 7 1530 0.988 2ng 
2 1639 1 1622 0.997 5ng 
2 1639 2 1639 1  5ng 
2 1639 3 1584 0.997 5ng 
2 1639 4 1549 0.997 5ng 
2 1639 5 1541 0.998 5ng 
2 1639 6 1543 0.998 5ng 
2 1639 7 1530 0.995 5ng 
3 1584 1 1622 0.992 10ng 
3 1584 2 1639 0.997 10ng 
3 1584 3 1584 1  10ng 
3 1584 4 1549 0.997 10ng 
3 1584 5 1541 0.995 10ng 
3 1584 6 1543 0.999 10ng 
3 1584 7 1530 0.999 10ng 
4 1549 1 1622 0.99 15ng 
4 1549 2 1639 0.997 15ng 
4 1549 3 1584 0.997 15ng 
4 1549 4 1549 1  15ng 
4 1549 5 1541 0.998 15ng 
4 1549 6 1543 0.998 15ng 
4 1549 7 1530 0.998 15ng 
5 1541 1 1622 0.993 30ng 
5 1541 2 1639 0.998 30ng 
5 1541 3 1584 0.995 30ng 
5 1541 4 1549 0.998 30ng 
5 1541 5 1541 1  30ng 
5 1541 6 1543 0.998 30ng 
5 1541 7 1530 0.995 30ng 
6 1543 1 1622 0.994 60ng 
6 1543 2 1639 0.998 60ng 
6 1543 3 1584 0.999 60ng 
6 1543 4 1549 0.998 60ng 
6 1543 5 1541 0.998 60ng 
6 1543 6 1543 1  60ng 
6 1543 7 1530 0.998 60ng 
7 1530 1 1622 0.988 100ng 
7 1530 2 1639 0.995 100ng 
7 1530 3 1584 0.999 100ng 
7 1530 4 1549 0.998 100ng 
7 1530 5 1541 0.995 100ng 
7 1530 6 1543 0.998 100ng 
7 1530 7 1530 1  100ng 

和代碼:

g <- ggplot(data = f, aes(x=xval, y=yval))+ 
    geom_point(data = f[(xtfrm(f$C)<xtfrm(f$V)),], colour = "darkblue", size = 1.5)+ 
    geom_smooth(data = f[(xtfrm(f$C)<xtfrm(f$V)),], aes(colour = "red"), method="lm", size = 0.1)+ 
    geom_text(data = f[(xtfrm(f$C)==xtfrm(f$V)),], aes(x = 4000, y = 4000, label = paste(V1)), size = 10, colour="red")+ 
    geom_tile(aes(fill=corr))+ 
    geom_text(data = f[(xtfrm(f$C)>xtfrm(f$V)), ], aes(x = 4000, y = 4000, label = corr), size = 10)+ 
    coord_cartesian(xlim=c(0,8000), ylim=c(0,8000))+ 
    facet_grid(V~C, space = "fixed") + 
    theme(panel.grid.major = element_blank(), strip.background = element_blank(), strip.text.y = element_blank(), strip.text.x = element_blank(), legend.position = "none") 

g 

enter image description here 不過,我有一個問題,對於提高圖形的外觀。我想用不同的背景顏色分隔圖的一部分,例如:散點圖的白色背景,變量名稱的灰色和相關係數的藍色。有誰知道該怎麼做?是否有可能或我必須分別創建圖形的每個部分?

+0

我不認爲你可以修改各個方面的主題選項(如背景顏色)。你可以默認爲白色背景,然後繪製'geom_rect'來欺騙它。你也可能想看看[GGally :: ggpairs'是怎麼做的](https://github.com/ggobi/ggally/blob/master/R/ggpairs.R)。 – Gregor

回答

6

這樣做的一種方法是分別創建每個繪圖,然後將它們放在一起。這允許您爲每個繪圖單獨設置所有繪圖元素,包括背景填充顏色等主題元素。

功能下面使用mapply創建的VC每個組合一個單獨的小區,並使用if語句來分配用於每個小區所需的背景顏色。

library(gridExtra) 

p.list = mapply(FUN=function(v,c) { 

    fvc = f[f$V==v & f$C==c, ] 

    g <- ggplot(data=fvc, aes(x=xval, y=yval))+ 
     coord_cartesian(xlim=c(0,8000), ylim=c(0,8000))+ 
     theme(axis.title=element_blank(), 
      axis.text=element_text(size=8)) 

    if (c == v) { 
     g = g + geom_text(aes(x = 4000, y = 4000, label=V1), 
         size = 6, colour="red") 
    } 

    if (c < v) { 
     g = g + geom_text(aes(x = 4000, y = 4000, label=corr), size = 6) + 
     theme(panel.background=element_rect(fill=hcl(180,100,60, alpha=0.3))) 
    } 

    if (c > v) { 
     g = g + geom_point(colour = "darkblue", size = 1.5) + 
       geom_smooth(aes(colour = "red"), method="lm", size = 0.1) + 
     theme(panel.background=element_blank(), 
       panel.grid.major=element_line(colour="grey80", size=0.3)) 
    } 

    if(v != 1) { 
     g = g + theme(axis.text.y=element_blank(), 
        axis.ticks.y=element_blank()) 
    } 

    if(c != max(f$C)) { 
     g = g + theme(axis.text.x=element_blank(), 
        axis.ticks.x=element_blank()) 
    } 

    return(g) }, 
    expand.grid(V=unique(f$V), C=unique(f$C))[[1]], 
    expand.grid(V=unique(f$V), C=unique(f$C))[[2]], SIMPLIFY=FALSE) 

現在奠定了所有的地塊一起,並添加x軸和y軸標籤:

grid.arrange(
    arrangeGrob(
    arrangeGrob(textGrob("Y vals", rot=90), 
       do.call(arrangeGrob, c(p.list, ncol=7)), widths=c(0.05,0.95)), 
    textGrob("X vals"), heights=c(0.95,0.05))) 

enter image description here

正如你可以在上面的情節看,有一個尚未解決的問題:由於軸標籤的原因,第一列和最後一行中的面板具有較小的繪圖區域。 This SO answer顯示如何均衡繪圖區域的大小,但必須分別列出每個繪圖對象。

幸運的是,我們沒有必要推出我們自己的函數來平衡大網格中所有地塊面積的大小,因爲cowplot包中的plot_grid可以做到這一點。不過,我認爲在運行plot_grid後,每張圖之間的差距太大。您可以在創建圖形時通過更改繪圖邊距來調整此值。要做到這一點,你可以按如下調整在return線的繪圖功能的情節邊距:

return(g + theme(plot.margin=unit(c(0,-0.15,0,-0.15), "lines"))) }, 

然後繪製使用plot_grid

library(cowplot) 

grid.arrange(
    arrangeGrob(
    arrangeGrob(textGrob("Y vals", rot=90), 
       do.call(plot_grid, c(p.list, align="hv")), widths=c(0.05,0.95)), 
    textGrob("X vals"), heights=c(0.95,0.05))) 

enter image description here

+0

太棒了!非常感謝!! – gosia