___ ____ ____ ____ ____ (R) /__ / ____/ / ____/ ___/ / /___/ / /___/ 16.0 Copyright 1985-2019 StataCorp LLC Statistics/Data Analysis StataCorp 4905 Lakeway Drive Special Edition College Station, Texas 77845 USA 800-STATA-PC http://www.stata.com 979-696-4600 stata@stata.com 979-696-4601 (fax) 10-user Stata network perpetual license: 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 5000; see help set_maxvar. running /home/krg/bin/profile.do ... Compile number 785 . do pidigits.do . /* NIST/ITL StRD > > Univariate Summary Statistics > > Difficulty=Lower N=5000 Observed > > File Name: pidigits.dat > > Dataset Name: Pi digits > > Description: This is a constructed/fabricated data set > to test accuracy in summary statistic calculations. > The numbers are the first 5000 digits of the > mathematical constant pi (= 3.1415926535897932384...). > > Stat Category: Univariate > > Reference: Mathematics of Computation. > January 1962, page 76. > > Data: Constructed Variable > --> 1 Response : y = pi digits > --> 0 Predictors > --> 5000 Observations > > Model: Lower Level of Difficulty > --> 2 Parameters : mu, sigma > --> 1 Response Variable : y > --> 0 Predictor Variables > > y = mu + e > > Certified Values > Sample Mean ybar: 4.53480000000000 > Sample Standard Deviation (denom. = n-1) s: 2.86733906028871 > Sample Autocorrelation Coefficient (lag 1) r(1): -0.00355099287237972 > > Number of Observations: 5000 > */ . . clear . . scalar N = 5000 . scalar mean = 4.53480000000000 . scalar sd = 2.86733906028871 . scalar rho1 = -0.00355099287237972 . . qui input byte y . . summarize y Variable | Obs Mean Std. Dev. Min Max -------------+--------------------------------------------------------- y | 5,000 4.5348 2.867339 0 9 . assert scalar(N) == r(N) . . lrecomp r(mean) mean sqrt(r(Var)) sd r(mean) sqrt(r(Var)) 15.2 . . gen time = _n . tsset time time variable: time, 1 to 5000 delta: 1 unit . corrgram y, lag(1) -1 0 1 -1 0 1 LAG AC PAC Q Prob>Q [Autocorrelation] [Partial Autocor] ------------------------------------------------------------------------------- 1 -0.0036 -0.0036 .06309 0.8017 | | . . lrecomp r(ac1) rho1 r(ac1) 14.9 . end of do-file