multtest bug fixed
0
0
Entering edit mode
kpollard ▴ 110
@kpollard-7578
Last seen 9.5 years ago
United States
Hello All, I've just found and fixed a bug in the boot.resample function of the multtest package. This function is called by the function MTP. So, if you have been using the MTP function with multtest version 1.5.x, you may have incorrect p-values (rawp and adjp). Other functions, such as mt.minP, mt.maxT and mt.rawp2adjp do not use this function and are OK. The problem was that when standardizing the test statistics null distribution, we used min where we needed pmin. The consequence was that each marginal null distribution was standardized by 1/variance of the gene with the largest variance (not its own variance). Hence, the marginal test statistics null distributions were too narrow and more rejections were made than should have been (ie: pvalues were smaller than they should have been). I have checked this correction into the cvs archive today and the fixed package should be available under developmental packages soon. Version 1.6 (for the upcoming BioC release) will be corrected also. Apologies for any inconvenience. Best, Katie
multtest multtest • 826 views
ADD COMMENT

Login before adding your answer.

Traffic: 582 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