これはggplot2 googleグループにクロスポストされています
私の状況は、私が 関数で作業中 であり、任意の数のプロットを出力することです(ユーザーが入力した入力データに応じて)。この関数はn個のプロットのリストを返します。これらのプロットを2 x 2のフォーメーションでレイアウトしたいと思います。私は次の同時問題に苦労しています:
私の現在の戦略ではgrid.arrange
gridExtra
パッケージから。これはおそらく最適ではないでしょう。特に、これが重要であるためまったく機能しません。以下に、3つのプロットを試した、コメント付きのサンプルコードを示します。
library(ggplot2)
library(gridExtra)
x <- qplot(mpg, disp, data = mtcars)
y <- qplot(hp, wt, data = mtcars)
z <- qplot(qsec, wt, data = mtcars)
# A normal, plain-jane call to grid.arrange is fine for displaying all my plots
grid.arrange(x, y, z)
# But, for my purposes, I need a 2 x 2 layout. So the command below works acceptably.
grid.arrange(x, y, z, nrow = 2, ncol = 2)
# The problem is that the function I'm developing outputs a LIST of an arbitrary
# number plots, and I'd like to be able to plot every plot in the list on a 2 x 2
# laid-out page. I can at least plot a list of plots by constructing a do.call()
# expression, below. (Note: it totally even surprises me that this do.call expression
# DOES work. I'm astounded.)
plot.list <- list(x, y, z)
do.call(grid.arrange, plot.list)
# But now I need 2 x 2 pages. No problem, right? Since do.call() is taking a list of
# arguments, I'll just add my grid.layout arguments to the list. Since grid.arrange is
# supposed to pass layout arguments along to grid.layout anyway, this should work.
args.list <- c(plot.list, "nrow = 2", "ncol = 2")
# Except that the line below is going to fail, producing an "input must be grobs!"
# error
do.call(grid.arrange, args.list)
私はやりたくないので、私は謙虚にハドリングし、私よりもはるかに賢いコミュニティの賢明なフィードバックを待ち望んでいます。これを必要以上に難しくしている場合は特に。
あなたはそこにほとんどいます!問題はそれです do.call
は、引数が名前付きlist
オブジェクトにあることを想定しています。それらをリストに入れましたが、名前付きリスト項目ではなく、文字列として。
私はこれがうまくいくと思う:
args.list <- c(plot.list, 2,2)
names(args.list) <- c("x", "y", "z", "nrow", "ncol")
ベンとジョシュアがコメントで指摘したように、リストを作成したときに名前を割り当てることもできました。
args.list <- c(plot.list,list(nrow=2,ncol=2))
または
args.list <- list(x=x, y=y, z=x, nrow=2, ncol=2)
これを試して、
require(ggplot2)
require(gridExtra)
plots <- lapply(1:11, function(.x) qplot(1:10,rnorm(10), main=paste("plot",.x)))
params <- list(nrow=2, ncol=2)
n <- with(params, nrow*ncol)
## add one page if division is not complete
pages <- length(plots) %/% n + as.logical(length(plots) %% n)
groups <- split(seq_along(plots),
gl(pages, n, length(plots)))
pl <-
lapply(names(groups), function(g)
{
do.call(arrangeGrob, c(plots[groups[[g]]], params,
list(main=paste("page", g, "of", pages))))
})
class(pl) <- c("arrangelist", "ggplot", class(pl))
print.arrangelist = function(x, ...) lapply(x, function(.x) {
if(dev.interactive()) dev.new() else grid.newpage()
grid.draw(.x)
}, ...)
## interactive use; open new devices
pl
## non-interactive use, multipage pdf
ggsave("multipage.pdf", pl)
私は少し遅れて答えていますが、R Graphics Cookbookで、multiplot
と呼ばれるカスタム関数を使用して非常に似たようなことを行うソリューションにつまずきました。おそらく、この質問を見つける他の人を助けるでしょう。解決策はこの質問に対する他の回答よりも新しい可能性があるため、私も回答を追加しています。
ここに現在の関数がありますが、上記のリンクを使用してください。作者はggplot2 0.9.3用に更新されており、再び変更される可能性があることを示しています。
# Multiple plot function
#
# ggplot objects can be passed in ..., or to plotlist (as a list of ggplot objects)
# - cols: Number of columns in layout
# - layout: A matrix specifying the layout. If present, 'cols' is ignored.
#
# If the layout is something like matrix(c(1,2,3,3), nrow=2, byrow=TRUE),
# then plot 1 will go in the upper left, 2 will go in the upper right, and
# 3 will go all the way across the bottom.
#
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))
}
}
}
プロットオブジェクトを作成します。
p1 <- ggplot(...)
p2 <- ggplot(...)
# etc.
そして、それらをmultiplot
に渡します:
multiplot(p1, p2, ..., cols = n)