[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