___  ____  ____  ____  ____ ®
 /__    /   ____/   /   ____/      18.0
___/   /   /___/   /   /___/       SE—Standard Edition

 Statistics and Data Science       Copyright 1985-2023 StataCorp LLC
                                   StataCorp
                                   4905 Lakeway Drive
                                   College Station, Texas 77845 USA
                                   800-STATA-PC        https://www.stata.com
                                   979-696-4600        stata@stata.com

Stata license: 10-user network perpetual
Serial number: 1
  Licensed to: Stata Developer
               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 5,000 but can be increased;
          see help set_maxvar.


Running /home/krg/bin/profile.do ...
Compile number 180110

. do agatom.do 

. /* NIST/ITL 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.638e-09          1   3.638e-09     15.95  0.0002
                         |
                 instrum |  3.638e-09          1   3.638e-09     15.95  0.0002
                         |
                Residual |  1.050e-08         46   2.282e-10  
              -----------+----------------------------------------------------
                   Total |  1.413e-08         47   3.007e-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