Constructs a covariance matrix and associated location using a variety
of (possibly robust) estimators. The returned object is suitable for use
by `plot.d.ellipse`

.

1 2 3 |

`x` |
An R |

`y` |
A numeric vector of the same length as |

`cov.method` |
A character value specifying the covariance method used. |

`scalefn` |
A function that computes univariate scale and (optionally) location estimates from a
numeric vector.
If provided, |

`locfn` |
A function that computes univariate location estimates from a numeric vector.
If used, |

`cov.control` |
A named list of arguments passed to the covariance calculation used. Note that this can
override |

`cov.dellipse`

is a wrapper for a range of covariance estimation methods found in
various packages. Its operation and defaults depend on the particular covariance
estimator specified by `cov.method`

. Details for each are as follows.

`spearman`

,`kendall`

By default, the median and mad are used as location and scale respectively, and the covariance is calculated from the product of scale estimates and the Spearman rank correlation or Kendall's tau respectively. If either

`scalefn`

or`locfn`

is supplied,`scalefn`

is used for scale estimation and`locfn`

for location. For both`spearman`

and`kendall`

,`scalefn`

is only used as a scale estimator and need not take a`mu.too`

argument.`pearson`

By default, the mean and sd are used as location and scale respectively, and the covariance is calculated from the product of scale estimates and the Pearson correlation. If either

`scalefn`

or`locfn`

is supplied,`scalefn`

is used for scale estimation and`locfn`

for location, making it possible (if not very sensible) to use a combination of robust scale or location functions with the Pearson correlation coefficient. For this case,`scalefn`

is only used as a scale estimator and need not take a`mu.too`

argument.`MCD`

,`mcd`

Both compute the Minimum Covariance Determinant (MCD) estimator, a robust multivariate location and scale estimate with a high breakdown point, via the 'Fast MCD' or 'Deterministic MCD' ("DetMcd") algorithm.

`"MCD"`

uses the implementation`covMcd`

in the robustbase package;`"mcd"`

uses`cov.mcd`

in the MASS package. Neither require or use`scalefn`

or`locfn`

. Note that these MCD implementations differ appreciably for small samples (at least to n=60).`MCD`

includes consistency and finite sample correction whereas`mcd`

apparently does not apply a finite sample correction. As a result, the`mcd`

scales can be considerably smaller for modest data set sizes.`OGK`

Computes the orthogonalized pairwise covariance matrix estimate described by Maronna and Zamar (2002), as implemented by the

`covOGK`

in the robustbase package. By default, scale and location use`scaleTau2`

from robustbase. Alternatives can be specified either by providing**both**`scalefn`

and`locfn`

or by including an argument`sigmamu`

in`cov.control`

, which is passed to`covOGK`

. See`covOGK`

for a description of`sigmamu`

. If`sigmamu`

is not present in`cov.control`

and both`scalefn`

and`locfn`

are provided, scale and location are constructed from`scalefn`

and`locfn`

. If only one of these is provided, a warning is issued and`]{scaleTau2}`

is used.`GK`

Computes a simple pairwise covariance estimate suggested by Gnanadesikan and Kettenring (1972), as implemented by the

`covGK`

in the robustbase package. By default, scale and location use`scaleTau2`

from robustbase. Alternatives can be specified either by providing`scalefn`

and`locfn`

or by including an argument`scalefn`

in`cov.control`

, which is passed to`covGK`

. See`covGK`

for a description of`scalefn`

. If`scalefn`

is not present in`cov.control`

, scale and location are constructed from`scalefn`

and`locfn`

. If`locfn`

is omitted,`scalefn`

is used if it takes an argument`mu.too`

and the median is used otherwise.`gk`

As

`GK`

, except that the variables are scaled to unit (robust) sd (using`scalefn`

) before calculating the covariance (which is then rescaled). This can prevent large scale differences from masking outliers in a variable with small scale.`mve`

Uses uses

`cov.mve`

in the MASS package, which is based on the location and covariance matrix for a minimum volume ellipsoid. The method neither requires nor uses`scalefn`

or`locfn`

.

An object of class `cov.dellipse`

, which is a list with (at least) components

- method
Character string describing method; identical to

`cov.method`

- cov
2x2 covariance matrix

- cor
2x2 correlation matrix

- center
vector (length 2) specifying centre of ellipse

- scale
vector, length 2, specifying scale estimates for each variable

- n.obs
number of points (rows) used in the covariance estimate

This list is intended to be consistent with that returned by `cov.wt`

.

Stephen L R Ellison

Maronna, R.A. and Zamar, R.H. (2002) Robust estimates of location and dispersion of
high-dimensional datasets; Technometrics **44(4)**, 307-317.

Gnanadesikan, R. and John R. Kettenring (1972) Robust estimates, residuals,
and outlier detection with multiresponse data. Biometrics **28**, 81-124

`cov.rob`

in MASS, `covMcd`

, `covOGK`

and
`covGK`

in robustbase.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ```
data(potassium)
cov.dellipse(potassium) #Defaults to Spearman rank correlation
#With different method
cov.dellipse(potassium, cov.method="OGK")
#Same as above but specifying control parameters
library(robustbase) #For scaleTau2
cov.dellipse(potassium, cov.method="OGK", cov.control=list(sigmamu=scaleTau2))
#With individually specified (mad) scale
cov.dellipse(potassium, cov.method="GK", scalefn=mad)
#Defaults to median for location because mad()
#does not accept a mu.too argument
cov.dellipse(potassium, cov.method="GK", scalefn=scaleTau2)
#Defaults to specified scalefn for location because scaleTau2
#accepts mu.too=TRUE
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.