* Read in the data;

PROC IMPORT DATAFILE = 'C:\Documents and Settings\richc\My Documents\STA580S11\FEV.xls'
     OUT = TwoSamp   DBMS = EXCEL REPLACE;
     SHEET = Sheet1;
     GETNAMES = YES;
RUN;


* Start writing to PDF;

ODS PDF FILE = 'C:\Documents and Settings\richc\My Documents\STA580S11\FEV.pdf';

 
* The following code will carry out inferential tasks concerning
* \mu_1, \mu_2, \sigma^2_1, and \sigma^2_2. ;

* Chapter 8 tasks: 
* 	Test  H_0: \mu_1 = \mu_2  vs.  H_1: \mu_1 \ne \mu_2
*		assuming normality and equal variances.
*	Test  H_0: \mu_1 = \mu_2  vs.  H_1: \mu_1 \ne \mu_2
*		assuming normality and unequal variances.
*	Test  H_0: \sigma^2_1 = \sigma^2_2  vs.  H_1: \sigma^2_1 \ne \sigma^2_2
*		assuming normality. ;

proc ttest data = TwoSamp h0=0 alpha = 0.05;
class smoke;
var FEV;
run;


* The following code will carry out inferential tasks concerning
* p_1 and p_2. ;

* Chapter 10 tasks: 
*	Test  H_0: p_1 = p_2  vs.  H_1: p_1 \ne p_2
*		with a large sample.
*	Test  H_0: p_1 = p_2  vs.  H_1: p_1 \ne p_2
*		with a small sample.
*	Construct a 95% confidence interval for  p_1 - p_2
*		with a large sample.

proc freq data=TwoSamp;
   tables smoke*FEVge3 / chisq riskdiff alpha=0.05 ;
run;


* Stop writing to PDF;

ODS PDF CLOSE;
RUN;