正在试图入门，求大神指导，怎么在R中求GAMMA函数的极大似然估计。
我写了这个之后提示错误：
ilprec<-scan("E:\\作业\\非参\\illinois601234.txt")
Read 227 items
> a=sum(ilprec)
> a
[1] 50.937
> b=prod(ilprec)
> b
[1] 5.974125e-293
> f<-function(params,ilprec){
+ lambda=params[1]
+ alpha=params[2]
+ n=227
+ L<-n*alpha*log(lambda)-n*log(Gamma(alpha))+(alpha-1)*log(b)-lambda*a
+ return(-L)}
> result=optim(c(1,2),f,ilprec)
提示说：Error in Gamma(alpha) :
  link "alpha" not available for gamma family; available links are ‘inverse’, ‘log’, ‘identity’
求各位指导~
小的感激不尽。。