Unveiling the Generalized Burr XII Distribution: Formulas & Fitting Codes

The Burr Type XII distribution, also known simply as the Burr XII or Singh-Maddala distribution, is a highly flexible continuous probability distribution renowned for its ability to model diverse phenomena, particularly in fields like reliability engineering, survival analysis, econometrics, and actuarial science. Its adaptability stems from its parameters controlling shape and tail behavior.

However, real-world data often exhibit complexities that even the standard Burr XII cannot fully capture. This has led to the development of various generalized Burr XII distributions. These extensions typically introduce additional parameters, enhancing the distribution's flexibility to model features like heavier tails, varying hazard rates, multimodality, or specific skewness and kurtosis patterns more accurately.

Highlights

• Enhanced Flexibility: Generalized Burr XII distributions add parameters to the standard model, allowing for better fits to complex, real-world datasets.
• Diverse Forms: "Generalized Burr XII" refers to a family of distributions (e.g., GBXII, Six-Parameter, Exponentiated, Weibull-Burr) rather than a single unique form, each with specific characteristics.
• Modern Modeling Tools: Bayesian inference using Stan and frequentist approaches using R packages like `gamlss` provide powerful methods for estimating the parameters of these complex distributions.

Understanding the Mathematical Foundation: PDF and CDF

To work with these distributions, understanding their Probability Density Function (PDF) and Cumulative Distribution Function (CDF) is essential. The PDF describes the likelihood of observing a specific value, while the CDF gives the probability of observing a value less than or equal to a specific point.

The Standard Burr XII Distribution (Baseline)

Before exploring the generalizations, let's define the standard two-parameter Burr XII distribution:

Probability Density Function (PDF)

Cumulative Distribution Function (CDF)

Here, \(c\) and \(k\) are shape parameters. Sometimes, a scale parameter \(\beta\) is included, replacing \(x\) with \(x/\beta\).

Key Generalized Burr XII Variants

Several generalizations exist. Here are some prominent examples found in the literature:

1. Generalized Burr XII (GBXII) Family

This introduces an additional parameter (\(\gamma\)) to modulate the exponent, increasing flexibility, particularly for lifetime data.

• CDF: \( F(x; \alpha, \beta, \gamma) = 1 - \left(1 + (x/\theta)^\alpha\right)^{-\beta \gamma} \) (Including scale \(\theta\))
• PDF: \( f(x; \alpha, \beta, \gamma) = \frac{\alpha \beta \gamma}{\theta} \left(\frac{x}{\theta}\right)^{\alpha-1} \left(1 + \left(\frac{x}{\theta}\right)^\alpha\right)^{-\beta \gamma - 1} \)

Parameters: \(\alpha, \beta, \gamma > 0\) (shape), \(\theta > 0\) (scale).

2. Six-Parameter Generalized Burr XII Distribution

This highly flexible version introduces parameters modifying the exponent based on \(x\), allowing for complex tail behaviors.

• CDF: \( F(x; a, b, c, d, e, f) = 1 - \left(1 + (a x^b)^c\right)^{-(d + e x^f)} \)
• PDF: \( f(x; ...) = \frac{d F(x)}{dx} \) (The derivative is complex but derived in sources like arXiv:0806.1579).

Parameters: \(a, b, c, d, e, f > 0\).

3. Exponentiated Generalized Burr XII (EGBXII) Distribution

This adds an exponentiation parameter (\(a\)) to the CDF of a generalized form, allowing for more varied hazard rate shapes.

• CDF: \( F(x; a, c, k, \theta) = \left(1 - \left(1 + \left(\frac{x}{\theta}\right)^c\right)^{-k}\right)^a \)
• PDF: \( f(x; a, c, k, \theta) = \frac{ack}{\theta} \left(\frac{x}{\theta}\right)^{c-1} \left(1 + \left(\frac{x}{\theta}\right)^c\right)^{-k-1} \left(1 - \left(1 + \left(\frac{x}{\theta}\right)^c\right)^{-k}\right)^{a-1} \)

Parameters: \(a, c, k > 0\) (shape), \(\theta > 0\) (scale).

4. Weibull-Burr XII Distribution

Combines features of Weibull and Burr XII distributions.

• CDF: \( F(x; \alpha, \beta, \gamma) = 1 - \exp\left(-\alpha \left(\frac{(x/\theta)^\beta}{(1+(x/\theta)^\beta)}\right)^\gamma\right) \) (Including scale \(\theta\))
• PDF: \( f(x; ...) = \frac{d F(x)}{dx} \) (Derivable from the CDF).

Parameters: \(\alpha, \beta, \gamma > 0\) (shape), \(\theta > 0\) (scale).

Visualizing the Burr XII Family Tree

The standard Burr XII distribution serves as a foundation for numerous extensions, each designed to capture different data characteristics. This mindmap illustrates the relationship between the standard form and some of its common generalizations.

Comparing Generalizations: A Radar Chart Perspective

Choosing the right generalized Burr XII distribution depends on the specific needs of the analysis and the characteristics of the data. This chart provides a qualitative comparison of some common generalizations across several dimensions (scaled 1-5, higher is better/more pronounced).

Note: This chart represents a qualitative assessment. Parameter Count is scaled inversely (lower count = higher score). Fitting Complexity is scaled inversely (less complex = higher score).

Fitting Generalized Burr XII Distributions: Stan Code

Stan is a powerful platform for statistical modeling and high-performance statistical computation. It is particularly well-suited for Bayesian inference using Markov Chain Monte Carlo (MCMC) methods. To fit a generalized Burr XII distribution in Stan, you typically define the log-probability density function (log PDF) within the model.

Example: Stan Code for the Six-Parameter Generalized Burr XII

This example demonstrates how to structure a Stan model for the complex six-parameter generalization. You would need to adapt the `log_lik` calculation based on the exact derivative of the six-parameter CDF if implementing directly, or use a simpler generalization like EGBXII (as shown in Answer B) or GBXII for easier implementation.

The following code structure is based on Answer C's example, illustrating the key components. The exact `log_lik` calculation for the six-parameter version requires careful derivation of its PDF.

// Stan model structure for a Generalized Burr XII Distribution
// (Illustrative - log_lik needs specific PDF for the chosen generalization)

functions {
  // Define the log PDF function if needed, e.g., for EGBXII
  // real EGBXII_lpdf(real x, real a, real c, real k, real theta) { ... }
}

data {
  int<lower=1> N;             // Number of observations
  vector<lower=0>[N] y;     // Data vector (must be positive)
  // Add prior hyperparameters if needed
  // vector[num_params] prior_means;
  // vector<lower=0>[num_params] prior_sds;
}

parameters {
  // Define parameters with constraints
  // Example for EGBXII:
  real<lower=0> a;
  real<lower=0> c;
  real<lower=0> k;
  real<lower=0> theta;
  // Example for Six-Parameter:
  // real<lower=0> param_a;
  // real<lower=0> param_b;
  // ... etc. for 6 parameters
}

// transformed parameters {
//   // Optional block for derived quantities or intermediate calculations
//   // Example: Calculate log-likelihood contributions for clarity
//   vector[N] log_lik;
//   for (i in 1:N) {
//      // Calculate log PDF for the chosen generalization
//      // log_lik[i] = EGBXII_lpdf(y[i] | a, c, k, theta); // Using function
//      // log_lik[i] = log(...) // Direct log PDF calculation
//   }
// }

model {
  // Priors: Assign prior distributions to parameters
  // Weakly informative priors are common
  a ~ normal(1, 5);      // Example prior for EGBXII 'a'
  c ~ normal(1, 5);      // Example prior for EGBXII 'c'
  k ~ normal(1, 5);      // Example prior for EGBXII 'k'
  theta ~ normal(mean(y), sd(y)); // Example prior for EGBXII 'theta'

  // Likelihood: Define how the data relates to the parameters
  // This uses the log PDF of the chosen generalized Burr XII distribution
  // If using a function defined in 'functions' block:
  // for (i in 1:N) {
  //   y[i] ~ EGBXII(a, c, k, theta); // Syntactic sugar for target += EGBXII_lpdf(...)
  // }
  // Or calculate directly using target += log_pdf(...)
  // Example for EGBXII directly in model block:
   for (i in 1:N) {
     target += log(a * c * k / theta) + (c - 1) * log(y[i] / theta)
               - (k + 1) * log(1 + pow(y[i] / theta, c))
               + (a - 1) * log1m_exp(-k * log(1 + pow(y[i] / theta, c))); // log1m_exp for numerical stability
   }

  // If using 'transformed parameters' block for log_lik:
  // target += sum(log_lik);
}

generated quantities {
  // Optional block for posterior predictive checks, simulated data, etc.
  // vector[N] y_rep;
  // for (n in 1:N) {
  //   y_rep[n] = generalized_burr_rng(...); // Requires RNG function
  // }
}

Key aspects of the Stan code:

• data block: Declares the data needed (observations y, sample size N).
• parameters block: Declares the parameters to be estimated (e.g., a, c, k, theta for EGBXII) with constraints (<lower=0>).
• model block: Defines the statistical model.
  • Priors: Specifies prior beliefs about the parameters (e.g., normal(1, 5)). Choosing appropriate priors is crucial in Bayesian analysis.
  • Likelihood: Specifies the probability of the observed data given the parameters, typically by adding the log-PDF contributions to the target variable. Using functions like log1m_exp can improve numerical stability.
• functions (optional): Can define reusable functions like the log PDF.
• transformed parameters (optional): For calculating intermediate quantities like the log-likelihood vector.
• generated quantities (optional): For simulations based on the fitted model (e.g., posterior predictive checks).

Visualizing Distribution Shapes

Generalized Burr XII distributions can adopt various shapes depending on their parameters. The image below, taken from a study involving a related Chen-Burr XII model (another extension), illustrates potential PDF and hazard rate function shapes that these flexible distributions can capture. Understanding these shapes helps in selecting an appropriate model for specific data characteristics like skewness, peakness, and failure rate patterns over time.

Image illustrating potential PDF and Hazard Rate Function (HRF) shapes from a study on a related Chen-Burr XII model. Source: MDPI Axioms.

These plots showcase how different parameter combinations can lead to distributions that are unimodal, decreasing, or exhibit complex hazard rate patterns (increasing, decreasing, bathtub, or inverted bathtub shapes), highlighting the value of generalized forms in accurately modeling diverse real-world phenomena.

Fitting Generalized Burr XII Distributions: R Code

R provides several ways to fit generalized Burr XII distributions.

1. Using RStan

The most flexible approach for Bayesian inference is to use the rstan package to interface with the Stan model defined above.

# Install and load rstan
# install.packages("rstan")
library(rstan)

# Assume 'stan_model_code' is a string containing the Stan code from the previous section
# Or load from a .stan file: stan_model_file <- "path/to/your/model.stan"

# Sample data (replace with your actual data)
set.seed(123)
# Example: Generate data from a known distribution for testing
# For instance, using EGBXII parameters: a=2, c=1.5, k=1, theta=10
# This requires a custom RNG or approximation; using gamma as placeholder
y_data <- rgamma(100, shape=2, rate=0.1) # Replace with actual data

# Prepare data list for Stan
stan_data <- list(
  N = length(y_data),
  y = y_data
  # Add prior hyperparameters if defined in Stan model
  # prior_means = c(...),
  # prior_sds = c(...)
)

# Compile and fit the Stan model
# Use model_code or file argument
fit <- stan(
  model_code = stan_model_code, # Or use file = stan_model_file
  data = stan_data,
  iter = 2000,          # Total iterations per chain
  warmup = 1000,        # Burn-in iterations
  chains = 4,           # Number of chains
  cores = parallel::detectCores() # Use multiple cores if available
  # init = list(list(a=1, c=1, k=1, theta=mean(y_data)), ...) # Optional initial values
)

# Check results
print(fit)
summary(fit)
plot(fit)
# Check convergence diagnostics (Rhat, effective sample size)

This workflow compiles the Stan code, runs the MCMC sampler, and provides posterior estimates for the parameters.

2. Using GAMLSS

For frequentist estimation, the Generalized Additive Models for Location, Scale and Shape (GAMLSS) framework in R is very powerful. While standard Burr XII distributions might be directly available (e.g., `BURR12` or `DBURR12` for discrete data in `gamlss.dist`), fitting *generalized* versions often requires defining a custom GAMLSS family.

This involves specifying the PDF/CDF, links functions for parameters, derivatives for fitting, initial value expressions, and the deviance increment function.

# Install and load gamlss
# install.packages("gamlss")
library(gamlss)
library(gamlss.dist) # May contain standard Burr or related families

# ---- Define Custom GAMLSS Family for a Generalized Burr XII ----
# (Conceptual Example - Requires implementing PDF, CDF, derivatives)

# Example: Defining functions for EGBXII (a, c, k, theta)
# Parameters mapped to GAMLSS: mu=theta, sigma=k, nu=c, tau=a

# PDF function
dEGBXII <- function(x, mu=1, sigma=1, nu=1, tau=1, log=FALSE) {
  # mu = theta, sigma = k, nu = c, tau = a
  if (any(x < 0)) stop("x must be positive")
  if (any(mu <= 0) || any(sigma <= 0) || any(nu <= 0) || any(tau <= 0)) stop("parameters must be positive")

  term1 <- (x/mu)^nu
  term2 <- 1 + term1
  log_pdf <- log(tau) + log(nu) + log(sigma) - log(mu) + (nu-1)*(log(x)-log(mu)) - (sigma+1)*log(term2) + (tau-1)*log(1 - term2^(-sigma))

  if (log == FALSE) return(exp(log_pdf)) else return(log_pdf)
}

# CDF function
pEGBXII <- function(q, mu=1, sigma=1, nu=1, tau=1, lower.tail = TRUE, log.p = FALSE) {
  # mu = theta, sigma = k, nu = c, tau = a
   if (any(q < 0)) stop("q must be positive")
   if (any(mu <= 0) || any(sigma <= 0) || any(nu <= 0) || any(tau <= 0)) stop("parameters must be positive")

   term1 <- (q/mu)^nu
   term2 <- 1 + term1
   cdf <- (1 - term2^(-sigma))^tau

   if (lower.tail == FALSE) cdf <- 1 - cdf
   if (log.p == TRUE) cdf <- log(cdf)
   return(cdf)
}

# Define the GAMLSS family structure
EGBXII <- function (mu.link="log", sigma.link="log", nu.link="log", tau.link="log") {
    mstats <- checklink("mu.link", "EGBXII", substitute(mu.link), c("log", "identity"))
    dstats <- checklink("sigma.link", "EGBXII", substitute(sigma.link), c("log", "identity"))
    vstats <- checklink("nu.link", "EGBXII", substitute(nu.link), c("log", "identity"))
    tstats <- checklink("tau.link", "EGBXII", substitute(tau.link), c("log", "identity"))

    structure(
        list(
            family = c("EGBXII", "Exponentiated Generalized Burr XII"),
            parameters = list(mu=TRUE, sigma=TRUE, nu=TRUE, tau=TRUE),
            nopar = 4,
            type = "Continuous",

            mu.link = as.character(substitute(mu.link)),
            sigma.link = as.character(substitute(sigma.link)),
            nu.link = as.character(substitute(nu.link)),
            tau.link = as.character(substitute(tau.link)),

            mu.linkfun = mstats$linkfun, sigma.linkfun = dstats$linkfun,
            nu.linkfun = vstats$linkfun, tau.linkfun = tstats$linkfun,

            mu.linkinv = mstats$linkinv, sigma.linkinv = dstats$linkinv,
            nu.linkinv = vstats$linkinv, tau.linkinv = tstats$linkinv,

            mu.dr = mstats$mu.eta, sigma.dr = dstats$mu.eta,
            nu.dr = vstats$mu.eta, tau.dr = tstats$mu.eta,

            # Need to define derivatives dldm, dldd, dldv, dldt
            # (Derivatives of log-likelihood w.r.t. parameters)
            # dldm = function(y, mu, sigma, nu, tau) { ... },
            # ... (similarly for dldsigma, dldnu, dldtau)
            # Need to define second derivatives for estimation algorithm
            # d2ldm2 = function(y, mu, sigma, nu, tau) { ... },
            # ... (all second derivatives: d2ldm2, d2lds2, d2ldn2, d2ldt2, and cross-derivatives)

            # Deviance increment function: -2 * log-likelihood
            G.dev.incr = function(y, mu, sigma, nu, tau, ...) -2 * dEGBXII(y, mu, sigma, nu, tau, log = TRUE),

            # Quantile residual function
            rqres = expression(rqres(pfun = "pEGBXII", type = "Continuous", y = y, mu = mu, sigma = sigma, nu = nu, tau = tau)),

            # Initial values for parameters on link scale
            mu.initial = expression(mu <- rep(log(mean(y)), length(y))),
            sigma.initial = expression(sigma <- rep(log(1), length(y))), # Start sigma (k) at 1
            nu.initial = expression(nu <- rep(log(1), length(y))),      # Start nu (c) at 1
            tau.initial = expression(tau <- rep(log(1), length(y))),       # Start tau (a) at 1

            mu.valid = function(mu) all(mu > 0),
            sigma.valid = function(sigma) all(sigma > 0),
            nu.valid = function(nu) all(nu > 0),
            tau.valid = function(tau) all(tau > 0),

            y.valid = function(y) all(y > 0)
        ),
        class = c("gamlss.family", "family")
    )
}

# Fit the model (requires full derivative definitions in the family function)
# Assuming y_data is your data
# fit_gamlss <- gamlss(y_data ~ 1, # Simple model for location parameter
#                      sigma.formula = ~ 1, # Simple model for scale/shape
#                      nu.formula = ~ 1,    # Simple model for shape
#                      tau.formula = ~ 1,   # Simple model for shape
#                      family = EGBXII(),   # Use the custom family
#                      data = your_data_frame) # Optional data frame
# summary(fit_gamlss)

Important Note: Creating a full custom GAMLSS family requires correctly implementing the first and second derivatives of the log-likelihood function with respect to each parameter. This can be mathematically intensive. The code above provides the structure but omits the complex derivative definitions.

3. Using `fitdistrplus` (for simpler cases/standard Burr XII)

The `fitdistrplus` package is excellent for fitting standard distributions. While not ideal for complex generalized forms, you can use it for the standard Burr XII or potentially simpler generalizations if you provide the density (`d...`) and optionally the cumulative (`p...`) functions.

# Install and load fitdistrplus
# install.packages("fitdistrplus")
library(fitdistrplus)

# Define density function for the standard Burr XII
# (parameters c and k, often called shape1 and shape2)
dburr12 <- function(x, shape1, shape2, log = FALSE) {
  # shape1 = c, shape2 = k
  if (any(shape1 <= 0) || any(shape2 <= 0)) stop("Shape parameters must be positive")
  log_pdf <- log(shape1) + log(shape2) + (shape1 - 1) * log(x) - (shape2 + 1) * log(1 + x^shape1)
  if (log) return(log_pdf) else return(exp(log_pdf))
}

# Define CDF function for the standard Burr XII
pburr12 <- function(q, shape1, shape2, lower.tail = TRUE, log.p = FALSE) {
  # shape1 = c, shape2 = k
  if (any(shape1 <= 0) || any(shape2 <= 0)) stop("Shape parameters must be positive")
  cdf <- 1 - (1 + q^shape1)^(-shape2)
  if (!lower.tail) cdf <- 1 - cdf
  if (log.p) return(log(cdf)) else return(cdf)
}

# Example data
y_data <- c(1.2, 2.5, 3.1, 4.0, 5.5, 7.1, 10.3) # Replace with your data

# Fit the standard Burr XII distribution
# Provide reasonable starting values
start_params <- list(shape1 = 1, shape2 = 1)
fit_fdp <- fitdist(y_data, "burr12", start = start_params)

# View results
summary(fit_fdp)
plot(fit_fdp)

This method is straightforward for standard distributions but less adaptable to the complex parameter structures of many generalized Burr XII variants compared to Stan or GAMLSS.

Getting Started with Stan in R

If you're new to Stan and R for Bayesian modeling, understanding the workflow and syntax is key. The following video provides an introduction to using Stan via R, covering basic syntax and model evaluation, which is essential for fitting distributions like the generalized Burr XII.

This video covers fundamental concepts like defining the Stan model, preparing data in R, running the sampler using `rstan`, and interpreting the output, including parameter estimates and convergence diagnostics.

Frequently Asked Questions (FAQ)

References

• A New Generalized Burr XII Distribution with Real Life Applications - ResearchGate
• A Six Parameter Generalized Burr XII Distribution - arXiv
• The Generalized Weibull-burr XII Distribution and Its Applications - JDS Online
• EXPONENTIATED GENERALIZED BURR XII DISTRIBUTION PROPERTIES AND APPLICATIONS - ResearchGate
• Stan Documentation - mc-stan.org
• RStan: the R interface to Stan - CRAN R Project
• DBURR12 function | R Documentation - R Project Search
• fitdistrplus: Help to Fit Distributions - CRAN R Project

Recommended Exploration

How do you select appropriate prior distributions for parameters when fitting models in Stan?

When should one use a Burr XII distribution versus a Generalized Extreme Value (GEV) distribution for modeling extremes?

What are specific examples of how generalized Burr XII distributions are applied in finance or insurance risk modeling?

How can you diagnose MCMC convergence issues when fitting complex models with Stan?

statisticshowto.com

Burr Distribution (Types III, XII & Others) - Statistics How To

pjsor.com

[PDF] A Generalization of Burr Type XII Distribution with Properties ...

mdpi.com

The Logistic Burr XII Distribution: Properties and Applications to Income Data

researchgate.net

(PDF) The extended Burr XII distribution: properties and applications

bookdown.org

14 Stan | Doing Bayesian Data Analysis in brms and the tidyverse

stackoverflow.com

r - Need to fix Stan code for the generalized pareto distribution to take in real arguments rather than vectors - Stack Overflow

scielo.org.co

Burr XII Distribution with Truncated and Censored Data ... - SciELO

arxiv.org

[0806.1579] On a six-parameter generalized Burr XII distribution

weirdfishes.blog

Fitting Bayesian Models using Stan and R

mdpi.com

The Exponentiated Burr XII Power Series Distribution: Properties and Applications

mc-stan.org

Fit a model with Stan — stan • rstan

researchgate.net

(PDF) A New Generalized Burr XII Distribution with Real ... - ResearchGate

pvamu.edu

PDF

ourcodingclub.github.io

Intro to Stan

docs.scipy.org

scipy.stats.burr12 — SciPy v1.15.1 Manual

researchgate.net

The pdf and cdf of the Burr Type XII distribution with skew (γ3)=3 and... | Download Scientific Diagram

pdfs.semanticscholar.org

PDF

cran.r-project.org

R-project

degruyterbrill.com

A New Method for Generalizing Burr and Related Distributions

en.wikipedia.org

Burr distribution - Wikipedia

search.r-project.org

R: Burr Type XII distribution

au.mathworks.com

Burr Type XII Distribution

researchgate.net

(Pdf) Exponentiated Generalized Burr Xii Distribution: Properties and ...

degruyter.com

A New Method for Generalizing Burr and Related Distributions

aimspress.com

The unit generalized log Burr XII distribution: properties and ...

cran.r-project.org

distributionsrd: Distribution Fitting and Evaluation

discourse.mc-stan.org

Model fitting in R stan

degruyter.com

A NEW METHOD FOR GENERALIZING BURR AND RELATED DISTRIBUTIONS

researchgate.net

Probability density function of log-Burr distribution for k =0 . 5,1,10

rdrr.io

DBURR12: The Discrete Burr type XII distribution for fitting a GAMLSS... in gamlss.dist: Distributions for Generalized Additive Models for Location Scale and Shape

epublications.marquette.edu

[PDF] The Unit Generalized Log Burr XII Distribution

jstatsoft.org

[PDF] Newdistns: An R Package for New Families of Distributions

opensiuc.lib.siu.edu

[PDF] On Simulating Univariate and Multivariate Burr Type III ... - OpenSIUC

reddit.com

r/statistics on Reddit: [Q] How to speed up model fitting with ...

cran.r-project.org

RStan: the R interface to Stan

epublications.marquette.edu

[PDF] The Burr XII-Burr XII Distribution: Mathematical Properties and ...

mathworks.com

Burr Type XII Distribution

r-bloggers.com

An introduction to Stan with R | R-bloggers

search.r-project.org

R: The Discrete Burr type XII distribution for fitting a GAMLSS...

stats.stackexchange.com

How do I fit a set of data to Burr Distribution in R? - Cross Validated

jds-online.org

The Generalized Weibull-burr Xii Distribution and Its Applications

geocities.ws

[PDF] A look at Burrs and related distributions.pdf -.:: GEOCITIES.ws ::.

researchgate.net

A Method for Simulating Burr Type III and Type XII Distributions ...