Skip to content

ParameterMixtureDistribution ​

A distribution whose own parameter θ\theta is itself random, drawn from priorprior — the marginal of dist(θ)dist(\theta) averaged over θ∼prior\theta \sim prior.

ParameterMixtureDistribution(dist(theta), Distributed(theta, prior))dist(theta) with theta itself drawn from prior.
Typed signatures
TypePackageApplies
(any, expression<Distributed>) -> distributionstatisticsalways

Domain: Statistics

Details
  • Two named conjugate closed forms: PoissonDistribution(theta) mixed over a GammaDistribution(shape, scale) prior on theta is exactly NegativeBinomialDistribution(shape, 1/(1+scale)) (the standard Gamma-Poisson identity — PDF/Mean/Variance all delegate to that equivalent NegativeBinomial). BinomialDistribution(n, theta) mixed over a BetaDistribution(a, b) prior is the Beta-Binomial: PDF(x)=(nx)B(x+a,n−x+b)/B(a,b)PDF(x) = \binom{n}{x} B(x+a,n-x+b)/B(a,b), Mean=na/(a+b)Mean = na/(a+b), Variance=nab(a+b+n)/((a+b)2(a+b+1))Variance = nab(a+b+n)/\left((a+b)^2(a+b+1)\right).
  • Anything else falls back to the LAW OF TOTAL EXPECTATION/VARIANCE — E[X]=Eθ[Mean(dist(θ))]E[X] = E_\theta[Mean(dist(\theta))], Var(X)=Eθ[Variance(dist(θ))]+Varθ[Mean(dist(θ))]Var(X) = E_\theta[Variance(dist(\theta))] + Var_\theta[Mean(dist(\theta))] — built entirely out of the (already fully extended) Mean/Variance/Expectation operators, so it composes with any inner distribution/prior pair those already answer. Resolves exactly when Mean(dist(theta)) (and, for Variance, Variance(dist(theta)) too) is linear or quadratic in theta — Expectation's own recognized shapes; anything past that stays unevaluated.
  • PDF has no generic fallback (only the two named closed forms above) — an unrecognized inner/prior pair's PDF stays unevaluated.

Examples

Basic #

matches PDF(NegativeBinomialDistribution(3, 1/3), 4) — p = 1/(1+scale) = 1/3.

live ↗
InPDF⁡(ParameterMixtureDistribution(PoissonDistribution⁡(θ),Distributed(θ,GammaDistribution(3,2))),4)\operatorname{PDF}(\mathrm{ParameterMixtureDistribution}(\operatorname{PoissonDistribution}(θ), \mathrm{Distributed}(θ, \mathrm{GammaDistribution}(3, 2))), 4)Out80729\frac{80}{729}

shape * scale = 3 * 2.

live ↗
Inmean(ParameterMixtureDistribution(PoissonDistribution⁡(θ),Distributed(θ,GammaDistribution(3,2))))\mathrm{mean}(\mathrm{ParameterMixtureDistribution}(\operatorname{PoissonDistribution}(θ), \mathrm{Distributed}(θ, \mathrm{GammaDistribution}(3, 2))))Out66

shapescale + shapescale^2 = 6 + 12.

live ↗
InVariance(ParameterMixtureDistribution(PoissonDistribution⁡(θ),Distributed(θ,GammaDistribution(3,2))))\mathrm{Variance}(\mathrm{ParameterMixtureDistribution}(\operatorname{PoissonDistribution}(θ), \mathrm{Distributed}(θ, \mathrm{GammaDistribution}(3, 2))))Out1818

n a/(a+b) = 10*2/5.

live ↗
Inmean(ParameterMixtureDistribution(BinomialDistribution⁡(10,θ),Distributed(θ,BetaDistribution(2,3))))\mathrm{mean}(\mathrm{ParameterMixtureDistribution}(\operatorname{BinomialDistribution}(10, θ), \mathrm{Distributed}(θ, \mathrm{BetaDistribution}(2, 3))))Out44
live ↗
InVariance(ParameterMixtureDistribution(BinomialDistribution⁡(10,θ),Distributed(θ,BetaDistribution(2,3))))\mathrm{Variance}(\mathrm{ParameterMixtureDistribution}(\operatorname{BinomialDistribution}(10, θ), \mathrm{Distributed}(θ, \mathrm{BetaDistribution}(2, 3))))Out66
live ↗
InPDF⁡(ParameterMixtureDistribution(BinomialDistribution⁡(4,θ),Distributed(θ,BetaDistribution(2,2))),2)\operatorname{PDF}(\mathrm{ParameterMixtureDistribution}(\operatorname{BinomialDistribution}(4, θ), \mathrm{Distributed}(θ, \mathrm{BetaDistribution}(2, 2))), 2)Out935\frac{9}{35}
About
Examples6

See also: NegativeBinomialDistribution, BetaDistribution, Expectation, Distributed