Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • glm edgeR issue

    Hi everybody
    I just performed a time series analysis on edgeR and I have a question.
    This is my data frame for the condition and time:
    sample condition time
    1 sample1 Placebo 1h
    2 sample2 Placebo 1h
    3 sample3 Placebo 2h
    4 sample4 Placebo 2h
    3 sample5 Placebo 4h
    4 sample6 Placebo 4h
    5 sample7 Drug 1h
    6 sample8 Drug 1h
    7 sample9 Drug 2h
    8 sample10 Drug 2h
    9 sample11 Drug 4h
    10 sample12 Drug 4h

    and my script to test interaction

    design<-model.matrix(~condition*time,data=meta)
    y<-DGEList(counts=countTable)
    y<-calcNormFactors(y,method="TMM")
    y<-estimateGLMCommonDisp(y,design)
    y<-estimateGLMTrendedDisp(y,design)
    y<-estimateGLMTagwiseDisp(y,design)
    fit<-glmFit(y,design)
    colnames(design)
    "(intercept)" "conditionDrug" "time2h" "time4h" "conditionDrug:time2h" "conditionDrug:time4h"
    lrt<-glmLRT(fit,coef=6)
    In this case with coef=6 I am able to detect the genes that respond differently to the drug, relative to placebo at 4h.
    But same gene that are modulated are not expressed either in placebo 4h and Drug 4h!!!
    Example(matrix of countTable)
    ID Placebo-1h Placebo-1h Placebo-2h Placebo-2h Placebo-4h Placebo-4h Drug-1h Drug-1h Drug-2h Drug-2h Drug-4h Drug-4h
    ENSMUSG000000xxxx 0 5 4 0 0 0 0 3 0 1 0 0
    And here the different comparison
    ID logFC logCPM LR PValue
    "ENSMUSG000000xxxx" -4.36 0.14 6.48 0.01 Placebo4-h
    "ENSMUSG000000xxxx" -2.19 0.14 1.21 0.27 Drug-4h
    "ENSMUSG000000xxxx" 2.16 0.14 100.14 1.41e-23 Contrast(Placebo-4h/Drug-4h)
    As you can see in the contrast Placebo-4h and Drug-4h I have a logFC of 2.16(that is the difference between Placebo-1h/Placebo-4h and Drug-1h/Drug-4h) and a PValue of 1.41e-23!!But I have no reads on this gene!
    How is possible?
    Could you help me to understand if I made same mistake in the analysis?

    Thank you

    Giuliano

  • #2
    No this isn't possible, so you have presumably made a mistake.

    The edgeR User's Guide asks you to send questions to the Bioconductor mailing list. If you write to that list, someone will help diagnose your problem.

    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
    50 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