我有一组数据，通过proc histogram 得出他符合gamma分布：
                              Threshold   Theta           0
                                 Scale       Sigma    4.529809
                                 Shape       Alpha    62.27631
                                 Mean                 282.0998
                                 Std Dev              35.74714
我想用蒙特卡罗的方法模拟输出：
1. 通过给定的gamma distribution 参数（如上），确定一个gamma分布，求出gamma分布的cdf，并通过这个函数求反函数。
2.产生随机数0－1之间，通过gamma函数的反函数， 得到模拟的值。

想请教有没有高人做过类似问题，sas有没有写好的gamma distribution 反函数的function。
我在gamma distribution求反函数的地方就卡壳了，目前找资料中。希望各位指点。