Question: DEseq package statistical test
gravatar for rohan bareja
5 months ago by
rohan bareja200
rohan bareja200 wrote:

Hi all,


I use DESeq2 package for differential expression analysis ,where df is the the count data.




##do differential expression
dds <- DESeqDataSetFromMatrix(countData = data1,colData = pdata,design = ~ condition )
#reference is group2,so up means up in group1 and down means down in group1.

dds$condition <- relevel(dds$condition, "group2")
res <- results(dds)

I wanted to know what test is used in this by default "Wald test or LRT test  or no test " because I am not specifying anything in my command.





ADD COMMENTlink modified 5 months ago by Steve Lianoglou12k • written 5 months ago by rohan bareja200
gravatar for Steve Lianoglou
5 months ago by
Steve Lianoglou12k wrote:

Take a look at the ?DESeq2 help file and take note of the function signature:

DESeq(object, test = c("Wald", "LRT"), fitType = c("parametric", "local",
  "mean"), betaPrior, full = design(object), reduced, quiet = FALSE,
  minReplicatesForReplace = 7, modelMatrixType, parallel = FALSE,
  BPPARAM = bpparam())

Look at the test parameter. What you see here is a common idiom in R programming. When a function parameter takes an array like this c("Wald", "LRT"), it almost always implies that it expects the user to provide one of the values given in the vector of choices. This is usually checked with the a call to the match.arg function inside the function definition. When a user neglects to provide a value for the parameter, the first element on the "choice" vector is chosen, which in this case is "Wald".


ADD COMMENTlink written 5 months ago by Steve Lianoglou12k
Please log in to add an answer.


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