---
title: Discrete Generalized Pareto Distribution
url: https://www.emergentmind.com/topics/discrete-generalized-pareto-distribution-dgpd
type: topic
---

# Discrete Generalized Pareto Distribution

The discrete generalized Pareto distribution (DGPD) is a theoretically grounded, parameterized family of probability distributions defined on the nonnegative integers, designed primarily to model the tail behavior of integer-valued data. It arises via a direct discretization of the continuous generalized Pareto distribution (GPD), inheriting its flexible tail properties and playing a central role in discrete extreme value theory, multivariate extremes, and advanced count-data modeling with heavy tails, threshold exceedances, and regression structures.

## 1. Definition, Parameterizations, and Theoretical Justification

The univariate DGPD is supported on $\{0, 1, 2, \ldots\}$ or $\{u, u+1, \ldots\}$ for some threshold or location parameter $u$. The standard definition uses the continuous GPD survival/cdf, yielding the following probability mass function for $k \geq 0$ (or $x \geq \mu$):

\[
p(k; \sigma, \xi) = \left(1 + \frac{\xi k}{\sigma}\right)^{-1/\xi} - \left(1 + \frac{\xi (k+1)}{\sigma}\right)^{-1/\xi}, \quad \sigma > 0, \xi \in \mathbb{R}
\]

For $\xi=0$, one obtains the discrete analogue of the exponential:
\[
p(k; \sigma, 0) = e^{-k/\sigma}(1-e^{-1/\sigma})
\]

The three-parameter location-scale-shape form is:
\[
p(x; \mu, \sigma, \xi) = \left[1 + \frac{\xi}{\sigma}(x - \mu)\right]^{-1/\xi} - \left[1 + \frac{\xi}{\sigma}(x - \mu + 1)\right]^{-1/\xi}, \quad x \geq \mu
\]
with $p(x) = 0$ for $x < \mu$ and $1+\xi(x-\mu)/\sigma > 0$ required [1312.4383].

The DGPD is the unique nondegenerate limit for discrete threshold exceedance distributions under natural domain-of-attraction and rounding conditions. If $X$ is discrete and $X = \lfloor Y \rfloor$ for $Y$ in the maximum domain of attraction of a continuous GPD, then for high thresholds $u$, the law of $X-u\,|\,X \geq u$ converges uniformly to DGPD as $u \to x_F$ [1707.05033, 2506.19361].

The tail behavior is governed by the shape parameter $\xi$:
- $\xi > 0$: heavy-tailed (power-law) decay
- $\xi = 0$: discrete exponential
- $\xi < 0$: finite upper endpoint

In the multivariate case, the DGPD generalizes to the MDGPD, which retains threshold stability and tail dependence properties via spectral or geometric constructions [2506.19361].

## 2. Properties, Moments, and Special Cases

The DGPD cdf, survival, and hazard functions are available in closed form:
- CDF: $F(x) = 1 - [1 + \lambda(x-\mu+1)]^{-\alpha}$
- Survival function: $\bar{F}(x) = [1 + \lambda(x-\mu)]^{-\alpha}$
- Hazard: $h(x) = 1 - \left(\frac{1+\lambda(x-\mu)}{1+\lambda(x-\mu+1)}\right)^{\alpha}$ [1312.4383]

The $r$th moment exists if $r < 1/\xi$ ($\xi>0$). In particular,
\[
E[X^r] = \sum_{x=\mu+1}^\infty \left[x^r - (x-1)^r\right] [1+\lambda(x-\mu)]^{-\alpha}
\]
The mean and variance converge only for suitable $\alpha, \lambda$ (e.g., mean exists for $\alpha > 1$).

Special cases include:
- Geometric distribution: DGPD with $\xi \to 0$
- Discrete Lomax: two-parameter DGPD with $\mu=0$

DGPD is tail-equivalent to both the continuous GPD (with a continuity correction) and generalized Zipf distribution for large scale $\sigma$, inheriting their tail indices [1707.05033]. It remains invariant under linear transformations and rounding operations, i.e., threshold stability.

## 3. Parameter Estimation and Inference Methods

Parameter estimation proceeds via two main methods:

- **Seed (μ-(μ+1) frequency) method**: For data $x_1,\ldots,x_n$, $\mu$ is set to the sample minimum. Relative frequencies at the smallest two counts are equated to the model pmf to solve numerically for $(\sigma_0, \xi_0)$ [2004.06150, 1312.4383].
- **Maximum Likelihood Estimation (MLE)**: The log-likelihood for observations $x_i$ is
\[
\ell(\mu, \sigma, \xi) = \sum_{i=1}^n \ln\left\{ [1+\xi(x_i-\mu)/\sigma]^{-1/\xi} - [1+\xi(x_i-\mu+1)/\sigma]^{-1/\xi} \right\}
\]
Partial derivatives yield nonlinear score equations, solved numerically (e.g., using R’s `optim()` with Simulated Annealing) [2004.06150].

Bootstrap methods are used to estimate standard errors and form confidence intervals: resample, re-fit DGPD, and use empirical quantiles of parameter estimates [2004.06150].

Goodness-of-fit is assessed via:
- Discrete Kolmogorov–Smirnov statistics (with parametric bootstrap for discrete null)
- Pearson’s chi-squared test with binned support and expected counts
- Discrete Q–Q plots [1312.4383, 1707.05033]

## 4. Flexible Extensions: Bulk, Tail, and Zero-Inflated Models

The DGPD is well-suited for modeling high-threshold exceedances, but becomes less reliable for lower thresholds or full-support modeling. Several flexible extensions have been proposed [2409.18719, 2210.15253]:
- **Discrete Extended GPD (DEGPD)**: Introduces a CDF transformation $G$:
\[
P(Y=y) = G(F_\mathrm{GPD}(y+1;\beta,\xi)) - G(F_\mathrm{GPD}(y;\beta,\xi)),
\]
where $G$ is usually $u^\kappa$, truncated Normal, or truncated Beta; $\kappa$ tunes the lower-tail or bulk behavior, and ordinary DGPD is recovered when $G(u) = u$.

- **Zero-Inflated DGPD (ZIDEGPD)**: For excess zeros,
\[
\Pr(Y=0) = \pi + (1-\pi) G(F_\mathrm{GPD}(1;\beta,\xi))
\]
with the bulk handled as above; $\pi$ is the structurally zero probability [2409.18719, 2210.15253].

- **Smooth Transition Models**: Bypass the need to preselect a threshold $u$ and allow parameters to vary as functions of covariates within a penalized likelihood regression (GAM) framework, using spline expansions and link functions for all DGPD (or DEGPD) parameters [2210.15253].

For bulk+tail and ZI models, maximum likelihood estimation is used, often with numerical optimization and BIC/AIC-based model selection; the DEGPD with $\kappa \approx 1$ indicates that classical DGPD is appropriate. Real data demonstrate that DEGPD and extensions outperform Poisson/NB and classical DGPD in bulk+tail or low-threshold scenarios [2409.18719, 2210.15253].

## 5. Multivariate Discrete Generalized Pareto Distributions

The MDGPD extends the DGPD to vectors of discrete exceedances over high thresholds [2506.19361]. The standard MDGPD is defined through the law of $\boldsymbol N$ with the following properties:
- $G = \max_i N_i$ is geometric,
- $\boldsymbol S = \boldsymbol N - G$ is independent and supported on $\max_i S_i = 0$.

Simulation draws on the generator representation: for i.i.d. $\Delta_i$, $G \sim \text{Geom}(1-e^{-1})$, $N_1 = G + \Delta_1 \mathbf{1}_{\Delta_1 < 0}$, $N_2 = G - \Delta_1 \mathbf{1}_{\Delta_1 \geq 0}$.

Bootstrap and likelihood-free inference in higher dimensions utilizes neural Bayes risk minimization: generate synthetic data under random parameter draws, train neural networks to predict parameter vectors, and minimize squared error over large simulation batches [2506.19361].

Multivariate DGPD preserves tail dependence, threshold invariance, and can accommodate flexible marginal and dependence structures, as illustrated in the modeling of spatial dry-spell extremes.

## 6. Applications and Empirical Performance

DGPD and its extensions have been empirically validated in diverse contexts:

- **Non-life insurance claims**: DGPD systematically outperforms negative binomial models by orders of magnitude in AIC/BIC across both yearly and aggregated datasets, accurately capturing overdispersion and heavy claim frequency tails [2004.06150].
- **Road accident blackspots**: DGPD and its discrete Lomax (DLomax) special case provide parsimonious parametric fits for counts per site; bootstrap-based discrete KS tests confirm adequacy in practical datasets [1312.4383].
- **Discrete extremes**: DGPD fits rare-event tails in Poisson simulations, word-frequency, tornado clusters, and multiple-birth data more accurately than continuous GPD (even with continuity corrections), demonstrating the necessity for discrete-tail modeling [1707.05033].
- **Flexible regression and full-support modeling**: ZIDEGPD–GAM and DEGPD outperform zero-inflated negative binomial GAMs in avalanche, doctor-visit, and complaint-count datasets, capturing both zero inflation and tail behavior robustly [2210.15253, 2409.18719].
- **Multivariate extreme value analysis**: MDGPD has been applied in spatial extremes such as Swiss dry-spell analysis, utilizing neural likelihood-free inference and exhibiting strong empirical fit and tail dependence estimation [2506.19361].

## 7. Practical Considerations, Limitations, and Outlook

The DGPD and its extensions offer several advantages:
- Asymptotically justified for integer exceedances,
- Closed-form pmfs and tractable parameterizations,
- Tail flexibility, threshold invariance, and multivariate analogues,
- Robust bootstrap-based confidence estimation and diagnostics.

Limitations include:
- Need for custom numerical optimization; no widely available off-the-shelf implementations exist (although some R packages provide basic versions for recent extensions),
- In small data or nearly homogeneous samples, the three-parameter model may offer minimal advantage,
- Model extensions introduce additional identifiability challenges when bulk/zero parameters are weakly identified  [2409.18719, 2210.15253].

A plausible implication is that, while DGPD is an essential tool for discrete extreme value modeling, extended models such as DEGPD and ZIDEGPD, as well as MDGPD for multivariate inference, are necessary when modeling the full distribution, handling lower thresholds, or addressing zero inflation. Robust threshold selection and diagnostics remain active areas of methodological research.

Source: https://www.emergentmind.com/topics/discrete-generalized-pareto-distribution-dgpd