gregor.gorjanc at bfro.uni-lj.si
2006-Aug-30 23:41 UTC
[Rd] Summary() fails after use of na.action="na.exclude" in lm() (PR#9191)
Hello! I have already reported on this behaviour at [1] and now I have checked the same example on r-devel and the problem is still there. Thomas Lumley inidicated problem [2]. Example: n <- 50 x <- runif(n=n) y1 <- 2 * x + rnorm(n=n) y2 <- 5 * x + rnorm(n=n) y2[sample(1:n, size=5)] <- NA y <- cbind(y1, y2) ## Goes ok here fit <- lm(y1 ~ 1, na.action="na.exclude") summary(fit) ## And here fit <- lm(y2 ~ 1, na.action="na.exclude") summary(fit) ## But fails if we use matrix of response variables fit <- lm(y ~ 1, na.action="na.exclude") summary(fit) Response y1 : Call: lm(formula = y1 ~ 1, na.action = "na.exclude") Residuals: Error in quantile.default(resid) : missing values and NaN's not allowed if 'na.rm' is FALSE [1]http://tolstoy.newcastle.edu.au/R/devel/06/07/6187.html [2]http://tolstoy.newcastle.edu.au/R/devel/06/07/6203.html --please do not edit the information below-- Version: platform = i686-pc-linux-gnu arch = i686 os = linux-gnu system = i686, linux-gnu status = Under development (unstable) major = 2 minor = 4.0 year = 2006 month = 08 day = 30 svn rev = 39022 language = R version.string = R version 2.4.0 Under development (unstable) (2006-08-30 r39022) Locale: LC_CTYPE=en_GB.UTF-8;LC_NUMERIC=C;LC_TIME=en_GB.UTF-8;LC_COLLATE=en_GB.UTF-8;LC_MONETARY=en_GB.UTF-8;LC_MESSAGES=en_GB.UTF-8;LC_PAPER=en_GB.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_GB.UTF-8;LC_IDENTIFICATION=C Search Path: .GlobalEnv, package:methods, package:stats, package:graphics, package:grDevices, package:utils, package:datasets, Autoloads, package:base -- Lep pozdrav / With regards, Gregor Gorjanc ---------------------------------------------------------------------- University of Ljubljana PhD student Biotechnical Faculty Zootechnical Department URI: http://www.bfro.uni-lj.si/MR/ggorjan Groblje 3 mail: gregor.gorjanc <at> bfro.uni-lj.si SI-1230 Domzale tel: +386 (0)1 72 17 861 Slovenia, Europe fax: +386 (0)1 72 17 888 ---------------------------------------------------------------------- "One must learn by doing the thing; for though you think you know it, you have no certainty until you try." Sophocles ~ 450 B.C.
Apparently Analagous Threads
- Bug?: summary() fails after use of na.action="na.exclude" in lm()
- Global setting for na.rm=TRUE
- Fwd: Extract just some fields from XML]
- [Fwd: Re: levels for list and data.frame]
- Re: [Rd] corrupt data frame: columns will be truncated or padded with NAs in: format.data.frame(x, digits = digits)