在R中,我准备了多个ggplot2
数字,这些数字保存到这样的变量中:
flashplot <- ggplot(flash, aes(x=flash$Year, y=flash$Proc))
+ stat_smooth(method="lm", formula = "flash$Proc ~ poly(flash$Year, 2, raw=TRUE)", n=14, level=0.9)
+ geom_point(shape=19) + ylab("Percentage of total existing versions")
+ xlab("Years of Existence") + scale_x_continuous(breaks = c(0,2,4,6,8,10,12,14,16,18,20))
绘制单个图形按预期工作,但我想使用多色绘制4个图形(2乘2),如in the R cookbook所述。
当我收到错误时,我尝试了食谱中的例子,我会这样做,但它给了我同样的错误:
> multiplot(p1, p2, p3, p4, cols=2)
Error in UseMethod("getModelInfo", model) :
no applicable method for 'getModelInfo' applied to an object of class "c('gg', 'ggplot')"
我在Windows 7下使用R 2.15.2,以及ggplot2
和coefplot
的最新软件包:
> sessionInfo()
R version 2.15.2 (2012-10-26)
Platform: i386-w64-mingw32/i386 (32-bit)
locale:
[1] LC_COLLATE=German_Austria.1252 LC_CTYPE=German_Austria.1252
[3] LC_MONETARY=German_Austria.1252 LC_NUMERIC=C
[5] LC_TIME=German_Austria.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] coefplot_1.1.8 ggplot2_0.9.3
loaded via a namespace (and not attached):
[1] colorspace_1.2-0 dichromat_1.2-4 digest_0.6.0 grid_2.15.2
[5] gtable_0.1.2 labeling_0.1 MASS_7.3-22 munsell_0.4
[9] plyr_1.8 proto_0.3-9.2 RColorBrewer_1.0-5 reshape2_1.2.2
[13] scales_0.2.3 stringr_0.6.2 useful_1.1.6
任何有助于我解决问题的提示都将受到高度赞赏。我错过了什么吗?我已经在网上搜索过,但没有找到任何帮助
- 根据要求,这是示例数据:
library(ggplot2)
# This example uses the ChickWeight dataset, which comes with ggplot2
# First plot
p1 <-
ggplot(ChickWeight, aes(x=Time, y=weight, colour=Diet, group=Chick)) +
geom_line() +
ggtitle("Growth curve for individual chicks")
# Second plot
p2 <-
ggplot(ChickWeight, aes(x=Time, y=weight, colour=Diet)) +
geom_point(alpha=.3) +
geom_smooth(alpha=.2, size=1) +
ggtitle("Fitted growth curve per diet")
# Third plot
p3 <-
ggplot(subset(ChickWeight, Time==21), aes(x=weight, colour=Diet)) +
geom_density() +
ggtitle("Final weight, by diet")
# Fourth plot
p4 <-
ggplot(subset(ChickWeight, Time==21), aes(x=weight, fill=Diet)) +
geom_histogram(colour="black", binwidth=50) +
facet_grid(Diet ~ .) +
ggtitle("Final weight, by diet") +
theme(legend.position="none") # No legend (redundant in this graph)
multiplot(p1, p2, p3, p4, cols=2)
答案 0 :(得分:7)
通过使用名为grid.arrange()
的内置函数,可以更轻松地执行相同操作。您需要安装名为gridExtra
的软件包才能使用它:
grid.arrange(p1,p2,p3,p4, ncol = 2, main = "Main title")
答案 1 :(得分:1)
问题解决了。我以为我应该使用multiplot
包中的预定义coefplot
函数,但事实上我必须使用自定义函数,该函数已在R cookbook page上定义。
multiplot <- function(..., plotlist=NULL, file, cols=1, layout=NULL) {
require(grid)
# Make a list from the ... arguments and plotlist
plots <- c(list(...), plotlist)
numPlots = length(plots)
# If layout is NULL, then use 'cols' to determine layout
if (is.null(layout)) {
# Make the panel
# ncol: Number of columns of plots
# nrow: Number of rows needed, calculated from # of cols
layout <- matrix(seq(1, cols * ceiling(numPlots/cols)),
ncol = cols, nrow = ceiling(numPlots/cols))
}
if (numPlots==1) {
print(plots[[1]])
} else {
# Set up the page
grid.newpage()
pushViewport(viewport(layout = grid.layout(nrow(layout), ncol(layout))))
# Make each plot, in the correct location
for (i in 1:numPlots) {
# Get the i,j matrix positions of the regions that contain this subplot
matchidx <- as.data.frame(which(layout == i, arr.ind = TRUE))
print(plots[[i]], vp = viewport(layout.pos.row = matchidx$row,
layout.pos.col = matchidx$col))
}
}
}