Hi guys. I have a question and I hope someone can help me out
Let X1,.....Xn be a random sample from gamma(a,b) with a known.
Find the best unbiased estimator of 1/b
Waiting for your response as soon as you can
Thanks in advance
OK, I'll work with that, BUT that is not how most people write the gamma density. I place beta in the denomiator.
The likelihood function is .
Thus with known we have suficient for .
And since we have UMVUE for .
I could have done this yesterday, but the mean of my gamma's is
and getting an unbiased estimator of is a lot harder in that case.
And I wasn't going to do this until I knew how you were writing your density.
It can be shown that Sum(Xi) is not only sufficient but complete for beta (as w(beta) of exponential function contains an open set), so try 1/sum(xi) to estimate 1/beta. Lehmann-Scheffe tells us that an unbiased estimator that is a function of a complete statistic is the best unbiased estimator.
It can be shown that, assuming iid xi, Sum(xi) is distributed as Gamma(n*alpha,beta).
Given the above, let Y=(1/sum(xi)). Y is distributed as an inverted gamma(n*alpha, 1/beta) with mean=(1/beta)/(n*alpha-1). Thus, E(1/sum(xi)) = E(Y) = (1/beta)/(n*alpha-1), which is obviously a biased estimator of 1/beta.
Now let T = (n*alpha - 1)/sum(xi) be the unbiased estimator of 1/beta which is a function of a complete statistic. Thus, T is the best unbiased estimator of 1/beta, but it can be shown that it does not attain the lower Cramer-Rao bound.
So we can conclude that 1/beta(^) is the best estimator. Is it right?