data gamma_sample;
   call streaminit(123);
   a=2.5;s=2;
   do i=1 to 100;
     x=s*RAND('GAMMA', a);
         output;
        end;
        keep x;
run;
proc nlmixed data=gamma_sample;
parms a=1 s=1;
p=PDF ('GAMMA', x, a , s);
loglik=log(p);
model x ~ general(loglik);
run;
这个程序的过程步看不懂啊？哪位大牛告诉我每步这是什么意思！！我本来是想用极大似然估计法估计求样本容量为50的100组数据的估计值的，数据来自伽马分布，小弟生成的程序如下data sample;
call streaminit(123);
a=2;
do i=1 to 50;
    do group=1 to 100;
           ran_gamma=RAND('GAMMA', a) ;
           output;
        end;
  end;
run;
proc.......
过程步就不会了，请问能用上面好心人给的程序计算我这生成数据的估计值吗（用极大似然估计），请大家帮帮忙！！！！非常着急