1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;7071 /* lynch.sas */7273 title 'Hovland and Sears Lynching data';7475 data lynch;76 infile '/home/brunner0/441s20/Lynch.data.txt' firstobs=6; /* Skipping the header */77 input year ayres cotton1 cotton2 Blynch Totlynch;78 label ayres = 'Ayres'' composite economic index'79 cotton1 = 'Per-acre value of cotton'80 cotton2 = 'Farm value of cotton'81 blynch = 'Negro lynchings'82 totlynch = 'Total lynchings';83 olynch = totlynch-blynch;84 label olynch = 'Non-Black lynchings';85 /* Created differenced and lagged variables after looking at86 autocorrelation and cross-correlation functions, etc. */87 dc1 = cotton1 - lag(cotton1);88 dc2 = cotton2 - lag(cotton2);89 dbl = blynch - lag(blynch);90NOTE: The infile '/home/brunner0/441s20/Lynch.data.txt' is:Filename=/home/brunner0/441s20/Lynch.data.txt,Owner Name=brunner0,Group Name=oda,Access Permission=-rw-r--r--,Last Modified=09Mar2020:12:40:13,File Size (bytes)=3083NOTE: 49 records were read from the infile '/home/brunner0/441s20/Lynch.data.txt'.The minimum record length was 54.The maximum record length was 55.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).1 at 87:20 1 at 88:20 1 at 89:19NOTE: The data set WORK.LYNCH has 49 observations and 10 variables.NOTE: DATA statement used (Total process time):real time 0.00 secondsuser cpu time 0.00 secondssystem cpu time 0.01 secondsmemory 773.78kOS Memory 33708.00kTimestamp 03/11/2020 03:19:00 PMStep Count 83 Switch Count 2Page Faults 0Page Reclaims 158Page Swaps 0Voluntary Context Switches 16Involuntary Context Switches 0Block Input Operations 0Block Output Operations 26491 proc corr;92 var year cotton1 -- olynch;93NOTE: PROCEDURE CORR used (Total process time):real time 0.09 secondsuser cpu time 0.09 secondssystem cpu time 0.00 secondsmemory 3083.93kOS Memory 34476.00kTimestamp 03/11/2020 03:19:00 PMStep Count 84 Switch Count 0Page Faults 0Page Reclaims 230Page Swaps 0Voluntary Context Switches 0Involuntary Context Switches 0Block Input Operations 0Block Output Operations 1694 proc reg plot=none;95 title2 'Naive regression';96 model blynch = cotton1;97NOTE: PROCEDURE REG used (Total process time):real time 0.04 secondsuser cpu time 0.05 secondssystem cpu time 0.00 secondsmemory 2518.65kOS Memory 36036.00kTimestamp 03/11/2020 03:19:00 PMStep Count 85 Switch Count 2Page Faults 0Page Reclaims 282Page Swaps 0Voluntary Context Switches 18Involuntary Context Switches 0Block Input Operations 0Block Output Operations 5698 proc reg plot=none;99 title2 'Naive regression, but request Durbin-Watson statistic';100 model blynch = cotton1 cotton2 / dw;101 output out=lynch2 residual=e;102NOTE: The data set WORK.LYNCH2 has 49 observations and 11 variables.NOTE: PROCEDURE REG used (Total process time):real time 0.05 secondsuser cpu time 0.05 secondssystem cpu time 0.00 secondsmemory 2668.18kOS Memory 36552.00kTimestamp 03/11/2020 03:19:00 PMStep Count 86 Switch Count 4Page Faults 0Page Reclaims 323Page Swaps 0Voluntary Context Switches 33Involuntary Context Switches 1Block Input Operations 0Block Output Operations 336103 proc sgplot;104 title2 'Residuals Against Time';105 series X = year Y = e;106107 proc autoreg;NOTE: PROCEDURE SGPLOT used (Total process time):real time 0.14 secondsuser cpu time 0.06 secondssystem cpu time 0.02 secondsmemory 12993.81kOS Memory 44972.00kTimestamp 03/11/2020 03:19:00 PMStep Count 87 Switch Count 2Page Faults 0Page Reclaims 3068Page Swaps 0Voluntary Context Switches 212Involuntary Context Switches 0Block Input Operations 0Block Output Operations 688NOTE: There were 49 observations read from the data set WORK.LYNCH2.108 title2 'Durbin-Watson test with proc autoreg';109 model blynch = cotton1 cotton2 / dwprob;110NOTE: The graph in the PDF(WEB) destination will be rendered as an image due to the use of a gradient contour.NOTE: The graph in the PDF(WEB) destination will be rendered as an image due to the use of a surface plot.NOTE: The graph in the PDF(WEB) destination will be rendered as an image due to the use of a bivariate histogram.NOTE: PROCEDURE AUTOREG used (Total process time):real time 0.39 secondsuser cpu time 0.24 secondssystem cpu time 0.02 secondsmemory 14324.56kOS Memory 58644.00kTimestamp 03/11/2020 03:19:01 PMStep Count 88 Switch Count 17Page Faults 0Page Reclaims 11034Page Swaps 0Voluntary Context Switches 653Involuntary Context Switches 1Block Input Operations 0Block Output Operations 1224111 proc autoreg;112 title2 'First-order autoregressive model with proc autoreg';113 model blynch = cotton1 cotton2 / nlag=1 method=ml;114 both: test cotton1=cotton2=0;115116 /* I don't think we need a higher order model, but let's try it and test. */117NOTE: The graph in the PDF(WEB) destination will be rendered as an image due to the use of a gradient contour.NOTE: The graph in the PDF(WEB) destination will be rendered as an image due to the use of a surface plot.NOTE: The graph in the PDF(WEB) destination will be rendered as an image due to the use of a bivariate histogram.NOTE: PROCEDURE AUTOREG used (Total process time):real time 0.52 secondsuser cpu time 0.32 secondssystem cpu time 0.03 secondsmemory 13698.90kOS Memory 58900.00kTimestamp 03/11/2020 03:19:01 PMStep Count 89 Switch Count 17Page Faults 0Page Reclaims 10515Page Swaps 0Voluntary Context Switches 831Involuntary Context Switches 1Block Input Operations 0Block Output Operations 1456118 proc autoreg plot=none;119 title2 'Sixth-order autoregressive model with cotton1 and cotton2';120 model blynch = cotton1 cotton2 / nlag=6 method=ml;121 both: test cotton1=cotton2=0;122NOTE: PROCEDURE AUTOREG used (Total process time):real time 0.28 secondsuser cpu time 0.17 secondssystem cpu time 0.02 secondsmemory 2892.37kOS Memory 49428.00kTimestamp 03/11/2020 03:19:01 PMStep Count 90 Switch Count 1Page Faults 0Page Reclaims 341Page Swaps 0Voluntary Context Switches 226Involuntary Context Switches 0Block Input Operations 0Block Output Operations 600123 proc autoreg;124 title2 'First-order autoregressive model with just cotton1';125 model blynch = cotton1 / nlag=1 method=ml;126127128 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;139