我试图通过以下方式在循环中创建一些图:我想要实现的是将图表彼此下方绘制成一个单独的图。为此,我已经阅读了我可以使用par(mfrow)
,但我不确定将它放在我的循环中的哪个位置。目前,创建的单个图表被分为两个,但只绘制了上半部分。
我知道我可以通过以下方式实现:
par(mfrow=2:1)
plot(...)
plot(...)
但是使用这个解决方案我每次都需要为每个绘图调用plot命令。最后,我希望有大约20个非常狭窄的地块堆叠在一起,我不想每次都打电话给他们。我怎么能独立完成这个?
以下是我目前使用的代码:
xy <- structure(list(NAME = structure(c(2L, 2L, 1L, 1L), .Label = c("CISCO", "JOHN"), class = "factor"), ID = c(41L, 41L, 57L, 57L), X_START_YEAR = c(1965L, 1932L, 1998L, 1956L), Y_START_VALUE = c(960L, -45L, 22L, -570L), X_END_YEAR = c(1968L, 1955L, 2002L, 1970L), Y_END_VALUE = c(960L, -45L, 22L, -570L), LC = structure(c(1L, 1L, 2L, 2L), .Label = c("CA", "US"), class = "factor")), .Names = c("NAME", "ID", "X_START_YEAR","Y_START_VALUE", "X_END_YEAR", "Y_END_VALUE", "LC"), class = "data.frame", row.names = c(NA,-4L))
ind <- split(xy,xy$ID)
# Plots
for (i in ind){
xx = unlist(i[,grep('X_',colnames(i))])
yy = unlist(i[,grep('Y_',colnames(i))])
fname <- paste0(i[1, 'ID'],'.png')
png(fname, width=1679, height=1165, res=150)
par(mar=c(6,8,6,5))
par(mfrow=2:1)
plot(xx,yy,type='n', xlab=NA, ylab="Value [mm]",ylim = range(c(yy,-.5,.5)))
i <- i[,-1]
segments(i[,2],i[,3],i[,4],i[,5],lwd=2)
abline(h=0)
dev.off()
}
编辑:到目前为止,我尝试使用一种方法,在按照此处所述(R: crop multiple pngs and combine them into a single plot)创建后,在一个图中添加生成的png,但这种方法不是因为裁剪和输出分辨率问题而直截了当...