[Rd] [Package car/data.ellipse]: confidence intervals off
by factor sqrt(2)??? (PR#2584)
John Fox
jfox@mcmaster.ca
Wed Feb 26 22:31:03 2003
Dear Volker,
If the data ellipse (or, in this case, circle) is scaled so that its
shadows (projections) on the axes each includes 68% of the data (that is of
the marginal distribution of each variable), then the ellipse will include
less than 68% of the data (i.e., of the joint distribution of the two
variables). Conversely, to include 68% of the data in the ellipse, the
shadows of the ellipse have to be larger.
Did I understand your point correctly?
John
At 09:40 PM 2/26/2003 +0100, volker.franz@tuebingen.mpg.de wrote:
>Full_Name: Volker Franz
>Version: Version 1.6.2 (2003-01-10)
>OS: Debian
>Submission from: (NULL) (192.124.28.104)
>
>
>Hi there,
>
>it seems to me that data.ellipse of package "car" (Version 1.0-1)
>produces confidence interval's which are too big. To see this, do:
>
>library(car)
>plot(c(-2,2),c(-2,2),pch=0)
>data.ellipse(rnorm(10000),rnorm(10000),levels=0.68,plot.points=F)
>abline(v=+1)
>abline(v=-1)
>abline(h=+1)
>abline(h=-1)
>
>To my knowledge, this should result in a circle with radius
>1. However, the circle is larger. It seems that the problem is due to
>an erroneous specification of the degrees of freedom and can be fixed
>with the following patch:
>
>======================================================================
>--- Ellipse.R Wed Feb 26 17:49:43 2003
>+++ Ellipse.orig Thu Sep 19 18:20:41 2002
>@@ -34,7 +34,7 @@
> stop("x and y must be vectors of the same length")
> if (plot.points & !add) plot(x, y, xlab=xlab, ylab=ylab, col=col,
> pch=pch,
>las=las, ...)
> if (plot.points & add) points(x, y, col=col, pch=pch, ...)
>- dfn<-1
>+ dfn<-2
> dfd<-length(x)-1
> if (robust) {
> require(MASS)
>======================================================================
>
>Or --- am I totally on the wrong track here?
>
>Best
>Volker
>
>______________________________________________
>R-devel@stat.math.ethz.ch mailing list
>http://www.stat.math.ethz.ch/mailman/listinfo/r-devel
____________________________
John Fox
Department of Sociology
McMaster University
email: jfox@mcmaster.ca
web: http://www.socsci.mcmaster.ca/jfox