如何做到循环画图

循环画图实例展示

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