Q2015J11

Q2015J11

data work1001;
	p=1/2;	/*確率*/
	m=36*p;  /*平均*/
	v=36*p*(1-p);  /*分散*/
	σ=sqrt(v);  /*標準偏差*/
	x=(24-0.5);
	p1=cdf('normal',x,m,σ);   /*正規分布の累積分布関数->CDF*/	
	p2=cdf('normal',25,m,σ);  /*(仮に半整数補正を行わなかった場合）正規分布の累積分布関数->CDF*/
	z=1-p1; /*問では超える確率としているので、正規分布の右側の面積の数値を求める*/
	z1=1-p2; /*(仮に半整数補正を行わなかった場合）*/
	n=36;
	d0=comb(n,0)*p**0*(1-p)**(n-0); /*0日の確率*/
	d1=comb(n,1)*p**1*(1-p)**(n-1); /*1日の確率*/
	d2=comb(n,2)*p**2*(1-p)**(n-2); /*2日の確率*/
	d3=comb(n,3)*p**3*(1-p)**(n-3); /*3日の確率*/
	d4=comb(n,4)*p**4*(1-p)**(n-4); /*4日の確率*/
	d5=comb(n,5)*p**5*(1-p)**(n-5); /*5日の確率*/
	d6=comb(n,6)*p**6*(1-p)**(n-6); /*6日の確率*/
	d7=comb(n,7)*p**7*(1-p)**(n-7); /*7日の確率*/
	d8=comb(n,8)*p**8*(1-p)**(n-8); /*8日の確率*/
	d9=comb(n,9)*p**9*(1-p)**(n-9); /*9日の確率*/
	d10=comb(n,10)*p**10*(1-p)**(n-10); /*10日の確率*/
	d11=comb(n,11)*p**11*(1-p)**(n-11); /*11日の確率*/
	d12=comb(n,12)*p**12*(1-p)**(n-12); /*12日の確率*/
	d13=comb(n,13)*p**13*(1-p)**(n-13); /*13日の確率*/
	d14=comb(n,14)*p**14*(1-p)**(n-14); /*14日の確率*/
	d15=comb(n,15)*p**15*(1-p)**(n-15); /*15日の確率*/
	d16=comb(n,16)*p**16*(1-p)**(n-16); /*16日の確率*/
	d17=comb(n,17)*p**17*(1-p)**(n-17); /*17日の確率*/
	d18=comb(n,18)*p**18*(1-p)**(n-18); /*18日の確率*/
	d19=comb(n,19)*p**19*(1-p)**(n-19); /*19日の確率*/
	d20=comb(n,20)*p**20*(1-p)**(n-20); /*20日の確率*/
	d21=comb(n,21)*p**21*(1-p)**(n-21); /*21日の確率*/
	d22=comb(n,22)*p**22*(1-p)**(n-22); /*22日の確率*/
	d23=comb(n,23)*p**23*(1-p)**(n-23); /*23日の確率*/
	d24=comb(n,24)*p**24*(1-p)**(n-24); /*24日の確率*/
	d25=comb(n,25)*p**25*(1-p)**(n-25); /*25日の確率*/
	d26=comb(n,26)*p**26*(1-p)**(n-26); /*26日の確率*/
	d27=comb(n,27)*p**27*(1-p)**(n-27); /*27日の確率*/
	d28=comb(n,28)*p**28*(1-p)**(n-28); /*28日の確率*/
	d29=comb(n,29)*p**29*(1-p)**(n-29); /*29日の確率*/
	d30=comb(n,30)*p**30*(1-p)**(n-30); /*30日の確率*/
	d31=comb(n,31)*p**31*(1-p)**(n-31); /*31日の確率*/
	d32=comb(n,32)*p**32*(1-p)**(n-32); /*32日の確率*/
	d33=comb(n,33)*p**33*(1-p)**(n-33); /*33日の確率*/
	d34=comb(n,34)*p**34*(1-p)**(n-34); /*34日の確率*/
	d35=comb(n,35)*p**35*(1-p)**(n-35); /*35日の確率*/
	d36=comb(n,36)*p**36*(1-p)**(n-36); /*36日の確率*/

	d37=d24+d25+d26+d27+d28+d29+d30+d31+d32+d33+d34+d35+d36;　/*実際に計算した場合*/
run;

proc print data=work1001;
run;