[R] Plotting Bivariate Normal Data
Liaw, Andy
andy_liaw at merck.com
Tue Oct 26 04:09:50 CEST 2004
> From: Berton Gunter
>
> Just a little addendum to Martin's comments below. It is well
> known that
> using LS centers and covariances for the M-distances is
> generally not a good
> way to do this, as these statistics, themselves, are
> distorted by the long
> "tails" (do > 1D distributions have "tails"?)
How about `long skirt' or some such? That seems descriptive, at least for
2D...
Andy
> so that the problems are
> hidden (see Brian Ripley's comments on the R-Help "robust
> regression with
> groups" thread from last week). Hence, one should use a
> resistant center
> (the medioid, say) and a resistant covariance matrix (e.g.,
> from cov.rob())
> to compute the M-distances.
>
> ... But then, this begs the question: Why do normality testing at all?
> (again, see BR's comments). Better to use robust/resistant statistical
> procedures for estimation from the beginning, though,
> unfortunately, this
> shatters the nice simple mathematical framework for inference.
>
> -- Bert Gunter
> Genentech Non-Clinical Statistics
> South San Francisco, CA
>
> "The business of the statistician is to catalyze the
> scientific learning
> process." - George E. P. Box
>
>
>
> > Since one of the more severe and common deviations from
> > normality is "long tailed"ness (in all it's vaguety), we have
> > been recommending to QQ-plot mahalanobis distances against chi
> > squared quantiles - even before looking at the univariate
> > QQ plots.
> >
> > Exactly for this reason, in R,
> > example(mahalanobis)
> > shows a version of how to do this!
> >
> > Martin Maechler, ETH Zurich
> >
> > ______________________________________________
> > R-help at stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide!
> > http://www.R-project.org/posting-guide.html
> >
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>
>
More information about the R-help
mailing list