Description Usage Arguments Details Value Author(s) References See Also Examples
View source: R/FRBhotellingS.R
Robust one-sample and two-sample Hotelling test using the S-estimator and the Fast and Robust Bootstrap.
1 2 3 4 5 6 7 |
formula |
an object of class |
data |
data frame from which variables specified in formula are to be taken. |
X |
a matrix or data-frame |
Y |
an optional matrix or data-frame in case of a two-sample test |
mu0 |
an optional vector of data values (or a single number which will be repeated p times) indicating the
true value of the mean (does not apply in case of the two-sample test). Default is the null vector |
R |
number of bootstrap samples. Default is |
bdp |
required breakdown point. Should have 0 < |
conf |
confidence level for the simultaneous confidence intervals. Default is |
method |
for the two-sample Hotelling test, indicates the way the common covariance matrix is estimated:
|
control |
a list with control parameters for tuning the computing algorithm, see |
na.action |
a function which indicates what should happen when the data contain NAs. Defaults to |
... |
allows for specifying control parameters directly instead of via |
The classical Hotelling test for testing if the mean equals a certain center or if two means
are equal is modified into a robust one through substitution of the empirical estimates
by the S-estimates of location and scatter. The S-estimator uses Tukey's biweight function where the constant is chosen
to obtain the desired breakdown point as specified by bdp
.
One-sample S-estimates are computed by a call to the implementation of the fast-S algorithm in the rrcov package of Todorov and Filzmoser (2009). For two-sample S-estimates an adaptation of the fast-S algorithm is used.
The tuning parameters of the algorithm can be changed via control
.
The fast and robust bootstrap is used to mimic the distribution of the test statistic under the null hypothesis. For instance, the 5% critical value for the test is given by the 95% quantile of the recalculated statistics.
Robust simultaneous confidence intervals for linear combinations of the mean (or difference in means) are developed similarly to the classical case
(Johnson and Wichern, 1988, page 239). The value CI
is a matrix with the confidence intervals for each element
of the mean (or difference in means), with level conf
. It consists of two rows, the first being the lower bound and the second the upper bound.
Note that these intervals are rather conservative in the sense that the simultaneous confidence level holds for all linear combinations
and here only p of these are considered (with p the dimension of the data).
For the two-sample Hotelling test we assume that the samples have an underlying distribution with the same covariance matrix.
This covariance matrix can be estimated in two different ways using the pooled covariance matrix or the two-sample
estimator of He and Fung (He and Fung 2000), and argument method
defaults to the second option.
For more details see Roelant et al. (2008).
In the two-sample version, the null hypothesis always states that the two means are equal. For the one-sample version, the default
null hypothesis is that the mean equals zero, but the hypothesized value can be changed and specified through argument mu0
.
Bootstrap samples are discarded if the fast and robust covariance estimate is not positive definite, such that the actual number
of recalculations used can be lower than R
. This number is returned as ROK
.
An object of class FRBhot
which extends class htest
and contains at least the following components:
statistic |
the value of the robust test statistic. |
pvalue |
p-value of the robust one or two-sample Hotelling test, determined by the fast and robust bootstrap |
estimate |
the estimated mean vector or vectors depending on whether it was a one-sample test or a two-sample test. |
alternative |
a character string describing the alternative hypothesis. |
method |
a character string indicating what type of Hotelling test was performed. |
data.name |
a character string giving the name(s) of the data. |
teststat.boot |
the bootstrap recalculated values of the robust test statistic. |
CI |
bootstrap simultaneous confidence intervals for each component of the center |
conf |
a copy of the |
Sigma |
covariance of one-sample or common covariance matrix in the case of two samples |
w |
implicit weights corresponding to the S-estimates (i.e. final weights in the RWLS procedure at the end of the fast-S algorithm) |
outFlag |
outlier flags: 1 if the robust distance of the observation exceeds the .975 quantile of (the square root of)
the chi-square distribution with degrees of freedom equal to the dimension of |
ROK |
number of bootstrap samples actually used (i.e. not discarded due to non-positive definite covariance |
Ella Roelant, Stefan Van Aelst and Gert Willems
X. He and W.K. Fung (2000) High breakdown estimation for multiple populations with applications to discriminant analysis. Journal of Multivariate Analysis, 72, 151–162.
R.A. Johnson, D.W. Wichern (1988) Applied Multivariate Statistical Analysis, 2nd Edition, Prentice-Hall.
E. Roelant, S. Van Aelst and G. Willems, (2008) Fast Bootstrap for Robust Hotelling Tests, COMPSTAT 2008: Proceedings in Computational Statistics (P. Brito, Ed.) Heidelberg: Physika-Verlag, 709–719.
M. Salibian-Barrera, S. Van Aelst and G. Willems (2008) Fast and robust bootstrap. Statistical Methods and Applications, 17, 41–71.
V. Todorov and P. Filzmoser (2009) An Object-Oriented Framework for Robust Multivariate Analysis. Journal of Statistical Software, 32, 1–47. URL http://www.jstatsoft.org/v32/i03/.
S. Van Aelst and G. Willems (2013). Fast and Robust Bootstrap for Multivariate Inference: The R Package FRB. Journal of Statistical Software, 53(3), 1–32. URL: http://www.jstatsoft.org/v53/i03/.
plot.FRBhot
,summary.FRBhot
, FRBhotellingMM
, Scontrol
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 39 40 41 42 43 44 45 | ## One sample robust Hotelling test
data(delivery)
delivery.x <- delivery[,1:2]
FRBhotellingS(delivery.x,R=199)
## One sample robust Hotelling test
data(ForgedBankNotes)
samplemean <- apply(ForgedBankNotes, 2, mean)
res = FRBhotellingS(ForgedBankNotes, mu0=samplemean,R=199)
res
# Note that the test rejects the hypothesis that the true mean equals the
# sample mean; this is due to outliers in the data (i.e. the robustly estimated
# mean apparently significantly differs from the non-robust sample mean.
# Graphical display of the results:
plot(res)
# It is clear from the (scaled) simultaneous confidence limits that the rejection
# of the hypothesis is due to the differences in variables Bottom and Diagonal
## Two sample robust Hotelling test
data(hemophilia)
grp <-as.factor(hemophilia[,3])
x <- hemophilia[which(grp==levels(grp)[1]),1:2]
y <- hemophilia[which(grp==levels(grp)[2]),1:2]
#using the pooled covariance matrix to estimate the common covariance matrix
## Not run: res = FRBhotellingS(x,y,method="pool")
#using the estimator of He and Fung to estimate the common covariance matrix
res = FRBhotellingS(x,y,method="HeFung",R=199)
# or using the formula interface
## Not run: res = FRBhotellingS(as.matrix(hemophilia[,-3])~hemophilia[,3],method="HeFung",R=99)
# From the confidence limits it can be seen that the significant difference
# is mainly caused by the AHFactivity variable. The graphical display helps too:
plot(res)
# the red line on the histogram indicates the test statistic value in the original
# sample (it is omitted if the statistic exceeds 100)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.