1

我有一个函数，它使用spatstat'scolourmap函数为绘图生成色标。该colourmap函数使用该image函数绘制色标。这是一个例子：

library(spatstat)

set.seed(1)
dat <- data.frame(x = 1:10, y = 10:1, z = rnorm(10,6,8))

plot_this <- function(x,y,z) {
colpal <- colorRampPalette(c("red", "blue"), bias=1)(20)
colmap <- colourmap(colpal, range = range(z))
layout(matrix(c(1,2), nrow = 1), widths = c(9,1))
par(mar = c(4, 4, 2.5, 1) + 0.1)
plot(x, y, type = "p", bg = colmap(z), pch = 21, cex = 3)
par(mar=c(1,0.5,3,2))
plot(colmap, vertical = TRUE)}

plot_this(dat$x, dat$y, dat$z)

在此处输入图像描述

我想并排绘制其中几个图，但我的尝试失败了：

par(mfcol = c(1,2))
plot_this(dat$x, dat$y, dat$z)
plot_this(dat$x, dat$y, dat$z)

dev.off()
layout(matrix(c(1,2), nrow = 1))
plot_this(dat$x, dat$y, dat$z)
plot_this(dat$x, dat$y, dat$z)

dev.off()
library(gridExtra)
grid.arrange(grob(plot_this(dat$x, dat$y, dat$z)), grob(plot_this(dat$x, dat$y, dat$z)), ncol = 2)

我找到了一个胎面，Paul Murrell 博士说那image.plot是不相容的layout。不过，我相信这个问题是可以解决的。如何改进代码以使用layout或进行绘图par？

我想做类似以下的事情，但要分别包含每个图的色阶：

plot_this2 <- function(x,y,z) {
colpal <- colorRampPalette(c("red", "blue"), bias=1)(20)
colmap <- colourmap(colpal, range = range(z))
plot(x, y, type = "p", bg = colmap(z), pch = 21, cex = 3)
}

layout(matrix(c(1,2), nrow = 1))
plot_this2(dat$x, dat$y, dat$z)
plot_this2(dat$x, dat$y, dat$z)

在此处输入图像描述

编辑： @IShouldBuyABoat 建议调查gridBase. 混合基本图形和 ggplot 图形有一个很好的 SO 答案，但我找不到一个有效的grid函数组合来制作绘图。这是其中一项试验：

library(gridBase)

plot_this <- function(x,y,z) {
plot.new()
colpal <- colorRampPalette(c("red", "blue"), bias=1)(20)
colmap <- colourmap(colpal, range = range(dat$z))
gl <- grid.layout(nrow=1, ncol=2, widths = c(9,1))
vp.1 <- viewport(layout.pos.col=1, layout.pos.row=1) 
vp.2 <- viewport(layout.pos.col=2, layout.pos.row=1) 
pushViewport(viewport(layout=gl))
pushViewport(vp.1)
par(new=TRUE, fig=gridFIG(), mar = c(4, 4, 2.5, 1) + 0.1)
plot(dat$x, dat$y, type = "p", bg = colmap(dat$z), pch = 21, cex = 3)
popViewport(1)
pushViewport(vp.2)
par(new=TRUE, fig=gridFIG(), mar = c(1,0.5,3,2))  
print(plot(colmap, vertical = TRUE))
popViewport(0)}

graphics.off()
plot.new()
gl <- grid.layout(ncol = 2, nrow = 1)
vp.1 <- viewport(layout.pos.col=1, layout.pos.row=1) 
vp.2 <- viewport(layout.pos.col=2, layout.pos.row=1) 
pushViewport(viewport(layout=gl))
pushViewport(vp.1)
plot_this(dat$x, dat$y, dat$z)
popViewport()
pushViewport(vp.2)
par(new=TRUE, fig=gridFIG())
plot_this(dat$x, dat$y, dat$z)
popViewport(1)
4

1 回答 1

0

主要问题是您一次只能使用一种布局控制机制。该命令与使用等layout操作布局不兼容。par(mfrow)

一个解决方案是从函数中删除layout调用plot_this。相反，在里面plot_this，首先创建一个基本的绘图坐标系

plot(0,0,xlim=c(0,1), ylim=c(0,1), type="n", axes=FALSE)

然后使用添加散点图points(...., add=TRUE) 并使用添加颜色图plot(..., add=TRUE, xlim, ylim)，后者是对spatstat函数的调用plot.colourmap。

于 2015-12-30T11:01:16.850 回答