Put everything you want to output into a single data.frame
, then use write.table()
with the sep="\t"
argument:
for(num in 0:10){
input<-paste("C:\\Users\\Desktop\\non-viewpoint 20-2-14\\HMEC(p9)\\Exp*_HMEC(p9)-IP",num,".txt",sep="")
files<-Sys.glob(input)
outfile<-paste("C:\\Users\\Desktop\\non-viewpoint 20-2-14\\exps_HMEC(p9)-IP",num,".txt",sep="")
exp1<-read.table(files[1],header=TRUE,sep="\t")
exp2<-read.table(files[2],header=TRUE,sep="\t")
exp3<-read.table(files[3],header=TRUE,sep="\t")
e1<-exp1[grep("9", exp1$Chromosome, invert=TRUE), ]
e2<-exp2[grep("9", exp2$Chromosome, invert=TRUE), ]
e3<-exp3[grep("9", exp3$Chromosome, invert=TRUE), ]
tableData <- cbind(e1$Probe, e2$Probe, e3$Probe)
write.table(tableData,file=outfile,sep="\t",row.names=FALSE,col.names=FALSE,quote=FALSE,append=TRUE)
}