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.......
过程步就不会了,请问能用上面好心人给的程序计算我这生成数据的估计值吗(用极大似然估计),请大家帮帮忙!!!!非常着急