Computes a robust multivariate location and scatter estimate with a high breakdown point, using one of the available estimators.
CovRobust(x, control, na.action = na.fail)An object derived from a CovRobust object, depending on the selected estimator.
a matrix or data frame.
a control object (S4) for one of the available control classes,
e.g. CovControlMcd-class, CovControlOgk-class,
CovControlSest-class, etc.,
containing estimation options. The class of this object defines
which estimator will be used. Alternatively a character string can be specified
which names the estimator - one of auto, sde, mcd, ogk, m, mve, sfast, surreal,
bisquare, rocke. If 'auto' is specified or the argument is missing, the
function will select the estimator (see below for details)
A function to specify the action to be taken if missing values are found. The default action is for the procedure to fail. An alternative is na.omit, which leads to rejection of cases with missing values on any required variable.
Valentin Todorov valentin.todorov@chello.at
This function simply calls the restimate method of the control object
control. If a character string naming an estimator is specified, a
new control object will be created and used (with default estimation options).
If this argument is missing or a character string
'auto' is specified, the function will select the robust estimator
according to the size of the dataset. If there are less than 1000
observations and less than 10 variables or less than 5000 observations
and less than 5 variables, Stahel-Donoho estimator will be used. Otherwise,
if there are less than 50000 observations either bisquare S-estimates
(for less than 10 variables) or Rocke type S-estimates (for 10 to 20 variables)
will be used. In both cases the S iteration starts at the initial MVE estimate.
And finally, if there are more than 50000 observations and/or more than 20 variables the Orthogonalized
Quadrant Correlation estimator (CovOgk with the corresponding parameters) is used.
Todorov V & Filzmoser P (2009), An Object Oriented Framework for Robust Multivariate Analysis. Journal of Statistical Software, 32(3), 1--47. tools:::Rd_expr_doi("10.18637/jss.v032.i03").
data(hbk)
hbk.x <- data.matrix(hbk[, 1:3])
CovRobust(hbk.x)
CovRobust(hbk.x, CovControlSest(method="bisquare"))
Run the code above in your browser using DataLab