如何存储ggplots列表以在多时隙中使用而不覆盖以前的图？

Question

我想使用从另一个函数创建的对象(cd下面的参数)在多时隙中绘制协方差/相关矩阵的一些热图.协方差矩阵存储在3维数组中,因此cd$covmat[,,i]调用第i个协方差矩阵.

最初我有一些问题与复制相同的情节.但是,我发现我遇到了环境问题.我已尝试解决这几种方法,下面的代码是最新的,但我无法弄清楚为什么它没有正确读取它.

这有什么特别的原因吗？我试过,包括和不包括环境参数(我希望不应该需要),我已经直接使用试过cd$covmat[,,i]在 aes()参数.

drawCovs<-function(cd,ncols){
    require(ggplot2)
    coords=expand.grid(x=1:cd$q,y=1:cd$q)
    climits = c(-1,1)*max(cd$covmat)
    cd$levels=c(cd$levels,"Total")
    covtext=ifelse(!(cd$use.cor),'Covariance','Correlation')
    plots=list()
    cmat=list()
    for (i in 1:(nlevels+1)){
        cmat[[i]]<-cd$covmat[,,i]
        .e<-environment
        plots[[i]]<-ggplot(environment=.e)+geom_tile(aes(x=coords$x,y=coords$y,
        fill=as.numeric(cmat[[i]]),color='white'))+
        scale_fill_gradient(covtext,low='darkblue',high='red',limits=climits)+ylab('')
        +xlab('')+guides(color='none')+scale_x_discrete(labels=cd$varnames,
        limits=1:cd$q, expand=c(0,0))+scale_y_discrete(labels=cd$varnames,
        limits=1:cd$q, expand=c(0,0))+theme(axis.text.x = element_text(angle = 90,
        hjust = 1))+labs(title=paste0(covtext,"s of data, ",cd$levels[i]))
    }

    multiplot(plotlist=plots,cols=ncols)
}

Answer 1

如果您最终尝试通过直接调用环境来修复问题,则可能会使代码过于复杂.这是一个简单的代码片段,可以作为您的函数的核心:

drawCovs <- function(cd, ncols) {
  require(ggplot2)
  require(reshape2)
  plots=list()
  cmat=list()
  for (i in 1:(length(cd$covmat))) {
    cmat[[i]] <- cd$covmat[[i]]
    plots[[i]] <- ggplot(melt(cmat), aes(x=Var1, y=Var2, fill=value)) + 
                  geom_tile(color='white')
  }  
  multiplot(plotlist=plots,cols=ncols)
}

cd <- list()
cd$covmat <- list(matrix(runif(25), 5), matrix(runif(25), 5))

drawCovs(cd, 1)