Proving Sufficient Estimators for Gamma Distribution | Homework Solution

  • Thread starter Thread starter cse63146
  • Start date Start date
  • Tags Tags
    Estimators
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
8 replies · 7K views
cse63146
Messages
435
Reaction score
0

Homework Statement



Show that the product of the sample observations is a sufficient statistic for theta if the random sample is taken from a gamma distribution with parameters alpha = theta and beta = 6.

Homework Equations





The Attempt at a Solution



So I need to make sure that [tex]Y = \prod X_i[/tex] is a sufficient estimator for theta, which is true if:

[tex]\frac{f(x_1; \theta)...f(x_n; \theta)}{Y = \prod X_i}[/tex] does not depend on theta.

But I keep getting getting 1 for the ratio. Am I correct?
 
Last edited:
Physics news on Phys.org
You take samples from populations, not other samples.
You neglected to state the type of distribution (normal, gamma, uniform, ?) being sampled; until you do that not much can be done.

When you do use the actual distribution, think about writing out the joint distribution - as you started in the numerator of your first post, and work on the factorization criterion for identifying sufficiency.
 
statdad said:
You take samples from populations, not other samples.
You neglected to state the type of distribution (normal, gamma, uniform, ?) being sampled; until you do that not much can be done.

When you do use the actual distribution, think about writing out the joint distribution - as you started in the numerator of your first post, and work on the factorization criterion for identifying sufficiency.

Sorry about the type, I meant to say a sample is taken from a gamma distribution with parameters alpha = theta and beta = 6;not from a sample.

The theoroem I'm trying to use is:

[tex]\frac{f(x_1; \theta)...(f(x_n \theta)}{f_{Y_1}[u_1(x_1,...,x_n); \theta]} = h(x_1,...x_n)[/tex]

and if h(x1,...,xn) doesn't depend on theta, then it's an sufficient estimator.

I get h(x1,...,xn)=1 and I'm wondering if this is correct.
 
Have you seen this result?

For a random sample from a distribution that has density [tex]f(x;\theta), \, \theta \in \Omega[/tex], a statistic
[tex] Y=u_1(X_1, X_2, \dots, X_n)[/tex]

is sufficient for [tex]\theta[/tex] if, and only if there are two nonnegative functions [tex]k_1, k_2[/tex] such that

[tex] \prod_{i=1}^n f(x_i ; \theta) = k_1[u_1(x_1, x_2, \dots, x_n);\theta] k_2(x_1, x_2, \dots, x_n)[/tex]

where [tex]k_2(x_1, x_2, \dots, x_n)[/tex] does not depend on [tex]\theta[/tex].

This is the factorization theorem I referred to in my earlier post.

To use this, write out the joint distribution and factor it so that the second factor does not depend on either unknown parameter. This is essentially a consequence of the approach you cited, but may be a bit easier to deal with.
 
That definition is in my textbook, but we haven't gotten to it yet and I'm unsure on how to use it; we've only done the definition I mentioned in my previous post, and I thought I could use it to solve the problem.

Guess I'll have to wait until tuesday to learn it. Thanks for your help.
 
May not need to wait. I'm away from my computer now, so typing is a pain: what do you get for the product of the densities? If you can, post it in its gory detail.
 
I was looking over my textbook and found an example that used it, and after a while I figured it out:

[tex]\frac{1}{\Gamma(\theta)^n 6^{\theta n}}\prod x_i^{\theta - 1}e^{-\Sigma x_i /6}<br /> = \frac{1}{\Gamma(\theta)^n 6^{\theta n}}\prod x_i^{\theta}e^{-\Sigma x_i /6} \frac{1}{\prod x_i}[/tex]

And since [tex]\frac{1}{\prod x_i}[/tex] doesn't depend on theta, Y is a sufficient estimator of theta.

Thanks for all your help once again.
 
Thanks for all your help once again