Derive the probability generating function (pfg) of a negative multinomial distribution with parameters $(k; p_{0}, p_{1}, ..., p_{r})$ where the k-th occurrence of the event with the probability $p_{0}$ stops the trials.
My approach: Find the pgf of the event that is stopped by the first occurrence of the event associated with $p_{0}$ then raise that expression to the k-th power.
This elementary event is a collection of multinomial event (excluding the stopping event) sequence of length $0 \rightarrow \infty$, followed by the stopping event:
$g_{elementary}(s_{1},...,s_{r}) = p_{0}(\sum_{j=0}^\infty (g_{multinomial}(s_{1},...,s_{r}))^j =\\ p_{0} / (1-g_{multinomial}(s_{1},...,s_{r}))$
where $g_{multinomial}(s_{1},...,s_{r})$ is the pgf of the multinomial event sequence of length $1$ with r possible outcomes, i.e. $\sum_{i=1}^r p_{i}s_{i}$
Raising $g_{elementary}(s_{1},...,s_{r})$ to the k-th power results in:
$g(s_{1},...,s_{r}) = p_{0}^k (1-\sum_{i=1}^r p_{i}s_{i})^{-k}$
While this result is correct, I am not sure about my reasoning.