[R] Results of MCD estimators in MASS and rrcov

rainer grohmann rainer.grohmann at gmx.net
Mon Jan 31 15:14:45 CET 2005


Thanks a lot!

Indeed, both implementations agree on the 'best' points. Your answer helped
me a great deal.

Rainer

> The two implementations use different consistency factors as well as
> different small sample correction factors.
> 
> 1. The search parts of both implementations produce the same result -
> compare rrcov.mcd$best and mass.mcd$best.
> 
> 2. The raw MCD covariance matrix is corrected as follows:
> 
> MASS:
>  - Rousseeuw and Leroy (1987), p.259 (eq. 1.26)
>  - Marazzi (1993) (or may be Rousseeuw and van Zomeren (1900) p.638 (eq 
> A.9)
> 
> rrcov:
>  - Croux and Haesbroeck (1999), Pison et.al. p. 337
>  - Pison et.al. (2002), p.338
> 
> 3. The reweighted (final) covariance matrix is corrected as follows:
> 
> MASS: no correction
> rrcov: Pison et.al. (2002) p. 339
> 
> This explains the different covariance matrices.
> As far as the location is concerned, in this particular case the raw MCD
> estimates in MASS identify one additional outlier - observation 53, which 
> is
> discarded from the computation of the reweighted estimates.
> Look at the following plots and judge yourself if this is an outlier or 
> not:
> 
>   covPlot(hbk, mcd=rrcov.mcd, which="distance", id.n=15)
>   covPlot(hbk, mcd=mass.mcd, which="distance", id.n=15)
> 
> valentin
> 

-- 
GMX im TV ... Die Gedanken sind frei ... Schon gesehen?
Jetzt Spot online ansehen: http://www.gmx.net/de/go/tv-spot




More information about the R-help mailing list