r - par(mfcol=c(2,1)) 和空白的第二行

Question

我想知道为什么下面的代码不会产生一个在另一个下面的两个图。

data(mtcars)
library(randomForest)

mtcars.rf <- randomForest(mpg ~ ., data=mtcars, ntree=1000, keep.forest=FALSE,
                          importance=TRUE)
png("rf1.png", width=6, height=6, units="in", res=100)
par(mfcol=c(2,1))
varImpPlot(mtcars.rf)
plot(mtcars.rf, log="y")
dev.off()

这只是产生



带有空白的第二行。

Answer 1

问题在于 varImpPlot 重新定义了绘图区域，然后将其重置为以前的值。这意味着它就像你调用 par(mfcol=c(2,1))后 varImpPlot线。如果您提取 varImpPlot 绘图数据，您可以自己绘制两个点图(您可以使用 layout 而不是 par 将绘图区域分成不同形状的区域):

data(mtcars)
library(randomForest)

mtcars.rf <- randomForest(mpg ~ ., data=mtcars, ntree=1000, keep.forest=FALSE,
                          importance=TRUE)
varImpData <- varImpPlot(mtcars.rf) # calculate this outside the plot

png("rf1.png", width=6, height=6, units="in", res=100)
layout(matrix(c(1,2,3,3), 2, 2, byrow = TRUE))
dotchart(varImpData[,c(1)])
dotchart(varImpData[,c(2)])
plot(mtcars.rf, log="y")
dev.off()