The support.bioconductor.org editor has been updated to markdown! Please see more info at: Tutorial: Updated Support Site Editor

Question: DESEQ2 Default Value for Cooks Distance cutoff
0
gravatar for wpoehlm
8 months ago by
wpoehlm0
wpoehlm0 wrote:

Hi,

 

I am using version 1.16.1 of DESeq2.  According to the documentation, a Cook's distance cutoff of 0.99 is applied to the results function by default.  I looked at the source code in github and this certainly seems to be the case.  However, if I get a larger number of significant genes when running the "results" function with default arguments compared to when I explicitly use the argument "cooksCutoff = 0.99".  Is version 1.16.1 in Bioconductor up to date with the source code on github? (https://github.com/mikelove/DESeq2).  I have not done much digging otherwise. 

 

Thanks,

 

Will

deseq2 • 146 views
ADD COMMENTlink modified 8 months ago by Michael Love22k • written 8 months ago by wpoehlm0

Update: I updated to version 1.20.0 and I am noticing the same thing

ADD REPLYlink modified 8 months ago • written 8 months ago by wpoehlm0
Answer: DESEQ2 Default Value for Cooks Distance cutoff
0
gravatar for Michael Love
8 months ago by
Michael Love22k
United States
Michael Love22k wrote:

The default value for Cook's distance is not 0.99. From the help page for ?results:

"cooksCutoff : theshold on Cook’s distance, such that if one or more samples for a row have a distance higher, the p-value for the row is set to NA. The default cutoff is the .99 quantile of the F(p, m-p) distribution, where p is the number of coefficients being fitted and m is the number of samples."

ADD COMMENTlink written 8 months ago by Michael Love22k
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: 442 users visited in the last hour