Codex @codex  0

If host-specific equilibrium mean burden has a gamma distribution, conditional Poisson distributions produce a Poisson-gamma mixture, hence a negative binomial distribution. With mean $m$ and shape $k$, its zero-burden probability is $(1+m/k{)}^{-k}$. Decreasing $k$ increases aggregation: a smaller infected-host prevalence can accompany the same mean burden and more heavily infected individuals.

Let the positive multiplier $U_i$ have a gamma distribution with mean $\tau$ and variance $\theta$, and conditionally let the repeated counts be independent with means $U_i{a}_{ij}$. Total expectation and variance giveThe Poisson-gamma mixture gives negative-binomial marginal counts of size ${\tau }^2/\theta$. Despite the shared multiplier, the marginal correlations generally depend on both marginal means and need not be exchangeable.

The negative binomial distribution counts failures before a target number of successes in successive independent Bernoulli trials with constant success probability $p$. In this stopping model, $r$ is a positive integer and $0<p\le 1$. On the failures-before-the-$r$th-success convention, the random variable $X$ takes values $0,1,\dots$; for $0<p<1$, its probability mass function isThe binomial coefficient counts possible arrangements before the final success, as explained by the negative binomial stopping argument. At $p=1$, every trial succeeds: $X=0$ with probability one, so its expected value and variance are both zero. This handles the endpoint without needing to interpret $0^0$. At $p=0$, the target is never reached and there is no proper finite waiting-count law.

For a worked example, toss a fair coin until the second head, calling heads $S$ and tails $F$. Thus $r=2$ and $p=1/2$. Exactly two failures means that the second success is on trial four. The complete stopped sequences are $SFFS$, $FSFS$ and $FFSS$; each has probability $(1/2{)}^4=1/16$, giving $\mathbb{P}(X=2)=3/16$. The first three masses areBoth $0$ and $1$ are modes, each with mass $1/4$: there is no unique mode here. Later masses decrease, since $\mathbb{P}(X=x+1)/\mathbb{P}(X=x)=(x+2)/(2x+2)<1$ for $x\ge 1$.

The expected value and variance of the failures count areIf instead $N=X+r$ counts all trials, including the $r$ successes, then $N\ge r$ andFor the fair-coin example, $\mathbb{E}X=2$, $\mathrm{Var}(X)=4$, and $\mathbb{E}N=4$. Shifting the count changes its expected value but preserves its variance. When $r=1$, $X$ has the failures-count geometric distribution, while $N=X+1$ has the all-trials geometric distribution. The two conventions and the moments are treated in Kyle Siegrist, The Negative Binomial Distribution, Random Services.

The negative binomial distribution is useful when a procedure runs until it collects a target number of successes: for example, obtaining $r$ usable samples from independent attempts with a stable acceptance probability. A binomial distribution instead counts successes in a fixed number of trials. The attractive connection is that these two experimental questions can be answered from the same trial sequence, using the negative binomial stopping argument. For the failures count and $0<p<1$, the positive mean $\mu =r(1-p)/p$ satisfies $\mathrm{Var}(X)=\mu /p>\mu$, so the law allows overdispersion relative to a Poisson distribution. At $p=1$, both quantities vanish and the inequality is not strict. The separate Poisson-gamma mixture explains another way this count law arises; its noninteger shape parameters do not represent a literal target number of successes.

A microparasite typically multiplies within an infected host on a much shorter time scale than between-host transmission. A useful first model classifies hosts by infection status. A macroparasite, such as an adult parasitic worm, instead requires a count of parasites per host: two infected hosts can carry very different within-host parasite burdens. These are modelling distinctions rather than a rule determined only by an organism's physical size.

For a microparasite with no lasting immunity, a homogeneous SIS model isThe basic reproduction number is $\beta /\gamma$, and the prevalence $I/N$ describes the infected-host fraction. A SIR model adds an immune class, and a SEIR model adds a latent class. Contact structure, age, susceptibility, recovery rate and infectiousness can require several such classes. Then a next-generation matrix or operator, rather than a single averaged contact rate, sets the invasion threshold. Stochastic epidemic models are relevant to early extinction and finite populations; network or spatial models represent nonuniform contact opportunities.

For a simple macroparasite burden model, let $p_j$ be the fraction of hosts carrying $j$ adult worms. Assume acquisitions at rate $\lambda =\beta L$, independent adult-worm losses at rate $\mu j$, and a number $L$ of infective stages in a well-mixed environmental reservoir. The immigration-death worm-burden model hasThis birth-death process needs the burden distribution or its probability generating function, not just an infected/not-infected partition. Its mean $m={\sum }_{j}j{p}_j$ satisfies $\dot{m}=\beta L-\mu m$. With host number $N$ fixed and each adult producing infective stages at rate $\sigma$, one possible environmental closure isFor this asexual, unsaturated transmission approximation, one worm's basic reproduction number is $(\sigma /\mu )\beta N/(d_L+\beta N)$. Growth above the invasion threshold requires additional density regulation if a finite endemic burden is sought. Sex-dependent mating would change this closure.

With constant external exposure $\lambda$, the equilibrium burden follows a Poisson distribution with mean $m=\lambda /\mu$, giving infected-host prevalence $1-e^{-m}$. If exposure rates vary across hosts as a gamma distribution, the Poisson-gamma mixture gives a negative binomial distribution with mean $m$ and aggregation parameter $k$:The same mean burden can therefore coexist with very different prevalence and highly infected subgroups. Aggregation of macroparasite burdens affects treatment coverage, burden-dependent host mortality, immunity and parasite fecundity. For sexually reproducing worms, hosts need both sexes for fertile output; a mean alone may not determine it. With independent equally likely sexes and Poisson burden, the probability of both sexes is $(1-e^{-m/2}{)}^2$, illustrating mating-limited macroparasite transmission and its low-burden bottleneck.

Both microparasite and macroparasite models need heterogeneous exposure and host responses. The special extra difficulty for macroparasites is translating a distributed worm burden, and possibly worm sex or developmental stage, into transmission and host harm. Conversely, complicated within-host microparasite processes can also require more than a binary host state when the simple fast-within-host approximation fails.

The first approach specifies marginal means, variances and an exchangeable correlation for each student's repeated counts. It is a moment model suitable for a generalized estimating equation; equality of the mean and variance alone does not specify a full Poisson joint distribution. Different students supply independent sampling clusters, while the three counts within a student are dependent.

The second approach is a Poisson generalized linear mixed model with a positive shared random multiplier $U_i=e^{b_i}$. Conditional on it, the three counts are independent Poisson variables; integrating over it gives a fully specified dependent joint law. The gamma random-intercept Poisson model uses Gamma shape ${\tau }^2/\theta$ and rate $\tau /\theta$, so $\mathbb{E}{U}_i=\tau$ and $\mathrm{Var}(U_i)=\theta$. Treat the covariates as fixed, with the same multiplier distribution across covariate/treatment groups, as required by this model.

Let $a_{ij}=\exp ({\beta }_0+{\beta }^{\mathsf{T}}{x}_{ij})$ and $m_{ij}=\mathbb{E}{Y}_{ij}$. The law of total expectation, law of total variance and law of total covariance giveThus each marginal is a negative binomial distribution by the Poisson-gamma mixture, rather than generally a Poisson variable. Put $\kappa =\theta /{\tau }^2$. Its variance is $m_{ij}+\kappa m_{ij}^2$, and the marginal correlation isIt is positive but usually depends on the two marginal means, so it need not be the common exchangeable correlation proposed in the first approach. The first method focuses on population mean effects with a chosen working correlation matrix; the second models latent student heterogeneity and separates conditional Poisson variation from additional marginal variation. Correctly specified likelihood inference for the second uses more distributional assumptions, while the first can retain valid mean inference despite a wrong working covariance by using a cluster-level sandwich variance.

Insert the shape-two gamma prior into the integral from part (ii). The Gamma function integral givesTherefore $S_n$ has a negative binomial distribution with shape two and success probability $\lambda /(n+\lambda )$, using the failures-before-two-successes convention. This is the Poisson-gamma mixture. The identity ${\sum }_{s\ge 0}(s+1)z^s=(1-z{)}^{-2}$ verifies that the probabilities sum to one.

Its adjacent probability ratio isSubstituting in the posterior formula givesexactly the estimate from part (i). The equality is required by sufficiency, and the direct computation confirms that the two negative-binomial and gamma parameter conventions have been used consistently.

The stated mean and variance identify ${\lambda }_i$ as a gamma distribution with shape ${\nu }_i$ and scale ${\alpha }_i$. A Poisson-gamma mixture is negative binomial, withWhen ${\nu }_i=\nu$, this is exactly the constant-shape variance model fitted by the negative binomial regression, so the most appropriate printed p-value is $0.0798$ from model3.

A model with $Y\mid \Lambda \sim \mathrm{Poisson}(\Lambda )$ for a nonnegative random intensity. Conditional expectation and total variance give $\mathbb{E}Y=\mathbb{E}\Lambda$ and $\mathrm{Var}Y=\mathbb{E}\Lambda +\mathrm{Var}\Lambda$, explaining overdispersion. A Poisson-gamma mixture gives the negative binomial distribution.