0

So I am still an R newb so be gentle ;) I basically should be able to figure this out but I think frustration has set in. I wish to simply write three columns to a file using write.table (easy enough) but I also wish that it is appended via "\t" separator so my resulting file is a tab delimited txt file. I have attached my code which works fine but they are appended via rows so I would love to be able to solve this. Any pointers would be awesome!

Cheers

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")
    e1<-exp1[grep("9", exp1$Chromosome, invert=TRUE), ]
    write.table(e1$Probe,row.names=FALSE,col.names=FALSE,quote=FALSE,append=TRUE,outfile)
    exp2<-read.table(files[2],header=TRUE,sep="\t")
    e2<-exp2[grep("9", exp2$Chromosome, invert=TRUE), ]
    write.table(e2$Probe,row.names=FALSE,col.names=FALSE,quote=FALSE,append=TRUE,outfile)
    exp3<-read.table(files[3],header=TRUE,sep="\t")
    e3<-exp3[grep("9", exp3$Chromosome, invert=TRUE), ]
    write.table(e3$Probe,row.names=FALSE,col.names=FALSE,quote=FALSE,append=TRUE,outfile)    

}
4

2 回答 2

1

将您想要输出的所有内容放入单个data.frame中，然后write.table()与sep="\t"参数一起使用：

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)


}
于 2014-02-20T14:27:11.537 回答
0

好的，非常感谢梅森·德卡米利斯！！！这是对我有用的解决方案......它可能看起来不漂亮，但它确实有用:) 再次感谢！

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")
    e1<-exp1[grep("9", exp1$Chromosome, invert=TRUE), ]

    exp2<-read.table(files[2],header=TRUE,sep="\t")
    e2<-exp2[grep("9", exp2$Chromosome, invert=TRUE), ]

    exp3<-read.table(files[3],header=TRUE,sep="\t")
    e3<-exp3[grep("9", exp3$Chromosome, invert=TRUE), ]


    e1$id <- 1:nrow(e1)
    e1<-cbind(e1,e1$id)
    e2$id <- 1:nrow(e2)
    e2<-cbind(e2,e2$id)
    e3$id <- 1:nrow(e3)
    e3<-cbind(e3,e3$id)
    tableData <- merge(e1, e2, by="id", all=TRUE)
    tableData <- merge(tableData, e3, by="id", all=TRUE)

    ire<-data.frame(tableData$Probe.x)
    eng<-data.frame(tableData$Probe.y)
    fra<-data.frame(tableData$Probe)

    jondata<-cbind(ire,eng,fra)
    write.table(jondata,row.names=FALSE,col.names=FALSE,quote=FALSE,append=TRUE,outfile,sep="\t")
}
于 2014-02-20T16:02:36.927 回答