___  ____  ____  ____  ____ (R)
 /__    /   ____/   /   ____/
___/   /   /___/   /   /___/   15.1   Copyright 1985-2017 StataCorp LLC
  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)

Single-user Stata perpetual license:
       Serial number:  15
         Licensed to:  Kreshna Gopal
                       StataCorp LLC

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.


running /home/krg/bin/profile.do ...
Compile number 629

. do sires.do 

. /* NIST StRD benchmark from http://www.nist.gov/itl/div898/strd/
> 
> ANOVA
> 
> Difficulty=Lower  n_i=5  k=5  Observed
> 
> Dataset Name:   Si_Resistivity     (Si_Resistivity.dat)
> 
> 
> Procedure:      Analysis of Variance
> 
> 
> Reference:      Ehrstein, James and Croarkin, M. Carroll.
>                 Unpublished NIST dataset.
> 
> 
> Data:           1 Factor
>                 5 Treatments
>                 5  Replicates/Cell
>                 25 Observations
>                 3 Constant Leading Digits
>                 Lower Level of Difficulty
>                 Observed Data
> 
> 
> Model:          6 Parameters (mu,tau_1, ... , tau_5)
>                 y_{ij} = mu + tau_i + epsilon_{ij}
> 
> 
> Certified Values:
> 
> Source of                  Sums of               Mean
> Variation          df      Squares              Squares             F Statist
> ic
> 
> Between Instrument  4 5.11462616000000E-02 1.27865654000000E-02
> 1.18046237440255E+00
> Within Instrument  20 2.16636560000000E-01 1.08318280000000E-02
> 
>                    Certified R-Squared 1.90999039051129E-01
> 
>                    Certified Residual
>                    Standard Deviation  1.04076068334656E-01
> */
. 
. clear

. 
. scalar N         = 25

. scalar df_r      = 20

. scalar df_m      = 4

. 
. scalar mss       = 5.11462616000000E-02

. scalar F         = 1.18046237440255E+00

. scalar rss       = 2.16636560000000E-01

. scalar r2        = 1.90999039051129E-01

. scalar rmse      = 1.04076068334656E-01

. 
. qui input byte instrum double resist

. 
. anova resist instrum

                         Number of obs =         25    R-squared     =  0.1910
                         Root MSE      =    .104076    Adj R-squared =  0.0292

                  Source | Partial SS         df         MS        F    Prob>F
              -----------+----------------------------------------------------
                   Model |  .05114626          4   .01278657      1.18  0.3494
                         |
                 instrum |  .05114626          4   .01278657      1.18  0.3494
                         |
                Residual |  .21663656         20   .01083183  
              -----------+----------------------------------------------------
                   Total |  .26778282         24   .01115762  

. 
. 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)                13.1
e(rmse)             13.4
e(r2)               13.2
e(mss)              14.0
e(rss)              13.1

. 
end of do-file