Skip to main content

pdpublishers.com

Article ID: PD2601208008

Views: 806
Volume 1 (2026)
Published 01 Oct 2026

Information and Uncertainty Analysis of the OBXII-IRD Using Seven Measures of Entropy with Application to Carbon Fibers Data

📚 Cited by: 0

⬇ Downloads: 37


Authors

1Ministry of Education, Diyala Education Directorate, Diyala, Iraq

2Department of Mathematics, College of Computer Science and Mathematics, Tikrit University, Tikrit, Iraq

3Department of Mathematics College of Education for Pure Sciences (Ibn AL-Haitham), University of Baghdad, Iraq

4Department of public Administration, Collage of Administration and Economics, Tikrit University, Tikrit, Iraq

5Department of Mathematics, College of Education for Pure Science, University of Tikrit, Iraq

Article History:

Received: 04 January, 2026

Accepted: 16 July, 2026

Revised: 14 July, 2026

Published: 01 October, 2026

ABSTRACT:

Introduction: This research studies the Odd Burr Inverse Rayleigh (OBXII-IRD) distribution, concentrating on parameter estimation methods and entropy measures.

Methodology: Three estimation methods: Maximum likelihood estimation, maximum product of spacings, and Kolmogorov were used and their outcomes were investigated using Monte Carlo simulation.

Results: The results demonstrated their ability to estimate parameters reasonably well. Seven entropy measures including Renyi, Shannon, Tsallis, Arimoto, Havrda-Charvat, Sharma-Mittal, and Sant’ Anna-Taneja were subject to distribution’s uncertainty. To show its practical perfomance, the OBXII-IRD was utilized to real-world data representing carbon fibers.

Conclusion: The results indicate that the OBXII-IRD is superior and fits the data better than other distributions. The statistical software R was used for all programming and computational tasks in this study.

Keywords: Inverse rayleigh distribution, non-bayesian estimation, MLE, information criterion, goodness-of-fit statistics, carbon fibers.

1. INTRODUCTION

Trayer [1] introduced the inverse Rayleigh distribution model for checking the reliability and survival data, as followed by Voda [2], whom investigated the distribution’s basic properties and acknowledged the potential for approximating the lifetime distribution of experimental units. Later, Graf [3] compared the parameter estimators for inverse Rayleigh and derived five measures of location the arithmetic, harmonic, and geometric means, as well as the mode and median.

Howlader et al., [4] determined Bayesian prediction limits for lifetime models based on the inverse Rayleigh distribution. Banerjee and Bohoni [5] introduced the exponentiated transformed inverse Rayleigh distribution and investigated various estimation methods. Khan and King [6] addressed the Transmuted modified inverse Rayleigh distribution and explored its statistical properties. Meanwhile, Goual and Yousef [7] proposed a modified inverse Rayleigh based on the Burr XII distribution, originally introduced by Cordeiro et al. [8] contributing to the study of its mathematical properties and associated regression model. Ali [9] study the inverse Rayleigh distribution mixture, determined its properties, and employed Bayesian estimation to estimate its parameters.

Fisher [10] and Rao [11] established the framework for weighted distributions, while Aydin [12] introduced the weighted inverse Rayleigh distribution, their comprehensive study of the mathematical properties of these distributions contributed to a better understanding of weighted distributions in statistical analysis. For further information on exponential and weighted distributions, as well as Bayesian estimation methods, please refer to the sources [13–17].

Cover and Thomas [18] provided a comprehensive foundation for information theory, building upon the work of Shannon [19], who established the concept of entropy to measure uncertainty. To explore different types of entropy, consult see the refer [20–25].

In this article, we extend the Odd Burr XII- inverse Rayleigh distribution, for which Khalaf et al., [26] established statistical characteristics and estimated parameters via the maximum likelihood estimation, by introducing the following novel contributions:

First, present three different classical methods for estimating the parameters of the OBXII-IR distribution, employing various sample sizes to ensure numerical consistency and methodological robustness.

Secondly, as our core theoretical contribution, we provide a thorough evaluation of uncertainty from an information-theoretic standpoint and compelling inaugural seven diverse generalized entropy measures, thereby providing a comprehensive and unprecedented characterization of model randomness from multiple mathematical viewpoints.

In the end, we investigate the practical implication of the model using real-world engineering data from materials science (carbon fibers) to present our developed estimation and entropy frameworks providing superior and highly competitive modeling in contrast to other leading distributions.

The rest of the paper is sectioned into 7 sections. Section 2 contains the mathematical part of the probability density function and the cumulative distribution function for the OBXII-IR distribution. Section 3 demonstrates the informational area by evaluating seven different measures of entropy to quantify the uncertainty in the model. Section 4 reviews the three estimation methods used to extract the parameters of the distribution, followed by Section 5, investigates the intensive simulation experiments to evaluate the efficiency of these estimators.  In Section 6, real world implication is discussed, and lastly in Section 7 concludes the study.

2. THE ODD BURR INVERSE RAYLEIGH (OBXII-IRD) DISTRIBUTION

Below are the cumulative distribution (cdf) and density function (pdf) of the OBXII-IR distribution Eq: (1 & 2):

where .

Fig. (1) depicts the probability density function (pdf) of the OBXII-IRD under different combinations of shape and scale parameters, reflecting the high structural flexibility of the model in modeling real data:

Fig. (1). Pdf plot for OBXII-IRD.

  • Effect of the β parameter (first line): When the parameters (λ = 2.5, ρ = 2) are held constant and the value of the scale parameter β is increased from 0.5 to 1.3, a clear shift of the function’s curve to the right is observed, along with an increase in the curve’s flattening and the extension of its tail.
  • In terms of the shape parameter  and its behavior as the parameter  increases, while keeping  constant.  At low values, we observe that the curve’s skewness is positive; as the value of beta increases, the function becomes steeper.
  • In terms of the shape parameter  and its behavior as the parameter  increases from 0.5 to 1.3, while keeping  constant observe the peak becoming sharper and the horizontal spread of the curve contracting, this indicates that  acts as a shape parameter, controlling the expansion and contraction of the function.

Fig. (2) illustrates the behavior of the cumulative distribution function (cdf) for the proposed OBXII-IRD. All curves reflect the standard mathematical properties of any continuous and true cumulative distribution function:

Fig. (2). Cdf plot for OBXII-IRD.

  • Checking the mathematical boundary conditions: It is clearly shown from all the graphs that the function starts precisely from zero at the minimum space , and increases monotonically, continuously and smoothly to reach the maximum integer one , thus proving the mathematical validity of the deduced equations.
  • The elasticity of S-shaped curves: The cumulative function takes on elastic forms, oscillating between sharp arrowhead curves and soft S-shaped curves, depending on the parameter values. We observe that as the values of the parameters β and λ increase, the function becomes more concave and the later it reaches one, which perfectly matches the heavy-tail behavior observed in the corresponding PDF graphs in Fig. (1).

3. VARIOUS ENTROPY MEASURES

As mentioned in the introduction, the entropy of the OBXII-IRD can be measured in various ways. This allows for the analysis and understanding of the uncertainty in the OBXII-IRD, and such analyses are important in machine learning and information theory. Due to the complexity of the OBXII-IRD function, the equations cannot be written in closed form; therefore, they are expressed as integrals and can be calculated numerically or using an Arithmetic programming language.

3.1. Rényi Entropy

The Renyi entropy is considered a generalization of Shannon entropy and serves as an important measure of uncertainty for probability distribution, the dispersion of data in OBXII-IRD can be described using the  order parameter. It can be expressed by the following equation Eq (3):

Substituting equation (2) into the integrand yields

3.2. Shannon Entropy

Shannon entropy is one of the most widely used measures of dispersion in information theory, and for the OBXII-IRD, it can use the following formula Eq (4):

Using (2), can be derived

3.3. Tsallis Entropy

The entropy Tsallis of the OBXII-IRD is calculated using the following formula Eq (5):

Substituting equation (2) into the integrand obtain:

3.4. Arimoto Entropy

Arimoto entropy prosses a highly flexible formula for handling complex probability distribution, therefore, it is considered a fundamental measure in information theory. It can be calculated for the OBXII-IRD using the following equation:

From (2) we obtain Eq (6):

3.5. Havrda and Charvat Entropy

The Havrda-Charvat measure was introduced as a generalized of Shannon entropy to account for changes in probability distribution. It is significant for describing in complex systems via the order parameter . For the OBXII-IRD. It is calculated using the following equation Eq (7):

provided by (2):

 

3.6. Sharma and Mittals Entropy

Sharma and Mittal entropy is one of the most comprehensive and flexible information measures, developed as a general framework combining the properties of Rennie and Tsalis entropy. Sharma and Mittal entropy is calculated using the following equation Eq (8):

From equation (2) we get:

 

3.7. Sant’ Anna and Taneja Entropy

The Sant’ Anna-Taneja entropy measure as it combines two entropy concepts to confirm a stabilize analysis subject to uncertainty. It is highly potent for measuring data sensitive to dispersion and described the OBXII-IRD by the following equation Eq (9):

By equation (2), we obtain:

Tables 1, 2, and 3 shows an in-depth evaluation of the seven entropy measures for the OBXII-IRD, presenting the degree of uncertainty as the distribution’s shape and scale parameters differs.

Table 1. Analysis of different entropy values for the OBXII-IRD assuming that ρ =1.1, β =2, and α=0.8.

λSETEREAEHCESMESTE
0.52.38603965.57764881.71876304.09029355.08953324.11281681.0219391
0.71.86711143.93679751.37036752.85084093.61784753.04441851.0174308
0.91.54192783.11397451.13589972.20918622.83531012.42924481.0123078
1.11.31382402.62111150.96360361.80824772.33818252.02101051.0068136
1.31.14238522.29258250.82967451.52999541.98913501.72621451.0011308
1.51.00723082.05737650.72143241.32344211.72776091.50102520.9953748
1.70.89735781.88020820.63137081.16275221.52301721.32194080.9895968
1.90.80527191.74156900.55478751.03335511.35721511.17519990.9838479
2.10.72674371.62985200.48848360.92634451.21945541.05208360.9781567
2.30.65862151.53769470.43028480.83597621.10266590.94685420.9725333
2.50.59875581.46020810.37857080.75835451.00198860.85552800.9669901
2.70.54552791.39400540.33218900.69074790.91402630.77525780.9615580
2.90.49775991.33671530.29021720.63114880.83627730.70394660.9562022
3.10.45454501.28564640.25196540.57807980.76688810.64000420.9509409
3.30.41516521.24220770.21268660.53042200.70443600.58223460.9457827
3.50.37906881.20266920.18453210.48730550.64781100.52966780.9407178
3.70.34578781.16715050.15456110.44803880.59614520.48154940.9357408
3.90.31496571.13054760.12666890.41206320.54873830.43725810.9308567
4.10.28629551.10584450.10060140.37893050.50501040.39630700.9260633
4.30.25951361.07914520.07616800.34829510.46451550.35827860.9213524
4.50.23443651.05462350.05318700.31982310.42683640.32281970.9167352
4.70.21084771.03199890.03149730.29328100.39166690.28965830.9121970
4.90.18860971.01104600.01098210.26845050.35891870.25852330.9077211

 

Table 2. Analysis of different entropy values for the OBXII-IRD assuming that λ=0.6, β=2.5, and α=0.8.

ΡSETEREAEHCESMESTE
0.53.61693829.80699012.28309509.166857910.7183977.13797151.0249472
0.72.95039237.76736422.04993075.83090677.08265555.40779421.0241710
0.92.51833606.23751931.83058284.39194025.44063824.4034151.0229331
1.12.20332915.13736301.63653993.55866354.46424823.72392381.0212614
1.31.95727144.33387141.46646162.99988143.78746233.22216291.0191607
1.51.75624553.73101751.31668042.59109083.30291472.83016761.0166272
1.71.58676863.26629621.18365652.27444672.91562252.51171631.0136615
1.91.44053942.89923461.06444622.01913862.60050082.24549411.0102648
2.11.31211222.60310820.95670681.80710032.33675292.01801541.0064353
2.31.19772432.35981140.85858131.62696212.11116571.82026771.0021786
2.51.09467812.15674770.76860061.47116301.91489331.64595780.9974959
2.71.00097451.98493680.68558691.33445521.74174801.49053990.9923917
2.90.91509171.83783390.60858771.21306251.58725611.35062960.9869705
3.10.83584981.71056630.53682441.10419041.44808641.22362550.9809393
3.30.76231191.59944710.46965711.00571881.32169771.10760240.9746031
3.50.69372461.50163340.40655320.91600271.20611451.00089110.9678695
3.70.62947241.41490420.34706130.83374721.09977270.90224620.9607459
3.90.56904721.33750160.29080320.75791531.00141280.81062500.9532396
4.10.51202431.26801520.23745200.68766340.91001030.72517210.9453595
4.30.45804511.20530440.18673110.62229780.82471960.64517430.9371134
4.50.40680451.14843380.13839840.56124290.74483490.57002970.9285121
4.70.35804081.09663170.09224350.50401320.66976170.49922780.9195655
4.90.31152621.04925520.04808070.45020200.59899840.43233390.9102819

 

Table 3. Analysis of different entropy values for the OBXII-IRD assuming that λ=0.9, ρ=1.1, and α=0.8.

βSETEREAEHCESMESTE
0.50.84877401.55698540.44275101.22128011.59773751.24422600.9689681
0.71.01701021.84225060.61098731.44556521.88253921.51683120.9852560
0.91.14266732.08891560.73664491.61935042.10160191.72668830.9944222
1.11.24300242.30938290.83698081.76208702.28051571.89800201.0002993
1.31.32653052.51056350.92050761.88367562.43222082.04326341.0043874
1.51.39808042.69677740.99205771.98986692.56420352.16964031.0073965
1.71.46066262.87093961.05463952.08431772.68120042.28166861.0097036
1.91.51627443.03512321.11025122.16950082.78640532.38240551.0115288
2.11.56631663.19087041.16029352.24716792.88207742.47401381.0130087
2.31.61180253.33936131.20577982.31861242.96987302.55808071.0142324
2.51.65349343.48152501.24747022.38481423.05104922.63580881.0152611
2.71.69197333.61810621.28595032.44653373.12657722.70812771.0161392
2.91.72770393.74971741.32168002.50437423.19722772.77577761.0168965
3.11.76104873.87686291.35502612.55882383.26362142.83935141.0175568
3.31.79230983.99996781.38628602.61028223.32626692.89933521.0181369
3.51.82172994.11939591.41570612.65908063.38558382.95613211.0186516
3.71.84951484.23545731.44349102.70549773.44192473.01008031.0191093
3.91.87583624.34842301.46981362.74976783.49558863.06146581.0195212
4.11.90084124.45857761.49481842.79209493.54682933.11053001.0198933
4.31.92465534.56597701.51863282.83265203.59586903.15748631.0202314
4.51.94738644.67095511.54136372.87159123.64289753.20251741.0205836
4.71.96912964.77362641.56310682.90904413.68808023.24578191.0208207
4.91.98996584.87413411.58394292.94512863.73156513.28741971.0210789

First: Table 1 shows the conduct of the seven entropy and its result of slowly increasing the value of  while keeping the parameters  fixed and .

  • From the Table 1 and Fig. (3), we inspect that all seven entropy measures decrease with increasing  For example, Reny entropy values decline abruptly from 1.1787 at  to 0.0109 at , and likewise, STE entropy decreases from 1.0219391 to 0.9077211.
  • From a mathematical point of view, we notice that the PDF constrict and reduces as the value of the parameter  This geometric outcome shows a decline in data dispersion, signifying a decrease in the system’s level of uncertainty.

Fig. (3). It illustrates the behavior of entropy measures and their changes at different values of the parameter λ.

Second: Table 2 depicts the behavior of the seven entropy and effect of gradually increasing the value of  while keeping the parameters  fixed and .

From the Table 2 and Fig. (4), we observe that all seven entropy measures decrease with increasing  value. For example, Shannon entropy values decline sharply from 3.6169382 at   to 0.3115262 at .

Fig. (4). It illustrates the behavior of entropy measures and their changes at different values of the parameter ρ.

From a mathematical view point, it’s been observed that the PDF constrict and cinch as the parameter  surge, resulting in data clustering around the model. This geometric inclination proposes a decrease in data dispersion, which indicates a lower degree of uncertainty in the system.

Third: Table 3 displays the behavior of the seven entropies and the impact of  slowly increasing the value of  while maintaining fixed values for  and α=0.8.

Table 3 and Fig. (5) show that when the  value increases, each of the seven entropy measures decrease. Arimoto entropy values, for instance, decrease dramatically from 1.2212801 at   =0.5 to 2.9451286 at   =4.9.

Fig. (5). It illustrates the behavior of entropy measures and their changes at different values of the parameter β.

Regarding the statistical interpretation, we observe that the situation for the parameter  beta differs from that of the parameters  here, increasing the values of  led to an increase in all entropy measures due to the heightened randomness of the system and the dispersion of the data.

Key Conclusions: All seven entropy measures showed complete consistency and logical agreement in their response to changes in the parameters of the OBXII-IRD, thus proving the mathematical stability of the model. Rényi entropy and Shannon entropy are the most sensitive and ideal in detecting and characterizing these dynamic shifts with extreme accuracy compared to the Sant’Anna-Taneja (ST) measure, which showed a slow response and very slight changes around the value (1).

4. NON-BAYESIAN ESTIMATION METHODS

This section reviews the classical inference methods for estimating the parameters of the OBXII-IRD, which include: Maximum Likelihood Estimation (MLE), Approximate Confidence Interval (ACI), Maximum Product space Estimation (MPS), and Kolmogorov Estimation (KE).

This section will identify the Inference Classica estimators for the OBXII-IRD parameter set .

4.1. MLE

Let   be a random sample of size n drawn from the OBXII-IRD; accordingly, the Likelihood Function denoted by L according to the sources [28–31], can be expressed as follows Eq: (10).

The derivation of the log-likelihood function for  is as follows Eq: (11).

 

The potential log-likelihood function is partially derived with respect to the coefficients  to obtain the maximum potential estimators Eq: (12 & 13).

 

Where  solve three equations at the same time to get an approximation for λ, ρ, and β. The integral of  becomes zero when  also becomes zero. Using statistical software like R, one may achieve the necessary results.

4.1.1. ACI

Here we examine the most common confidence intervals; the asymptotic properties of MLEs were used to calculate the range of confidence intervals at a significance level of 100(1 – α) % for the unknown parameters λ, ρ, and β. If the regularity conditions are met, the MLEs follow a normal asymptotic distribution with a mean of (λ, ρ, β) and a variance-covariance matrix represented by , i.e., that:

Assuming the existence of a variance and covariance matrix, denoted by:

The asymptotic variances for λ, ρ, and β.are determined by the elements on the main diagonal of the  matrix. Accordingly, confidence intervals at a significance level of 100(1 – α) % for λ, ρ, and β can be calculated using normal approximation according to the following equations:

4.2. MPSE

It minimizes the MPS function to estimate λ, σ, and β for further details on MPS, see to [32].

From equation (1) we get Eq: (15).

4.3. KE

The Kolmogorov estimation is used in the process of estimating the parameters for the OBXII-IRD. The object of this estimating method is to find the minimum value of the following equation [33].

From equation (1) we get Eq: (16).

 

5. SIMULATION

To assess the statistical effectiveness and numerical efficiency of the different estimation methods proposed in Section (5) for the OBXII-IRD, a detailed Monte Carlo simulation study was introduced. This study was formulated to ensure the accuracy and reproducibility of the results by integrating the following methodological aspects:

  • Number of Replications: N = 1000 independent random samples were originated for each simulated scenario to maintain the stability of the estimates and minimize random error.
  • Testing five different sample sizes, namely:30, 60, 120, 200, and 300.
  • True parameter values:

Group 1: (λ= 0.7, ρ=3.2, β=1.1).

Group 2: (λ= 0.7, ρ=3, β=0.8).

Group 3:( λ= 0.8, ρ=3.5, β=1.2).

To examine the properties of the estimators, three cases were selected with different values of skewness, kurtosis, and measures of dispersion. This variation ensures an assessment of the robustness of the estimation methods for the OBXII-IRD.

5.1. Statistical Evaluation Metrics

1. The numerical average (Mean): 

Mean= 

2.   Bias: Measures the range to which the estimated value deviates from the true value:

Bias = 

3. Root Mean Square Error (RMSE):

RMSE  = 

4. Convergent Confidence Intervals (ACI): Confidence intervals were calculated at a significance level of 95% using the observed information matrix.

5. Average interval length (AIL): Measures how narrow and useful the confidence interval is (the shorter the interval, the more accurate the estimate):

AIL= 

6. Probability of Convergence (CP): Measures the proportion of times (out of 1000 iterations) that the confidence interval successfully contained the true value of the parameter θ.

where : is actual value of the parameter, : estimated value of the parameter, and i: is iteration number.  

Based on this methodology, estimation methods were tested to determine the optimal and most appropriate option for modeling real data, and the numerical results and trends will be reviewed in detail in the following tables.

Our modeling findings are sensibly reported in Tables 4–6. Additionally, we have skillfully visualized the numerical data in these tables. Lower positions indicate more strength and dominance in this situation. These classifications measure the performance of each estimation method. The Table provide a comprehensive overview of the effectiveness and strengths of the estimation methods associated with the various approaches.

Table.4. Simulation results for the OBXII-IRD parameter estimators at λ= 0.7, ρ=3.2, β=1.1.

N–E. P.MLEMPSEKE–MLEMPSEKE
30Mean1.6228321.6700510.874213ACI0.201573
5.997310
0.110347
6.593772
0.094540
1.653886
3.6181423.9794773.3325890.342760
8.576670
0.014837
10.70739
1.536750
5.128428
1.7902212.1078401.1424150.290987
4.814944
0.104590
5.623921
0.843442
1.441388
RMSE2.7810513.0726860.460859AIL5.9973106.5937721.559345
2.6783733.5776650.9764398.57667010.7073933.591677
1.8910072.3513470.1663414.1849445.6239210.597945
Bias0.9228320.9700510.174213CP90.9091.7095.90
0.4181420.7794770.13258996.1096.3095.40
0.6902211.0078400.0424159289.5095.10
60Mean1.0345051.0469230.803711ACI0.257890
3.734370
0.270166
3.948830
0.250642
1.356780
3.4665173.5878033.2924580.595910
6.337125
0.065259
7.110348
1.931589
4.653327
1.3224671.4757251.1277530.390142
3.067305
0.230197
3.821698
0.877966
1.377539
RMSE1.6618491.6276030.338568AIL3.7343703.9488301.106138
1.6834501.8960670.7582105.7412157.0450892.721737
1.0586831.4084520.1323933.0673053.8216980.499572
Bias0.3345050.3469230.103711CP97.4097.4094.40
0.2665170.3878030.09245895.5097.4093.80
0.2224670.3757250.02775396.2093.9095.10
120Mean0.8296460.8687090.761395ACI0.280124
2.459609
0.301275
2.663477
0.323185
1.199605
3.2689953.3359763.2504791.963295
4.574695
1.939246
4.732705
2.092412
4.408546
1.1531481.1841951.1193030.416327
1.889969
0.294774
2.073615
0.953626
1.284980
RMSE0.8655120.9682700.247901AIL2.4596092.6634770.876419
0.7322880.8176220.5995002.6113992.7934582.316133
0.3876010.5435030.1016191.4736421.7788400.331353
Bias0.1296460.1687090.061395CP99.4098.7094.60
0.0689950.1359760.05047993.2092.6094.50
0.0531480.0841950.01930399.4099.1091.10
200Mean0.7459580.7617700.742843ACI0.302968
1.188948
0.323754
1.199786
0.342542
1.143143
3.2304263.2706643.2300802.393511
4.067340
2.378038
4.163291
2.324641
4.135519
1.1179211.1250031.1136850.946407
1.289435
0.930309
1.319697
0.953471
1.273899
RMSE0.2339470.2593180.210393AIL0.8859790.8760310.800601
0.5013490.5332720.4910251.6738291.7852531.810878
0.0984190.1061750.0856120.3430270.3893870.320427
Bias0.0459580.0617700.042843CP95.4094.4094.80
0.0304260.0706640.03008091.5091.7093.30
 0.0179210.0250030.0136859494.6094
300Mean 0.7215510.7311620.734181ACI0.416339
1.026762
0.422068
1.040255
0.422579
1.045784
 3.2371583.2660843.2271922.546719
3.927596
2.456371
4.075798
2.409439
4.044945
 1.1078731.1124591.1106350.973505
1.242241
0.969974
1.254944
0.961063
1.260207
RMSE 0.1638440.1781280.188929AIL0.6104220.6181870.623205
 0.4068830.4299470.4359411.3808771.6194271.635505
 0.0754970.0770470.0783930.2687360.2849690.299144
Bias 0.0215510.0311620.034181CP94.1092.7093.40
0.0271580.0660840.0271929194.2094.30
0.0078730.0124590.0106359393.9095.10

 

Table 5. Simulation results for the OBXII-IRD parameter estimators at λ= 0.7, ρ=3, β=0.8.

N–Est. Par.MLEMPSEKE–MLEMPSEKE
30Mean1.7700231.9850290.879650ACI

0.083715

6.393954

0.011204

7.455226

0.116109

1.643192

 3.4060143.8322283.109433

0.290347

8.174990

0.702398

11.80781

1.351819

4.867046

 1.2900371.5821520.834830

0.184204

3.359213

0.021457

4.181349

0.604868

1.064792

RMSE 2.9723873.4847610.454819AIL6.3939547.4552261.527083
 2.5765904.2224740.9531068.17499011.807813.515226
 1.3002421.7544220.1235253.3592134.1813490.459923
Bias 1.0700231.2850290.179650CP90.2088.9095.20
 0.4060140.8322280.10943396.209795.10
 0.4900370.7821520.03483092.1089.8095.30
60Mean 1.0521901.2102510.808010ACI

0.170390

3.669367

0.063788

4.955028

0.239749

1.376271

 3.2445483.3509943.085509

0.488223

6.000874

1.012487

6.748188

1.730357

4.440661

 0.9855031.1189700.823305

0.279934

2.328251

0.060237

2.974303

0.624842

1.021767

RMSE 1.6168052.1145710.348314AIL3.6693674.9550281.136522
 1.6146711.8240290.7542535.5126506.7481882.710304
 0.8178391.1199060.1054602.3282512.9743030.396924
Bias 0.3521900.5102510.108010CP96.7095.6095.10
 0.2445480.3509940.08550996.7097.4093.70
 0.1855030.3189700.02330595.609395
120Mean 0.8221770.8786480.767653ACI

0.240321

2.156483

0.098732

2.620398

0.339574

1.195732

 3.0634603.1209983.022349

1.813244

4.313676

1.765536

4.476459

1.987373

4.057325

 0.8393400.8761770.817348

0.362339

1.316344

0.132213

1.620140

0.688134

0.946562

RMSE0.7110990.9423810.244185AIL2.1564832.6203980.886157
0.7009310.7917160.5343442.5004232.7109222.069952
 0.2451640.4555440.0797220.9540041.4879270.298427
Bias0.1221770.1786480.067536CP9998.8093.60
 0.0634600.1209980.02234993.1092.5094.20
 0.0393410.0761770.01734899.1098.8091.60
200Mean 0.7484940.7722090.748512ACI

0.288772

1.208216

0.119433

1.424986

0.313089

1.183936

 3.0274103.0630393.026606

2.222836

3.831984

2.198354

3.927724

2.135948

3.917264

 0.8149260.8225050.812455

0.676776

0.953076

0.617802

1.027207

0.677128

0.947781

RMSE 0.2429410.3822210.229250AIL0.9194441.3055520.870847
 0.4818670.5158900.4828361.6091471.7293691.781316
 0.0793650.1108030.0724610.2763010.4094050.270653
Bias 0.0484940.0722090.048512CP95.8097.2095.40
 0.0274100.0630390.02660691.4091.7093.30
 0.0149260.0225050.01245594.1097.3095
300Mean 0.7183250.7284460.735705ACI

0.410950

1.025700

0.389956

1.066936

0.426551

1.044810

 3.0527003.0764893.016812

2.275297

3.825242

2.311706

3.841272

2.276289

3.757335

 0.8048280.8086250.808429

0.706218

0.903438

0.711219

0.906031

0.704036

0.91282

RMSE 0.1652970.1785200.177168AIL0.6147500.6769800.61830
 0.4053350.4283020.4008031.5499441.5295651.481045
 0.0556010.0590300.0573540.1972200.1948110.208786
Bias 0.0183250.0284460.035750CP94.7095.7092.60
 0.0502700.0764890.01681293.9092.2092.80
 0.0048280.0086250.0084299392.4093.30

 

Table 6. Simulation results for the OBXII-IRD parameter estimators at λ= 0.8, ρ=3.5, β=1.2.

N–E. P.MLEMPSEKE–MLEMPSEKE
30Mean 1.7745071.4880120.992217ACI

0.096783

6.804978

0.037489

5.886295

0.140345

1.844090

3.8705764.3140093.668577

0.991023

8.932019

0.078649

12.29458

1.755870

5.581285

 1.9437622.2989311.236282

0.192674

5.138058

0.078136

6.052449

0.946624

1.525941

RMSE 3.1703282.6936620.504147AIL6.8049785.8862951.703744
 2.7257544.2214731.0440508.93201912.294583.825415
 2.0024882.5200890.1599965.1380586.0524490.579316
Bias0.9745070.6880120.192217CP92.8095.6093.40
 0.3705760.8140090.16857796.3096.8094.60
 0.7437621.0989310.03628291.8089.3093.10
60Mean 1.0753681.0945460.918114ACI

0.101289

3.488727

0.060127

3.716324

0.257503

1.578725

 3.7227283.8482823.608937

1.269901

6.175555

0.294581

7.401983

2.123780

5.094093

 1.4537701.6373791.225184

0.360297

3.414366

0.172037

4.307581

0.975343

1.475025

RMSE 1.4809081.4665630.402672AIL3.4887273.7163241.321221
 1.4376621.9045170.8284764.9056537.1074012.970313
 1.1903971.6057320.1319063.4143664.3075810.499681
Bias 0.2753680.2945460.118111CP97.909894.20
 0.2227280.3482820.10893792.2097.5093.80
 0.2537700.4373790.02518496.2093.5094.40
120Mean0.9096200.9353910.879270ACI

0.240256

2.030847

0.096246

1.774536

0.369911

1.376030

3.5622733.6235133.556115

2.202556

4.921989

2.175559

5.071467

2.319159

4.793072

 1.2572991.2897671.219103

0.435934

2.078663

0.352258

2.227275

1.051502

1.386704

RMSE 0.5987740.4658950.285214AIL2.0308471.6782901.006118
 0.7617400.8448810.6405282.7194322.8959072.473913
 0.4318120.5730480.1027191.6427281.8750170.335201
Bias0.1096200.1353910.072970CP98.2095.1093.90
 0.0622730.1235130.05611593.5092.9094.10
 0.0572990.0897670.01910399.5099.2091.30
200Mean 0.8597970.8841020.853816ACI

0.323385

1.396209

0.348237

1.419968

0.392118

1.315515

 3.2525443.5608583.520731

2.653212

4.397277

2.633077

4.488639

2.593466

4.447996

 1.2189631.2273731.214774

1.043974

1.393952

1.028040

1.426706

1.052755

1.375799

RMSE 0.2841270.3193880.243597AIL1.0728231.0717300.923396
 0.5220340.5527470.5023451.7440651.8555611.854298
 0.1005390.1091370.0863890.3499780.3986660.323044
Bias 0.0597970.0841020.053816CP95.7094.8094.40
 0.03252440.0608580.04073191.4091.9092.50
 0.0189630.0273730.01427794.109593.50
300Mean 0.8279010.8432250.832136ACI

0.470419

01.185382

0.476654

1.209796

0.487849

1.176423

 3.5340593.5589833.536905

2.818367

4.249750

2.721176

4.396790

2.718194

4.355417

 1.2082781.2139081.207418

1.073331

1.343225

1.069948

1.357867

1.063687

1.351149

RMSE 0.1922710.2124380.207802AIL0.7149630.7331420.688574
 0.4213800.4435170.4371011.4313821.6756141.637222
 0.0728500.0780680.0750030.2698930.2879190.287461
Bias0.0279010.0432250.032136CP94.4093.3092.20
 0.0340590.0589830.03680591.3094.2094.20
 0.0082780.0139080.00741893.1093.6094.90

Tables 4–6 provide a summary of the simulation results illustrating capability performance. The main conclusions can be summarized as follows:

7. All estimation converge toward the true parameter values as the value of n increases, confirming their consistency with respect to the parameters.

8. A decrease in bias value was observed for all methods as the sample size increased, indicating improved accuracy in estimating the parameters of the OBXII-IRD.

9. The results Table 7 and Fig. (6) demonstrated all methods achieved satisfactory performance; however, the best performing method was the KE method, which exhibited high accuracy and minimized the gap between theoretical and empirical distribution making it an effective alternative to MLE method.

Table 7. Results of the descriptive analysis of the data used.

DatasetNMeanSDMedian madMin Max Rang CSCKSE
Data1001.660.61.540.40.925.314.393.1314.080.06

Note: mad: Median Absolute Deviation; CS: Coefficient of Skewness; CK: Coefficient of Kurtosis; SE: Standard Error of the mean

Fig. (6). Display data and its representation.

6. APPLICATION

To demonstrate the applicability and high suitability of the OBXII-IRD, we compared it with three competing distributions (such as, IR, HLIR [34], and OFIR [35]). based on real carbon fiber data. Table 8 presents the maximum likelihood estimates (MLEs) for the model parameters, along with the quality of fit test statistics Cramér-von Mises )W(, Anderson-Darling (A) and Kolmogorov-Smirnov (K) with the probability value (P-value). Table 9 summarizes the integrated statistical trade-off criteria, which include the negative probability function value (-LL), Akaike Information Criteria, Consistent AIC, Bayesian Information Criteria, Hanan and Quinn.

Table 8. MLEs and differentiation criteria for various models.

ModelMLEsWAKSP-value
OBXII-IR

:8.981019

:0.138496

:2.544796

0.0556220.4273780.0598030.866877
IR:2.2021590.0861960.7053850.325146 
HLIR

:4.827274

:3.921237

0.1398561.1019670.0805990.534419
OFIR

:1.296087

:1.551353

0.1965411.3026360.1171630.128402

Table 9. Statistical preference criteria and quality of fit for the data set.

Model-LLAICCAICBICHQIC
OBXII-IR51.83538109.6708109.9208117.4863112.8338
IR89.86446181.7289181.7697184.3341182.7833
HLIR57.17079118.3416118.4653123.5519120.4503
OFIR57.37916118.7583118.8820123.9680120.8670

Based on the results in Tables 8 and 9, the superiority of the OBXII-IRD over the comparative distribution can be explained by the following reasons:

  • Table 8 shows that the OBXII-IRD recorded the smallest distance between the empirical and theoretical distributions in the KS specifically 0.059803, which is the lowest value among the competing distributions. This superiority is directly reflected in the P-value 0.866877, which exceeded the standard significance level (0.05), indicating a failure to reject the null hypothesis. Consequently, the proposed distribution offers a better fit for the data and outperforms the comparative distributions.
  • The OBXII-IRD also yields the lowest statistical values W= 0.055622 and A= 0.427378. this superiority minimizes the squared differences between the empirical and cdf thereby enhancing the model’s reliability at the extremes of the data.
  • Table 9 demonstrates the efficiency of the OBXII-IRd based on information criteria, it achieved the lowest values for (-LL=51.8353, AIC=109.67, BIC=117.48, and HQIC=112.83). consequently, the OBXII-IRD is considered the most suitable fit for the data compared to the other distributions evaluated, as the distribution yielding the lowest values for these criteria is deemed superior.

The study included breaking stress data for 100 carbon fibers, obtained from reference [36]

(0.92, 0.928, 0.997, 0.9971, 1.061, 1.117,1.162, 1.183, 1.187, 1.192, 1.196, 1.213, 1.215, 1.2199, 1.22, 1.224, 1.225, 1.228, 1.237, 1.24, 1.244, 1.259,1.261, 1.263, 1.276, 1.31, 1.321, 1.329, 1.331, 1.337,1.351, 1.359, 1.388, 1.408, 1.449, 1.4497, 1.45, 1.459, 1.471, 1.475, 1.477, 1.48, 1.489, 1.501, 1.507, 1.515,1.53, 1.5304, 1.533, 1.544, 1.5443, 1.552, 1.556, 1.562, 1.566, 1.585, 1.586, 1.599, 1.602, 1.614, 1.616, 1.617, 1.628, 1.684, 1.711, 1.718, 1.733, 1.738, 1.743, 1.759, 1.777, 1.794, 1.799, 1.806, 1.814, 1.816, 1.828, 1.83, 1.884, 1.892, 1.944, 1.972, 1.984, 1.987, 2.02, 2.0304, 2.029, 2.035, 2.037, 2.043, 2.046, 2.059, 2.111, 2.165, 2.686, 2.778, 2.972, 3.504, 3.863, 5.306).

To enhance the numerical results listed in Tables 8 and 9, the model fit quality was graphically represented across (Fig. 7–9) using real experimental data for carbon fibers, and the geometric reading of these drawings shows the following:

Fig. (7). Graphical comparison of pdf for different models with the experimental distribution of data.

Fig. (8). Graphical comparison of the cdf for different models.

Fig. (9). Shows the P–P Plot for the OBXII-IRD and the comparative models, to show how well the proposed distribution matches the data.

Fig. (7) shows the pdf plotted for the various distributions overlaid on the bar graph (histogram) of the actual experimental data distribution. It is observed that the OBXII-IRD curve (represented by its corresponding color) is the only distribution that correctly tracks the rise and fall of the bar graphs, covering the distribution peak with exceptional accuracy and providing a flexible fit to the rightward skew of the data.

Fig. (8) depicts the cdf for the OBXII-IRD and the comparative distributions. A high degree of geometric agreement is noticed between the OBXII-IRD curve and the data curve. This agreement decressed progressively for the IR, OFIR, and HLIR distributions, which explains why the proposed distribution achieved the lowest KS statistic value.

Fig. (9) presents the p-p plot considered one of the accurate tools for assessing goodness of fit, which illustrates the relationship between the cumulative probabilities and its theoretical counterpart. The figure clearly demonstrates the probability points for the OBXII-IRD constant matching with the reference line. This shows a lack of dispersion or deviation, therefore presenting that the OBXII-IRD is more efficient than the other distributions.

CONCLUSION

This study provided an in-depth analysis of the efficiency of different estimation methods for the OBXII-IRD, where the results showed the consistency of all estimation methods, with the estimated values converging from the true values ​​of the parameters as the sample size increased, while the probability of coverage (CP) remained stable. As for the comparison between non-Bayesian methods, the KE method has consistently proven its superiority, followed by the MLE method.

Regarding entropy measures, seven distinct entropy measures were derived and analyzed to assess the degree of uncertainty associated with the OBXII-IRD. The Reny entropy measure proved to be the most effective at representing the behavior of the OBXII-IRD compared to the other measures. The analysis revealed that increasing the values of the parameters  resulted in lower entropy values, indicating reduced uncertainty. Conversely, increasing the value of the parameter  led to a higher level of dispersion. A practical application using carbon fiber data analyzed via the MLE demonstrated that the OBXII-IRD provided a better fit than IR, HLIR, and OFIR distributions, a conclusion supported by a p-value of 0.8668.

LIST OF ABBREVIATIONS

ACI

=

Approximate Confidence Interval

KE

=

Kolmogorov Estimation

MLE

=

Maximum Likelihood Estimation

MPS

=

Maximum Product space Estimation

OBXII-IRD

=

Odd Burr Inverse Rayleigh

AUTHORS’ CONTRIBUTION

A.A.K. contributed to the study concept and design and writing of the paper. M.A.K. contributed to the study concept and design and writing of the paper. A.M.S. contributed to data collection and writing of the paper. A.T.M. contributed to data collection, review, and editing. Q.N.H. contributed to data analysis or interpretation, review, and editing. M.K.A.H. contributed to data analysis or interpretation, review, and editing. All authors reviewed and approved the final version of the manuscript.

ETHICAL APPROVAL & INFORMED CONSENT

The following study does not involve any human participants, clinical trials, or animal tests. Moreover, ethical approval and informed consent were not required for this research. Section 6 represents previously published physical measurements of carbon fiber materials science tensile strength, which are publicly available.

AVAILABILITY OF DATA AND MATERIALS

The data used in this study are publicly available and were obtained from previously published sources. The dataset analyzed in this work is cited appropriately in the manuscript. The simulation results and computational procedures can be reproduced using the methods and details provided in the article.

FUNDING

This research received no financial support from any public, commercial, or not-for-profit funding agency. 

CONFLICT OF INTEREST

The authors declare that they have no competing interests that may influence this study.

ACKNOWLEDGEMENTS

None.

DECLARATION OF AI

The authors declare that this research project did not use any artificial intelligence tools or technologies in the process of preparing, analyzing, or writing the study and the authors take full responsibility for the content and originality of this work.

REFERENCES

[1]Treyer VN. Inverse Rayleigh (IR) model. Proc. USSR Acad. Sci. 1964.

[2] Voda VG. On the inverse Rayleigh distributed random variable. Rep Statis App Res JUSE. 1972;19(4):13–21.

[3] Gharraph MK. Comparison of estimators of location measures of an inverse Rayleigh distribution. Egypt Stat J. 1993; 37(2): 295–309.
https://doi.org/10.21608/esju.1993.426926

[4] Howlader HA, Hossain AM, Makhnin O. Bayesian prediction bounds for Rayleigh and inverse Rayleigh lifetime models. J Appl Stat Sci. 2009; 17(1).
https://doi.org/10.47013/19.2.3

[5] Banerjee P, Bhunia S. Exponential Transformed Inverse Rayleigh distribution: Statistical properties and different methods of estimation. Austrian J Stat. 2022; 51(4): 60–75.
http://dx.doi.org/10.17713/ajs.v51i4.1338

[6] Khan MS, King R. Transmuted modified Inverse Rayleigh distribution. Austrian J Stat. 2015; 44(3): 17–29.
https://doi.org/10.17713/ajs.v44i3.21

[7] Goual H, Yousof HM. Validation of Burr XII inverse Rayleigh model via a modified chi-squared goodness-of-fit test. J Appl Stat. 2020; 47(3): 393–423.
https://doi.org/10.1080/02664763.2019.1639642

[8] Cordeiro GM, Yousof HM, Ramires TG, Ortega EMM. The Burr XII System of densities: properties, regression model and applications. J Stat Comput Simul. 2018;88(3):432–456.
https://doi.org/10.1080/00949655.2017.1392524

[9] Ali S. Mixture of the inverse Rayleigh distribution: Properties and estimation in a Bayesian framework. Appl Math Model. 2015; 39(2): 515–530.
http://dx.doi.org/10.1016/j.apm.2014.05.039

[10] Fisher RA. The effect of methods of ascertainment upon the estimation of frequencies. Ann Eugen. 1934; 6(1): 13–25.
https://doi.org/10.1111/j.1469-1809.1934.tb02105

[11] Rao CR. On discrete distributions arising out of methods of ascertainment. Sankhyā: Indian J Stat Ser A. 1965; 27: 311–324.
https://www.jstor.org/stable/25049375

[12] Aydın D. The new weighted inverse Rayleigh distribution and its application. Facta Univ Ser Math Inform. 2019; 34(3): 511-523.
https://doi.org/10.22190/FUMI1903511A

[13] Rao GS, Mbwambo S. Exponentiated inverse Rayleigh distribution and an application to coating weights of iron sheets data. J Probab Stat. 2019; 2019: 1–13.
https://doi.org/10.1155/2019/7519429

[14] Kamnge JS, Chacko M. Half logistic exponentiated inverse Rayleigh distribution: Properties and application to life time data. PLoS One. 2025; 20(1): e0310681.
https://doi.org/10.1371/journal.pone.0310681

[15] Yanuar F, Iqbal M, Devianto D, Zetra A, Asdi Y, Ilahi R, & Sani, RF. Bayesian estimation under different loss functions for the case of inverse Rayleigh distribution. Kuwait J Sci. 2025; 52(1): 100343.
https://doi.org/10.1016/j.kjs.2024.100343

[16] Asif M, Raftab M, Usman M, Khan AA. Analysis of parameters of the exponentiated inverse Rayleigh distribution under the Bayesian framework. Kuwait J Sci. 2025; 52(3): 100424.
https://doi.org/10.1016/j.kjs.2025.100424

[17] Alghamdi AS. Marshall-Olkin Exponentiated Inverse Rayleigh Distribution Using Bayesian and Non-Bayesian Estimation Methods. Symmetry. 2025; 17(5): 707.
https://doi.org/10.3390/sym17050707

[18] Cover TM, Thomas JA. Elements of Information Theory. 2nd ed. New York: John Wiley & Sons; 2006.
https://doi.org/10.1002/047174882X

[19] Shannon CE. A mathematical theory of communication. Bell Syst Tech J. 1948; 27(3): 379–423.
https://doi.org/10.1002/j.1538-7305.1948.tb01338.x

[20] Rényi A. On measures of entropy and information. In: Proceedings of the Fourth Berkeley Symposium on Mathematical Statistics and Probability. 1960; 1: 547–61. Available From: https://www.scirp.org/reference/referencespapers?referenceid=2095576

[21] Tsallis C. Possible generalization of Boltzmann-Gibbs statistics. J Stat Phys. 1988; 52(1–2): 479–87.
http://dx.doi.org/10.1007/bf01016429

[22] Arimoto S. Information-theoretical considerations on estimation problems. Inf Contr. 1971; 19(3): 181–94.
http://dx.doi.org/10.1016/s0019-9958(71)90065-9

[23] Havrda J, Charvát F. Quantification method of classification processes. Concept of structural a-entropy. Kybernetika. 1967; 3(1): 30–5.
http://dml.cz/dmlcz/125526

[24] Gupta HC, Sharma BD. On non-additive measures of inaccuracy. Czechoslov Math J. 1976; 26(4): 584–95.
http://dx.doi.org/10.21136/CMJ.1976.101429

[25] Sant’Anna AP, Taneja IJ. Trigonometric entropies, Jensen difference divergence measures, and error bounds. Inf Sci. 1985; 35(2): 145–56.
https://doi.org/10.1016/0020-0255(85)90046-5

[26] Khalaf AA, Khaleel MA, Jawa TM, Sayed-Ahmed N, Tolba AH. A Novel Extension of the Inverse Rayleigh Distribution: Theory, Simulation, and Real-World Application. Appl Math Inf Sci. 2025; 19(2): 467–488.
http://dx.doi.org/10.18576/amis/190220

[27] Khalaf AA, El-Saeed AR, Khaleel MA, Tolba AH. The four-parameter odd generalized Rayleigh Lomax distribution: Theory, simulation, and applications. Symmetry. 2026; 18(2): 244.
http://dx.doi.org/10.3390/sym18020244

[28] Abdulqader Salih S, Mahmood Ali D, Khaled Kolaib H, Abdulrahman Khalaf A. Truncated Exponentiated Ailamujia Exponential distribution: Properties and applications. Baghdad Sci J. 2026; 23(1): 271–90. Available from: http://dx.doi.org/10.21123/2411-7986.5183

[29] Khalaf AA, Haleel MA. The new strange Generalized Rayleigh family: Characteristics and applications to COVID-19 data. Iraqi J Comput Sci Math. 2024; 5(3).
http://dx.doi.org/10.52866/ijcsm.2024.05.03.005

[30] Khalaf AA, Khaleel MA. The Odd Burr XII Exponential distribution: Properties and applications. In: AIP Conference Proceedings. 2025; 3264(1).
https://doi.org/10.1063/5.0258451

[31] Khalaf AA, Khaleel MA, Tolba AH, Ahmed NS. Classical Inference for the Five Parameter Exponentiated Weibull Distribution: Properties and Applications in Health and Reliability. 2025; 14: 469–492.
http://dx.doi.org/10.18576/jsap/140310

[32] Gemeay A M, Sapkota L P, Tashkandy Y A, Bakr M E, Balogun OS, Hussam E. New bounded probability model: Properties, estimation, and applications. Heliyon. 2024; 10 (23).‏
https://doi.org/10.1016/j.heliyon.2024.e38965

[33] Leão J, Saulo H, Bourguignon M, Cintra RJ, Rêgo LC, Cordeiro GM. On some properties of the beta inverse Rayleigh distribution. arXiv. 2022.
http://arxiv.org/abs/2206.01229

[34] Almarashi AM, Badr MM, Elgarhy M, Jamal F, Chesneau C. Statistical inference of the half-logistic inverse Rayleigh distribution. Entropy. 2020; 22(4): 449.
http://dx.doi.org/10.3390/e22040449

[35] Elgarhy M, Alrajhi S. The odd Fréchet inverse Rayleigh distribution: statistical properties and applications. J Nonlinear Sci Appl. 2018; 12(5): 291–299.
http://dx.doi.org/10.22436/jnsa.012.05.03

[36] Nichols MD, Padgett WJ. A bootstrap control chart for Weibull percentiles. Quality and reliability engineering international. 2006; 22: 141–51.
https://doi.org/10.1002/qre.691

Insert math as
Block
Inline
Additional settings
Formula color
Text color
#333333
Type math using LaTeX
Preview
\({}\)
Nothing to preview
Insert