我昨天问了这个关于在对象中存储图的问题。我尝试实施第一种方法(意识到我没有qplot()
在原始问题中指定我正在使用)并注意到它没有按预期工作。
library(ggplot2) # add ggplot2
string = "C:/example.pdf" # Setup pdf
pdf(string,height=6,width=9)
x_range <- range(1,50) # Specify Range
# Create a list to hold the plot objects.
pltList <- list()
pltList[]
for(i in 1 : 16){
# Organise data
y = (1:50) * i * 1000 # Get y col
x = (1:50) # get x col
y = log(y) # Use natural log
# Regression
lm.0 = lm(formula = y ~ x) # make linear model
inter = summary(lm.0)$coefficients[1,1] # Get intercept
slop = summary(lm.0)$coefficients[2,1] # Get slope
# Make plot name
pltName <- paste( 'a', i, sep = '' )
# make plot object
p <- qplot(
x, y,
xlab = "Radius [km]",
ylab = "Services [log]",
xlim = x_range,
main = paste("Sample",i)
) + geom_abline(intercept = inter, slope = slop, colour = "red", size = 1)
print(p)
pltList[[pltName]] = p
}
# close the PDF file
dev.off()
在这种情况下,我使用了示例编号,因此如果它只是被复制,代码就会运行。我确实花了几个小时对此感到困惑,但我无法弄清楚出了什么问题。它可以毫无问题地编写第一组 pdf,因此我有 16 个带有正确绘图的 pdf。
然后当我使用这段代码时:
string = "C:/test_tabloid.pdf"
pdf(string, height = 11, width = 17)
grid.newpage()
pushViewport( viewport( layout = grid.layout(3, 3) ) )
vplayout <- function(x, y){viewport(layout.pos.row = x, layout.pos.col = y)}
counter = 1
# Page 1
for (i in 1:3){
for (j in 1:3){
pltName <- paste( 'a', counter, sep = '' )
print( pltList[[pltName]], vp = vplayout(i,j) )
counter = counter + 1
}
}
dev.off()
我得到的结果是abline
每个图表上的最后一个线性模型线 ( ),但数据没有改变。当我检查我的情节列表时,似乎所有这些都被最近的情节覆盖(abline
对象除外)。
一个不太重要的次要问题是如何生成一个多页 pdf,每页上有几个图,但我的代码的主要目标是将图存储在一个列表中,以便以后访问。