Hi all, I have build the following data frame head(href) chr tx_start tx_end g_id strand cds_start cds_end exon_count 1 chr1 8384389 8404227 NM_001080397 + 8384389 8404073 8 2 chr1 16767166 16786584 NM_001145277 + 16767256 16785491 7 3 chr1 16767166 16786584 NM_001145278 + 16767256 16785385 8 4 chr1 16767166 16786584 NM_018090 + 16767256 16785385 8 5 chr1 48998526 50489626 NM_032785 - 48999844 50489468 14 6 chr1 33546713 33585995 NM_052998 + 33547850 33585783 12
Now when I'm trying to export it into a bed file I did: output_href<- write.table(new_CTTS, file="new_href.bed", quote=FALSE, sep = "\t", na="NA", row.names=FALSE, col.names=TRUE) The thing is when I look at the file I only have: head(test) chr1 X564620 X564649 X564644 X565645 X94 X. X10 1 chr1 565369 565404 565371 566372 217 + 8 2 chr1 565463 565541 565480 566481 1214 + 15 3 chr1 565653 565697 565662 566663 1031 + 28 4 chr1 565861 565922 565883 566884 316 + 12 5 chr1 566537 566573 566564 567565 119 + 11 6 chr1 567535 567579 567562 568563 2085 + 39 but what I want is to include the corresponding gene ID as well in a last column. How can I do that? Best, Nanami [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.