Hi, I am new to using DE Seq. I have 2 files, output from htseq, containing the name of the genes with the # of reads next to them. One file has reads before induction of rna interference (uninduced) and another after (induced). So I want to compare the read count of each gene in both files and I believe I can do that with DE-Seq. When I read the manual, it only mentions inputting one file name in the command line in R. How do I input both files? This is my first time using DE Seq. Thank you1
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Just merge the two files. The easiest way for you to do that would be to open both in Excel and just copy and paste the second column from one of the files into the empty third column of the other. If you end up doing this sort of thing often then you can either write a small program to do that sort of thing or just do it in R.
Don't forget to remove the "no_feature..." and similar lines from the end.
-
Or you could run this R script to merge your HTSeq-count output files. Probably a little quicker then Excel if you have a lot of files. Anybody know how to do this with a BASH script??
Code:#!/usr/bin/R #multimergRNA.R #Ethan Ford #for use with ezDESeq.sh #Put all your output files from HTseq-count in a directory #change directories to that directory in BASH. #Run script from dirctory with your files by typing in the command line (in BASH not in R) $ Rscript path/to/this/scirpt/multimergeRNA.R #outputs a file called countstable.txt in the directory you ran the script from #The column with your gene names must have the header "gene" #To modify for other uses change: # 1)'by="gene"' to the name of the column used as a reference for merging. # 2)"countstable.txt" to the name of your desired outputfile multimerge = function(mypath){ filenames=list.files(path=mypath, full.names=TRUE) datalist = lapply(filenames, function(x){read.delim(file=x,header=T)}) Reduce(function(x,y) {merge(x,y,all=TRUE,by="gene",sort=FALSE)}, datalist) } merged <- multimerge(WD) write.table(merged, "countstable.txt", sep="\t", quote=FALSE, row.names=FALSE, na="0")
--------------
Ethan
Comment
-
@AdrianJ217: They should be sorted the same to begin with. The DESeq vignette expects a file where the first column is the gene ID (or whatever the counts were done by) and subsequent columns are other samples. So, if you just copied both columns in then things would get more complicated.
@ETHANol: I don't have a bash script, but I did write a program in C to do that (it also changes the label on each column so it's easier to remember what is what). I'm sure I can send you a copy if you'd like, it makes making processing pipelines a bit easier.
Comment
-
dpryan, if you could send me a copy of that I'd really appreciate it. My email is:
[email protected]
Thank you so much!
Comment
-
Anybody know how to do this with a BASH script??
This of course assumes that the genes are always in the same order (which they are from HTSeq) but you may get into trouble if you have, e.g., HTSeq output using different versions of GTF files (like different ENSEMBL releases), when you will lose some genes from one of the files. That's why I think it's still better to join in R, Python etc where it is easier to check that nothing gets lost.
Comment
Latest Articles
Collapse
-
by seqadmin
The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...-
Channel: Articles
05-06-2024, 07:48 AM -
-
by seqadmin
The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...-
Channel: Articles
04-22-2024, 07:01 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Yesterday, 06:35 AM
|
0 responses
14 views
0 likes
|
Last Post
by seqadmin
Yesterday, 06:35 AM
|
||
Started by seqadmin, 05-09-2024, 02:46 PM
|
0 responses
18 views
0 likes
|
Last Post
by seqadmin
05-09-2024, 02:46 PM
|
||
Started by seqadmin, 05-07-2024, 06:57 AM
|
0 responses
17 views
0 likes
|
Last Post
by seqadmin
05-07-2024, 06:57 AM
|
||
Started by seqadmin, 05-06-2024, 07:17 AM
|
0 responses
19 views
0 likes
|
Last Post
by seqadmin
05-06-2024, 07:17 AM
|
Comment