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
Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...-
Channel: Articles
03-22-2024, 06:39 AM -
-
by seqadmin
The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.
Avian Conservation
Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...-
Channel: Articles
03-08-2024, 10:41 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 06:37 PM
|
0 responses
7 views
0 likes
|
Last Post
by seqadmin
Today, 06:37 PM
|
||
Started by seqadmin, Today, 06:07 PM
|
0 responses
7 views
0 likes
|
Last Post
by seqadmin
Today, 06:07 PM
|
||
Started by seqadmin, 03-22-2024, 10:03 AM
|
0 responses
49 views
0 likes
|
Last Post
by seqadmin
03-22-2024, 10:03 AM
|
||
Started by seqadmin, 03-21-2024, 07:32 AM
|
0 responses
66 views
0 likes
|
Last Post
by seqadmin
03-21-2024, 07:32 AM
|
Comment