Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Cuffdiff - ouput problem

    I am attempting to use cuffdiff to do a comparison analysis of various yeast transcriptome samples. I have processed the reads via Tophat and compiled transcripts via Cufflinks. I would like to use Cuffdiff to generate a comparison between samples so that I can start compiling which genes are differentially expressed. I have attempted to run Cuffdiff several times, but I keep getting the same "segmentation fault" error. Here is an example of my input:

    $ cuffdiff ../ethanol/transcripts.gtf ~/RNA_Seq/Yeast/reads_for_comp/tophat_output/ethanol/accepted_hits_e.sam ~/RNA_Seq/Yeast/reads_for_comp/tophat_output/glucose/accepted_hits_g.sam

    I quickly get the following output:

    Segmentation fault (core dumped)

    Does anyone have any experience with this type of error arising from Cuffdiff, and any suggestions as to how I can correct this?

    Thanks in advance for your help.

  • #2
    I don't have a solution, but I do have the same problem..though not always.

    I used tophat and cufflinks with no problems on a total of 7 different samples. running cuffdiff on two particular samples threw the segmentation fault back at me, but worked fine for other samples. Other than the source organism (it worked for mammal, failed for insect), there was no difference in the way that these samples were prepared or analysed up to that point.

    Comment


    • #3
      cuffdiff v0.9.3 seg. fault

      Hello everyone,
      I've found another way (maybe linked with your problem) to get the "Segmentation fault":
      Code:
      cuffdiff  -o ./ input.gtf input1.bam input2.bam
      works fine, but
      Code:
      cuffdiff  -o ./another_directory input.gtf input1.bam input2.bam
      causes the "Segmentation fault"

      The "another_directory is visible and I've the rigth to write in... But I can't reproduce it on any directories (sometimes it works) :-?
      I'm using the current version, v0.9.3
      Anyone has the same problem ?
      Last edited by nicolallias; 02-25-2011, 06:28 AM.

      Comment


      • #4
        Cuffdiff is a great tool but I stopped using it -> too many bugs

        Comment


        • #5
          Originally posted by pasta View Post
          Cuffdiff is a great tool but I stopped using it -> too many bugs
          Do you know of a good alternative?

          Comment


          • #6
            Originally posted by nicolallias View Post
            Hello everyone,
            I've found another way (maybe linked with your problem) to get the "Segmentation fault":
            Code:
            cuffdiff  -o ./ input.gtf input1.bam input2.bam
            works fine, but
            Code:
            cuffdiff  -o ./another_directory input.gtf input1.bam input2.bam
            causes the "Segmentation fault"

            The "another_directory is visible and I've the rigth to write in... But I can't reproduce it on any directories (sometimes it works) :-?
            I'm using the current version, v0.9.3
            Anyone has the same problem ?
            This is actually exactly the same problem I was facing until I read your post. It appears that specifying an outdir other than the current directory caused a segmentation fault. So strange.

            Comment


            • #7
              Dont specify output directory

              Originally posted by fongchun View Post
              This is actually exactly the same problem I was facing until I read your post. It appears that specifying an outdir other than the current directory caused a segmentation fault. So strange.
              I had the segmentation fault despite specifying the same outdir as the current directory. However, if you DO NOT specify the output directory, the segmentation fault goes away. Cuffdiff is working for me now. Strange.

              Comment


              • #8
                Originally posted by Aurelien Mazurie View Post
                Do you know of a good alternative?
                Two alternatives exist:

                1. Use HTSeq followed by DESeq : you must note though that there is a significant difference between how cufflinks/cuffdiff and HTSeq/DESeq look at mapping for quantification, there is an enlightening discussion here cufflinks/cuffdiff - HTSeq/DESeq debate

                2. There is RSEQTools from the Gerstein lab, this I have not used personally, but the sequencing center attached to where I work base their RNA-Seq analysis pipeline on RSEQTools.

                Vineeth

                Comment


                • #9
                  Addendum

                  Originally posted by vineeth_s View Post
                  Two alternatives exist:

                  there is an enlightening discussion here cufflinks/cuffdiff - HTSeq/DESeq debate
                  As that thread has quite many extraneous posts, follow the conversation between Lior Pachter (lpachter) and Simon Anders

                  Comment


                  • #10
                    Originally posted by flobpf View Post
                    I had the segmentation fault despite specifying the same outdir as the current directory. However, if you DO NOT specify the output directory, the segmentation fault goes away. Cuffdiff is working for me now. Strange.
                    You are right. I specified -o ./cuffdiff_output and it gave me a segmentation fault halfway into the execution. It seems that the only surefire way of not crashing is to NOT specify any output directory.

                    Comment


                    • #11
                      Hi All,
                      If you don't know already, you can submit your cuffdiff problems to [email protected]

                      Comment


                      • #12
                        Hi All,

                        Besides poisson200's suggestion, please make sure you are using the latest version to see if it helps. I know this package is under active development.

                        To me personally, different programs generate different numbers (F/R PKM), which must be unified down the road. It is not good for the science.

                        Douglas

                        Comment


                        • #13
                          Originally posted by flobpf View Post
                          I had the segmentation fault despite specifying the same outdir as the current directory. However, if you DO NOT specify the output directory, the segmentation fault goes away. Cuffdiff is working for me now. Strange.
                          I'm glad you all figured it out. I had the same problem and now it works with "-o ./" setting..Thanks!!!!

                          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, Yesterday, 06:37 PM
                          0 responses
                          10 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, Yesterday, 06:07 PM
                          0 responses
                          9 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
                          67 views
                          0 likes
                          Last Post seqadmin  
                          Working...
                          X