___ ____ ____ ____ ____ 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 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 | .051146262 4 .012786565 1.18 0.3494 | instrum | .051146262 4 .012786565 1.18 0.3494 | Residual | .21663656 20 .010831828 -----------+---------------------------------------------------- Total | .267782822 24 .011157618 . . 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