Corefw: Fully Efficient Robust Correlation Estimation

Description Usage Arguments Details Value Author(s) References Examples

Description

Calculates a robust and asymptotically fully efficient correlation matrix, see Gervini (2003) for details.

Usage

1
Corefw(data, boundq = 0.975, startestimator = "wMCD")

Arguments

data

matrix with the observations in rows.

boundq

quantile bound for comparison of theoretical and empirical distribution function, see details.

startestimator

start estimator to be used: "rMCD", "wMCD" and "S"-estimator are possible, see details.

Details

This implementation calculates the (asymptotically) fully efficient scatter estimator proposed by Gervini (2003). Based on a an initial scatter estimator which can be determined by startestimator, residuals and their cumulative distribution function are calculated. This empirical distribution is compared with the theoretical one (a chi square distribution). Only if the empirical distribution function lays under the theoretic one, from the boundq quantile onwards, observations are identified as outliers, the number depends on the distance between the distribution functions. Based on all observations which are not marked as outliers, the usual correlation is returned.

This procedure has asymptotically the same efficiency as the usual empirical correlation and retains the breakdown point of the initial scale estimator.

Value

Numeric correlation matrix.

Author(s)

Alexander Dürre

References

Gervini, D. (2003): A robust and efficient adaptive reweighted estimator of multivariate location and scatter, Journal of multivariate analysis, vol 84, 116–144, doi: 10.1016/S0047-259X(02)00018-0.

Examples

1
2
3
4
5
set.seed(1066)
x <- rnorm(100)
set.seed(1485)
X <- cbind(x, x + rnorm(100), x^2)
Corefw(X)


Search within the robts package
Search all R packages, documentation and source code

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

Please suggest features or report bugs with the GitHub issue tracker.

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