title 'Group Comparison ';
   data graze;
      length GrazeType $ 10; 
      input GrazeType $ WtGain @@; 
      datalines;
   controlled  45   controlled  62  
   controlled  96   controlled 128  
   controlled 120   controlled  99  
   controlled  28   controlled  50  
   controlled 109   controlled 115  
   controlled  39   controlled  96  
   controlled  87   controlled 100  
   controlled  76   controlled  80  
   continuous  94   continuous  12 
   continuous  26   continuous  89     
   continuous  88   continuous  96 
   continuous  85   continuous 130          
   continuous  75   continuous  54 
   continuous 112   continuous  69     
   continuous 104   continuous  95 
   continuous  53   continuous  21 
   ;
   run;

 proc sort  data=graze;
      by GrazeType;
	  run;
  proc print data=graze;
   run;


  /* Normality Check*/
   title 'Normality Check';
   proc univariate data=graze normality plot;
    by Grazetype; 
   var WtGain;
   run;

   title '2-sample t-test';
     proc ttest data=graze;
      class GrazeType;
      var WtGain;
   run;


   /*  Shoe Material Data */

   title 'Paired Comparison (paired t-test)';
   data sole;
      input MaterialA MaterialB @@;
      datalines;
13.2	14
8.2	8.8
10.9	11.2
14.3	14.2
10.7	11.8
6.6	6.4
9.5	9.8
10.8	11.3
8.8	9.3
13.3	13.6
   ;
   run;

   proc ttest data=sole;
      paired MaterialA*MaterialB;
   run;

   /* Fuel Additive Test */

title "Fuel Additive Test" ;

PROC IMPORT OUT= WORK.raw 
            DATAFILE= "E:\STAT481-2014\FuelAdditive.txt"
            DBMS=DLM REPLACE;
            GETNAMES=yes;
RUN;

proc print data=raw;
run;

data additive;
set raw;
Ratio=Control/Test;
run;

proc sort data=additive;
by Laps;
run;

/* Normality Check*/
   title 'Normality Check';
   proc univariate data=additive normal;
   by Laps;
   var Ratio;
   QQPlot / Normal(mu=est sigma=est);
   run;

title '2-sample t-test';
     proc ttest data=additive;
      class Laps;
      var Ratio;
   run;

 /* Chisquare Test for 2x2 Contingency Table */

   data twowaycont;
   length performance $ 4 motivation $5; 
   input performance motivation Count @@;
   datalines;
high money 49 
high satis 51
marg money 53
marg satis 47
;

proc sort data=twowaycont;
   by performance;
run;

proc freq data=twowaycont;
    tables performance * motivation / chisq ;
    weight Count;
    title 'Job Performance vs Job Motivation';
run;