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

10-user Stata network perpetual license:
       Serial number:  1
         Licensed to:  Stata Developer
                       StataCorp LP

Notes:
      1.  Stata is running in batch mode.
      2.  Unicode is supported; see help unicode_advice.
      3.  Maximum number of variables is set to 5000; see help set_maxvar.

. do sl4.do 

. /* NIST StRD benchmark from http://www.nist.gov/itl/div898/strd/
> 
> ANOVA
> 
> Difficulty=Average  n_i=21  k=9  Generated
> 
> Dataset Name:   Simon-Lesage4   (Simon-Lesage4.dat)
> 
> 
> Procedure:      Analysis of Variance
> 
> 
> Reference:      Simon, Stephen D. and Lesage, James P. (1989).
>                 "Assessing the Accuracy of ANOVA Calculations in
>                 Statistical Software".
>                 Computational Statistics & Data Analysis, 8, pp. 325-332.
> 
> 
> Data:           1 Factor
>                 9 Treatments
>                 21 Replicates/Cell
>                 189 Observations
>                 7 Constant Leading Digits
>                 Average Level of Difficulty
>                 Generated Data
> 
> 
> Model:          10 Parameters (mu,tau_1, ... , tau_9)
>                 y_{ij} = mu + tau_i + epsilon_{ij}
> 
> 
> Certified Values:
> 
> Source of                  Sums of               Mean
> Variation          df      Squares              Squares             F Statist
> ic
> 
> Between Treatment   8 1.68000000000000E+00 2.10000000000000E-01
> 2.10000000000000E+01
> Within Treatment  180 1.80000000000000E+00 1.00000000000000E-02
> 
>                   Certified R-Squared 4.82758620689655E-01
> 
>                   Certified Residual
>                   Standard Deviation  1.00000000000000E-01
> */
. 
. clear

. 
. scalar N         = 189

. scalar df_r      = 180

. scalar df_m      = 8

. 
. scalar mss       = 1.68

. scalar F         = 21

. scalar rss       = 1.8

. scalar r2        = 4.82758620689655E-01

. scalar rmse      = 0.1

. 
. qui input byte treat double resp

. 
. anova resp treat

                         Number of obs =        189    R-squared     =  0.4828
                         Root MSE      =         .1    Adj R-squared =  0.4598

                  Source | Partial SS         df         MS        F    Prob>F
              -----------+----------------------------------------------------
                   Model |       1.68          8         .21     21.00  0.0000
                         |
                   treat |       1.68          8         .21     21.00  0.0000
                         |
                Residual |        1.8        180         .01  
              -----------+----------------------------------------------------
                   Total |       3.48        188   .01851064  

. 
. 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.4
e(rmse)             10.6
e(r2)               10.7
e(mss)              10.1
e(rss)              10.3

. 
end of do-file