6
$\begingroup$

I have same data and I would like to choose a model for it. To start with I fit an exponential distribution and a gamma distribution. Now I wanted to do a simple likelihood ratio test . However, I am told that to do this properly the two models have to be nested (which they are) and the parameter space of one has to be in the interior of the other, not on the boundary. This second condition doesn't seem to hold unfortunately.

What can I do?

$\endgroup$

1 Answer 1

4
$\begingroup$

Fortunately, you're mistaken.

The shape parameter for a gamma ($\alpha$, say) has to be $\ge 0$.

http://en.wikipedia.org/wiki/Gamma_distribution

The exponential has $\alpha=1$.

http://en.wikipedia.org/wiki/Gamma_distribution#Others

So the exponential is not at the boundary and you should be able to apply a likelihood ratio test without difficulty.

(I would say, however, that hypothesis tests are not necessarily a good approach to model selection.)

$\endgroup$
8
  • $\begingroup$ I am very pleased I am mistaken :) Thank you. Maybe "model selection" has a technical meaning that is not exactly what I am looking for. In my case I am trying to find a good model that accurately describes the data I have. Is there a better way to do this? $\endgroup$
    – Simd
    Aug 31, 2014 at 7:35
  • $\begingroup$ What do you mean by "a good model" and why would a hypothesis test be a good way to achieve it? $\endgroup$
    – Glen_b
    Aug 31, 2014 at 8:11
  • $\begingroup$ I am in no way a statistician so I may not be using the right words. I want a model for my data so I can, for example, make predictions or look for anomalies. I just mean a good model in the normal non-formal sense of fitting the data reasonably well. I was using the hypothesis test to see if one model was at least better than another. Any suggestions for better ways are gratefully received. $\endgroup$
    – Simd
    Aug 31, 2014 at 8:35
  • 1
    $\begingroup$ Yes, that's fine. "Fitting the data reasonably well" doesn't really change much with sample size, but our ability to detect even trivial deviations from a specific model (like an exponential model, against a gamma model) does change with sample size. At larger sample sizes, you become practically certain to reject an exponential model against a more general alternative like the gamma (and then reject a gamma against a generalized gamma, and on an on in turn), but this doesn't imply that an exponential model would be perfectly fine for everything you want to use it for. ...(ctd) $\endgroup$
    – Glen_b
    Aug 31, 2014 at 9:02
  • $\begingroup$ (ctd)... Your data will almost certainly be neither exponential nor gamma (George Box's epigram is relevant); the relevant question is not whether you can detect a difference (as sample size grows, this becomes a certainty) - it's how much that matters to your inference. That's not a question answered by hypothesis tests; it's better addressed by diagnostic tools like Q-Q plots and simulations. $\endgroup$
    – Glen_b
    Aug 31, 2014 at 9:06

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.