option ps=70 ls=100 nodate nonumber;
data bin;
underb26=probbnml(0.75 ,36 ,26);
overb27=1-underb26;
LABEL
underb26=’P(X<=26)'
overb27=’P(X>=27)’;
run;
proc print data=bin noobs label;
run;

data T;
undert2=probt(2*(sqrt(40)), 39);
undert1=probt(-2*(sqrt(40)), 39);
bt12=undert2-undert1;
label
undert2=’P(X-mu<=2sigma)'
undert1=’P(X-mu<=-2sigma)'
bt12=’P(-2sigma<=X-mu<=2sigma)';
run;
proc print data=T noobs label;
run;


data F;
underf2=probf(2*3/2, 30-1, 40-1);
underf1=probf(3/2, 30-1, 40-1);
bf12=underf2-underf1;
label
underf2=’P(SA/SB<=2)'
underf1=’P(SA/SB<=1)'
bf12=’P(1<=SA/SB<=2)';
run;
proc print data=f noobs label;
run;

data univ;
input score @@;
cards;
3.4 3.7 3.2 2.4 2.7 1.6 3.1
1.3 3.9 3.1 3.3 1.6 3.8 2.3
2.2 2.4 2.0 3.6 4.4 4.0 2.9
1.8 4.0 4.6 1.7 3.0 2.8 2.5
;
run;
proc univariate data=univ plot;
var score;
run;