如何用SAS拟合Beta分布

小弟想请教诸位大神一个问题，假如我有100个数据，如何在SAS里对这100个数据用Beta分布进行拟合呢，并且求出Beta分布的分位点。请诸位大神不吝赐教，谢谢！！！

1. /*贝塔分布概率密度曲线*/
   
2. data beta;
   
3. do x=0 to 1 by 0.01;
   
4. y1=x**(0.5-1)*(1-x)**(0.5-1)/(gamma(0.5)*gamma(0.5)/gamma(1));
   
5. y2=x**(5-1)*(1-x)**(1-1)/(gamma(5)*gamma(1)/gamma(6));
   
6. y3=x**(1-1)*(1-x)**(3-1)/(gamma(1)*gamma(3)/gamma(4));
   
7. y4=x**(2-1)*(1-x)**(2-1)/(gamma(2)*gamma(2)/gamma(4));
   
8. y5=x**(2-1)*(1-x)**(5-1)/(gamma(2)*gamma(5)/gamma(7));
   
9. output;
   
10. end;
    
11. run;
    
12. proc sgplot data=beta;
    
13. series x=x y=y1 / lineattrs=(color='#cf86fb' pattern=solid thickness=2);
    
14. series x=x y=y2 / lineattrs=(color='#00b8e5' pattern=solid thickness=2);
    
15. series x=x y=y3 / lineattrs=(color='#00bc57' pattern=solid thickness=2);
    
16. series x=x y=y4 / lineattrs=(color='#c7b13f' pattern=solid thickness=2);
    
17. series x=x y=y5 / lineattrs=(color='#CD0000' pattern=solid thickness=2);
    
18. yaxis ranges=(0-2.5) label='概率密度';
    
19. xaxis label='随机变量';
    
20. title "贝塔分布概率密度曲线";
    
21. run;
    
22. 
    
23. /*贝塔分布累积分布函数曲线*/
    
24. data beta1;
    
25. set beta;
    
26. total1+y1;
    
27. total2+y2;
    
28. total3+y3;
    
29. total4+y4;
    
30. total5+y5;
    
31. run;
    
32. proc sgplot data=beta1;
    
33. series x=x y=total1 / lineattrs=(color='#cf86fb' pattern=solid thickness=2);
    
34. series x=x y=total2 / lineattrs=(color='#00b8e5' pattern=solid thickness=2);
    
35. series x=x y=total3 / lineattrs=(color='#00bc57' pattern=solid thickness=2);
    
36. series x=x y=total4 / lineattrs=(color='#c7b13f' pattern=solid thickness=2);
    
37. series x=x y=total5 / lineattrs=(color='#CD0000' pattern=solid thickness=2);
    
38. title "贝塔分布累积分布函数曲线";
    
39. yaxis label='累积概率';
    
40. xaxis label='随机变量';
    
41. run;

复制代码