EXST7034 Multiple Regression Geaghan Chapter 11 Bootstrapping (Toluca example) Page 1 Toluca Company Example (Problem from Neter, Kutner, Nachtsheim & Wasserman 1996,1.21) A particular part needed for refigeration equipment replacement parts are produced periodically. Production requires a setup plus the time needed to produce the part. The company wants to study costs by examining 25 production runs over 3 years. 1 dm'log;clear;output;clear'; 2 ************************************************************************; 3 *** EXST7034 Example 1 using PC-SAS - Toluca Company Example ***; 4 *** Problem from Neter, Kutner, Nachtsheim & Wasserman 1996,1.21 ***; 5 ************************************************************************; 6 *** A particular part needed for refigeration equipment replacement ***; 7 *** parts are produced periodically. Production requires a setup ***; 8 *** plus the time needed to produce the part. The company wants to ***; 9 *** study costs by examining 25 production runs over 3 years. ***; 10 ************************************************************************; 11 12 ODS HTML style=minimal rs=none 13 body='C:\Geaghan\EXST\EXST7034New\Fall2002\SAS\03b-Toluca-Bootstrap.html' ; NOTE: Writing HTML Body file: C:\Geaghan\EXST\EXST7034New\Fall2002\SAS\03b-Toluca-Bootstrap.html 14 15 OPTIONS LS=155 PS=256 NOCENTER NODATE NONUMBER nolabel; 16 DATA ONE; INFILE CARDS MISSOVER; 17 TITLE1 'EXST7034 - Chapter 3 examples : Toluca example'; 18 LABEL X = 'Lot size' 19 Y = 'work hours'; 20 INPUT X Y; 21 CARDS; NOTE: The data set WORK.ONE has 25 observations and 2 variables. NOTE: DATA statement used (Total process time): real time 0.01 seconds cpu time 0.02 seconds 47 ; 48 PROC SORT DATA=ONE; BY X Y; run; NOTE: There were 25 observations read from the data set WORK.ONE. NOTE: The data set WORK.ONE has 25 observations and 2 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.02 seconds cpu time 0.02 seconds 49 proc print DATA=ONE; run; NOTE: There were 25 observations read from the data set WORK.ONE. NOTE: The PROCEDURE PRINT printed page 1. NOTE: PROCEDURE PRINT used (Total process time): real time 0.10 seconds cpu time 0.02 seconds EXST7034 - Chapter 3 examples : Toluca example Obs X Y 1 20 113 2 30 121 3 30 212 4 30 273 5 40 160 6 40 244 7 50 157 8 50 221 9 50 268 10 60 224 11 70 252 12 70 323 13 70 361 14 80 342 15 80 352 16 80 399 17 90 376 18 90 377 19 90 389 20 90 468 21 100 353 22 100 420 23 110 421 24 110 435 25 120 546EXST7034 Multiple Regression Geaghan Chapter 11 Bootstrapping (Toluca example) Page 2 51 DATA BOOTS; 52 DO Run = 1 TO 100 BY 1; 53 DO rep = 1 TO 25 BY 1; 54 get = INT(RANUNI(7882271)*25+1); 55 SET ONE POINT=get NOBS=subjects; OUTPUT; 56 END; 57 END; 58 stop; 59 RUN; NOTE: The data set WORK.BOOTS has 2500 observations and 4 variables. NOTE: DATA statement used (Total process time): real time 0.00 seconds cpu time 0.01 seconds 60 proc print data=boots; run; NOTE: There were 2500 observations read from the data set WORK.BOOTS. NOTE: The PROCEDURE PRINT printed pages 2-11. NOTE: PROCEDURE PRINT used (Total process time): real time 0.70 seconds cpu time 0.55 seconds 61 62 PROC SORT DATA=BOOTS; BY run rep; run; NOTE: There were 2500 observations read from the data set WORK.BOOTS. NOTE: The data set WORK.BOOTS has 2500 observations and 4 variables. NOTE: PROCEDURE SORT used (Total process time): real time 0.00 seconds cpu time 0.01 seconds EXST7034 - Chapter 3 examples : Toluca example Obs Run rep X Y 1 1 1 30 212 2 1 2 30 273 3 1 3 80 342 4 1 4 70 252 5 1 5 50 157 6 1 6 90 389 7 1 7 120 546 8 1 8 50 157 9 1 9 20 113 10 1 10 80 352 11 1 11 20 113 12 1 12 70 252 13 1 13 40 244 14 1 14 30 121 15 1 15 60 224 16 1 16 60 224 17 1 17 120 546 18 1 18 50 268 19 1 19 50 268 20 1 20 90 376 21 1 21 20 113 22 1 22 90 376 23 1 23 110 421 24 1 24 30 212 25 1 25 50 157 26 2 1 50 268 27 2 2 30 121 28 2 3 50 221 29 2 4 70 323 30 2 5 80 342 31 2 6 40 244 32 2 7 50 221 33 2 8 110 421 34 2 9 20 113 35 2 10 80 352 36 2 11 20 113 37 2 12 70 252 38 2 13 90 468 39 2 14 30 212 40 2 15 110 435 41 2 16 100 420 42 2 17 90 468 43 2 18 30 121 44 2 19 90 376 45 2 20 30 212 46 2 21 70 252 47 2 22 50 221 48 2 23 70 252 49 2 24 30 273 50 2 25 90 468 51 3 1 70 361 52 3 2 110 421 53 3 3 70 252 54 3 4 30 273 55 3 5 70 361 56 3 6 110 435 57 3 7 60 224 58 3 8 110 421 59 3 9 120 546 60 3 10 70 323 61 3 11 100 353 62 3 12 70 252 63 3 13 90 376 64 3 14 70 252 65 3 15 30 273 66 3 16 50 157 67 3 17