Question: 51 views DESeq2 RNAseq differential expression analysis
0
gravatar for adeler001
12 weeks ago by
adeler0010
adeler0010 wrote:

Hello Im trying to do an RNA seq differential expression analysis using DESeq2 by doing a likelihood ratio test The DESeq2 manual states " The likelihood ratio test can be performed by ...providing a reduced design formula, e.g. one in which a number of terms from design(dds) are removed"

They then shows this as an example script you would use to this : dds <- DESeq(dds, test="LRT", reduced=~1)

1) how do you determine how many terms are removed from design(dds) ? 2) Is the number of terms removed always 1 like they show in the example script above ?

ADD COMMENTlink modified 12 weeks ago by Michael Love26k • written 12 weeks ago by adeler0010
Answer: 51 views DESeq2 RNAseq differential expression analysis
0
gravatar for Michael Love
12 weeks ago by
Michael Love26k
United States
Michael Love26k wrote:

You can remove as many coefficients as you want. This is up to you as the analyst.

You can either provide a reduced design formula, as above, or you can provide model matrices directly:

dds <- DESeq(dds, test="LRT", full=mm, reduced=mm0)

Where mm and mm0 are created by model.matrix().

ADD COMMENTlink written 12 weeks ago by Michael Love26k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 174 users visited in the last hour