1 The SAS System 09:12 Wednesday, March 24, 2004 NOTE: Copyright (c) 1999 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) Proprietary Software Version 8 (TS M0) Licensed to UNIVERSITY OF FLORIDA STATISTICS DEPARTMENT, Site 0007029011. NOTE: This session is executing on the SunOS 5.6 platform. This message is contained in the SAS news file, and is presented upon initialization. Edit the files "news" in the "misc/base" directory to display site-specific news and information in the program log. The command line option "-nonews" will prevent this display. NOTE: SAS initialization used: real time 0.08 seconds cpu time 0.04 seconds 1 /* Report on the Best Method of Proportioning the Excise upon Spirituous Liquors 2 by Charles Blagden 3 Philosophical Transactions of the Royal Society of London, Vol 80 (1790),pp321-345 4 */ 5 6 options nodate nonumber ps=54 ls=80; 7 8 data onea; 9 do temp=30 to 100 by 5; 10 do dilute=0 to 50 by 5; 11 input weight @@; 12 weight=weight+2000; 13 dilute0=0; 14 if dilute=0 then dilute0=1; 15 output; 16 end; end; 17 cards; NOTE: SAS went to a new line when INPUT statement reached past the end of a line. NOTE: The data set WORK.ONEA has 165 observations and 4 variables. NOTE: DATA statement used: real time 0.01 seconds cpu time 0.02 seconds 33 ; 34 run; 35 36 37 data oneb; 38 do temp=30 to 100 by 5; 39 do dilute=55 to 100 by 5; The SAS System 40 input weight @@; 41 weight=weight+2000; 42 dilute0=0; 43 output; 44 end; end; 45 cards; NOTE: SAS went to a new line when INPUT statement reached past the end of a line. NOTE: The data set WORK.ONEB has 150 observations and 4 variables. NOTE: DATA statement used: real time 0.00 seconds cpu time 0.01 seconds 61 ; 62 run; 63 64 data one; 65 set onea oneb; 66 run; NOTE: There were 165 observations read from the dataset WORK.ONEA. NOTE: There were 150 observations read from the dataset WORK.ONEB. NOTE: The data set WORK.ONE has 315 observations and 4 variables. NOTE: DATA statement used: real time 0.00 seconds cpu time 0.01 seconds 67 68 proc glm; 69 model weight=temp dilute temp*dilute temp*temp dilute*dilute; 70 output out=sprtout p=yhat r=e; 71 run; 72 NOTE: There were 315 observations read from the dataset WORK.ONE. NOTE: The data set WORK.SPRTOUT has 315 observations and 6 variables. NOTE: The PROCEDURE GLM printed pages 1-2. NOTE: PROCEDURE GLM used: real time 0.06 seconds cpu time 0.04 seconds 73 proc univariate normal; 74 var e; 75 run; NOTE: There were 315 observations read from the dataset WORK.SPRTOUT. NOTE: The PROCEDURE UNIVARIATE printed pages 3-4. The SAS System NOTE: PROCEDURE UNIVARIATE used: real time 0.02 seconds cpu time 0.01 seconds 76 77 proc print; 78 run; NOTE: There were 315 observations read from the dataset WORK.SPRTOUT. NOTE: The PROCEDURE PRINT printed pages 5-11. NOTE: PROCEDURE PRINT used: real time 0.01 seconds cpu time 0.01 seconds 79 80 proc glm data=one; 81 model weight=temp dilute temp*dilute temp*temp dilute*dilute dilute0; 82 output out=sprtout0 p=yhat0 r=e0; 83 run; 84 NOTE: There were 315 observations read from the dataset WORK.ONE. NOTE: The data set WORK.SPRTOUT0 has 315 observations and 6 variables. NOTE: The PROCEDURE GLM printed pages 12-13. NOTE: PROCEDURE GLM used: real time 0.04 seconds cpu time 0.03 seconds 85 proc univariate normal; 86 var e0; 87 run; NOTE: There were 315 observations read from the dataset WORK.SPRTOUT0. NOTE: The PROCEDURE UNIVARIATE printed pages 14-15. NOTE: PROCEDURE UNIVARIATE used: real time 0.02 seconds cpu time 0.01 seconds 88 89 proc print; 90 quit; NOTE: There were 315 observations read from the dataset WORK.SPRTOUT0. NOTE: The PROCEDURE PRINT printed pages 16-22. NOTE: PROCEDURE PRINT used: real time 0.01 seconds cpu time 0.02 seconds The SAS System 91 /* 92 symbol1 c=black v=:; 93 proc gplot; 94 plot e*yhat / vref=0 frame; 95 run; 96 97 proc gchart; 98 vbar e / space=0; 99 run; 100 101 quit; 102 quit; NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414 NOTE: The SAS System used: real time 0.81 seconds cpu time 0.21 seconds