IHJPAS. 36 (3) 2023 437 This work is licensed under a Creative Commons Attribution 4.0 International License *Corresponding Author: kehindebashiru@uniosun.edu.ng Abstract A new family of distribution named Double-Exponential-X family is proposed. The proposed family is generated from the double exponential distribution. The forms of the probability densities and hazard functions of two distinct subfamilies of the proposed family are examined and reported. General properties such as moment, survival, order statistics, probability weighted moments and quartile functions of the models are investigated. A sub family of the developed family of double –Exponential-X family of the distribution known as double-Exponential-Pareto distribution was used to fit a real life data on the use of antiretroviral drugs. Montecarlo simulation of Efficacy of antiretroviral drugs is conducted to evaluate the performance of the model. The models were tested using some models diagnostic tests and it was revealed that the proposed model was better than the ones proposed before it from the same family and also, the stochastic dominance method was used to affirm the best antiretroviral drugs used in the study. Keywords: Double Exponential Distribution, Hazard Function, Molecular Simulation, Stochastic Dominance. doi.org/10.30526/36.3.3377 Article history: Received 25 January 2023, Accepted 14 March 2023, Published in July 2023. Ibn Al-Haitham Journal for Pure and Applied Sciences Journal homepage: jih.uobaghdad.edu.iq Double-Exponential-X Family of Distributions: Properties and Applications Kehinde Adekunle Bashiru* Department of Statistics, Osun State University, Osogbo, Nigeria Taiwo Adetola Ojurongbe Department of Statistics, Osun State University, Osogbo, Nigeria Lawal Sola Department of Statistics, Osun State College of Technology, Esa Oke, Nigeria Nureni .Olawale Adeboye Department of Statistics, Osun State University, Osogbo, Nigeria Habeeb Abiodun Afolabi Department of Statistics, Osun State University, Osogbo, Nigeria https://creativecommons.org/licenses/by/4.0/ mailto:%20kehindebashiru@uniosun.edu.ng mailto:kehindebashiru@uniosun.edu.ng mailto:taiwo.ojurongbe@uniosun.edu.ng mailto:slawal@gmail.com mailto:nureni.adeboye@uniosun.edu.ng mailto:habeeb.afolabi@uniosun.edu.ng IHJPAS. 36 (3) 2023 438 1. Introduction In the past few decades, a variety of classical or general distributions have been widely used to structure data in the fields of engineering, actuarial, environmental science, biology, demography, economy, insurance, and finance. Despite this, there continue to be problems in each of the aforementioned fields. Naturally, this may be attributed to the fact that, as a result of the growth of science and technology, the underlying distributions are no longer able to suit the data set. Because of this, it is necessary to create new kinds of extended distributions that can handle these fresh difficulties. As a result, numerous methods have been proposed and researched to develop or introduce new families of distribution by including a parameter or several parameters in the baseline random variable model while providing a great degree of flexibility in actual data modeling, including beta-G [1,2], Gamma-G type-1 [3], Mc-G [4], Log-Gamma-G Type-2 [5], Gamma-G type-3 [6], Weibull-X family of distributions [7], odd generalized exponential-G [8,9], type-1 half-logistic family [10], Kumaraswamy-Weibull-generated family [11], new Weibull- G family [12], generalized transmuted- G [13] Taking 𝑣(𝑑) to be the probability density function (PDF) of a random variable, say, which belongs to m and n, that is T∈ [π‘š, 𝑛] for -∞ ≀ m < n <∞ and also, let π‘Š[𝐹(π‘₯, πœ‰)] be a function of cumulative distribution function (CDF) of a random variable, say X which depends on the parameter πœ‰ which satisfies the conditions: i π‘Š[𝐹(π‘₯, πœ‰)] is bounded between m and n ii. π‘Š[𝐹(π‘₯, πœ‰)] is differentiable and monotonically increasing iii. π‘Š[𝐹(π‘₯, πœ‰)] β†’m as xβ†’-∞ and π‘Š[𝐹(π‘₯, πœ‰)]β†’n as xβ†’βˆž The CDF of the T-X family of distribution as recently used by [14] as 𝐺(π‘₯) = ∫ 𝑣(𝑑)𝑑𝑑 𝑀(𝐹(π‘₯;πœ‰)) π‘š , π‘₯ ∈ 𝑅 (1) And the PDF corresponding to (1) is given as: 𝑔(π‘₯) = { 𝑑 𝑑π‘₯ 𝑀[𝐹(π‘₯; πœ‰)]} 𝑉{𝑀[𝐹(π‘₯: πœ‰)]} π‘₯ ∈ 𝑅 Several new families of distributions have been postulated using the concept of T-X family as proposed [14]. Different functions of π‘Š[𝐹(𝑋; ∈)] have been employed in the literature for the T- X family which are revealed in Table1. This research aims to introduce a new family of continuous distributions, which is named as Double Exponential-X (β€œDE-X” for short) family. This family's two distinct sub-models are looked into and discussed. Both real life and simulated data are used to illustrate the essence of the new model in the area of applied science. This paper's remaining sections are divided as follows: In the section titled "Double-Exponential X family," we define the DE-X family of distributions, provide a mathematical derivation of some of the family's fundamental statistical properties, including the moment, survival function, order statistics, probability weighted moment (PWM), and quartile function, and provide a maximum likelihood estimation method for determining the model parameters. The DE-Pareto Distribution sub-model constructed from this family was fitted to real- life data of viral load (Copies/ml) of some patients receiving antiretroviral medication. In order to calculate the average binding affinity IHJPAS. 36 (3) 2023 439 of Abacavir and Neviranpine as they docked (using molecular docking), simulated data (from simulation of molecular dynamics) was also used. Table 1: A few T-X Family Members Currently in Existence 𝑾[𝑭(𝒙; 𝝃)] Range of 𝑿 Members of T – X family 𝐹(π‘₯, πœ‰) (0,1) Beta – G [1] βˆ’ log[𝐹(π‘₯; πœ‰)] (0, ∞) Gamma-GType-22) [15] βˆ’ log[1 βˆ’ 𝐹(π‘₯; πœ‰)] (0, ∞) Gamma-GType-1[3] 𝐹(π‘₯, πœ‰) 1 βˆ’ 𝐹(π‘₯, πœ‰) (0, ∞) Gamma-GType-3 [6] βˆ’ log[1 βˆ’ 𝐹(π‘₯, πœ‰)] (0, ∞) Exponentiated T – X [14] log { 𝐹(π‘₯, πœ‰) 1 βˆ’ 𝐹(π‘₯, πœ‰) } (βˆ’βˆž, ∞) Logistics – G [16] log[βˆ’ log{1 βˆ’ 𝐹(π‘₯, πœ‰)}] (0, ∞) The logistic – X family [8,9] 2. Materials and Methods 2.0 The Double Exponential-X Family Given the PDF of the double exponential distribution 𝑉(π‘₯; 𝑏; Β΅) = 1 2𝑏 𝑒 |π‘₯βˆ’Β΅| 𝑏 ( 21) (2) And the corresponding CDF with high reliability on the support is given by: 𝐹(π‘₯; Β΅; 𝑏) = 1 2 𝑒 βˆ’(Β΅βˆ’π‘₯) 𝑏 ifπ‘₯ ≀ Β΅ (3) 𝐹(π‘₯; Β΅; 𝑏) = 1 βˆ’ 1 2 𝑒 βˆ’(π‘₯βˆ’Β΅) 𝑏 ifπ‘₯ β‰₯ Β΅ (4) Then, if π‘₯ β‰₯ 0, then the CDF changes to 𝐹(π‘₯; 𝑏) = 1 βˆ’ 1 2 𝑒 βˆ’ π‘₯ 𝑏 (5) And its corresponding PDF is given as: 𝑉(π‘₯; 𝑏) = 1 2𝑏 𝑒 βˆ’π‘₯ 𝑏 (6) If π‘₯ follows (6) and setting π‘Š[𝐹(π‘₯, πœ‰)]= 1 1βˆ’πΉ(π‘₯,πœ‰) in (1), we define the CDF of the Double Exponential-X family by G(π‘₯, 𝑏, πœ‰) = 1 2 [1 βˆ’ 𝑒 βˆ’ 1 𝑏 (π‘Š[𝐹(π‘₯,πœ‰)]) ] = 1 2 βˆ’ 1 2 𝑒 βˆ’ 1 𝑏 (π‘Š[𝐹(π‘₯,πœ‰)]) , b, πœ‰ > 0 , π‘₯ ∈ 𝑅 (7) Where the baseline distribution's CDF, which is dependent on the vector parameter, is denoted by 𝐹(π‘₯, πœ‰), πœ‰ and π‘Š[𝐹(π‘₯, πœ‰)]= 1 1βˆ’πΉ(π‘₯,πœ‰) is a special transformation. The corresponding PDF is obtained via the first derivative of the CDF, G(π‘₯, 𝑏, πœ‰) and it is shown by: 𝑔(π‘₯, 𝑏, πœ‰) = 1 2𝑏 [ 𝑑 𝑑π‘₯ (π‘Š[𝐹(π‘₯, πœ‰)])] 𝑒 βˆ’{ 1 𝑏 (π‘Š[𝐹(π‘₯,πœ‰)])} , 𝑏, πœ‰ > 0 (8) 𝑔(π‘₯, 𝑏, πœ‰)= 𝑔(π‘₯, 𝑏, πœ‰) = 1 2𝑏 [(π‘Š[𝐹(π‘₯, πœ‰)])β€²]𝑒 βˆ’{ 1 𝑏 (π‘Š[𝐹(π‘₯,πœ‰)])} , where (π‘Š[𝐹(π‘₯, πœ‰)])β€² = 𝑑 𝑑π‘₯ (π‘Š[𝐹(π‘₯, πœ‰)]) and 𝑏, πœ‰ > 0 IHJPAS. 36 (3) 2023 440 The motivations behind using DE-X family in the real life situations include: ο‚· It is a popular technique for inserting extra parameter(s) to produce an expanded baseline distribution, ο‚· To enhance the properties of the traditional distributions, ο‚· To continuously offer better fits compared to other produced distributions with the same or more parameters, 2.1 Special Sub-Models Most often for one reason or another, extended versions of distributions are introduced for the following reasons; 1. An alternate model to the existing distribution that has been successfully applied in the past 2. A flexible statistical model that provides a good match to the data 3. A model with closed CDF, Survival Function (SF), Hazard Rate Function( HRF) forms and a strong propensity to lessen estimation challenges. Thus, we introduce two unique sub-models of the DE-X family that have at least one of the aforementioned characteristics. 2.1.1 The DE-Pareto (DE-P) distribution The PDF and the CDF of the Pareto random variable are f(π‘₯, ⍬, 𝛼) = 𝛼(𝛼⍬) π‘₯⍬+1 and F(π‘₯, 𝛼, πœƒ) = 1 βˆ’ ( 𝛼 π‘₯ ) ⍬ , π‘₯>0, 𝛼, πœƒ>0, . Then the CDF and the PDF of the DE-P model are given by (9) and (11) respectively. 𝐺(π‘₯, 𝑏, 𝛼, πœƒ) = 1 2 βˆ’ 1 2 𝑒 βˆ’ 1 𝑏 ( π‘₯ 𝛼 ) πœƒ , b, πœƒ, 𝛼>0 (9) Hence, 𝐺(π‘₯, 𝑏, 𝛼, πœƒ) = 1 2 (1 βˆ’ 𝑒 βˆ’ 1 𝑏 ( π‘₯ 𝛼 ) πœƒ )b, πœƒ, 𝛼 > 0(10) (10) Hence,𝑔(π‘₯, 𝑏, 𝛼, πœƒ)= ⍬ 2𝛼𝑏 ( π‘₯ 𝛼 ) β¬βˆ’1 𝑒 βˆ’ 1 𝑏 ( π‘₯ 𝛼 ) ⍬ , π‘₯, b,πœƒ, Ξ± > 0 (11) Plots of the DE-P density and Hazard rate function for selected parameter values are presented in Fig.1 and Fig.2 respectively below: Figure.1:PDF plot at fixed b =1 and increasing 𝛼 π‘Žπ‘›π‘‘ πœƒ IHJPAS. 36 (3) 2023 441 Figure.2:Hazard plot at fixed b =1 and increasing 𝛼 π‘Žπ‘›π‘‘ πœƒ 2.1.2 The DE-ParetoII (DE-P2) Given the PDF and CDF of the Pareto type-2 distribution𝑓(π‘₯, π‘Ž, 𝑏) = 𝑏 π‘Ž ( π‘Ž π‘₯βˆ’π‘Ž ) 𝑏+1 , b,a, π‘₯ > 0 and 𝐹(π‘₯; π‘Ž, 𝑏) = 1 βˆ’ ( π‘Ž π‘₯βˆ’π‘Ž ) 𝑏 , x, b,a,> 0. Then the CDF and PDF of DE-PII model are given as 𝐺(π‘₯, 𝑏, π‘Ž, πœ‰) = 1 2 [1 βˆ’ 𝑒 βˆ’ 1 𝑏 ( π‘₯βˆ’π‘Ž π‘Ž ) πœƒ ], a,b,πœƒ > 0; π‘₯ > 0 0 (12) For the PDF, 𝑔(π‘₯, a, b, πœ‰) = πœƒ 2π‘π‘Ž ( π‘₯βˆ’π‘Ž π‘Ž ) πœƒβˆ’1 𝑒 βˆ’ 1 𝑏 ( π‘₯βˆ’π‘Ž π‘Ž ) πœƒ , a, b, πœƒ > 0 ; π‘₯> 0 (13) Fig.3 and Fig.4 give the plots of the PDF and HRF of DE-PII distribution for selected parameter values. Figure.3 The pdf of double exponential-ParetoII distribution at different values of parameters IHJPAS. 36 (3) 2023 442 Figure.4 showing the hazard of double exponential-ParetoII distribution at different values of parameters 2.2 Basic Mathematical Properties of DE-Pareto Here, we provide the basic properties of DE-X family of distributions. 2.2.1 Moments of DE-Pareto Assuming X follows DE-X family with density 𝑔(π‘₯, a, b, πœ‰), then its π‘Ÿπ‘‘β„Ž moment is Β΅r 1=∫ xr ∞ βˆ’βˆž g(x, a, b, ΞΎ)dx E(Xr) = ∫ xr ∞ 0 ⍬ 2Ξ±b ( x Ξ± ) β¬βˆ’1 e βˆ’ 1 b ( x Ξ± ) ⍬ dx (14) E(Xr) = 1 2 (Ξ± √b ΞΈ ) r β”Œ ( r ΞΈ + 1) E(X)= Ξ± √b ΞΈ 2 β”Œ ( 1 ΞΈ + 1); r =1 (15) Variance: Var(X)= (Ξ± √b ΞΈ ) 2 [ 1 2 β”Œ ( 2 ΞΈ + 1) βˆ’ ( 1 2 β”Œ ( 1 ΞΈ + 1)) 2 ] (16) 2.2.2 Survival Function of DE-Pareto The survival function of DE-X family is given as 𝑆(π‘₯) = 1 βˆ’ 𝐺(π‘₯, 𝑏, 𝛼, πœƒ) (17) 𝑆(π‘₯) = 1 2 (1 + 𝑒 βˆ’ 1 𝑏 ( π‘₯ 𝛼 ) πœƒ ) (18) 2.2.3 Probability Weighted Moments of DE-Pareto The Probability Weighted Moment (PWM) of DE-X random variable is π‘€π‘Ÿπ‘ (π‘₯) = ∫ π‘₯ π‘ βˆž βˆ’βˆž 𝑔(π‘₯, 𝑏, π‘Ž, πœ‰)𝐺(π‘₯, 𝑏, π‘Ž, πœ‰))𝑠𝑑π‘₯ (19) 2.2.4 Hazard Function of DE-Pareto 𝐻(π‘₯; 𝛼, 𝑏, πœƒ) = 𝑔(π‘₯) 𝑆(π‘₯) = ⍬ 2𝛼𝑏 ( π‘₯ 𝛼 ) β¬βˆ’1 𝑒 βˆ’ 1 𝑏 ( π‘₯ 𝛼 ) ⍬ 1 2 (1+𝑒 βˆ’ 1 𝑏 ( π‘₯ 𝛼 ) πœƒ ) = ⍬ 𝛼𝑏 ( π‘₯ 𝛼 ) β¬βˆ’1 𝑒 βˆ’ 1 𝑏 ( π‘₯ 𝛼 ) ⍬ 1+𝑒 βˆ’ 1 𝑏 ( π‘₯ 𝛼 ) πœƒ (20) 2.2.5 Renyi Entropy of DE-Pareto π›Ώπ‘Ÿ = 1 1βˆ’π‘Ÿ π‘™π‘œπ‘” ∫ 𝑓 π‘Ÿ (π‘₯)𝑑π‘₯ ∞ 0 = 1 1βˆ’π‘Ÿ π‘™π‘œπ‘” ∫ ( ⍬ 2𝛼𝑏 ( π‘₯ 𝛼 ) β¬βˆ’1 𝑒 βˆ’ 1 𝑏 ( π‘₯ 𝛼 ) ⍬ ) 𝒓 𝑑π‘₯ ∞ 𝟎 (21) IHJPAS. 36 (3) 2023 443 π›Ώπ‘Ÿ = 1 1βˆ’π‘Ÿ π‘™π‘œπ‘” ( π‘πœƒ 2π‘Ÿπ›Όπ‘ √ π‘Ÿ 𝑏 πœƒ ) π‘Ÿ ( 𝛼 πœƒ √ 𝑏 π‘Ÿ πœƒ ) β”Œ (π‘Ÿ βˆ’ π‘Ÿ πœƒ + 1 πœƒ ) (22) 2.3. Estimation We describe the maximum likelihood method of estimation in this section for the unknown parameters of the DE-X family. The method of Maximum Likelihood estimation is adjudged to be the best method of estimating parameters of distributions. 2.3.1 Maximum Likelihood Estimation DE-Pareto In this, let 𝑋1, 𝑋2, … , π‘‹π‘˜ be the observed values for DE-X distribution with parameters𝑏, π‘Ž, πœ‰. The total log-likelihood function corresponding to (11) is presented as Log L(π‘₯, 𝑏, π‘Ž, πœ‰) = ∏ ⍬ 2𝛼𝑏 ( π‘₯ 𝛼 ) β¬βˆ’1 𝑒 βˆ’ 1 𝑏 ( π‘₯ 𝛼 ) ⍬ 𝑛 𝑖=1 (23) 3. Results 3.1 Fitting the Model to Molecular Dynamics Simulated Data The two approved HIV drugs were obtained from The Drug bank archive online, modeled using (Spartan’14 version 1.18 by wave function Inc) i. Geometry optimization and calculation of the two drugs were done using density functional theory (DFT) with 6-31G* basis set to obtain the molecular parameters responsible for the cytotoxicity of the compounds under investigation and then saved as pdb files for simulation study ii. The protein receptors responsible for HIV(that is, the HIV-causing organisms) used were gotten from protein data bank (PDB) [18] iii. with no missing residues and a single strand were obtained by deleting multiple ligands, non- protein and other protein parts from the pdb files. Autoduck tools version 1.5.6 was used to locate the binding site and to add hydrogen to the prepared protein receptor and finally converted to a protein data bank format (pdbqt) for the molecular simulation. Autodockvina were used for the simulation process and the simulation was repeated for seven different in order to obtain the most stable conformation of the drug-receptor complex and post docking analysis was done using EduPymolversion 1.7.4.4 and discovery studio 4.1 visualizer. Figure 5a: Fitting the Models to Neviranpine data (Obtained via Molecular Simulation) IHJPAS. 36 (3) 2023 444 From the figure 5a above, we can say that the selected distributions approximately follow the distribution of the data, thus, making the distributions suitable for data collected when using Neviranpine to dock HIV protein receptor Figure 5b: Fitting the Models to Abacavir data (Obtained via Molecular Simulation) From the Fig. 5b above, we can say that the selected distributions approximately follow the distribution of the data, thus, making the distributions suitable for data collected when using Abacavir to dock HIV protein receptor Table1. Fitting cells table with the including numbers Model Parameter Criteria Mean Affinity EPD οΏ½Μ‚οΏ½ = 0.029978(0.034) οΏ½Μ‚οΏ½ = 0.9980(0.367) οΏ½Μ‚οΏ½ = 0.996486(0.648) AIC=303.405 BIC=303.6433 CAIC=309.405 1.663157(0.902308) DEPD οΏ½Μ‚οΏ½=1.978e+06(1.639e- 02) οΏ½Μ‚οΏ½ = 1.503e- 04(5.307e-05) οΏ½Μ‚οΏ½=9.024e-01(2.211e- 04) AIC =-27203950 BIC= -27203950 CAIC=27203944 0.4511979(0.225) PD οΏ½Μ‚οΏ½= 8.089e+03 (7.567e-02) οΏ½Μ‚οΏ½ =(3.372e+03) 9.268e-02 AIC=-810083 BIC=-810082.8 CAIC=-810080.6 3.3724(0.41697) From table 2 above, we can say that the parameter of the model is stable and the Double exponential Pareto distribution (DEPD) appears to be the best of all the competing model (model with the lowest information criteria is the best of all the models). The 45.1198% of the binding affinity is attributed to Neviranpine. IHJPAS. 36 (3) 2023 445 Fig 4: Investigating Stochastic Dominance From fig 4 above, we can say that Neviranpine is first order (Stochastically) dominated Abacavir. That is, the mean effect of Nevirapine has dominated that of Abacavir when it comes to the viral load measurement (copies/ml) from the patients subjected to Anti-retroviral therapy. 3D Structures of Nevirapine 3D of Abacavir Figure 5: The Ligands (Nevirapine and Abacavir) used as the drug in Anti-retroveral therapy IHJPAS. 36 (3) 2023 446 3D Receptor Responsible for HIV (Protein responsible for HIV) Interaction of Abacavir and 3PLT Interaction of Nevirapine and 3LPT Figure 6. The Docking Result of Abacavir and Nevirapine The Abacavir and Nevirapine were used as the Ligands to dock HIV protein receptor, the fig 6 showed interaction between the ligands and HIV protein receptor. IHJPAS. 36 (3) 2023 447 Table 5: The Molecular Properties of the Drug used in this study (Molecular Docking) Descriptors Nevirapine Abacavir Molecular Formula C15H14N4O C14H18N6O Molecular Weight 266.30amu 286.34amu Area 281.25A2 305.46A2 Volume 269.04A3 285.66A3 Polarizability 62.14 63.34 Hydrogen Bond Donor (HBD) 1 2 Hydrogen Bond Acceptor (HBA) 5 5 Polar Surface Area (PSA) 36.84 67.47 Ovality 1.40 1.46 Log P 0.13 0.16 Dipole Movement 2.33 1.93D HOMO -5.80Ev -5.28eV LUMO -1.32eV -0.16eV BG 4.58eV 5.12eV 4. Discussion From Table5, the properties of the drugs used in determining the toxicity are clearly displayed, as behavior, and the general effect of each drug on the body system. Parameters such as HBA, HBD, PSA, LogP and polarizability are used to determine if a particular compound is a good drug candidate. This reveals the active properties of the drugs used in determining the cytotoxicity of each drug towards the receptor responsible for HIV. Failure of most drugs in clinical experiments is a result of there poor ADME (Absorption, Distribution, Metabolism and Excretion) properties, hence for a ligand to be chosen as a drug candidate with good ADME properties HBD must be less than five, HBA must be less than ten and its molecular weight in g/mol must be below 500, any ligand that violates this set of rules might result to the problem if ingested [19]. Polar Surface Area (PSA) gives information about the ability of a drug-like compound to penetrate into the cell membrane. It has been established that ligands with high PSA value greater than 140Γ…2 will have poor penetrating power into the cell membrane. On the other hand, molecules with PSA value below 60Γ…2 will have high penetrating power [20]. Hence, we are expecting high penetrating power into the cell membrane from Abacavir and Nevirapine, and these are used to determine if a particular compound is good drug candidate. The docking result used for stochastic dominance in order to determine the most effective drugs that perform actively in curing the diseases shows better binding affinity and inhibition constant from Nevirapine as compared to Abacavir and this is an indication that Nevirapine drug is more potent for curing the disease(HIV). 5. Conclusion A new family of distribution (Double Exponential-X family) has been introduced. Two distinct sub-models of the distribution were discussed. It is shown that the two special models' densities can generally be skewed to the left or right, symmetrical, reverse j-shaped, and can also have both growing and decreasing, as well as, of course, unimodal. The maximum likelihood estimation approach is used to estimate the parameters. A molecular simulation study that compares the newly IHJPAS. 36 (3) 2023 448 built model's performance to previous well-known distributions is presented. There are derived and explained several fundamental mathematical characteristics of the new family. Finally, real- life applications of the model to data sets on the effectiveness of various antiretroviral medications were provided, and they compare favorably to other well-known distributions in terms of data fit, and the efficacy of the pharmaceuticals utilized was compared using the stochastic dominance approach. References 1. Eugene, N.; Lee, C; Famoye, F. Beta-normal distribution and its applications. Communications in Statistic: Theory and methods, 2002; 31(4), 497-512. 2. Jones, M. C. Families of distributions arising from distributions of orderstatistics. Test, (2004) 13: 1 – 42. 3. Zografos, K. ; Balakrishnan, N. On families of beta- and generalized gamma-generated distributions and associated inference. Statistical Methodology. 2009, 6, 344–362. 4. Alexander C; Cordeiro G.M; Ortega E.M.M.; Sara-bia J.M. Generalized beta-generated distributions. Computational Statistics and Data Analysis., 2012, 56:1880 – 1897. 5. Amini, M.; Mir Mostafaee; S.M. T.K.; Ahmadi, J. (2012). Log-gamma-generated families of distributions. Statistics, iFirst.doi:10.1008/02331888.2012.748775 6. Torabi, H. ; Montazari, N.H.. The gamma-uniform distribution and its application. Kybernetika 2012, 48:16-30 7. Alzaatreh, A.; Famoye, F. ; Lee, C. The gamma-normal distribution:properties and applications. Computational Statistics and Data Analysis. 2014, 69, 67–80. 8. Tahir, M.H.; Cordeiro, G.M.; Alizadeh, M.; Mansoor, M.; Zubair, M.; Hamedani, G.G. The odd generalized exponential family of distributions with applications. Journal of Statistical Distributions and Applications, 2015, 2(1), 1–28. 9. Tahir, M. H.; Cordeiro, G. M.; Alzaatreh, A.; Mansoor, M.; Zubair, M. (2015). The Logistic- X Family of Distributions and its Applications. 10. Cordeiro, G.M., Tahir, M.H., Zubair, M., Mansoor, M., Alizadeh, M. and Hamedani, G.G. A new Weibull-G family of distributions. Hacettepe Journal of Mathematics and Statistics, 2016, 45(2), 629–647. 11. Hassan A.S. ;Elgarhy Kumaraswamy Weibull- Generated Family of Distributions with Applications. Advances and Applications in Statistics, 2016, Volume 48, Number 3, 2016:205- 239 12. Nofal M.Z.; Ahmed Z. A.; Haitham M.Y. Gauss M.C. Communications in Statistic-Theory and Methods. Volume 46, 2017-Issue 8 13. Nofal, Z.M.; Afify, A.Z.; Yousof, H.M.; Cordeiro, G.M. The generalized transmuted-G family of distributions. Communication in Statistics-Theory Methods. 2017, 46, 4119–4136 14. Alzaatreh, A.; Famoye, F. ;Lee, C. A new method for generating families ofcontinuous distributions. Metron. 2013, 71, 63–79. 15. Ristic, M.M.; Balakrishnan, N. The gamma-exponentiated exponentialdistribution. Journal of Statistical Computation and Simulation, 2012, 82(8),1191-1206 16. Torabi, H. and Montazari, N.H. The logistic-uniform distribution and itsapplication. Communications in Statistics–Simulation and Computation. 2014, 43, 2551–2569. IHJPAS. 36 (3) 2023 449 17. Alzaghal, A.; Lee, C. ; Famoye, F. Exponentiated T-X family of distributions with some applications. International Journal of Probability and Statistics.2013,. 2, 31–49. 18. Protein data bank (2015) www.pdb.org 19. Syed, J.; Palani, S.; Clarke, S.T.; Asad, Z.; Bottrill, A.R.; Jones, A.M.E., Sampath, K.,and Balasubraman, M.K. Expanding the Zebrafish Genetic Code through Site-Specific Introduction of Azido-Iysine, Bicyclononyne-Iysine, and Diazirine-Iysine.International Journal of Molecular Sciences. 20(10). 20. Adejoro, I.A; Waheed, S.O ; Adeboye, O.O. Molecular Docking Studies of Lonchocarpus cyanescens Triterpenoids as inhibitors for Malaria. Journal of Physical Chemistry & Biophysics, 2016, 6,213 21. Hardle, W. & Simar, L. Applied Multivariate Statistical Analysis. Springer. 2015. http://www.pdb.org/