[R] Wrong results from anova

William Dunlap wdunlap at tibco.com
Tue Apr 14 19:40:32 CEST 2015


It is a printing problem - the default number of digits in print.summary.aov
is max(3L, getOption("digits") - 3L).   Set options(digits=7) instead of
your current 6 (?) or try print(summary(fit.Y), digits=7) to see more
digits.


Bill Dunlap
TIBCO Software
wdunlap tibco.com

On Tue, Apr 14, 2015 at 10:23 AM, li li <hannah.hlx at gmail.com> wrote:

> Hi all,
>   I have following data. When I perform an anova, the residual sum of
> square returns to be zero.
> But this is wrong, since we can hand calculate the RSS to be around 0.0308.
>   Did anyone come across the same problem before? Any suggestions?
>   Thanks.
>     Hanna
>
> > ydata
>         Y sample.Y
> 1   0.477      0.5
> 2   0.477      0.5
> 3   0.478      0.5
> 4  27.320       27
> 5  27.420       27
> 6  27.300       27
> 7  29.440       29
> 8  29.620       29
> 9  29.610       29
> 10 35.840       35
> 11 35.900       35
> 12 35.850       35
> > fit.Y <- aov(Y~sample.Y, data=ydata)
> > summary(fit.Y)
>             Df Sum Sq Mean Sq F value Pr(>F)
> sample.Y     3   2203   734.2  190706 <2e-16 ***
> Residuals    8      0     0.0
> ---
> Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.

	[[alternative HTML version deleted]]



More information about the R-help mailing list