把之前的巩固和前进一些:
args<-commandArgs(T)
dat<-read.table(args[1],header = T, sep="\t",row.names=1)
for(i in 1:(ncol(dat)-1))
{
for(j in (i+1):ncol(dat))
{
a = as.numeric(dat[,i])
b = as.numeric(dat[,j])
c = a + b
a = a[c>0]
b = b[c>0]
cor <- round(cor(a, b, method = "pearson"), 4)
x <- sub('_[0-9]$', '', colnames(dat[i]))
y <- sub('_[0-9]$', '', colnames(dat[j]))
fp <- paste(args[2],"/", x,"_vs_", y, ".pearson.pdf", sep = "")
pdf(file=fp)
plot(a, b, log = "xy", xlab=paste("FPKM of ",x,sep=""),ylab=paste("FPKM of ",y,sep=""),col=rgb(255,0,0,max=255,alpha=100), pch=20, cex=0.5)
legend("topleft", legend = paste("Pearson correlation: ",cor,sep=""), col =rgb(255,0,0,max=255,alpha=255), pch=20)
dev.off()
}
}