dm 'log;clear;output;clear'; 
ods html close; ods html style=minimal body='Assign06.html'; 
ods graphics on; 
options nodate nocenter pageno=1 ls=90 ps=56;
ODS listing;   

***EXST 7005***************************;
*** James P Geaghan                 ***;
*** Assignment 06                   ***;
***  Multiple Observation input and ***; 
***  Confidence Interval Example    ***;
**********************************************************************;
*** Examples are from your text book                               ***;
*** Freund, Rudolph J., William J. Wilson and Donna L. Mohr. 2010. ***;
*** Statistical Methods, Academic Press (ELSIVIER), N.Y            ***;
**********************************************************************;

title1 "Assignment 06"; 

data STEEL (keep = CLEAN TemPr020 TemPr093 TemPr165 TemPr238 TemPr310) 
     uni (keep = CLEAN TemPrID TemPr Elast) ; 
   INFILE 'steel.csv' dlm=',' dsd missover firstobs=2;
	  input CLEAN TemPr020 TemPr093 TemPr165 TemPr238 TemPr310;
	  output steel;
      TemPrID = 'A'; TemPr = 0.20; Elast = TemPr020; output uni;
      TemPrID = 'B'; TemPr = 0.93; Elast = TemPr093; output uni;
      TemPrID = 'C'; TemPr = 1.65; Elast = TemPr165; output uni;
      TemPrID = 'D'; TemPr = 2.38; Elast = TemPr238; output uni;
      TemPrID = 'E'; TemPr = 3.10; Elast = TemPr310; output uni;
datalines; run;
;
run;

title2 "Listing of Multiple observations per line"; 
proc print data=STEEL; run;
*proc print data=uni; run;

title2 "Proc means with class statement"; 
proc means data=uni n mean stderr; 
   class clean TemPr; 
   var Elast; 
run;

title2 "Proc means with BY statement and output"; 
proc sort data=uni; by clean TemPr TemPrID; run;
proc means data=uni noprint; by clean TemPr TemPrID; 
   var Elast; 
   output out = next1 n=n mean=mean var=var std=std stderr=stderror;
run;

title2 "Listing of means by CLEAN and TemPr"; 
proc print data=next1; 
run;

options ls=90 ps=48;
title2 "Plot of mean on CLEAN with TemPr identifier"; 
proc plot data=next1; 
   plot mean*clean=TemPrID;
run;

PROC UNIVARIATE data=uni CIBASIC;
   VAR elast;
   ods exclude BasicMeasures ExtremeObs ExtremeValues Modes
               MissingValues Quantiles TestsForLocation; 
RUN;