Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • independet filtering and experimental design in DESeq

    Hi everybody,

    I know this problem to be discussed quite a lot. I read the posts here and here (and also the papers mentioned in them).

    I have two questions concerning my experiment. One is about the experimental design, the second about how to set the filtering. I think they are both somehow connected, so I would like to place them in one post.

    In my experiment we have three conditions (ctrl, KO1 and KO2) and three separate cell types ( I, P, and NP).
    I would like to understand better how to analyse the data in one go.

    The aim of the experiment is not only to compare the ctrl vs. KO1 and/or KO2, but also to analyse the efficiency of cellular processes by comparing NP vs. P in ctrl and/or KO1 and KO2.

    I ran the analysis once with all genes (without any filtering at all, first!). I compared the ctrl vs. KO1 and KO2. It was interesting to see, that in all the comparisons of ctrl vs. KO1 I get a long list of significantly deregulated genes (FDR=0.1%), but in the comparison ctrl vs. KO2 I get only 2-5 genes.
    So I thought a good explanation for that will be filtering the low-count genes. In search of a good cutoff I tried the genefilter package and got the following rank plot:
    Click image for larger version

Name:	rank_scatterplot.png
Views:	1
Size:	95.1 KB
ID:	308315

    Q1: I was wondering if cutting the data set at 0.57 is a good decision.

    Than I looked for a FDR value and did the rejection plot, to see how many genes I am left with, with each of the different FDR values.
    Click image for larger version

Name:	rejection_plot.png
Views:	1
Size:	67.2 KB
ID:	308316
    It was interessting to see, that from 0%-50% they are all overlap each other.

    Q2: Does that mean, that there is no difference between ϑ=0.5 and ϑ=0.1?

    pair-wise vs. multifactor design:

    I read the DESeq manual and ran the analysis as described here:
    Code:
    pd <- read.delim2("../phenoData.txt", sep="\t",quote="", row.names=1)
    
    featureCountTable = read.table( "countTable.txt", header=TRUE, row.names=1, quote="")
    
    conditions = factor(pd$comparison) # I have [COLOR="Red"]nine conditions[/COLOR] are ctrl_I, ctrl_NP, ctrl_P, KO1_I, KO1,NP, KO1_P KO2_I, KO2_NP and KO2_P
    
    cds = newCountDataSet( featureCountTable, conditions )
    
    cds = estimateSizeFactors( cds )
    normResults <- counts( cds, normalized=TRUE ) 
    
    #Variance estimation
    cds = estimateDispersions(cds)
    
    # I than ran for each comparison a binomial test
    res_I_ctrl_KO1 = nbinomTest( cds, "ctrl_I", "KO1_I" )
    res_P_ctrl_KO1 = nbinomTest( cds, "ctrl_P", "KO1_P" )
    ...
    I was wondering if DESeq can work this way or if I need to run a multi-factor design such as

    Code:
    fit1 = fitNbinomGLMs( cdsFullDataSet, count ~ libType + condition )
    fit0 = fitNbinomGLMs( cdsFullDataSet, count ~ libType )
    whereas libType will be the ctrl, KO1 and KO2 and condition will be I, NP and P.

    It will be great if I can get some help.

    thanks a lot

    Assa

Latest Articles

Collapse

  • seqadmin
    Recent Advances in Sequencing Analysis Tools
    by seqadmin


    The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
    05-06-2024, 07:48 AM
  • seqadmin
    Essential Discoveries and Tools in Epitranscriptomics
    by seqadmin




    The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
    04-22-2024, 07:01 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 05-14-2024, 07:03 AM
0 responses
19 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-10-2024, 06:35 AM
0 responses
44 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-09-2024, 02:46 PM
0 responses
54 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-07-2024, 06:57 AM
0 responses
42 views
0 likes
Last Post seqadmin  
Working...
X