# plots
pdf("plot_corr.pdf",width=7,height=3.5)
require(corrplot)
par(mfcol=c(1,2))
a=read.table("gg_d")
b=read.table("gg_wHWEd")
colnames(a)=paste("var",c("A","D","AA","AD","DD","_e"),sep="")
colnames(b)=colnames(a)
corrplot.mixed(cor(a),col=gray.colors(10))
corrplot.mixed(cor(b),col=gray.colors(10))
dev.off()