NOTE: Copyright (c) 1999-2001 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) Proprietary Software Release 8.2 (TS2M0) Licensed to GEORGE WASHINGTON UNIVERSITY-CAMPUSWIDE-T/R, Site 0001180006. NOTE: This session is executing on the WIN_ME platform. NOTE: SAS initialization used: real time 0.76 seconds 1 ****** ST210B rev 2003; 2 3 options formdlim='-'; 4 5 filename sb_lst 'c:\sasjobs\stat210\jl2003\sb_lst.pdf'; 6 ods printer pdf file=sb_lst; NOTE: Writing ODS PRINTER output to DISK destination "c:\sasjobs\stat210\jl2003\sb_lst.pdf", printer "PDF". 7 8 TITLE1 'FITNESS DATA EXAMPLE FROM SAS MANUALS'; 9 10 DATA FITNESS; 11 INPUT AGE WEIGHT OXY RUNTIME RSTPULSE RUNPULSE MAXPULSE; 12 CARDS; NOTE: The data set WORK.FITNESS has 31 observations and 7 variables. NOTE: DATA statement used: real time 0.11 seconds 44 ; 45 46 PROC REG DATA=FITNESS; 47 MODEL OXY=RUNTIME AGE WEIGHT RUNPULSE MAXPULSE RSTPULSE 48 / SELECTION=BACKWARD SLS=0.05; 49 TITLE2 'BACKWARDS ELIMINATION VARIABLE SELECTION'; 50 NOTE: 31 observations read. NOTE: 31 observations used in computations. NOTE: PROCEDURE REG used: real time 0.76 seconds 51 PROC REG DATA=FITNESS; 52 MODEL OXY=RUNTIME AGE WEIGHT RUNPULSE MAXPULSE RSTPULSE 53 / SELECTION=STEPWISE SLE=0.05 SLS=0.05; 54 TITLE2 'STEPWISE VARIABLE SELECTION'; 55 NOTE: 31 observations read. NOTE: 31 observations used in computations. NOTE: PROCEDURE REG used: real time 0.17 seconds 56 PROC REG DATA=FITNESS; 57 MODEL OXY=RUNTIME AGE WEIGHT RUNPULSE MAXPULSE RSTPULSE 58 / SELECTION=RSQUARE ADJRSQ CP MSE SSE B; 59 TITLE2 'ALL VARIABLE SUBSETS REGRESSIONS'; 60 NOTE: 31 observations read. NOTE: 31 observations used in computations. NOTE: PROCEDURE REG used: real time 0.21 seconds 61 PROC REG; 62 MODEL OXY=RUNTIME AGE WEIGHT RUNPULSE MAXPULSE/scorr2 influence; 63 TITLE2 'REDUCED MODEL EXCLUDING RSTPULSE, ALL OBSERVATIONS'; 64 NOTE: 31 observations read. NOTE: 31 observations used in computations. NOTE: PROCEDURE REG used: real time 0.38 seconds 65 DATA SPLIT; SET FITNESS; 66 U=UNIFORM(44217); 67 BY=.; AY=.; 68 IF U<=0.5 THEN BY=OXY; 69 IF U> 0.5 THEN AY=OXY; 70 NOTE: There were 31 observations read from the data set WORK.FITNESS. NOTE: The data set WORK.SPLIT has 31 observations and 10 variables. NOTE: DATA statement used: real time 0.05 seconds 71 PROC PRINT; VAR OXY AY BY RUNTIME AGE; 72 title1 'cross validation data sets'; 73 NOTE: There were 31 observations read from the data set WORK.SPLIT. NOTE: PROCEDURE PRINT used: real time 0.04 seconds 74 PROC REG; 75 MODEL AY=RUNTIME AGE WEIGHT RUNPULSE MAXPULSE /scorr2; 76 OUTPUT OUT=B PREDICTED=PBY; 77 TITLE2 'REDUCED MODEL, A TO B SPLIT-HALF VALIDATION'; NOTE: 31 observations read. NOTE: 12 observations have missing values. NOTE: 19 observations used in computations. NOTE: The data set WORK.B has 31 observations and 11 variables. NOTE: PROCEDURE REG used: real time 0.17 seconds 78 DATA C; SET B; 79 RBY=BY-PBY; NOTE: Missing values were generated as a result of performing an operation on missing values. Each place is given by: (Number of times) at (Line):(Column). 19 at 79:11 NOTE: There were 31 observations read from the data set WORK.B. NOTE: The data set WORK.C has 31 observations and 12 variables. NOTE: DATA statement used: real time 0.10 seconds 80 proc print; var by pby rby; NOTE: There were 31 observations read from the data set WORK.C. NOTE: PROCEDURE PRINT used: real time 0.00 seconds 81 PROC UNIVARIATE FREQ NORMAL PLOT; VAR BY RBY; 82 NOTE: PROCEDURE UNIVARIATE used: real time 0.27 seconds 83 PROC REG; 84 MODEL BY=RUNTIME AGE WEIGHT RUNPULSE MAXPULSE /scorr2; 85 OUTPUT OUT=A PREDICTED=PAY; 86 TITLE2 'REDUCED MODEL, B TO A SPLIT-HALF VALIDATION'; NOTE: 31 observations read. NOTE: 19 observations have missing values. NOTE: 12 observations used in computations. NOTE: The data set WORK.A has 31 observations and 13 variables. NOTE: PROCEDURE REG used: real time 0.22 seconds 87 DATA C; SET A; 88 RAY=AY-PAY; NOTE: Missing values were generated as a result of performing an operation on missing values. Each place is given by: (Number of times) at (Line):(Column). 12 at 88:11 NOTE: There were 31 observations read from the data set WORK.A. NOTE: The data set WORK.C has 31 observations and 14 variables. NOTE: DATA statement used: real time 0.10 seconds 89 proc print; var ay pay ray; NOTE: There were 31 observations read from the data set WORK.C. NOTE: PROCEDURE PRINT used: real time 0.00 seconds 90 PROC UNIVARIATE FREQ NORMAL PLOT; VAR AY RAY; 91 NOTE: PROCEDURE UNIVARIATE used: real time 0.26 seconds 92 PROC STEPWISE; 93 MODEL AY=RUNTIME AGE WEIGHT RUNPULSE MAXPULSE RSTPULSE 94 / BACKWARD SLS=0.10; 95 TITLE2 'STEPWISE WITHIN A-HALF OF OBSERVATIONS'; 96 NOTE: 31 observations read. NOTE: 12 observations have missing values. NOTE: 19 observations used in computations. NOTE: PROCEDURE STEPWISE used: real time 0.22 seconds 97 PROC STEPWISE; 98 MODEL BY=RUNTIME AGE WEIGHT RUNPULSE MAXPULSE RSTPULSE 99 / BACKWARD SLS=0.10; 100 TITLE2 'STEPWISE WITHIN B-HALF OF OBSERVATIONS'; 101 102 RUN; NOTE: 31 observations read. NOTE: 19 observations have missing values. NOTE: 12 observations used in computations. NOTE: PROCEDURE STEPWISE used: real time 0.26 seconds 103 104 ods printer close; NOTE: ODS PRINTER printed 36 pages to c:\sasjobs\stat210\jl2003\sb_lst.pdf. 105 RUN;