I am facing a problem when trying to calculate the distribution of the sum of iid Beta-Negative-Binomial random variables or for that matter if only parameter $r$ is different.
To get a hint to how they might be distributed I calculated the characteristic function of the BNB distribution:

$$\varphi_X(t)=\int_0^1 \! \sum_{x=0}^{\infty} \binom{x+r-1}{x}p^r(1-p)^x\frac{p^{\alpha-1}(1-p)^{\beta-1}}{B(\alpha,\beta)}e^{itx} \, dp$$
where $B$ is the Beta-function. One can now factor out the terms that do not contain $x$ and use the generating function $\sum_{n=0}^\infty \binom{n+k}{k}x^n=\frac{1}{(1-x)^{k+1}}$ to get the following expression for the characteristic function:

$$\varphi_X(t)=\int_0^1 \! \frac{p^{\alpha-1}(1-p)^{\beta-1}}{B(\alpha,\beta)}\frac{p^r}{(1-(1-p)e^{it})^r} \, dp.$$
The characteristic function of the sum of two independent variables $X$ and $Y$ is in this case the product of their corresponding characteristic functions.

QUOTE: The characteristic function of the sum of two variables $X$ and $Y$ is defined as the product of their corresponding characteristic functions. END OF QUOTE. That's not true. If they're independent, then the charateristic function of their sum can be shown to be the product of their characteristic functions, but it is not defined to be that. Rather, that is a demonstrable fact, not a definition. The definition is still the same as for characteristic functions in general: it is $t\mapsto \mathbb E(e^{t(X+Y)})$.
–
Michael HardyDec 4 '12 at 17:34

1

Some of your characteristic functions $\varphi$ are such that $\varphi(0)\ne1$. This is odd.
–
DidDec 4 '12 at 18:52

@MichaelHardy : That is obviously correct. I will correct that in my question.
–
MarkDec 4 '12 at 19:26

@did: I will add the pmf of the BNB distribution too so it will be more obvious why i tried it this way. Trying to find the characteristic function when using the textbook definition of the pmf has proven to be too hard for me.
–
MarkDec 4 '12 at 19:29

1

The trouble is not with what you tried or did not try but with the fact that characteristic functions are defined as $\varphi_X(t)=\mathbb E(\mathrm e^{itX})$. Hence $\varphi_X(0)=1$ for every random variable $X$. If $\psi(0)\ne1$, then $\psi\ne\varphi_X$ for every $X$. Thus, there is a problem with your formulas, which I suggest to get rid of before going any further.
–
DidDec 4 '12 at 23:23