aov bug? (PR#1930)
iwhite@staffmail.ed.ac.uk
iwhite@staffmail.ed.ac.uk
Thu, 22 Aug 2002 21:46:00 +0200 (MET DST)
R : Copyright 2001, The R Development Core Team
Version 1.4.0 (2001-12-19)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type `license()' or `licence()' for distribution details.
R is a collaborative project with many contributors.
Type `contributors()' for more information.
Type `demo()' for some demos, `help()' for on-line help, or
`help.start()' for a HTML browser interface to help.
Type `q()' to quit R.
> invisible(options(echo = TRUE))
> R.version
_
platform i686-pc-linux-gnu
arch i686
os linux-gnu
system i686, linux-gnu
status
major 1
minor 4.0
year 2001
month 12
day 19
language R
> r <- c(10,23,23,26,17,5,53,55,32,46,10,8,10,8,23,0,3,22,15,32,3)
> n <- c(39,62,81,51,39,6,74,72,51,79,13,16,30,28,45,4,12,41,30,51,7)
> trt <- factor(rep(1:4,c(5,6,5,5)))
> Y <- r/n
> z <- aov(Y ~ trt,weights=n)
> summary(z)
Df Sum Sq Mean Sq F value Pr(>F)
trt 3 16.0626 5.3542 12.330 0.0001577 ***
Residuals 17 7.3820 0.4342
---
Signif. codes: 0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1
> z
Call:
aov(formula = Y ~ trt, weights = n)
Terms:
trt Residuals
Sum of Squares 16.062629 0.395016
Deg. of Freedom 3 17
Residual standard error: 0.1524343
Estimated effects may be unbalanced
>
> # BUG? summary(z) correctly reports the weighted residual sum of squares,
> # but the output from z gives an unweighted residual sum of squares.
> # At the least the discrepancy between the two forms of output is
> # confusing.
> proc.time()
[1] 1.60 0.08 1.65 0.00 0.00
>
======================================
I.White
ICAPB, University of Edinburgh
Ashworth Laboratories, West Mains Road
Edinburgh EH9 3JT
Fax: 0131 667 3210 Tel: 0131 650 5490
E-mail: iwhite@staffmail.ed.ac.uk
======================================
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._