[R] R hang-up using lm

Martin Maechler maechler at stat.math.ethz.ch
Wed Apr 28 10:23:53 CEST 2004


>>>>> "kjetil" == kjetil  <kjetil at acelerate.com>
>>>>>     on Tue, 27 Apr 2004 19:19:59 -0400 writes:

    kjetil> On 27 Apr 2004 at 16:46, Raubertas, Richard wrote:
    >> Within the last few weeks, someone else reported a similar
    >> problem when using the results of tapply in a call to rlm().
    >> Note that the result of tapply is a 1D array, and it appears
    >> there is a general problem with using such a thing on the
    >> RHS in formula-based modeling functions:
    >> 
    >> set.seed(3)
    >> yy <- rnorm(20)
    >> gg <- rep(1:10, 2)
    >> y <- tapply(yy, gg, median)
    >> x <- 1:10
    >> z <- lm(y ~ x)  # OK
    >> z <- lm(x ~ y)  # crashes R
    >> 
    >> (R 1.8.1 on Windows XP Pro)
    >> 

    kjetil> What exactly do you mean by "crashes R"

    kjetil> Doing this in R1.9.0, windows XP pro, there is no indication of 
    kjetil> problems.

nor is there with 1.9.0 or R-patched on Linux,
nor with R 1.8.1 on Linux.

no warning, no error, no problem at all.
Is it really the above (reproducible, thank you!) example
that crashes your R 1.8.1 ?

Martin




More information about the R-help mailing list