data catch;
input tow Dwt Dno Mkwt MKno;
Diff = Dno-MKno;
if diff < -1000 then group =1;
if -1000 <= diff < -500 then group=2;
if -500 <= diff < 0 then group=3;
if 0 <= diff < 500 then group=4;
if 500 <= diff < 1000 then group=5;
if  diff  >=  1000 then group=6;
cards;
1	677.9	564.9		1744.9	1109.5
2	1811.2	1059.5		1335.5	883.5
3	1238.2	1116.8		1501.3	1061.2
4	1583.1	1309.7		702.3	655.5
5	515.6	334.2		256.3	110.2
6	1013.2	460.6		406.5	192.3
7	1271.8	496.8		913.9	348.6
8	1255.5	479.2		1370.0	576.8
9	1412.0	574.4		767.8	343.8									
10	678.9	511.5		1229.8	585.6									
11	1928.5	846.7		881.8	454.1									
12	2493.0	989.7		2491.4	782.3									
13	1679.4	489.1		1191.5	525.7									
14	526.5	300.9		1087.5	1258.6									
15	717.5	701.3		630.5	437.0									
goptions target=winprtc rotate=landscape lfactor=5;
proc gplot ;
plot MKno*Dno /haxis=axis1 vaxis=axis2;
symbol1 v=dot i=none;
axis1 label = (f=swiss h=1.5 'Drake Catch (no)') value=(f=swiss h=1.25) order = (0 to 3000 by 500);
axis2 label = (f=swiss h=1.5 a=90 r= 0 'Mary K Catch (no) ') value=(f=swiss h=1.25) order = (0 to 3000 by 500);


Proc Univariate;
var Dno;
Proc Univariate;
Var MKno;
proc univariate normal plot;
var diff;
proc print;
var Dno MKno diff;
Proc sort;
by group;
Proc Univariate noprint;
by group;
var Diff;
output out=new n=n;
proc print data=new;

run;