See below,

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;

太感谢您了，我知道您是先生成100个来自伽玛的随机数，但是下面的过程步，我就解读的不太明白，请您稍微给出每部的解释，真是麻烦您了，我查资料也没弄明白