Entering edit mode
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