Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • DE Seq

    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

  • #2
    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.

    Comment


    • #3
      Thank you

      Great, thank you so much, and then can you tell me what commands I need to type into R in order to run it using my file?

      Comment


      • #4
        Save it as a CSV file and then follow the DESeq vignette, which is very well written.

        Comment


        • #5
          dpryan, thank you so much. Just one more question. You said to copy the 2nd column of one to the empty 3rd column of the other, but why wouldnt I copy both columns of one to the other, to have 4 columns. Won't they mismatch if I just copy and paste the 2nd column?

          Comment


          • #6
            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


            • #7
              @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


              • #8
                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


                • #9
                  dpryan, thanks!!! but in that I don't know any C at all, I don't think I could work it into my data analysis pipeline.
                  --------------
                  Ethan

                  Comment


                  • #10
                    Anybody know how to do this with a BASH script??
                    You should be able to just loop over all files and use the 'join' command.

                    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

                    • seqadmin
                      Strategies for Sequencing Challenging Samples
                      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...
                      03-22-2024, 06:39 AM
                    • seqadmin
                      Techniques and Challenges in Conservation Genomics
                      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...
                      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 seqadmin  
                    Started by seqadmin, Today, 06:07 PM
                    0 responses
                    7 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-22-2024, 10:03 AM
                    0 responses
                    49 views
                    0 likes
                    Last Post seqadmin  
                    Started by seqadmin, 03-21-2024, 07:32 AM
                    0 responses
                    66 views
                    0 likes
                    Last Post seqadmin  
                    Working...
                    X