我想使用从另一个函数创建的对象(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)
}