Fully Efficient Robust Correlation Estimation
Calculates a robust and asymptotically fully efficient correlation matrix, see Gervini (2003) for details.
Corefw(data, boundq = 0.975, startestimator = "wMCD")
matrix with the observations in rows.
quantile-bound for comparism of theoretical and empirical distribution function, see details.
start estimator to be used:
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 chisquare distribution). Only if the empirical distribution function lays under the theoretic one, from the
boundq-quantile onwords, 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.
Numeric correlation matrix.
Gervini, D. (2003): A robust and efficient adaptive reweighted estimator of multivariate location and scatter, Journal of multivariate analysis, vol 84, 116–144.
1 2 3 4 5
Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.