___ ____ ____ ____ ____ 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 wampler3.do . /* NIST StRD benchmark from http://www.nist.gov/itl/div898/strd/ > > Linear Regression > > Difficulty=Higher Polynomial k=6 N=21 Generated > > Dataset Name: Wampler-3 (wampler3.dat) > > Procedure: Linear Least Squares Regression > > Reference: Wampler, R. H. (1970). > A Report of the Accuracy of Some Widely-Used Least > Squares Computer Programs. > Journal of the American Statistical Association, 65, pp. 549-5 > 65. > > Data: 1 Response Variable (y) > 1 Predictor Variable (x) > 21 Observations > Higher Level of Difficulty > Generated Data > > Model: Polynomial Class > 6 Parameters (B0,B1,...,B5) > > y = B0 + B1*x + B2*(x**2) + B3*(x**3)+ B4*(x**4) + B5*(x**5) > > > Certified Regression Statistics > > Standard Deviation > Parameter Estimate of Estimate > > B0 1.00000000000000 2152.32624678170 > B1 1.00000000000000 2363.55173469681 > B2 1.00000000000000 779.343524331583 > B3 1.00000000000000 101.475507550350 > B4 1.00000000000000 5.64566512170752 > B5 1.00000000000000 0.112324854679312 > > Residual > Standard Deviation 2360.14502379268 > > R-Squared 0.999995559025820 > > > Certified Analysis of Variance Table > > Source of Degrees of Sums of Mean > Variation Freedom Squares Squares F Statistic > > Regression 5 18814317208116.7 3762863441623.33 675524.458240122 > Residual 15 83554268.0000000 5570284.53333333 > */ . . clear . . scalar N = 21 . scalar df_r = 15 . scalar df_m = 5 . . scalar rmse = 2360.14502379268 . scalar r2 = 0.999995559025820 . scalar mss = 18814317208116.7 . scalar F = 675524.458240122 . scalar rss = 83554268.0000000 . . scalar b_cons = 1 . scalar se_cons = 2152.32624678170 . scalar bx1 = 1 . scalar sex1 = 2363.55173469681 . scalar bx2 = 1 . scalar sex2 = 779.343524331583 . scalar bx3 = 1 . scalar sex3 = 101.475507550350 . scalar bx4 = 1 . scalar sex4 = 5.64566512170752 . scalar bx5 = 1 . scalar sex5 = 0.112324854679312 . . qui input long y byte x1 . . gen int x2 = x1*x1 . gen long x3 = x1*x2 . gen long x4 = x1*x3 . gen long x5 = x1*x4 . . reg y x1-x5 Source | SS df MS Number of obs = 21 -------------+------------------------------ F( 5, 15) = . Model | 1.8814e+13 5 3.7629e+12 Prob > F = 0.0000 Residual | 83554268 15 5570284.53 R-squared = 1.0000 -------------+------------------------------ Adj R-squared = 1.0000 Total | 1.8814e+13 20 9.4072e+11 Root MSE = 2360.1 ------------------------------------------------------------------------------ y | Coef. Std. Err. t P>|t| [95% Conf. Interval] -------------+---------------------------------------------------------------- x1 | .9999997 2363.552 0.00 1.000 -5036.791 5038.791 x2 | 1 779.3435 0.00 0.999 -1660.131 1662.131 x3 | 1 101.4755 0.01 0.992 -215.2899 217.2899 x4 | 1 5.645665 0.18 0.862 -11.03345 13.03345 x5 | 1 .1123249 8.90 0.000 .7605852 1.239415 _cons | 1 2152.326 0.00 1.000 -4586.575 4588.575 ------------------------------------------------------------------------------ . di "R-squared = " %20.15f e(r2) R-squared = 0.999995559025820 . . assert N == e(N) . assert df_r == e(df_r) . assert df_m == e(df_m) . . lrecomp _b[_cons] b_cons _b[x1] bx1 _b[x2] bx2 /* > */ _b[x3] bx3 _b[x4] bx4 _b[x5] bx5 () /* > */ _se[_cons] se_cons _se[x1] sex1 _se[x2] sex2 /* > */ _se[x3] sex3 _se[x4] sex4 _se[x5] sex5 () /* > */ e(rmse) rmse e(r2) r2 e(mss) mss e(F) F e(rss) rss _b[_cons] 6.8 _b[x1] 6.5 _b[x2] 6.9 _b[x3] 7.7 _b[x4] 9.0 _b[x5] 10.7 ------------------------- min 6.5 _se[_cons] 11.4 _se[x1] 10.9 _se[x2] 10.8 _se[x3] 10.8 _se[x4] 10.8 _se[x5] 10.8 ------------------------- min 10.8 e(rmse) 14.3 e(r2) 16.0 e(mss) 14.7 e(F) 14.1 e(rss) 14.1 . . end of do-file