资源描述:
《二项分布与Poisson分布实验课程序》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、二项分布程序1:二项分布的概率%letpai=0.8;%letn=3;dataa;dox=0to&n;p=pdf('binomial',x,&pai,&n);output;end;procprint;sump;run;程序2:二项分布的左侧累计概率dataa;pai=0.8;n=3;k=2;p_low=cdf('binomial',k,pai,n);procprint;run;程序3:二项分布的右侧累计概率dataa;pai=0.8;n=3;k=2;p_up=1-cdf('binomial',k-1,pai,n);procprint;run;程序4:二项分布的图形%letpai=
2、0.5;%letn=10;title'Binomialdistributionwithpai=’&pai‘andn=’&n;dataa;dox=0to&n;p=pdf('binomial',x,&pai,&n);output;end;procprint;sump;run;symboli=needlecolor=red;procgplot;plotp*x;run;程序5:二项分布的均数和方差%letpai=0.5;%letn=10;title‘MeanandVarianceofBinomialdistributionwithpai=’&pai‘andn=’&n';dataa;dox
3、=0to&n;p=pdf('binomial',x,&pai,&n);xx=(x-&n*&pai)*(x-&n*&pai);output;end;procmeansmean;varxxx;weightp;run;程序6:二项分布的正态近似%letpai=0.3;%letn=100;title'Binomialdistributionwithpai=’&pai‘andn=’&n;dataa;mean=&n*&pai;sigma=sqrt(&n*&pai*(1-&pai));dox=0ton;p1=pdf('binomial',x,&pai,&n);p2=pdf('normal',x
4、,mean,sigma);output;end;run;symbol1i=needlecolor=red;symbol2i=splinecolor=red;symbol3i=splinecolor=bluew=2;procgplot;plotp1*xp1*xp2*x/overlayhaxis=10to50by5;run;程序7:估计总体率的可信区间(查表法)datalimits;x=8;n=20;p=x/n;doCIlower=0.001topby0.000001;p1=1-cdf('binomial',x-1,CIlower,n);ifabs(p1-0.025)<0.00001
5、thengotoOK1;end;OK1:PUTxnpCIlowerP1;doCIupper=pto0.999by0.000001;p2=CDF('binomial',x,CIupper,n);ifabs(P2-0.025)<0.00001thengotoOK2;end;OK2:PUTxnpCIupperP2;CIlower=CIlower*100;CIupper=CIupper*100;Procprint;varxnpP1P2CIlowerCIupper;run;程序8:单样本的假设检验(直接计算概率法)%letpai=0.65;%letn=10;%letk=8;Title'Bi
6、nomialdistributionwithpai=’&pai‘andn=’&n;dataa;dox=0to&n;p=pdf('binomial',x,&pai,&n);output;end;run;symboli=needlecolor=red;procgplot;plotp*x;run;datab;seta;ifx>=&k;procprint;sump;run;Poisson分布程序1:Poisson分布的概率、均数和方差、图形%letmiu=5;title'Poissondistributionwithmiu=’&miu;dataa;dox=0to5*&miu;p=pdf(
7、'poisson',x,&miu);xx=(x-&miu)*(x-&miu);output;end;run;procprint;varxp;sump;run;procmeansmean;varxxx;weightp;run;symboli=needle;procgplot;plotp*x;run;程序2:Poisson分布的正态近似%letmiu=20;title'Poissondistributionwithmiu=’&miu;dataa;mean=&miu;sigma=sqr