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

****************************************************************************;
*** Source: 01/10/04, http://taylor.ifas.ufl.edu/marine/Thesis.pdf       ***; 
*** Sharks of the Inshore Coastal Atlantic Waters of Brevard County,     ***;
*** Florida, with Emphasis on the Spinner Shark, Carcharhinus Brevipinna ***;
*** by CRAIG W. AUBREY, M S Thesis, College of Arts and Sciences         ***;
*** University of Central Florida, Orlando, Florida, 2001                ***;
****************************************************************************;

title1 "EXST 7005 Assignment 4"; 
title2 "James P Geaghan, Section 0"; 
title3 'Atlantic Sharpnose Shark catch on the Atlantic coast of Brevard County, Fla. (1995 - 1998)';

***********************;
*** Assignment 04   ***;
*** James P Geaghan ***;
*** Section 0       ***;
***********************;

PROC IMPORT OUT= WORK.Sharks95_97 
            DATAFILE= "Sharks (1995-97).xls" 
            DBMS=EXCEL REPLACE;
     RANGE="Sheet1$"; 
     GETNAMES=YES;
     MIXED=YES; 
     SCANTEXT=YES;
     USEDATE=YES;
     SCANTIME=YES;
RUN;

*Proc contents data=Sharks95_97; run;

data Sharks98; length wt 8;
   INFILE 'Sharks (1998).csv' dlm=',' dsd missover firstobs=2;
	  input Shark_No Year Month Day Sex $ PCL FL STL Wt;
datalines; run;
;
run;

*Proc contents data=Sharks98; run;

data combined; set Sharks95_97 Sharks98; 
      label PCL = 'Precaudal Length (cm)'
            FL  = 'Fork length (cm)'
            STL = 'Stretched Total Length (cm)'
            Wt  = 'Weight (kg)';
run;

*Proc contents data=combined; run;

*proc print data=combined; run;

proc plot data=Combined; plot stl * fl = sex; 
   Title4 'Scatter plot';
run; 

proc freq data=Combined; 
   Title4 'Two-way frequencies of age and sex';
   table month * sex;
run;

proc chart data=Combined; 
   Title4 'Histogram (horizontal) of total lengths';
   Title5 'Sex specified as subgroup';
   hbar stl / midpoints=35 to 95 by 10 subgroup=sex;
run;

proc chart data=Combined; where sex eq 'Female'; 
   Title4 'Histogram (horizontal) of total lengths';
   Title5 'WHERE statement used to plot females only';
   hbar stl / midpoints=35 to 95 by 10 subgroup=sex;
run;

proc sort data=combined; by sex; run;
proc means data=combined; by sex; var stl; 
   Title4 'Proc MEANS by sex';
run;

proc means data=combined; class sex; 
   var stl; 
   Title4 'Proc MEANS classed by sex';
run;

proc means data=combined noprint; class sex; 
   var stl; 
   Title4 'Proc MEANS OUTPUT classed by sex';
   output out=next01 n=n mean=mean range=range median=median;
run;

proc print data=next01; run;

proc means data=combined noprint; by sex; 
   var stl; 
   Title4 'Proc MEANS OUTPUT classed by sex';
   output out=next01 n=n mean=mean range=range median=median;
run;

proc print data=next01; run;