binomRDtest {MCPAN} | R Documentation |
Simultaneous test for contrasts of independent binomial proportions (in a oneway layout)
Description
P-value of maximum test and adjusted p-values for M contrasts of I groups in a one-way layout. Tests are performed for contrasts of proportions, which can be interpreted as differences of (weighted averages of) proportions.
Usage
binomRDtest(x, ...)
## Default S3 method:
binomRDtest(x, n, names=NULL,
type="Dunnett", cmat=NULL, method="Wald",
alternative="two.sided", dist="MVN", ...)
## S3 method for class 'formula'
binomRDtest(formula, data,
type="Dunnett", cmat=NULL, method="Wald",
alternative="two.sided", dist="MVN", ...)
## S3 method for class 'table'
binomRDtest(x, type="Dunnett",
cmat=NULL, method="Wald", alternative="two.sided",
dist="MVN", ...)
## S3 method for class 'matrix'
binomRDtest(x, type="Dunnett",
cmat=NULL, method="Wald", alternative="two.sided",
dist="MVN", ...)
Arguments
x |
a numeric vector, giving the number of successes in I independent samples, or an object of class "table", representing the 2xk-table, or an object of class "matrix", representing the 2xk-table |
n |
a numerioc vector, giving the number of trials (i.e. the sample size) in each of the I groups |
names |
an optional character vector, giving the names of the groups in x, n; if not specified, possibly availbale names of x are taken as group names |
formula |
a two-sided formula of the style 'response ~ treatment', where 'response' should be a categorical variable with two levels, while treatment should be a factor specifying the treatment levels |
data |
a data.frame, containing the variables specified in formula |
type |
a character string specifying the contrast type |
cmat |
an optional user defined contrast matrix of dimension MxI |
method |
a single charcter string, specifying the method for adjustment, with options: "Wald" (Maximum likelihood estimators), "ADD1" (add1-adjustment on the raw proportion estimates) "ADD2" (add2-adjustment on proportion estimates following Agresti Caffo (2000)) |
alternative |
a character string specifying the direction of the alternative hypothesis |
dist |
a character string, where "MVN" invokes the computation of p-values using the multivariate normal distribution, and "N" invokes use p-value computation using the univariate normal distribution |
... |
arguments to be passed to binomest, currently only success labelling the event which should be considered as success |
Details
For usage, see the examples.
Value
An object of class "binomRDtest", a list containing:
teststat |
a numeric vector of teststatistics of length M |
pval |
a single numeric p-value, the p-value of the maximum test (minimum p-value) |
p.val.adj |
a vector of length M, the adjusted p-values of the single contrasts |
dist |
character string indicating whether the multivariate normal or normal distribution was used for computation of p-values |
alternative |
a single character vector, as the input |
x |
the observed number of successes in the treatment groups |
n |
the number of trials in the treatment groups |
p |
the estimated proportions in the treatment groups |
success |
a character string labelling the event considered as success |
method |
as input, a character string |
cmat |
used contrast matrix |
Note
Note, that all implemented methods are approximate only. The size of the test might seriously deviate from the nominal level for small sample sizes and extreme success probabilities. See the simulation results in Sill (2007) for details.
References
Statistical procedures and characterization of coverage probabilities are described in: Sill, M. (2007): Approximate simultaneous confidence intervals for multiple comparisons of binomial proportions. Master thesis, Institute of Biostatistics, Leibniz University Hannover.
See Also
Examples
ntrials <- c(40,20,20,20)
xsuccesses <- c(1,2,2,4)
names(xsuccesses) <- LETTERS[1:4]
binomRDtest(x=xsuccesses, n=ntrials, method="ADD1",
type="Dunnett")
binomRDtest(x=xsuccesses, n=ntrials, method="ADD1",
type="Williams", alternative="greater")
binomRDtest(x=xsuccesses, n=ntrials, method="ADD2",
type="Williams", alternative="greater")