Question: DEseq package statistical test
gravatar for rohan bareja
15 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 15 months ago by Steve Lianoglou12k • written 15 months ago by rohan bareja200
gravatar for Steve Lianoglou
15 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 15 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: 181 users visited in the last hour