我一直在学习ggplot2，并希望将其用于我的所有R绘图。但是，我还没有找到一种方法来制作与传统等高线图类似的等高线图，例如可以使用grid:filled.contour()获得的等高线图。例如:
#define data
x<-seq(1,11,1)
y<-seq(1,11,1)
xyz.func<-function(x,y) {-10.4+6.53*x+6.53*y-0.167*x^2-0.167*y^2+0.0500*x*y}
#contour plot using lattice graphics and R Color Brewer
library(lattice) #for filled.contour()
library(RColorBrewer) #for brewer.pal()
z.lattice<-outer(x,y,xyz.func)
filled.contour(x,y,z.lattice,nlevels=6,col=brewer.pal(6,"YlOrRd"))
这给了我一个不错的轮廓图。
现在，让我们在ggplot2中尝试相同的操作。根据我已阅读的所有内容(尤其是Drawing labels on flat section of contour lines in ggplot2)，我能想到的最好的结果是:
#contour plot using ggplot2
library(ggplot2)
library(reshape2) #for melt()
z.molten<-melt(z.lattice)
names(z.molten) <- c("x", "y", "z")
v<-ggplot(z.molten, aes(x,y,z=z))+
geom_tile(aes(fill=z))+
stat_contour(bins=6,aes(x,y,z=z), color="black", size=0.6)+
scale_fill_gradientn(colours=brewer.pal(6,"YlOrRd"))
v
该图具有与fill.contour()相同的基本思想，但是彩色图块与轮廓线的吻合性很差。
我也没有成功更改磁贴的大小。
关于如何使ggplot2的输出更接近fill.contour()的输出的任何建议？
最佳答案
您的问题的实质似乎是如何在ggplot中生成具有离散填充轮廓的等高线图，而不是使用常规geom_tile(...)方法获得的连续等高线。这是一种方法。
x<-seq(1,11,.03) # note finer grid
y<-seq(1,11,.03)
xyz.func<-function(x,y) {-10.4+6.53*x+6.53*y-0.167*x^2-0.167*y^2+0.0500*x*y}
gg <- expand.grid(x=x,y=y)
gg$z <- with(gg,xyz.func(x,y)) # need long format for ggplot
library(ggplot2)
library(RColorBrewer) #for brewer.pal()
brks <- cut(gg$z,breaks=seq(0,100,len=6))
brks <- gsub(","," - ",brks,fixed=TRUE)
gg$brks <- gsub("\\(|\\]","",brks) # reformat guide labels
ggplot(gg,aes(x,y)) +
geom_tile(aes(fill=brks))+
scale_fill_manual("Z",values=brewer.pal(6,"YlOrRd"))+
scale_x_continuous(expand=c(0,0))+
scale_y_continuous(expand=c(0,0))+
coord_fixed()
例如，使用
scale_x_continuos(...)只是为了消除ggplot围绕轴限制放置的多余空间；对大多数事情都很好，但会分散等高线图。 coord_fixed(...)的使用只是将宽高比设置为1:1。这些是可选的。