Approximately unbiased tests of regions using multistep-multiscale bootstrap resampling

Abstract

Approximately unbiased tests based on bootstrap probabilities are considered for the exponential family of distributions with unknown expectation parameter vector, where the null hypothesis is represented as an arbitrary-shaped region with smooth boundaries. This problem has been discussed previously in Efron and Tibshirani [Ann. Statist. 26 (1998) 1687-1718, and a corrected p-value with second-order asymptotic accuracy is calculated by the two-level bootstrap of Efron, Halloran and Holmes [Proc. Natl. Acad. Sci. U.S. A. 93 (1996) 13429-13434] based on the ABC bias correction of Efron [J. Amer. Statist. Assoc. 82 (1987) 171-185]. Our argument is an extension of their asymptotic theory, where the geometry, such as the signed distance and the curvature of the boundary, plays an important role. We give another calculation of the corrected p-value without finding the "nearest point" on the boundary to the observation, which is required in the two-level bootstrap and is an implementational burden in complicated problems. The key idea is to alter the sample size of the replicated dataset from that of the observed dataset. The frequency of the replicates falling in the region is counted for several sample sizes, and then the p-value is calculated by looking at the change in the frequencies along the changing sample sizes. This is the multiscale bootstrap of Shimodaira [Systematic Biology 51 (2002) 492-508], which is third-order accurate for the multivariate normal model. Here we introduce a newly devised multistep-multiscale bootstrap, calculating a third-order accurate pvalue for the exponential family of distributions. In fact, our p-value is asymptotically equivalent to those obtained by the double bootstrap of Hall [The Bootstrap and Edgeworth Expansion (1992) Springer, New York] and the modified signed likelihood ratio of Barndorff-Nielsen [Biometrika 73 (1986) 307-322] ignoring O(n −3/2 ) terms, yet the computation is less demanding and free from model specification.

New phylogenetic hypotheses for the core Chlorophyta based on chloroplast sequence data

“…The analysis settings used in the constrained analyses were identical to those in the unconstrained analyses (GTR+I+ model, partitioned by codon position). The best tree from each constrained analysis was compared to the ML tree from the corresponding unconstrained analysis using the Approximately Unbiased (AU) test in Consel (Shimodaira, 2004) with sitewise likelihoods calculated by RAxML. …”

A Human Ubiquitin Conjugating Enzyme (E2)-HECT E3 Ligase Structure-function Screen

“…Agglomerative hierarchical clustering of the similarity distance matrix was performed with a euclidean metric and average linkage using the pvclust library (42,43), which performs average-linkage hierarchical clustering while assessing its uncertainty through bootstrap resampling and providing "approximately unbiased" (AU) probability values (p values) and "bootstrap probability" (BP) values. Our analysis was performed using 1,000,000 bootstrap replications.…”

The Cytosolic Protein Response as a Subcomponent of the Wider Heat Shock Response inArabidopsis

“…Cluster analysis of our microarray data with public data of responses following heat shock (Kilian et al, 2007;Gene Expression Omnibus [GEO] accession number GSE5628; Control. GEO accession number GSE5620) was computed by multiscale bootstrap resampling (10,000 iterations; Shimodaira, 2004). The distance matrix was evaluated using the Pearson correlation coefficient, and clusters were created using the complete linkage method.…”