CovNASde | R Documentation |
Compute a robust estimate of location and scale using the Stahel-Donoho projection based estimator
CovNASde(x, nsamp, maxres, tune = 0.95, eps = 0.5, prob = 0.99,
impMeth = c("norm" , "seq", "rseq"), seed = NULL, trace = FALSE, control)
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 |
impMeth |
select imputation method to use - choose one of "norm" , "seq" or "rseq". The default is "norm" |
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 CovNASde
which is a subclass of the
virtual class CovNARobust
.
Valentin Todorov valentin.todorov@chello.at
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. <doi:10.18637/jss.v032.i03>.
data(bush10)
CovNASde(bush10)
## the following four statements are equivalent
c0 <- CovNASde(bush10)
c1 <- CovNASde(bush10, nsamp=2000)
c2 <- CovNASde(bush10, control = CovControlSde(nsamp=2000))
c3 <- CovNASde(bush10, control = new("CovControlSde", nsamp=2000))
## direct specification overrides control one:
c4 <- CovNASde(bush10, nsamp=100,
control = CovControlSde(nsamp=2000))
c1
summary(c1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.