Description Usage Arguments Value Note Author(s) References Examples
Compute a robust estimate of location and scale using the Stahel-Donoho projection based estimator
1 |
x |
a matrix or data frame. |
nsamp |
a positive integer giving the number of resamples required;
|
maxres |
a positive integer specifying the maximum number of
resamples to be performed including those that are discarded due to linearly
dependent subsamples. If |
tune |
a numeric value between 0 and 1 giving the fraction of the data to receive non-zero weight.
Defaults to |
prob |
a numeric value between 0 and 1 specifying the probability of high breakdown point;
used to compute |
eps |
a numeric value between 0 and 0.5 specifying the breakdown point; used to compute
|
seed |
starting value for random generator. Default is |
trace |
whether to print intermediate results. Default is |
control |
a control object (S4) of class |
An S4 object of class CovSde-class
which is a subclass of the
virtual class CovRobust-class
.
The Fortran code for the Stahel-Donoho method was taken almost with no changes from
package robust
which in turn has it from the Insightful Robust Library
(thanks to by Kjell Konis).
Valentin Todorov valentin.todorov@chello.at and Kjell Konis kjell.konis@epfl.ch
R. A. Maronna and V.J. Yohai (1995) The Behavior of the Stahel-Donoho Robust Multivariate Estimator. Journal of the American Statistical Association 90 (429), 330–341.
R. A. Maronna, D. Martin and V. Yohai (2006). Robust Statistics: Theory and Methods. Wiley, New York.
Todorov V & Filzmoser P (2009), An Object Oriented Framework for Robust Multivariate Analysis. Journal of Statistical Software, 32(3), 1–47. URL http://www.jstatsoft.org/v32/i03/.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | data(hbk)
hbk.x <- data.matrix(hbk[, 1:3])
CovSde(hbk.x)
## the following four statements are equivalent
c0 <- CovSde(hbk.x)
c1 <- CovSde(hbk.x, nsamp=2000)
c2 <- CovSde(hbk.x, control = CovControlSde(nsamp=2000))
c3 <- CovSde(hbk.x, control = new("CovControlSde", nsamp=2000))
## direct specification overrides control one:
c4 <- CovSde(hbk.x, nsamp=100,
control = CovControlSde(nsamp=2000))
c1
summary(c1)
plot(c1)
## Use the function CovRobust() - if no estimation method is
## specified, for small data sets CovSde() will be called
cr <- CovRobust(hbk.x)
cr
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.