Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Log ratios and/or ratios of the log's

    Hello all.
    I have a very basic question about the analisys of my RNA-seq data...can some one help me out?

    Usually we work with the log ratio of RPKM’s, log2(X1 RPKM/ X2 RPKM).
    I would like to know if instead of using log ratios I can aply log to all my RPKM data sets and from then use this log2(rpkm) values for the rest of the analysis. So, basically istead of doing log ratios, do the ratio between log2(RPKM). Example log2(X1 RPKM)= 5 and log2(X2 RPKM) = 5 so ratio=1 meaning equal expression. This is completely wrong and I have to use the log ratios or can I do this way also?


    Thanks

  • #2
    log(a / b) = log(a) - log(b),

    so you can do what you're proposing, but you have to substract the two logs, not divide them.

    Comment


    • #3
      Originally posted by ffinkernagel View Post
      log(a / b) = log(a) - log(b),

      so you can do what you're proposing, but you have to substract the two logs, not divide them.

      I don't understand how that changes the problem. Log(0) is undefined. In that case should I just use 0? If so how will the log(25) - 0 results be different from a result like log(25) - log(x!=0).

      Comment


      • #4
        This mabye a dumb question and I needs me some splainin', but ...

        I know folks always did log transforms on old school array data (because it better reflected the known inputs in spike in data).

        But ... why do people do log transforms on RPKMs from RNA seq data ? Isn't 2X reads on a gene really mean 2X expression?

        Comment


        • #5
          Originally posted by Richard Finney View Post
          This mabye a dumb question and I needs me some splainin', but ...

          I know folks always did log transforms on old school array data (because it better reflected the known inputs in spike in data).

          But ... why do people do log transforms on RPKMs from RNA seq data ? Isn't 2X reads on a gene really mean 2X expression?

          I'm not doing a typical experiment. I'm using RNA-Seq to predict transcriptome secondary structure in bacteria by the PARS method outlined in the 2010 nature paper that can be found here: http://genie.weizmann.ac.il/pubs/PARS10/index.html

          I run two protocols, one with a digestion that cuts at single stranded positions and one with a digestion that cuts at double stranded positions.

          To say if any position is in a secondary structure requires that I have counts of the how many reads start at each position, so it isn't measured in RPKM or any other measure usually used for differential expression. Now I have cleaned, aligned, and produced files containing counts for number of reads starting at each position.

          What I measure is the log(protocol1/protocol2) counts at each position to determine if there is secondary structure at each point. However when there is a case that is log(25 single strand/0 double strand) i cant accurately compute a score for that position.

          I still don't have a workable solution to this. I've been in the literature and found that there are some statistical methods that can apply to DE on a per gene basis that I may be able to apply to a per position basis if i tweak the statistics a bit, but I'd rather go with something that has been peer reviewed and tested if ya know what I mean.

          Comment


          • #6
            What I plan to try, is to get distributions of read counts for both protocols and use that to compute if a 0 is significant or not.

            Comment

            Latest Articles

            Collapse

            • 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
            • seqadmin
              The Impact of AI in Genomic Medicine
              by seqadmin



              Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
              02-26-2024, 02:07 PM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 03-14-2024, 06:13 AM
            0 responses
            33 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-08-2024, 08:03 AM
            0 responses
            72 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-07-2024, 08:13 AM
            0 responses
            81 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-06-2024, 09:51 AM
            0 responses
            68 views
            0 likes
            Last Post seqadmin  
            Working...
            X