Seqanswers Leaderboard Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • jwalker_tgi
    Junior Member
    • Mar 2012
    • 1

    Cuffdiff DE significance of zero FPKM values

    Hi-

    I was using cuffdiff 1.3.0 to compare four BAMs aligned with TopHat 1.4.1. I have noticed that zero FPKM values are never considered significant when compared with large FPKM values. Here is one example that I am concerned about:

    XLOC_017834 XLOC_017834 Olfr1307 2:111784672-111785611 FF FMC OK 0.257545 46.5301 7.4972 -4.50075 6.77156e-06 0.000658507 yes
    XLOC_017834 XLOC_017834 Olfr1307 2:111784672-111785611 FFC FMC OK 0 46.5301 1.79769e+308 1.79769e+308 0.0789649 0.365761 no
    XLOC_017834 XLOC_017834 Olfr1307 2:111784672-111785611 FM FMC OK 0.0380371 46.5301 10.2565 -5.37179 7.79602e-08 2.39776e-05 yes

    If line wrapping destroys the above genes_exp.diff output, here is a brief summary: FMC has an FPKM of 46.5301 and when compared to FF (FPKM=0.257545) and FM (FPKM=0.0380371) the q-value is below the FDR of 0.05. Appropriately the significance value is labeled as 'yes'. However, FFC has an FPKM of zero and this differential test is not significant. The goal of the experiment was to identify uniquely expressed genes only identified in the FMC data set. As such an infinite log(fold_change) seems more significant than the other comparisons. Can anyone explain this?
  • kenietz
    Member
    • Nov 2011
    • 86

    #2
    Hi,
    that is what is bugging my brain as well these days

    Maybe the math is not working properly with that zero value. But thats a guess only. Came to the idea to replace all the zeros with something much smaller than the minimum FPKM found in the GTF files. If the minimum FPKM is 0.01 i will try to replace all complete zeros with 0.00001 for example and then do the DE and see whats the difference.

    Comment

    • IBseq
      Member
      • Jul 2012
      • 56

      #3
      hi did you find out an answer to your question?
      I do have the same prob and cnt find a solution,

      best,
      ib

      Comment

      • bioliyezhang
        Member
        • Mar 2011
        • 19

        #4
        Hi,

        about this point.
        I do see significant ones in my results, shown in here.

        I am using Tophat 2.0 and Cufflink 2.0.

        Comment

        • IBseq
          Member
          • Jul 2012
          • 56

          #5
          but did u have replicates?i have one to one comparison, no replicates (biological or technical)

          Comment

          • bioliyezhang
            Member
            • Mar 2011
            • 19

            #6
            I do have replicates(3 in total).

            Therefore I checked two other samples, which I only have one replicate.
            The result is the same as yours, any zero value on one sample is not considered as significant. So I guess they do not considered one zero and one non-zero to be significant in their algorithms.

            Comment

            • IBseq
              Member
              • Jul 2012
              • 56

              #7
              Thanks....how would you approach this issue?I cannot sequence replicates thus I have to stick with these results, but then THEY ARE NOT RELIABLE...??
              I have use a recent program (recently published) called GFOLD. This, ranks all the DE but does not tell you if it is significant, but at lwast I can see the assigned values. Still, cant figure out of to test significance..

              ib

              Comment

              • bioliyezhang
                Member
                • Mar 2011
                • 19

                #8
                Hi,

                I am not a real expert on statistics...
                so I still think some of them may be reliable, however, they might not pass the statistics test.
                I haven't tried to tackle that yet, but I am planning to maybe combine some other similar control sample? In this way, we may have > 1 control sample and 1 test sample, which may give us more confidence in DE detection. This is my plan but I haven't tested it out yet. Just my 2 cents.

                liye

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Pathogen Surveillance with Advanced Genomic Tools
                  by seqadmin




                  The COVID-19 pandemic highlighted the need for proactive pathogen surveillance systems. As ongoing threats like avian influenza and newly emerging infections continue to pose risks, researchers are working to improve how quickly and accurately pathogens can be identified and tracked. In a recent SEQanswers webinar, two experts discussed how next-generation sequencing (NGS) and machine learning are shaping efforts to monitor viral variation and trace the origins of infectious...
                  03-24-2025, 11:48 AM
                • seqadmin
                  New Genomics Tools and Methods Shared at AGBT 2025
                  by seqadmin


                  This year’s Advances in Genome Biology and Technology (AGBT) General Meeting commemorated the 25th anniversary of the event at its original venue on Marco Island, Florida. While this year’s event didn’t include high-profile musical performances, the industry announcements and cutting-edge research still drew the attention of leading scientists.

                  The Headliner
                  The biggest announcement was Roche stepping back into the sequencing platform market. In the years since...
                  03-03-2025, 01:39 PM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, 03-20-2025, 05:03 AM
                0 responses
                49 views
                0 reactions
                Last Post seqadmin  
                Started by seqadmin, 03-19-2025, 07:27 AM
                0 responses
                57 views
                0 reactions
                Last Post seqadmin  
                Started by seqadmin, 03-18-2025, 12:50 PM
                0 responses
                50 views
                0 reactions
                Last Post seqadmin  
                Started by seqadmin, 03-03-2025, 01:15 PM
                0 responses
                201 views
                0 reactions
                Last Post seqadmin  
                Working...