[R] Single Factor Anova

[christopher ciotti]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Hello all -

As I progress in R I am trying to automate functions I would have 
normally farmed out to Excel, SPSS or Statistica.  Single factor anova 
is one of them.  For example, a dataset from NIST StRD 
(http://www.itl.nist.gov/div898/strd/anova/AtmWtAg.html) has two groups:

1             2

107.8681568   107.8681079
107.8681465   107.8681344
107.8681572   107.8681513
107.8681785   107.8681197
107.8681446   107.8681604
107.8681903   107.8681385
107.8681526   107.8681642
107.8681494   107.8681365
107.8681616   107.8681151
107.8681587   107.8681082
107.8681519   107.8681517
107.8681486   107.8681448
107.8681419   107.8681198
107.8681569   107.8681482
107.8681508   107.8681334
107.8681672   107.8681609
107.8681385   107.8681101
107.8681518   107.8681512
107.8681662   107.8681469
107.8681424   107.8681360
107.8681360   107.8681254
107.8681333   107.8681261
107.8681610   107.8681450
107.8681477   107.8681368


The certified values are (I hope the spaces remain intact):

Sums of                    Mean               
Variation          df      Squares              Squares             F 
Statistic

Between Instrument  1 3.63834187500000E-09 3.63834187500000E-09 
1.59467335677930E+01
Within Instrument  46 1.04951729166667E-08 2.28155932971014E-10


Any assistance you can offer is greatly appreciated.  Thanks.


- --
chris ciotti (chris_ciotti at yahoo.com)
PGP ID: 0xE94BB3B7


-----BEGIN PGP SIGNATURE-----
Version: PGP 8.0

iQA/AwUBQG16A1kgIqbpS7O3EQI+0wCbBfmexypwUit+JQfHx/ePbq3csyMAoJp+
bvKJsTwB/sjUAf4sXWpWQieb
=w7Bq
-----END PGP SIGNATURE-----