___  ____  ____  ____  ____ (R)
 /__    /   ____/   /   ____/
___/   /   /___/   /   /___/   11.0   Copyright 1984-2009
  Statistics/Data Analysis            StataCorp
                                      4905 Lakeway Drive
                                      College Station, Texas 77845 USA
                                      800-STATA-PC        http://www.stata.com
                                      979-696-4600        stata@stata.com
                                      979-696-4601 (fax)

3-user Stata network perpetual license:
       Serial number:  999
         Licensed to:  Brian P. Poi
                       StataCorp LP

Notes:
      1.  (-m# option or -set memory-) 10.00 MB allocated to data
      2.  Command line editing disabled
      3.  Stata running in batch mode


running /home/bpp/bin/profile.do ...

. do agatom.do 

. /* NIST StRD benchmark from http://www.nist.gov/itl/div898/strd/
> 
> ANOVA
> 
> Difficulty=Average  n_i=24  k=2  Generated
> 
> Dataset Name:   Atomic Weight of Silver   (agwt.dat)
> 
> 
> Procedure:      Analysis of Variance
> 
> 
> Reference:      Powell, L.J., Murphy, T.J. and Gramlich, J.W. (1982).
>                 "The Absolute Isotopic Abundance & Atomic Weight
>                 of a Reference Sample of Silver".
>                 NBS Journal of Research, 87, pp. 9-19.
> 
> 
> Data:           1 Factor
>                 2 Treatments
>                 24 Replicates/Cell
>                 48 Observations
>                 7 Constant Leading Digits
>                 Average Level of Difficulty
>                 Observed Data
> 
> 
> Model:          3 Parameters (mu, tau_1, tau_2)
>                 y_{ij} = mu + tau_i + epsilon_{ij}
> 
> Certified Values:
> 
> Source of                  Sums of               Mean
> Variation          df      Squares              Squares             F Statist
> ic
> 
> 
> Between Instrument  1 3.63834187500000E-09 3.63834187500000E-09
> 1.59467335677930E+01
> Within Instrument  46 1.04951729166667E-08 2.28155932971014E-10
> 
>                    Certified R-Squared 2.57426544538321E-01
> 
>                    Certified Residual
>                    Standard Deviation  1.51048314446410E-05
> */
. 
. clear

. 
. scalar N         = 48

. scalar df_r      = 46

. scalar df_m      = 1

. 
. scalar mss       = 3.63834187500000E-09

. scalar F         = 1.59467335677930E+01

. scalar rss       = 1.04951729166667E-08

. scalar r2        = 2.57426544538321E-01

. scalar rmse      = 1.51048314446410E-05

. 
. qui input byte instrum double agwt

. 
. anova agwt instrum

                           Number of obs =      48     R-squared     =  0.2574
                           Root MSE      = .000015     Adj R-squared =  0.2413

                  Source |  Partial SS    df       MS           F     Prob > F
              -----------+----------------------------------------------------
                   Model |  3.6383e-09     1  3.6383e-09      15.95     0.0002
                         |
                 instrum |  3.6383e-09     1  3.6383e-09      15.95     0.0002
                         |
                Residual |  1.0495e-08    46  2.2816e-10   
              -----------+----------------------------------------------------
                   Total |  1.4134e-08    47  3.0071e-10   

. 
. assert N    == e(N)

. assert df_r == e(df_r)

. assert df_m == e(df_m)

. 
. lrecomp e(F) F e(rmse) rmse e(r2) r2 e(mss) mss /*
> */ e(rss) rss

e(F)                10.2
e(rmse)             11.2
e(r2)               10.3
e(mss)              10.2
e(rss)              10.9

. 
end of do-file