OPTIONS LS=64 PS=90 NONUMBER NODATE;
/*
FILENAME TEMP 'FITDIST.PS';
GOPTIONS DEVICE=PS
GSFNAME=TEMP GSFLEN=132
VSIZE = 5 GSFMODE=APPEND;
*/
DATA in; INPUT y @@; CARDS;
16.1507 17.6973 16.6672 20.5913 16.6893 19.8868 18.8799 
22.0835 13.7296 17.9421 15.8432 17.9852 19.3464 13.4110 
20.0157 17.7576 16.7895 20.3240 17.4261 22.1498 18.5270 
17.3749 17.3625 18.3860 16.8319 18.5711 18.0273 19.1030 
16.0399 16.1083 20.0836 18.1312 18.4339 16.9291 18.8201 
15.5290 17.8874 19.4222 16.6248 18.1778 18.0600 14.9758 
17.4088 18.7244 26.0930 14.0824 20.8100 14.2168 17.5969 
15.1043 16.4125 13.5466 22.2426 16.9249 17.8638 16.1002 
18.1734 14.7194 23.0044 17.3785 16.7997 14.7654 20.9882 
14.9023 14.5986 19.8888 21.6152 22.0045 15.7987 16.3589 
16.7591 19.7733 17.5538 16.0035 18.6896 16.5447 16.8620 
15.3421 20.4876 19.7069 15.5354 18.5424 18.5762 17.3273 
17.6809 21.0111 19.3912 17.4237 18.8838 17.4451 15.4780 
17.7717 20.8119 16.8480 16.5819 19.8021 19.6766 18.9843 
18.8448 18.3172 19.3955 18.1663 15.7902 19.1205 13.8259 
19.1926 18.5689 18.9382 16.2522 18.5664 15.3433 17.9058 
17.6603 15.2576 16.8374 19.1917 17.4229 14.4113 16.2594 
24.3697 17.7119 15.5733 20.4894 23.1986 15.1375 19.1289 
19.6562 19.4416 18.0469 16.3447 16.9120 13.0575 15.7251 
18.9501 15.3099 22.1498 22.9714 19.1280 18.6097 22.0950 
19.5696 15.2504 18.8854 18.4809 19.1902 17.5347 22.2945 
20.1644 19.5403 14.4368

SYMBOL1 VALUE=DOT WIDTH=3 L=1;

TITLE 'FITTING A DISTRIBUTION TO A DATA SET';

PROC CAPABILITY DATA=in;
     VAR y;

******************************************;
*** HISTOGRAMS -- THE OPTION SPECIFIES ***;
***  THE DISTRIBUTION TO SUPERIMPOSE   ***;
******************************************;
     HISTOGRAM y / BETA(THETA=10 SCALE=20);
     HISTOGRAM y / EXPONENTIAL(THETA=10);
     HISTOGRAM y / GAMMA(THETA=10);
     HISTOGRAM y / LOGNORMAL(THETA=10);
     HISTOGRAM y / NORMAL;
     HISTOGRAM y / WEIBULL(THETA=10);

     HISTOGRAM y / MIDPOINTS = 13 TO 26 BY 1
                   BETA(THETA=10 SCALE=20 L=1)
                   GAMMA(THETA=10 L=4)
                   LOGNORMAL(THETA=10 L=8);

*****************************************;
*** CDF PLOTS -- THE OPTION SPECIFIES ***;
***  THE DISTRIBUTION TO SUPERIMPOSE  ***;
*****************************************;
     CDFPLOT    y / BETA(THETA=10 SCALE=20);
     CDFPLOT    y / GAMMA(THETA=10);
     CDFPLOT    y / WEIBULL(THETA=10);

***********************************************;
*** QQ AND PP PLOTS -- THE OPTION SPECIFIES ***;
***     THE DISTRIBUTION TO SUPERIMPOSE     ***;
***********************************************;
     QQPLOT y / NORMAL SQUARE;
     PPPLOT y / NORMAL SQUARE;

RUN;