Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • DESeq2 fold change does not match normalized counts data

    I am trying to understand how DESeq2 calculates Log2foldchange, because it does not match my calculations of fold change from my normalized counts.

    My script:
    CountTable_POS = read.table("DIF_pos_counts", header=TRUE, sep="")
    Design = data.frame( row.names = colnames( CountTable_POS ),
    condition = c( "untreated", "untreated", "untreated",
    "treated", "treated", "treated" ),
    type = c( "single-end", "single-end", "single-end",
    "single-end", "single-end","single-end" ) )
    dds_pos<-DESeqDataSetFromMatrix(countData = CountTable_POS,
    colData = Design,
    design = ~ condition)
    dds_pos<-DESeq(dds_pos, fitType="local")
    norm<-counts(dds_pos, normalized=T)
    write.csv(as.data.frame(norm), file="norm_counts.csv")
    res_pos<-results(dds_pos, contrast=c("condition", "treated", "untreated"))
    write.csv(as.data.frame(res_pos), file="DESeq2_out.csv")

    RESULTS
    from norm_counts.csv:
    Locus Control1 Control2 Control3 Treated1 Treated2 Treated3 Mean_Control Mean_Treated Basemean Fold change
    2R_188654 3.8127 5.7220 5.7191 1.9351 4.9898 2.9518 5.0846 3.2922 4.1884 0.6475
    2R_188736 3.8127 3.8147 6.6723 1.9351 4.9898 5.9036 4.7666 4.2762 4.5214 0.8971

    From DESeq2_out:
    Locus Basemean Log2foldchange P value Fold change
    2R_188654 4.1884 -0.3108 0.4093 0.8062
    2R_188736 4.5214 -0.0819 0.8279 0.9448

    As you can see, the basemeans (mean of control and treated counts) match- which tells me there is no internal error- but the fold changes do not. Why is this???

  • #2
    Fold-changes are shrunken, so they won't match what you calculate by hand (the shrinkage produces more reliable results). If you want the number that will more or less match what you're calculating them have it output the MLE results (use "addMLE=T" in results()).

    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