/********************************************
 THIS COMMAND FILE SHOWS EXAMPLES OF USING
 SIMPLE DESCRIPTIVE STATISTICS IN SAS.

 FILENAME:  DESCRIPT.SAS
*********************************************/

OPTIONS FORMCHAR="|----|+|---+=|-/\<>*";
title;

data pulse;
   infile "pulse.dat";
   input pulse1 pulse2 ran smokes sex height weight activity;   
run;


PROC IMPORT OUT= WORK.PULSE 
            DATAFILE= "PULSE.XLS" 
            DBMS=EXCEL REPLACE;
     SHEET="pulse$"; 
     GETNAMES=YES;
     MIXED=NO;
     SCANTEXT=YES;
     USEDATE=YES;
     SCANTIME=YES;
RUN;

/*PROC PRINT - DIFFERENT OPTIONS*/

proc print;
run;

proc print data=pulse;
run;
proc print data=pulse(obs=6);
run;

proc print data=pulse(firstobs=82 obs=85);
run;

proc print data=pulse;
  var ran height pulse1;
run;

proc print data=pulse;
  var sex -- activity;
run;

proc print data=pulse label;
  var sex -- activity;
run;

proc print data = pulse label noobs;
   var sex -- activity;
run;



/*PROC CONTENTS*/

proc contents data=pulse;
run;

proc contents data=pulse varnum;
run;



/*PROC MEANS*/


proc means data=pulse;
run;

proc means data = pulse;
  var height weight pulse1;
run;

proc means data = pulse;
  var ran -- activity;
run;

proc means data = pulse;
    class ran;
run;

proc means data = pulse n mean min max;
    class ran activity;
run;

proc means data = pulse n mean min max skewness kurtosis;
  var pulse1 pulse2;
run;

proc means data = pulse n mean clm;
  var pulse1 pulse2;
run;

proc means data = pulse n mean clm alpha = .01;
  var pulse1 pulse2;
run;

/*PROC FREQ*/

proc freq data = pulse;
   tables ran activity smokes;
run;

proc freq data = pulse;
   tables activity  /chisq testp = (.2,.5,.3);
run;

proc freq data = pulse;
   tables sex * activity;
run;

proc freq data = pulse;
   tables sex * activity / expected norow nocol nopercent chisq;
run;

/*PROC UNIVARIATE*/

proc univariate data = pulse;
   var pulse1;
  run;

proc univariate data = pulse;
   var  pulse2;
   id ran;
run;

proc univariate data = pulse plot;
   var pulse1;
   histogram;
run;

proc univariate data = pulse normal;
   var pulse1;
   qqplot / normal(mu=est sigma=est);
run;

proc univariate data=pulse;
  class ran;
  var pulse2;
  histogram;
run;