[R] Gamma MLE

[rach.s]

Hello,

I was asked to try the following code on R,

gamma.mles
function (xx,shape0,rate0)
{
n<- length(xx)
xbar<- mean(xx)
logxbar<- mean(log(xx))
theta<-c(shape0,rate0)
repeat {
theta0<- theta
shape<- theta0[1]
rate<- theta0[2]
S<- n*matrix(c(log(rate)-digamma(shape)+logxbar,shape/rate-xbar),ncol=1)
I<- n*matrix(c(trigamma(shape),-1/rate,-1/rate,shape/rate^2),ncol=2)
theta<- theta0 + solve(I) %*% S
if(max(abs(theta-theta0)) < 1e-08)
break
}
list(estimates=theta, infmat=I)
}

However, this appears: Error: object "gamma.mles" not found

I tried looking in the packages for gamma.mles, but I couldn't find it
anywhere. Can someone tell me where can I load it?

Thanks
-- 
View this message in context: http://www.nabble.com/Gamma-MLE-tf4122769.html#a11724823
Sent from the R help mailing list archive at Nabble.com.