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
                            Quality Control Essentials for Next-Generation Sequencing Workflows
                            by seqadmin




                            Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

                            Nucleic Acid Quality Control
                            Preparing for NGS starts with isolating the...
                            02-10-2025, 01:58 PM
                          • seqadmin
                            An Introduction to the Technologies Transforming Precision Medicine
                            by seqadmin


                            In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
                            01-27-2025, 07:46 AM

                          ad_right_rmr

                          Collapse

                          News

                          Collapse

                          Topics Statistics Last Post
                          Started by seqadmin, 02-07-2025, 09:30 AM
                          0 responses
                          64 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 02-05-2025, 10:34 AM
                          0 responses
                          99 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 02-03-2025, 09:07 AM
                          0 responses
                          78 views
                          0 likes
                          Last Post seqadmin  
                          Started by seqadmin, 01-31-2025, 08:31 AM
                          0 responses
                          44 views
                          0 likes
                          Last Post seqadmin  
                          Working...
                          X