循环画图实例展示
m(list = ls())
data = read.delim("E:\\Rwork\\0228.txt",skip = 3,header = T)
cancer.types = c("BLCA","KIRC","LUSC","PRAD")
pdf("E:\\Rwork\\5.12.pdf",width = 2.5*2,height = 2.5*2)
par(mfrow = c(2,2),mgp = c(2,1,1),mar=c(2,3,2,2))
for(c in cancer.types){
dat1 = data[data$Cancer.type == c,]
dat2 = read.delim(paste0("E:\\Rwork\\",c,".txt"),header = TRUE)
## select common samples between two data
isc = intersect(substr(dat1$Sample.ID,1,16),substr(dat2$SampleName,1,16))
idx = match(isc,substr(dat1$Sample.ID,1,16))
idy = match(isc,substr(dat2$SampleName,1,16))
data1 = dat1[idx,]
data2 = dat2[idy,]
data0 = cbind(data1$ESTIMATE,data2$Purity_InfiniumPurify)
rownames(data0) = isc
data.plot = na.omit(data0)
plot(data.plot[,1],data.plot[,2],pch = 20,cex = 0.5,col = "red",xlab