Discrete Joint Random Variables in Fréchet-Weibull Distribution: A Comprehensive Mathematical Framework with Simulations, Goodness-of-Fit Analysis, and Informed Decision-Making

Abstract

This paper introduces a novel four-parameter discrete bivariate distribution, termed the bivariate discretized Fréchet–Weibull distribution (BDFWD), with marginals derived from the discretized Fréchet–Weibull distribution. Several statistical and reliability properties are thoroughly examined, including the joint cumulative distribution function, joint probability mass function, joint survival function, bivariate hazard rate function, and bivariate reversed hazard rate function, all presented in straightforward forms. Additionally, properties such as moments and their related concepts, the stress–strength model, total positivity of order 2, positive quadrant dependence, and the median are examined. The BDFWD is capable of modeling asymmetric dispersion data across various forms of hazard rate shapes and kurtosis. Following the introduction of the mathematical and statistical frameworks of the BDFWD, the maximum likelihood estimation approach is employed to estimate the model parameters. A simulation study is also conducted to investigate the behavior of the generated estimators. To demonstrate the capability and flexibility of the BDFWD, three distinct datasets are analyzed from various fields, including football score records, recurrence times to infection for kidney dialysis patients, and student marks from two internal examination statistical papers. The study confirms that the BDFWD outperforms competitive distributions in terms of efficiency across various discrete data applications.