___ ____ ____ ____ ____ ® /__ / ____/ / ____/ 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 sl7.do . /* NIST/ITL StRD benchmark > > ANOVA > > Difficulty=Higher n_i=21 k=9 Generated > > Dataset Name: Simon-Lesage7 (Simon-Lesage7.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 > 13 Constant Leading Digits > Higher 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 Statisti > c > > 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 = .100003 Adj R-squared = 0.4598 Source | Partial SS df MS F Prob>F -----------+---------------------------------------------------- Model | 1.6801549 8 .21001936 21.00 0.0000 | treat | 1.6801575 8 .21001969 21.00 0.0000 | Residual | 1.8000994 180 .01000055 -----------+---------------------------------------------------- Total | 3.4802542 188 .01851199 . . 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) 4.4 e(rmse) 4.6 e(r2) 4.7 e(mss) 4.0 e(rss) 4.3 . end of do-file