[Rd] problem with plot.lm?
ligges at statistik.tu-dortmund.de
Tue Aug 25 18:26:27 CEST 2009
yes, this looks like a bug. Will investiagte.
Benilton Carvalho wrote:
> I'm bringing this issue to r-devel as my idea of sending it to r-help
> appears to be wrong.
> The following is reproducible in R-patched and R-devel (also in older
> An outlier is inserted in group E and the last plot by plot.lm suggests
> that the point belongs to D.
> It appears to me that:
> * (facval, rsp) and
> * x$xlevels[[1L]][order(sapply(split(yh,mf[,1]), mean))]
> are sorted differently.
> a <- rnorm(50,0,2)
> x <- rep(seq(10, 50, by=10), each=10)
> y <- x+a
> f <- rep(c("e","a","c","b","d"), rep(10,5))
> f <- as.factor(f)
> y <- 200
> mod1 <- aov(y ~ f)
> My latest attempt was on:
> > sessionInfo()
> R version 2.10.0 Under development (unstable) (2009-08-20 r49335)
>  LC_CTYPE=en_US.iso885915 LC_NUMERIC=C
>  LC_TIME=en_US.iso885915 LC_COLLATE=en_US.iso885915
>  LC_MONETARY=C LC_MESSAGES=en_US.iso885915
>  LC_PAPER=en_US.iso885915 LC_NAME=C
>  LC_ADDRESS=C LC_TELEPHONE=C
>  LC_MEASUREMENT=en_US.iso885915 LC_IDENTIFICATION=C
> attached base packages:
>  stats graphics grDevices utils datasets methods base
> thank you,
> R-devel at r-project.org mailing list
More information about the R-devel