___ ____ ____ ____ ____ (R)
/__ / ____/ / ____/
___/ / /___/ / /___/ 11.0 Copyright 1984-2009
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 network perpetual license:
Serial number: 999
Licensed to: Brian P. Poi
StataCorp LP
Notes:
1. (-m# option or -set memory-) 10.00 MB allocated to data
2. Command line editing disabled
3. Stata running in batch mode
running /home/bpp/bin/profile.do ...
. do sl7.do
. /* NIST StRD benchmark from http://www.nist.gov/itl/div898/strd/
>
> 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.68015611 8 .210019514 21.00 0.0000
|
treat | 1.68015627 8 .210019534 21.00 0.0000
|
Residual | 1.80009812 180 .010000545
-----------+----------------------------------------------------
Total | 3.48025423 188 .018511991
.
. 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