5
$\begingroup$

I'm dealing with some biological data (read counts from gene expression measurements). For each individual we have 2 measurements, which correspond to two different version of the same gene. These measurements are discrete count data, and generally we would model them with with a negative binomial distribution, since they contain plenty of noise in a addition to the poisson distributed shot noise.

What I"m interested in is the ratio between counts. I'm given to understand that if we consider two poisson distributed random variables, and condition on their total counts, then the count for one in particular will be binomial distributed (correct me if I'm wrong). My question is this - what if we replace the poisson distributions here with negative binomial distributions? What is the distribution of the counts for each one, conditional on the total? Is it a beta binomial?

$\endgroup$
8
  • 1
    $\begingroup$ Do the two negative binomial distributions have the same parameter values? $\endgroup$
    – jbowman
    Sep 26, 2016 at 13:22
  • $\begingroup$ No, I don't have any reason to think they do, although they're likely to be quite similiar. $\endgroup$ Sep 26, 2016 at 15:01
  • 2
    $\begingroup$ Can these counts be 0? $\endgroup$
    – Glen_b
    Sep 27, 2016 at 4:43
  • $\begingroup$ The total must be at least 1, and will probably get thresholded higher than that at something like 20 or 50. $\endgroup$ Sep 27, 2016 at 15:30
  • $\begingroup$ ttu-ir.tdl.org/ttu-ir/bitstream/handle/2346/59954/… $\endgroup$ Oct 2, 2016 at 19:07

1 Answer 1

3
$\begingroup$

I don't have a formal proof, but I think your intuition is right and that it will be beta-binomial, under the assumption that the two Negative Binomials have the same "p" parameter. A hand-wavy explanation would be that Negative Binomial is just a Poisson with a gamma prior. If you have two Poissons and condition on the sum you get a Binomial. If you have two gammas with the same scale parameter and divide by the total you get a Beta.

More generally, if you have many Poissons and condition on the sum you get a Multinomial. If you have many Gammas with same scale and divide by the total you get a Dirichlet. So this suggests that if you take a bunch of Negative Binomials that have the same "p" parameter but possibly different "r" parameters and condition on their sum, it is likely to follow a Dirichlet-Multinomial distribution.

$\endgroup$
1

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

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