DESEQ2 Default Value for Cooks Distance cutoff
1
0
Entering edit mode
wpoehlm • 0
@wpoehlm-16199
Last seen 5.8 years ago

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 • 1.1k views
ADD COMMENT
0
Entering edit mode

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

ADD REPLY
0
Entering edit mode
@mikelove
Last seen 4 hours ago
United States

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 COMMENT

Login before adding your answer.

Traffic: 886 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6