Evaluating the Discrete Generalized Rayleigh Distribution: Statistical Inferences and Applications to Real Data Analysis

: Various discrete lifetime distributions have been observed in real data analysis. Numerous discrete models have been derived from a continuous distribution using the survival discretization method, owing to its simplicity and appealing formulation. This study focuses on the discrete analog of the newly generalized Rayleigh distribution. Both classical and Bayesian statistical inferences are performed to evaluate the efficacy of the new discrete model, particularly in terms of relative bias, mean square error, and coverage probability. Additionally, the study explores different important submodels and limiting behavior for the new discrete distribution. Various statistical functions have been examined, including moments, stress–strength, mean residual lifetime, mean past time, and order statistics. Finally, two real data examples are employed to evaluate the new discrete model. Simulations and numerical analyses play a pivotal role in facilitating statistical estimation and data modeling. The study concludes that the discrete generalized Rayleigh distribution presents a notably appealing alternative to other competing discrete distributions.


Introduction
As each day passes, the volume of data in our world increases exponentially, necessitating the development of new statistical distributions to better characterize the features of many phenomena and experiments.While a great deal of lifetime data appear to be continuous, they are originally discrete.This discrepancy ensures the need for more appropriate methods to generate discrete distributions that more accurately represent the data in the experiment.Discrete distributions are frequently employed in statistical modeling for several reasons.
Discrete distributions are used to model data that can only take on a finite or countably infinite number of values, such as counts, proportions, and binary outcomes, for example, the number of customers in a store, the number of heads in a coin flip, or the number of defective items in a production line.Discrete distributions are often easy to understand and interpret as they model data that take on a limited number of values.The probability mass function (pmf ) or probability generating function (pgf ) of a discrete distribution is a simple function that provides the probability of each possible outcome.Also, many discrete distributions have closed-form expressions for their pmf or pgf, which makes it easy to work with them mathematically.This allows for efficient computation of probabilities and moments without the need for integration.Furthermore, discrete distributions can be used to model a wide variety of real-world phenomena, such as the distribution of species in an ecosystem, the distribution of genetic variations in a population, or the distribution of traffic on a road network.
The characterization of continuous random variables can be performed either by their probability density function, cumulative distribution function, moments, momentgenerating function, hazard rate functions, or others.Different discretization methods appeared in the literature to create an appropriate discrete distribution based on the underlying continuous model.
By deriving discrete analogs or counterparts of well-known continuous distributions, statisticians can better tailor their models to the specific nature of the data.Usually, creating a discrete analog from a continuous distribution is based on the principle of preserving one or more characteristic properties of the continuous one.Consequently, different ways to discretize a continuous distribution appear in the literature depending on the property the researcher intends to preserve; for example, Lai [10] used the survival and the hazard rate preservation methods to create discrete distributions from different continuous ones.Haj Ahmad and Almetwally [11] used the survival, hazard rate, and probability distribution function preservation methods to discretize the generalized Pareto distribution.
The benefit of using the survival discretization method is that it can maintain the statistical properties of the original distribution, including median and percentiles, in addition to the overall shape of the distribution.A drawback of this method is that it can be computationally intensive and may require numerical methods for complex distributions.
For the hazard preservation method, the main benefit is that it preserves the hazard function of the continuous distribution.This is important in applications like reliability analysis where the failure rate is a key parameter.On the other hand, mathematical complexity can be viewed using this method, especially for continuous distributions with nonlinear hazard functions.This complexity can increase computational time and resource requirements.Another drawback of the hazard preservation method is that it only preserves the hazard function, but other characteristics of the distribution (like mean, variance, or skewness) may not be as accurately retained.For more details about other discretization methods and their properties, one may refer to [12,13] who provided a review of several discretization methods.
From the previous research work, it is evident that the results look appealing and motivational to continue creating new discrete distributions to model new data.
In the present study, we obtain a discrete analog of the continuous new generalized Rayleigh distribution (NGRD) using the survival discretization method that depends on the survival function.See for example [6,7], in which the survival discretization approach was used to obtain the discrete normal and discrete Rayleigh distributions, respectively.Using the same approach, more discrete distributions have been introduced and studied; see for example [14][15][16][17][18][19][20][21][22][23].
Still, there is an enduring need to create and develop more discrete models and to generate new ones because of modeling and fitting real data, which appear and spread constantly in human life.The efficiency in discretization methods refers to the ability of a method to produce accurate and useful discretized versions of continuous data with minimal loss of information.Also, discrete distributions derived from continuous ones can inherit their flexibility and adaptability.This allows statisticians to model a wider range of data characteristics, such as skewness or kurtosis, which might be difficult with standard discrete distributions.In statistical methodology, continuous distributions may have characteristics that are missing in the discrete space; hence, creating discrete analogs can fill these gaps, providing new tools for data analysis.
The suggested discrete model with three parameters offers an immense degree of fitness to skewed, symmetric, monotone, and inverse J-shape types of data.Therefore, some statistical functions and properties are achieved, in addition to observing the submodels and limiting behavior of the proposed discrete distribution.Examining statistical inference is crucial; therefore, point and interval estimation for the unknown parameters using the maximum likelihood estimator and the Bayesian method is performed.
Simulation analysis via numerical techniques such as Monte Carlo simulation is employed to evaluate the estimators using the maximum likelihood and the Bayesian estimation methods to compare the performance of these two methods.The efficiency is assessed using the relative bias, the mean squared error, and the coverage probability of the confidence intervals.Two real datasets are analyzed to emphasize the empirical validation of the new model, where several goodness-of-fit measures are employed.The first example is related to the industrial field, where several strikes that occurred in coal mining in the UK were recorded over four weeks.Modeling and predicting the number of strikes will save human lives and money.The second example is related to the number of fires that occurred in Greece's forests in the year 1998 during the summer months.The main purposes of this study are, first, to introduce new discrete analogs of the continuous NGRD and evaluate some of its important statistical functions, second, to perform the inferential statistics related to the new distributions' parameters and compare the results, and, third, to assess the efficiency of the new discrete distribution by modeling real data examples and comparing the goodness-of-fit measures with other discrete distributions that were studied earlier in the literature.
The originality of this work emanates from the basis of exploring the creation of a new discrete analog from less commonly used continuous distributions and investigating their properties, potential applications, and how they compare to existing distributions.Also, it focuses on specific application areas, such as the industrial, engineering, and reliability fields.The remaining parts of this work are organized as follows: Section 2 describes the new generalized Rayleigh distribution.The discretization methods are presented in Section 3, along with some statistical functions.In Section 4, the maximum likelihood and the Bayesian inference are presented.In Section 5, simulation analysis and the tabulated results are carried out.Some real data examples are provided in Section 6.Finally, conclusions are remarked on in Section 7.

Model Description
The Rayleigh distribution (RD) is a continuous distribution that has much practical importance; hence, many of its statistical characteristics, inference, and reliability analysis have been studied by several authors, and numerous extended forms of the Rayleigh distribution have been proposed.For example, Ref. [24] applied the inverse Rayleigh to the failure times data.Ref. [25] introduced the transmuted Rayleigh and used it to model the amount of nicotine in the blood.In [26], the authors studied the beta-generalized Rayleigh distribution and its application.Ref. [27] obtained the transmuted inverse Rayleigh distribution to lifetime data.Ref. [28] obtained a new modified Rayleigh distribution named the Kumaraswamy generalized Rayleigh distribution with application to real data.For more information, refer to [29][30][31].In this work, we are interested in studying a new form of the Rayleigh distribution called a new generalized Rayleigh distribution (NGR), which was first introduced by Shen et al. [32].It has three parameters and it was shown that the NGR is suitable for modeling large data values rather than small data values.However, as a continuous distribution, it is restricted from describing discrete data forms.Discretizing the NGR distribution is our goal; therefore, it yields a subsequent distribution that accommodates the countable data while retaining the influential tail modeling characteristics of the NGR.In this study, we carry out a discrete version of the NGR and use it to model real data.
The probability density function (pdf ) and the survival function (S) of the continuous NGR are provided respectively as: and in which the parameters α > 1, β > 0, θ > 0.
The hazard rate function is (3) To identify submodels or special distributions that arise from this general form, we can consider different values or limits of parameters α, β, and θ.Here are some special cases: 1.
Standard Rayleigh Distribution: When θ = 1 and β = 1, the term (1 − e −x 2 ) simplifies to the CDF of the standard Rayleigh distribution.This is observed if the parameter α also approaches infinity, which simplifies the formula to 1 − e −x 2 , the CDF of the standard Rayleigh distribution.

2.
Exponential Distribution: If β approaches infinity, the term Λ = (1 − e −θx 2 ) β can approach an exponential-like behavior for small values of x, depending on how θ is defined.

3.
Modified Rayleigh Distribution: For specific fixed values of α and β, you can obtain various forms of modified Rayleigh distributions, where the behavior is measured by the degree of skewness and kurtosis determined by these parameters.

4.
Weibull-like Distribution: By interchanging between θ and β, especially when β is not equal to 1, the distribution can possess Weibull-like properties.
By discretizing the continuous range of x, discrete versions of this distribution can be derived, which may be useful for certain types of count data or integer-valued measurements.
Since our goal in this work is to define a new discrete NGR distribution, we generate a discrete analog based on the survival discretization method, which is denoted by DNGR.The pmf and cumulative distribution function (CDF) are obtained.Furthermore, the moments, stress-strength function, the mean residual, and mean past lifetimes, order statistics, and L-moments are obtained.All these statistical functions are used for studying the features of the DNGR.

The Discrete New Generalized Rayleigh Distribution
Roy and Gupta [3,4] defined the probability mass function (pmf ) for a discrete distribution using the survival function and it was expressed as follows: where S(x) is the survival function provided by Equation ( 2); hence, the pmf of the discrete analog of NGR distribution, namely DNGR, is written as where . The CDF of the DNGR distribution can be written as The quantile function with given values of parameters as α, β, and θ of the DNGR distribution is provided by The hazard rate function (HRF) of the DNGR distribution is provided by We also observe that the reversed hazard rate function for the DNGR of this distribution is provided by In Figure 1, the bar charts represent each parameter α, θ, and β that has a specific role in the behavior of the pmf, and their effects are observable when we fix one and vary the others.An explanation of the effect of each parameter based on the plots is as follows: 1. Effect of α when θ and β are changeable: When α is fixed, the variations in θ and β create different trends in the probability values.Higher values of θ tend to stretch the curve horizontally, meaning, for a given α, as θ increases, the decrease in probability values with increasing k is less steep.Higher values of β tend to amplify the curve vertically, making the probability have fewer values for higher k values.The reaction between θ and β at a fixed α demonstrates that θ affects the spread of the distribution, while β affects the sharpness of the probability decrease.2. Effect of θ when α and β are changeable: With θ fixed, the changes in α and β show distinct patterns.An increase in α generally results in higher probability values across all k.This is because a higher α relative to Λ(k; θ, β) increases the numerator and decreases the denominator of the function P(X = k), resulting in a larger overall value.The effect of β at a fixed θ is similar to its effect when α is fixed; it controls the sharpness of the decrease in probability values.
Higher β values cause a quicker decline in probability as k increases.

Fixing β and varying α and θ, we can see that
As α increases, for a fixed β, the overall probability values increase, similar to when θ is fixed.The role of θ here is nuanced; for lower values of k, the impact of changing θ is minimal, but, as k increases, higher θ values preserve higher probabilities, indicating a wider spread in the distribution.From the above explanations, it is clear that α primarily scales the probability values, β determines the rate at which the probability values decline as k increases (sharpness of the distribution), and θ controls the spread or dispersion of the distribution across k values.The combination of these three parameters can thus shape the function's distribution in various ways, and each has a distinctive role in the form of the probability curve.In Figure 2, the plots represent the HRF of the DNGR distribution for various combinations of parameters α, θ, and β.Each subplot corresponds to a different set of these parameters.The values of k range from 1 to 10.The curves are increasing for different values of the parameters; we can realize the effect of increasing the parameter θ while keeping other parameters fixed by going steeply to the left.For the effect of β, assuming other parameters are fixed, it can be figured by increasing in a slower mode when k takes small values.Finally, the effect of increasing the values of α while fixing the remaining parameters is going to the left more steeply.
From the above limiting behavior of the DNGR, some submodels and special cases can be derived, such as 1.
Discrete standard Rayleigh Distribution: When θ = 1 and β = 1, and α approaches infinity, the pmf simplifies to (1 ), which represents the pmf of the discrete Rayleigh distribution created from applying the survival discretization method.

2.
Discrete Exponential-like Distribution: For large values of β and specific values of θ, the DNGR distribution might possess characteristics similar to an exponential distribution for smaller values of k, where the exponential decay behavior is more evident, since the term Λ = (1 − e −θk 2 ) β has a decaying form and can be considered as exponential-like function.

3.
Discrete Uniform Distribution: If the parameters α, β, and θ are chosen such that the pmf becomes constant for all k within a certain range, the DNGR distribution could approximate a discrete uniform distribution.4.
Geometric-like Distribution: By adjusting θ and β, you might be able to create a distribution that behaves like the geometric distribution, especially if the probability of larger k values decays like the geometric series.
These possible submodels and special cases demonstrate the versatility and adaptability of the DNGR distribution.The ability to derive such a variety of distributions from a single distribution highlights the potential utility of the DNGR distribution in modeling a wide range of discrete data scenarios.Each submodel or special case would be suited to different types of data and could provide unique insights depending on the context of the analysis.

Moments
Assume a non-negative random variable k ∼ DNGR(α, β, θ).The s th moment, say ψ ' s , can be expressed as follows and then It is impossible to write an exact form of the s th moment; hence, R programming with version (4.3.0) is helpful, and the moment is evaluated numerically.Equation ( 11) is convergent for α > 1, β > 0 , and θ > 0 .
Table 1 explores some functions like the minimum, mean, variance, maximum, skewness (SK), and kurtosis (Kt) for different values of α, β, and θ.In addition, the DNGR distribution is appropriate for modeling both over-and under-dispersed data since, in this model, the variance can be smaller than the mean, which is the case with some standard classical discrete distributions, in addition to the positive and negative skewness values, which show that this distribution can be skewed to the right or left.Also, a very small skew value that tends to zero indicates a symmetry possible curve for the pmf.A higher kurtosis means more of the variation is due to infrequent extreme deviations as opposed to frequent modestly sized deviations.By varying θ, α, and β, one can realize the distribution changes.For instance, with θ = 0.8 and α = 0.5, β changing from 0.84 to 2.73 drastically increases the kurtosis, indicating a heavier tail.

Stress-Strength Analysis
The stress-strength (reliability) analysis is an important tool in mechanical design.The idea relies on the probability of failure that is obtained from the probability of r exceeding r * .Assume that both r and r * are in the positive domain.The expected reliability (R * ) can be calculated by in which K r ∼ DNGR(α 1 , β 1 , θ 1 ) and K r * ∼ DNGR(α 2 , β 2 , θ 2 ), and then R * can be ex- pressed as follows where We cannot obtain a closed form for the above equation; consequently, simulation analysis is utilized to obtain a numerical solution.In Section 6, numerical analysis is performed to obtain the value of the stress-strength function under two real data applications.

The Mean Residual and the Mean Past Lifetimes
In reliability and survival analysis, many lifetime measures have been discussed in the literature.They were defined to study the aging behavior of the experimental units.One of these measures is the mean residual lifetime (MRL), which is a helpful tool in determining burn-in and maintenance policies.For discrete distributions, the MRL is defined as follows: where 0 < l < ∞.
If the random variable k follows the DNGR distribution with parameters α, β, and θ, which is denoted by k ∼ DNGR(α, β, θ), then the MRL is expressed as The mean past lifetime (MPL) is another important measure in reliability analysis.The MPL measures the time elapsed after the failure of K units given that the system has failed sometime earlier to i.In the discrete case, the MPL is defined as follows where ζ * (0) = 0; see [33].Then,

Order Statistics
Let K 1 , K 2 , . . ., K n be a random sample with the DNGR distribution and K 1:n , K 2:n , . . ., K n:n denote the corresponding order statistics.Then, the CDF of i th order statistics at the value k can be written as follows By using the negative binomial theorem, then Therefore, Consequently, the pmf of the i th order statistic under the DNGR can be derived and expressed as follows So, the υ th moments of k i:n can be written as follows

Estimation
Two estimation methods are considered in this work: frequentist maximum likelihood estimation (MLE) and the Bayesian estimation method.Simulation analysis and numerical techniques are performed in Section 5 to assess the performance of these estimation methods.

Maximum Likelihood Estimation
In this section, we use the maximum likelihood estimation MLE method to estimate the unknown parameters of the DNGR distributions.To evaluate the required estimators, numerical techniques are used, such as the well-known Newton-Raphson technique.
Let x 1 , . . ., x n be a random sample following the DNGR, and then, from pmf in Equation (5), the log-likelihood function is written as The MLEs for α, β, and θ are obtained by finding the partial derivatives of ℓ(α, β, θ) for α, β, and θ, then equating the three equations to zero and solving numerically. and . To solve the system of nonlinear Equations ( 21)-( 23), only numerical methods are helpful.Many numerical techniques were used in the literature; here, we use the Newton-Raphson method, and all results are illustrated in Section 5.

Bayesian Inference
The Bayesian estimation method is used in this section to estimate the unknown parameters of the DNGR.The basic assumption of the Bayesian method is that the model parameters are considered random variables that follow a distribution known as the prior distribution.Since prior information is usually only available, we must specify a suitable prior option.We choose the gamma conjugate prior distribution for the parameters α, β, and θ.It is defined by assuming gamma distributions for α, β, and θ.
Therefore, the prior distributions for α, β, and θ can be written as where a 1 , a 2 , a 3 , b 1 , b 2 , and b 3 are nonnegative hyper parameters of the assumed distributions.Hence, the joint prior for α, β, and θ is The joint posterior of α, β, and θ given the data is defined as where L(x | α, β, θ) is the likelihood function of the DNGRD and The DNGRD parameters are estimated using a symmetric squared error (SE) loss function.A simulation study is used to investigate the performance of the estimators using the aforementioned loss function.As criteria for the superiority of the estimation methods, the bias, the mean square error (MSE), the average length (AL) of the confidence intervals, and the coverage probability (CP) are computed.
Under the SE loss function, Bayes estimation for the parameters α, β, and θ is defined as the mean or expected value regarding the joint posterior, provided as To evaluate the expected values and triple integration in Equations ( 26)-( 28) , numerical methods are required.We choose to use the Markov chain Monte Carlo (MCMC) technique via the Gibbs sampling method and by developing appropriate R code.The joint posterior density is Bayes estimation for parameters α, β, and θ under SE loss function is performed respectively using Equations ( 26)-( 28) and the posterior density Equation (29).
The estimators are numerically evaluated simulations using R codes under the SE loss function, and their results are summarized and presented in Tables 2 and 3.

Simulation Analysis
In this section, Monte Carlo simulations are performed to assess the effectiveness of the suggested estimators for the parameters α, θ, and β that were established in Section 4.1.We will sum up by providing the simulation scenario.The findings of the simulation are then offered for debate.

Simulation Scenario
In this subsection, several Monte Carlo simulation studies are carried out to assess the effectiveness of the acquired maximum likelihood estimates and Bayesian estimation of α, θ, and β.Now, we suggest the following steps to gather a sample from the DNGR model:

•
Set α, θ, and β to their actual values as shown: In Table 2: In Table 3: Specific values for n (total test units) should be determined as 30, 70, 100, 200.• Generate a uniform random variable within the range of 0 to 1. Utilize the quantile function described in Equation ( 7) to produce a random sample from the DNGR(α, θ, β) distribution.Afterwards, round the quantity of samples to the nearest whole number.The relative bias (RB), mean squared error (MSE), average lower, average upper, and coverage probability (CP) of the parameter are specifically determined for each group (n, or actual value of the parameter).For more details about comparing interval estimates, we discuss using the CP requirement in our evaluations.R 4.2.2 programming language is used to carry out all numerical analyses.In Tables 2 and 3, respectively, all numerical findings for α, β, and θ are obtained and presented.

Simulation Conclusion
The performance of the suggested point and interval estimate algorithms is the main topic of this subsection.We can infer the following facts from Tables 2 and 3:

•
The acquired estimates of the unknown parameters α, θ, and β generally perform well in terms of lowest MSE, RB, and difference between upper and lower values with CP.

•
The MSE, RB, and CI of α, θ, and β tend to decline as n rises.This result supports the associated estimates' consistency property of DNGR distribution when the necessary sample size is raised.

•
As the true value of β increases, for each setting, the MSE, RB, and CL measures of unknown parameters α and β decrease, while they increase regarding unknown parameter θ.

•
The MSE, RB, and CL measures of all unknown parameters α, θ, and β increase for each set as the true value of θ grows.

•
For CI of Bayesian, the credible interval decreases when the sample size increases.

•
Almost always, and regardless of sample size, Bayesian estimation based on the SE loss function yields the minimal RB and MSE values.

Real Data Examples
This section presents the analysis of two applications using different real datasets.The main goals of this section are

•
Examine the usefulness and applicability of the proposed model to real phenomena; • Show the applicability of the inferential results to a real practical situation; • Evaluate whether the proposed model is a better choice than the other seven models.
Based on the first and second datasets, the DNGRD probability model is contrasted and compared with the other seven competing models to show the reliability and superiority of the proposed model, including Poisson, binomial, geometric, discrete Burr (DB) by [37], discrete Marshall-Olkin Lomax (DMOL) by [38], new discrete Lindley (NDL) by [39], and discrete odd perks exponential (DOPE) by [40] distributions.To specify the best model, several criteria are used, namely: Akaike (AIC = 2p − 2 l), where p is the length of the model parameter and l is log-likelihood value, consistent Akaike (CAIC = −2 l + 2np n−p−1 ), Bayesian (BIC = −2 l + p ln(n)), and Hannan-Quinn (HQIC = −2 l + p ln(ln(n)) information criteria.Along with these, the X 2 − square statistic and its p-value are taken into account.If a probability model distribution has the highest p-value and the lowest values for all other metrics, it is obvious that it will provide the best fit for a particular collection of data.The maximum likelihood estimates (with their standard errors (St.Es)), as well as the fitted model selection criteria, are shown in Tables 5 and 6 using the R programming language and the 'bbmle' package in R program with version number (4.3.0), that was recommended.
To compare the performance and efficiency of the DNGR distribution with other distributions listed in Table 5, using measures of goodness of fit and p-values, we can proceed as follows: 1-DNGR versus DOPE: DNGR shows a better fit with a lower AIC, CAIC, BIC, and HQIC.The chi-squared value is lower for DNGR, indicating a better fit.DNGR has a higher p-value, suggesting a better fit to the data than DOPE.
2-DNGR versus Binomial/Poisson: DNGR has a higher p-value than both binomial and Poisson distributions, indicating a more suitable model for the data.The information criteria (AIC, CAIC, HQIC) for DNGR are lower compared to binomial and Poisson, suggesting a better fit than binomial and Poisson.
3-DNGR versus DMOL: DNGR and DMOL have comparable p-values, but DNGR shows better performance in terms of information criteria.
4-DNGR versus DB: DNGR has a higher p-value than DB, indicating that DNGR shows slightly better performance in terms of information criteria.
5-DNGR versus Geometric/NDL: DNGR outperforms both geometric and NDL distributions in terms of p-value, indicating a significantly better fit.DNGR has lower information criteria values, further suggesting its superiority in model fitting.Overall, DNGR appears to offer a more efficient and suitable fit for Data I compared to the other listed distributions.6 that the DNGR distribution is the best distribution among all the examined models in terms of the Pvalue, whereas Figures 6-8 confirm these results for Data II. Figure 3 confirms the results of MLE fitting and demonstrates the existence, uniqueness, and maximum point of likelihood value of the likelihood estimates for Data I. Figure 4 regarding associated empirical CDF and estimated CDF plot illustrates the connection between observed cumulative probability and observation through a visual plot and also Q-Q plot for Data I. Figure 5 highlights estimated frequency by using PMF for each comparative model for Data I. Table 7 indicates survival and hazard rate functions for DNGR distribution with different values of Data I, noting that the survival value decreased when the values of Data I increased, while the hazard rate value increased when the values x of Data I increased.
Figure 6 confirms the results of MLE fitting and demonstrates the existence, uniqueness, and maximum point of likelihood value of the likelihood estimates for Data II. Figure 7 regarding associated empirical CDF and estimated CDF plot illustrates the connection between observed cumulative probability and observation through a visual plot and Q-Q plot for Data II. Figure 8 highlights estimated frequency by using PMF for each comparative model for Data II.Table 8 indicates survival and hazard rate functions for DNGR distribution with different values of Data II, noting that the survival values decreased when the values of Data II increased, while the hazard rate value increased when the values x of Data II increased.

Conclusions
In this study, the authors successfully developed a novel discrete analog from the continuous generalized Rayleigh distribution denoted by DNGR through the application of the survival discretization method.Several key attributes of the DNGR model were studied, such as its unimodal probability mass function, which exhibits varying degrees of symmetry and skewness based on parameter selection.Comprehensive statistical measures for DNGR were derived, including moments, stress-strength function, momentgenerating function, and mean residual and mean past lifetimes.The potential submodels and special cases derived from the DNGR demonstrate the versatility and adaptability of the DNGR distribution, which can be suitable for modeling different types of data and could provide unique insights depending on the context of the analysis.Furthermore, the practical applicability of the work was enhanced by conducting detailed simulation analyses and presenting the results in tabular form.Point and interval estimation using both the maximum likelihood and the Bayesian methods were obtained, supplementing these with simulation analyses executed using R code.This was complemented by a numerical analysis aimed at evaluating the estimation methods for DNGR's unknown parameters and assessing the efficiency of these methods.A significant aspect of their contribution is the application of the DNGR model to real-world data.Two real data examples were selected, one from the industrial sector concerning UK coal mining strikes and another focusing on environmental issues related to fires in Greece.Their analysis revealed that the DNGR model outperformed seven competitive discrete distributions in various goodness-of-fit measures, demonstrating its superior ability to model the given datasets effectively.This finding was further illustrated through detailed tables and figures showcasing the properties and efficiency of the new model.As a pathway for future work, the authors suggest exploring alternative discretization methods to assess their performance and applicability to a broader range of real-life data scenarios.

Figure 2 .
Figure 2. The HRF of the DNGR distribution.

Figure 3 .
Figure 3. Likelihood profile (blue line) with the maximum likelihood estimation (red dot): Data I.

FFigure 4 .
Figure 4.Estimated CDF and Q-Q plot of DNGR by using MLE: Data I.

Figure 8 .
Figure 8. Likelihood profile (blue line) with the maximum likelihood estimation (red dot): Data II.
To our knowledge, no previous work has studied the discrete new generalized Rayleigh distribution and employed it to model real-life data examples from different scientific fields.

Table 1 .
Summary of descriptive statistics for the DNGR distribution.
• Compute the MLEs and 100(1 − γ)% via 'maxLik' package in R program with version number (4.3.0), with Fisher information matrix (Hessian matrix).• Use 'coda' package in R program with version number (4.3.0), to obtain the Bayes' inferences by running the MCMC sampler 12,000 times and 2000 is burn-in.• Repeat the above steps 5000 times.•

Table 4 .
Data I: The number of strikes and their frequency that occurred in the UK coal mining industry.

Table 5 .
MLE estimates and different measures of fit for Data I.

Table 6 .
MLE estimates and different measures for Data II.these results for Data I (the black point refer to data; the pink point refer to DNGR distribution).Additionally, it is evident from Data II in Table

Table 7 .
Survival and hazard rate functions for DNGR distribution with different values of Data I.

Table 8 .
Survival and hazard rate functions for DNGR distribution with different values of Data II.