1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
\name{hommel.test}
\alias{hommel.test}
\title{Test Multiple Comparisons Using Hommel's Method}
\description{Given a set of p-values and a test level, returns vector of test results for each hypothesis.}
\usage{hommel.test(p, alpha=0.05)
}
\arguments{
\item{p}{numeric vector of p-values}
\item{alpha}{numeric value, desired significance level}
}
\value{logical vector indicating whether each hypothesis is accepted}
\details{
This function implements the multiple testing procedure of Hommel (1988).
Hommel's method is also implemented as an adjusted p-value method in the function \code{p.adjust} but the accept/reject approach used here is faster.
}
\references{
Hommel, G. (1988). A stagewise rejective multiple test procedure based on a modified Bonferroni test. \emph{Biometrika}, \bold{75}, 383-386.
Shaffer, J. P. (1995). Multiple hypothesis testing. \emph{Annual Review of Psychology} \bold{46}, 561-576. (An excellent review of the area.)
}
\author{Gordon Smyth}
\seealso{
\code{\link[stats:p.adjust]{p.adjust}}
}
\examples{
p <- sort(runif(100))[1:10]
cbind(p,p.adjust(p,"hommel"),hommel.test(p))
}
\keyword{htest}
|
