___ ____ ____ ____ ____ tm /__ / ____/ / ____/ ___/ / /___/ / /___/ 10.0 Copyright 1984-2007 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 for Linux64 (network) perpetual license: Serial number: 999 Licensed to: Brian P. Poi, PhD StataCorp LP Notes: 1. (-m# option or -set memory-) 1.00 MB allocated to data 2. Command line editing disabled 3. Stata running in batch mode running /home/bpp/bin/profile.do ... . do sl8.do . /* NIST StRD benchmark from http://www.nist.gov/itl/div898/strd/ > > ANOVA > > Difficulty=Higher n_i=201 k=9 Generated > > Dataset Name: Simon-Lesage8 (Simon-Lesage8.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 > 201 Replicates/Cell > 1809 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 Statis > tic > > Between Treatment 8 1.60800000000000E+01 2.01000000000000E+00 > 2.01000000000000E+02 > Within Treatment 1800 1.80000000000000E+01 1.00000000000000E-02 > > Certified R-Squared 4.71830985915493E-01 > > Certified Residual > Standard Deviation 1.00000000000000E-01 > */ . . clear . . scalar N = 1809 . scalar df_r = 1800 . scalar df_m = 8 . . scalar mss = 16.08 . scalar F = 201 . scalar rss = 18 . scalar r2 = 4.71830985915493E-01 . scalar rmse = 0.1 . . qui input byte treat double resp . . anova resp treat Number of obs = 1809 R-squared = 0.4718 Root MSE = .100003 Adj R-squared = 0.4695 Source | Partial SS df MS F Prob > F -----------+---------------------------------------------------- Model | 16.0819143 8 2.01023929 201.01 0.0000 | treat | 16.0819143 8 2.01023929 201.01 0.0000 | Residual | 18.0009782 1800 .010000543 -----------+---------------------------------------------------- Total | 34.0828925 1808 .018851157 . . 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.2 e(rmse) 4.6 e(r2) 4.5 e(mss) 3.9 e(rss) 4.3 . end of do-file