New diverse types of soliton solutions to the Radhakrishnan-Kundu-Lakshmanan equation

Emad H. M. Zahran; Omar Abu Arqub; Ahmet Bekir; Marwan Abukhaled; Emad H. M. Zahran; Omar Abu Arqub; Ahmet Bekir; Marwan Abukhaled

doi:10.3934/math.2023450

AIMS Mathematics

2023, Volume 8, Issue 4: 8985-9008. doi: 10.3934/math.2023450

Previous Article Next Article

Research article Special Issues

New diverse types of soliton solutions to the Radhakrishnan-Kundu-Lakshmanan equation

1.
Department of Mathematical and Physical Engineering, Benha University, Faculty of Engineering, Shubra, Egypt
2.
Department of Mathematics, Faculty of Science, Al-Balqa Applied University, Jordan
3.
Neighbourhood of Akcaglan, Imarli Street, 28/4, 26030, Eskisehir, Turkey
4.
Department of Mathematics and Statistics, American University of Sharjah, United Arab Emirates

Received: 28 October 2022 Revised: 11 January 2023 Accepted: 18 January 2023 Published: 10 February 2023
MSC : 35C07, 35Q51, 83C15

The main purpose of this study was to produce abundant new types of soliton solutions for the Radhakrishnan-Kundu-Lakshmanan equation that represents unstable optical solitons that emerge from optical propagations through the use of birefringent fibers. These new types of soliton solutions have behaviors that are bright, dark, W-shaped, M-shaped, periodic trigonometric, and hyperbolic and were not realized before by any other method. These new forms have been detected by using four different techniques, which are, the extended simple equation method, the Paul-Painlevé approach method, the Ricatti-Bernoulli-sub ODE, and the solitary wave ansatz method. These new solitons will be arranged to create a soliton catalog with new impressive behaviors and they will contribute to future studies not only for this model but also for the optical propagations through birefringent fiber.

Keywords:

Citation: Emad H. M. Zahran, Omar Abu Arqub, Ahmet Bekir, Marwan Abukhaled. New diverse types of soliton solutions to the Radhakrishnan-Kundu-Lakshmanan equation[J]. AIMS Mathematics, 2023, 8(4): 8985-9008. doi: 10.3934/math.2023450

Related Papers:

[1]	Wael W. Mohammed, Kalpasree Sharma, Partha Jyoti Hazarika, G. G. Hamedani, Mohamed S. Eliwa, Mahmoud El-Morshedy . Zero-inflated discrete Lindley distribution: Statistical and reliability properties, estimation techniques, and goodness-of-fit analysis. AIMS Mathematics, 2025, 10(5): 11382-11410. doi: 10.3934/math.2025518
[2]	Alanazi Talal Abdulrahman, Khudhayr A. Rashedi, Tariq S. Alshammari, Eslam Hussam, Amirah Saeed Alharthi, Ramlah H Albayyat . A new extension of the Rayleigh distribution: Methodology, classical, and Bayes estimation, with application to industrial data. AIMS Mathematics, 2025, 10(2): 3710-3733. doi: 10.3934/math.2025172
[3]	Khaled M. Alqahtani, Mahmoud El-Morshedy, Hend S. Shahen, Mohamed S. Eliwa . A discrete extension of the Burr-Hatke distribution: Generalized hypergeometric functions, different inference techniques, simulation ranking with modeling and analysis of sustainable count data. AIMS Mathematics, 2024, 9(4): 9394-9418. doi: 10.3934/math.2024458
[4]	Ehab M. Almetwally, Ahlam H. Tolba, Dina A. Ramadan . Bayesian and non-Bayesian estimations for a flexible reduced logarithmic-inverse Lomax distribution under progressive hybrid type-Ⅰ censored data with a head and neck cancer application. AIMS Mathematics, 2025, 10(4): 9171-9201. doi: 10.3934/math.2025422
[5]	Mohamed S. Eliwa, Essam A. Ahmed . Reliability analysis of constant partially accelerated life tests under progressive first failure type-II censored data from Lomax model: EM and MCMC algorithms. AIMS Mathematics, 2023, 8(1): 29-60. doi: 10.3934/math.2023002
[6]	Mohamed S. Algolam, Mohamed S. Eliwa, Mohamed El-Dawoody, Mahmoud El-Morshedy . A discrete extension of the Xgamma random variable: mathematical framework, estimation methods, simulation ranking, and applications to radiation biology and industrial engineering data. AIMS Mathematics, 2025, 10(3): 6069-6101. doi: 10.3934/math.2025277
[7]	Juxia Xiao, Ping Yu, Zhongzhan Zhang . Weighted composite asymmetric Huber estimation for partial functional linear models. AIMS Mathematics, 2022, 7(5): 7657-7684. doi: 10.3934/math.2022430
[8]	Mohamed Ahmed Mosilhy . Discrete Erlang-2 distribution and its application to leukemia and COVID-19. AIMS Mathematics, 2023, 8(5): 10266-10282. doi: 10.3934/math.2023520
[9]	Nora Nader, Dina A. Ramadan, Hanan Haj Ahmad, M. A. El-Damcese, B. S. El-Desouky . Optimizing analgesic pain relief time analysis through Bayesian and non-Bayesian approaches to new right truncated Fréchet-inverted Weibull distribution. AIMS Mathematics, 2023, 8(12): 31217-31245. doi: 10.3934/math.20231598
[10]	Mustafa M. Hasaballah, Yusra A. Tashkandy, Oluwafemi Samson Balogun, M. E. Bakr . Reliability analysis for two populations Nadarajah-Haghighi distribution under Joint progressive type-II censoring. AIMS Mathematics, 2024, 9(4): 10333-10352. doi: 10.3934/math.2024505

Abstract

1. Introduction

The count data sets emerge in various fields like the yearly number of destructive earthquakes, number of patients of a specific disease in a hospital ward, failure of machines, number of patients due to coronavirus, number of monthly traffic accidents, hourly bacterial growth, and so on. Various discrete probability models have been utilized to model these kinds of data sets. Poisson and negative binomial distributions are frequently for modeling count observations. On the other hand, in the advanced scientific eon, the data generated from different fields is getting complex day by day, however, existing discrete models do not provide an efficient fit.

Discretization of continuous distribution can be applied by using different approaches (survival discretization-mixed-Poisson-infinite series). The most widely used technique is the survival discretization approach by [1]. One of the important virtues of this methodology is that the generated discrete model retains the same functional form of the survival function as that of its continuous counterpart. Due to this feature, many survival characteristics of the distribution remain unchanged. The discretization approach to any continuous model depends on the domain of the random variable $X$ . In literature, several discrete distributions are introduced and studied based on the discretization survival function. See for example, discrete Rayleigh [2], discrete a mixture of gamma and exponential [3], discrete Burr and Pareto [4], discrete Rayleigh generator [5], discrete Lindley [6], discrete Burr-XII [7], discrete Burr III [8], two-parameter discrete Lindley [9], discrete log-logistic [10], discrete Gompertz [11], discrete alpha power inverse Lomax [12], discrete Poisson and Ailamujia [13], discrete Half-Logistic [14], discrete Marshall-Olkin Weibull [15], discrete Gompertz-G family [16], discrete Burr-Hatke [17], three-parameter discrete Lindley [18], new discrete Lindley [19], exponentiated discrete Lindley [20], discrete generalized Lindley [21], discrete Gumble [22], discrete inverted Topp-Leone [23], and references cited therein.

Although various distributions are available in literate to analyze count observations, there is still a need to introduce a more flexible and suitable distribution under different conditions. The fundamental purpose of this paper is to propose discrete Ramos-Louzada distribution, which is a one-parameter lifetime distribution introduced by [24]. The proposed one-parameter distribution herein has distinctive properties which makes it among the best choice for modeling over-dispersed and positively skewed data with leptokurtic-shaped. A continuous random variable $X$ is said to have Ramos-Louzada distribution if its probability density function (pdf) can be written as

$g\left(x;\lambda \right) = \frac{1}{{\lambda }^{2}\left(\lambda -1\right)}\left({\lambda }^{2}-2\lambda +x\right){e}^{-\frac{x}{\lambda }}; x\ge 0, \lambda \ge 2,$

(1)

where λ is the shape parameter. The corresponding survival function (sf) to Eq (1) can be formulated as

$G\left(x;\lambda \right) = \frac{{\lambda }^{2}-\lambda +x}{\lambda \left(\lambda -1\right)}{e}^{-\frac{x}{\lambda }}; x\ge 0, \lambda \ge 2.$

(2)

In this article, the discrete version of Ramos and Louzada distribution is proposed and studied in detail. The following are some interesting features of the proposed distribution: Its statistical and reliability characteristics can be expressed as closed forms. Its failure rate is showing an increasing pattern. The suggested distribution evaluated time and count data sets more effectively than competing distributions. As a result, we feel that the proposed model is the greatest option for attracting a wider range of applications and industries.

The rest of the study is organized as follows: In Section 2, we introduce a new distribution using survival discretization methodology. Different mathematical properties are derived in Section 3. Parameter estimation and simulation study are presented in Section 4. Four data sets are utilized to show the flexibility of the proposed model in Section 5. Finally, Section 6 provides some conclusions.

2. Synthesis of discrete Ramos-Louzada distribution

Let $Y$ be a continuous random variable with sf $G\left(y;\boldsymbol{\eta }\right),$ then the pmf of the discrete random variable $X = \left\lfloor Y \right\rfloor$ can be expressed as

$Pr\left(X = x;\boldsymbol{\eta }\right) = G\left(x;\boldsymbol{\eta }\right)-G\left(x+1;\boldsymbol{\eta }\right); x\in {\mathbb{Z}}_{+},$

where $\left\lfloor \cdot \right\rfloor$ denotes the floor function, which returns the highest integer value smaller or equal than its argument, and $\boldsymbol{\eta }$ is a parameter vector $1\times k$ . If the random variable $Y$ have Ramos-Louzada (RL) distribution, then the pmf of discrete RL (DRL) distribution can be written as

$\mathrm{P}\mathrm{r}\left(X = x;\lambda \right) = p\left(x;\lambda \right) = \frac{{e}^{-\frac{x}{\lambda }}}{\lambda \left(\lambda -1\right)}\left[({\lambda }^{2}-\lambda +x)(1-{e}^{-\frac{1}{\lambda }}){- e}^{-\frac{1}{\lambda }}\right]; x = \mathrm{0, 1}, 2, \dots ,$

(3)

where $\lambda \ge 2$ is the shape parameter. The pmf $p\left(x+1;\lambda \right)$ can be expressed as a weight from the pmf $p\left(x;\lambda \right)$ as follows

$p\left(x+1;\lambda \right) = \frac{{e}^{-\frac{1}{\lambda }}\left[\left({\lambda }^{2}-\lambda +x+1\right)\left(1-{e}^{-\frac{1}{\lambda }}\right){-e}^{-\frac{1}{\lambda }}\right]}{\left[\left({\lambda }^{2}-\lambda +x\right)\left(1-{e}^{-\frac{1}{\lambda }}\right){-e}^{-\frac{1}{\lambda }}\right]}p\left(x;\lambda \right).$

Figure 1 illustrates some pmf plots of the DRL models based on different values of the model parameter $\lambda .$ It is found that the pmf can be used a positively skewed data with a uni-modal shape.

Figure 1. The pmf plots of the DRL distribution.

DownLoad: Full-Size Img PowerPoint

Based on $F\left(x;\lambda \right) = \mathrm{P}\mathrm{r}\left(X\le x;\lambda \right) = 1-G\left(x;\lambda \right)+\mathrm{P}\mathrm{r}\left(X = x;\lambda \right)$ , the cumulative distribution function (cdf) of the DRL distribution can be formulated as

$F\left(x;\lambda \right) = 1-\frac{{\lambda }^{2}-\lambda +x+1}{\lambda \left(\lambda -1\right)}{e}^{-\frac{x+1}{\lambda }}; x = \mathrm{0, 1}, 2, \dots .$

(4)

The corresponding sf to Eq (4) can be expressed as

$S\left(x;\lambda \right) = \frac{{\lambda }^{2}-\lambda +x}{\lambda \left(\lambda -1\right)}{e}^{-\frac{x}{\lambda }}; x = \mathrm{0, 1}, 2, \dots .$

(5)

The hazard rate function (hrf) of the DRL model is given by

$h\left(x;\lambda \right) = \frac{\left({\lambda }^{2}-\lambda +x\right)\left(1-{e}^{-\frac{1}{\lambda }}\right)-{e}^{-\frac{1}{\lambda }}}{\left({\lambda }^{2}-\lambda +x\right)}; x = \mathrm{0, 1}, 2, \dots ,$

(6)

where $h\left(x;\lambda \right) = \frac{\mathrm{P}\mathrm{r}\left(X = x;\lambda \right)}{1-F\left(x-1;\lambda \right)}.$ Mathematically, the shape of the hrf of the DRL model is always increasing, which makes it an effective statistical tool for modeling data, especially in the engineering and medical fields. Figure 2 shows the hrf plots of the new discrete model based on various values of the model parameter.

Figure 2. The hrf plots of the DRL distribution.

DownLoad: Full-Size Img PowerPoint

The reversed hazard rate function (rhrf) and the second rate of failure are given as

$\stackrel{˘}{r} = \frac{{e}^{-\frac{x}{\lambda }}\left[({\lambda }^{2}-\lambda +x)(1-{e}^{-\frac{1}{\lambda }}){-e}^{-\frac{1}{\lambda }}\right]}{\lambda \left(\lambda -1\right)-\left({\lambda }^{2}-\lambda +x+1\right){e}^{-\frac{x+1}{\lambda }}}; x = \mathrm{0, 1}, 2, \dots$

(7)

and

${r}^{*}\left(x\right) = log \left[\frac{\left({\lambda }^{2}-\lambda +x\right){e}^{\frac{1}{\lambda }}}{{\lambda }^{2}-\lambda +x+1}\right]; x = \mathrm{0, 1}, 2, \dots ,$

(8)

where $\stackrel{˘}{r} = \frac{\mathrm{P}\mathrm{r}\left(X = x;\lambda \right)}{F\left(x;\lambda \right)}$ and ${r}^{*}\left(x\right) = \mathrm{log}\left[\frac{S(x;\lambda )}{S(x+1;\lambda )}\right]$ . Mathematically, and after simple algebra steps, it is found that the shape of the rhrf of the DRL model is decreasing only. Figure 3 shows some rhrf plots of the proposed model based on specific parameter values.

Figure 3. The rhrf plots of the DRL distribution.

DownLoad: Full-Size Img PowerPoint

3. Distributional properties

In this Section, the probability generating function (pgf) as well as its rth moment are investigated. Assume the random variable $X$ have a DRL model, then the pgf can be expressed as

${W}_{X}\left(z\right) = \sum\limits_{x = 0}^{\infty }{z}^{x}\mathrm{P}\mathrm{r}\left(X = x;\lambda \right) = 1+\left(z-1\right)\sum\limits_{x = 1}^{\infty }{z}^{x-1}S\left(x;\lambda \right) \\ = 1+{e}^{-\frac{1}{\lambda }}\left(z-1\right)\left[\frac{1}{\left(1-z{e}^{-\frac{1}{\lambda }}\right)}+\frac{1}{\lambda \left(\lambda -1\right){\left(1-z{e}^{-\frac{1}{\lambda }}\right)}^{2}}\right].$

(9)

On replacing z by ${e}^{z}$ in Eq (9), the moment generating function (mgf) can be written as

${M}_{X}\left(z\right) = 1+{e}^{-\frac{1}{\lambda }}\left({e}^{z}-1\right)\left[\frac{1}{\left(1-{e}^{z}{e}^{-\frac{1}{\lambda }}\right)}+\frac{1}{\lambda \left(\lambda -1\right){\left(1-{e}^{z}{e}^{-\frac{1}{\lambda }}\right)}^{2}}\right].$

(10)

The first four moments around the origin $\left({\mu \text{'}}_{1}, {\mu \text{'}}_{2}, {\mu \text{'}}_{3}, {\mu \text{'}}_{4}\right)$ can be written as

${\mu \text{'}}_{1} = \frac{{\lambda -\lambda }^{2}+\left({\lambda }^{2}-\lambda +1\right){e}^{\frac{1}{\lambda }}}{\lambda \left(\lambda -1\right){\left({e}^{\frac{1}{\lambda }}-1\right)}^{2}},$

${\mu \text{'}}_{2} = \frac{\left(\lambda \left(\lambda -1\right)+1\right){e}^{\frac{2}{\lambda }}+{3e}^{\frac{1}{\lambda }}-\lambda \left(\lambda -1\right)}{\lambda \left(\lambda -1\right){\left({e}^{\frac{1}{\lambda }}-1\right)}^{3}},$

${{\mu }^{\text{'}}}_{3} = \frac{\left[\lambda \left(\lambda -1\right)+1\right]{e}^{\frac{3}{\lambda }}+\left[3\lambda \left(\lambda -1\right)+10\right]{e}^{\frac{2}{\lambda }}-\left[3\lambda \left(\lambda -1\right)-7\right]{e}^{\frac{1}{\lambda }}-\lambda \left(\lambda -1\right)}{\lambda \left(\lambda -1\right){\left({e}^{\frac{1}{\lambda }}-1\right)}^{4}}$

and

${\mu \text{'}}_{4} = \frac{\left[\lambda \left(\lambda -1\right)+1\right]{e}^{\frac{4}{\lambda }}+\left[10\lambda \left(\lambda -1\right)+25\right]{e}^{\frac{3}{\lambda }}+67{e}^{\frac{2}{\lambda }}+\left[10\lambda \left(\lambda -1\right)-3\right]{e}^{\frac{1}{\lambda }}-\lambda \left(\lambda -1\right)}{\lambda \left(\lambda -1\right){\left({e}^{\frac{1}{\lambda }}-1\right)}^{5}}.$

Based on the rth moments, the variance can be expressed as

${\sigma }^{2} = \frac{\left[{\lambda }^{4}-2{\lambda }^{3}+2{\lambda }^{2}-\lambda \right]{e}^{\frac{3}{\lambda }}-\left[2{\lambda }^{4}-4{\lambda }^{3}+2{\lambda }^{2}-1\right]{e}^{\frac{2}{\lambda }}+\left[{\lambda }^{4}-2{\lambda }^{3}+\lambda \right]{e}^{\frac{1}{\lambda }}}{{\lambda }^{2}{\left(\lambda -1\right)}^{2}{\left({e}^{\frac{1}{\lambda }}-1\right)}^{4}}.$

(11)

The dispersion index (di) is defined by variance to mean ratio. The di indicates that the reported model is suitable for under-, equi- or over-dispersed data sets. Using the derived moments, the coefficients skewness and kurtosis can be listed in closed forms. Some numerical computations for mean, variance, di, skewness, and kurtosis based on DRL parameters are listed in Table 1.

Table 1. Mean, variance, di, skewness, and kurtosis of the DRL distribution.

	$\lambda$
Measure	2.0	2.5	3.0	3.5	4.0	4.5	5.0	5.5	8.0
Mean	3.500	3.678	4.014	4.414	4.847	5.299	5.762	6.234	8.652
Variance	8.079	11.80	15.71	20.03	24.82	30.10	35.87	42.13	80.90
di	2.308	3.207	3.913	4.538	5.121	5.680	6.224	6.757	9.351
Skewness	1.395	1.519	1.629	1.708	1.764	1.806	1.837	1.862	1.928
Kurtosis	5.940	6.353	6.821	7.199	7.495	7.727	7.910	8.057	8.484

| Show Table

DownLoad: CSV

According to Table 1, it is noted that the DRL model can be used effectively to model overdispersion data as di is greater than one, which makes it a proper probability tool to discuss actuarial data. Moreover, the new discrete probabilistic model can be utilized to analyze positively skewed data with leptokurtic-shaped.

4. Estimation methods

In this section, six estimation methods are used to estimate the unknown parameter of DRL distribution. The considered estimation methods are maximum likelihood estimation (mle), method of moments (mom), least-squares estimation (lse), Anderson-Darling estimation (ade), Cramer von-Misses estimation (cvme), and maximum product of spacing estimator (mpse).

4.1. Maximum likelihood estimation

Assume a random sample ${x}_{1}, {x}_{2}, \dots , {x}_{n}$ from the DRL model, then the log-likelihood function can be expressed as

$L\left(\lambda |x\right) = -\frac{1}{\lambda }\sum _{i = 1}^{n}{x}_{i}+\sum _{i = 1}^{n}\mathrm{ln}\left[\left({\lambda }^{2}-\lambda +{x}_{i}\right)\left(1-{e}^{-\frac{1}{\lambda }}\right){-e}^{-\frac{1}{\lambda }}\right]-n\mathrm{ln}\lambda -n\mathrm{ln}\left(\lambda -1\right).$

(12)

Differentiating the Eq (12) with respect to the parameter $\lambda$ , we get the non-linear equation as follows

$\frac{\partial L\left(\lambda |x\right)}{\partial \lambda } = \frac{1}{{\lambda }^{2}}\sum _{i = 1}^{n}{x}_{i}+\sum _{i = 1}^{n}\frac{\left[\left(2\lambda -1\right)\left(1-{e}^{-\frac{1}{\lambda }}\right)-\left(1-\frac{1}{\lambda }+\frac{1}{{\lambda }^{2}}+\frac{{x}_{i}}{{\lambda }^{2}}\right){e}^{-\frac{1}{\lambda }}\right]}{\left[\left({\lambda }^{2}-\lambda +{x}_{i}\right)\left(1-{e}^{-\frac{1}{\lambda }}\right){-e}^{-\frac{1}{\lambda }}\right]}-\frac{n}{\lambda }-\frac{n}{\lambda -1},$

(13)

the exact solution of Eq (13) is not easy, so we will maximize it by using optimization approaches, for example, the Newton-Raphson approach using R software.

4.2. Method of moment estimation

Based on the mom definition, we must equate the sample mean to the corresponding population mean, and then solve the non-linear equation for the parameter

$\frac{{\lambda -\lambda }^{2}+\left({\lambda }^{2}-\lambda +1\right){e}^{\frac{1}{\lambda }}}{\lambda \left(\lambda -1\right){\left({e}^{\frac{1}{\lambda }}-1\right)}^{2}} = \frac{1}{n}\sum _{i = 1}^{n}{x}_{i}.$

(14)

To solve Eq (14), the uniroot function should be utilized.

4.3. Least squares estimation

To estimate the parameter minimizing the sum of squares of residuals, a standard approach like the lse should be used. For the estimation of the parameter of DRL distribution, the lse can be obtained by minimizing

$\mathrm{l}\mathrm{s}\mathrm{e}\left(\lambda \right) = \sum _{i = 1}^{n}{\left[1-\frac{\left({\lambda }^{2}-\lambda +{\mathrm{x}}_{i:n}\right)}{\lambda \left(\lambda -1\right)}{e}^{-\frac{{\mathrm{x}}_{i:n}+1}{\lambda }}-\frac{i}{n+2}\right]}^{2},$

(15)

with respect to the parameter $\lambda .$

4.4. Anderson-Darling estimation

The ade of the parameter can be derived by minimizing the following equation

$\mathrm{a}\mathrm{d}\mathrm{e}\left(\lambda \right) = -n-\frac{1}{n}\sum _{i = 1}^{n}{\left({2}_{i}-1\right)\left[\mathrm{log}\left(1-\frac{\left({\lambda }^{2}-\lambda +{\mathrm{x}}_{i:n}\right)}{\lambda \left(\lambda -1\right)}{e}^{-\frac{{\mathrm{x}}_{i:n}+1}{\lambda }}\right)+\mathrm{log}\left(\frac{\left({\lambda }^{2}-\lambda +{\mathrm{x}}_{i:n}\right)}{\lambda \left(\lambda -1\right)}{e}^{-\frac{{\mathrm{x}}_{i:n}+1}{\lambda }}\right)\right]}^{2},$

(16)

with respect to the parameter $\lambda .$

4.5. Cramer von-Misses estimation

The cvme is an estimation method. This method is derived as the difference between the empirical cdf and fitted cdf where

$\mathrm{c}\mathrm{v}\mathrm{m}\mathrm{e}\left(\lambda \right) = \frac{1}{12n}+\sum _{i = 1}^{n}{\left[1-\frac{\left({\lambda }^{2}-\lambda +{\mathrm{x}}_{i:n}\right)}{\lambda \left(\lambda -1\right)}{e}^{-\frac{{\mathrm{x}}_{i:n+1}}{\lambda }}-\frac{2i-1}{2n}\right]}^{2}.$

(17)

4.6. Maximum product of spacing estimator

For $u = \mathrm{1, 2}, 3, \dots , h+1$ , assume ${D}_{u}\left(\lambda \right) = F\left({x}_{\left(u\right)}|\lambda \right)-F\left({x}_{\left(u-1\right)}|\lambda \right),$ be the uniform spacings of a random sample from the DRL model, where $F\left({x}_{\left(0\right)}|\lambda \right) = 0$ , $F\left({x}_{(h+1)}|\lambda \right) = 1$ and ${\sum }_{r = 1}^{h+1}{D}_{u}\left(\lambda \right) = 1$ . The mpse of the parameter $\lambda$ , say $\widehat{\lambda }$ , can be estimated by maximizing the geometric mean of the spacings

$\mathrm{m}\mathrm{p}\mathrm{s}\mathrm{e}\left(\lambda \right) = {\left[\prod _{u = 1}^{h+1}{D}_{u}\left(\lambda \right)\right]}^{\frac{1}{h+1}},$

(18)

with respect to the parameter $\lambda$ .

5. Simulation

In this section, we discussed the results of the simulation study to compare the estimation performance of the proposed estimators based on the DRL model. The performance of considered estimators is evaluated via absolute biases and mean square errors. We simulate $\mathrm{10,000}$ random samples from the DRL distribution using the following sample sizes $n = 10, 20, 50,100,200$ , and $500$ . The results calculated for parameter values $\lambda = \left(2.0, 2.5, 3.0, 4.0, 5.0, 6.0, 8.0, 10.0\right)$ using R Software. The absolute bias and mse for the parameter $\lambda$ are presented in Tables 2 and 3.

Table 2. Simulation results of the DRL distribution for different parameter values.

Para.	n	Bias						mse
$\lambda$		mle	mom	ade	cvme	lse	mpse	mle	mom	ade	cvme	lse	mpse
2.0	10	0.099	0.072	0.878	0.606	0.956	0.801	0.762	1.496	2.897	2.059	3.034	1.797
	20	0.034	0.035	0.486	0.259	0.444	0.435	0.307	1.103	1.408	0.749	1.212	0.626
	50	0.017	0.160	0.122	0.037	0.065	0.172	0.055	0.811	0.304	0.086	0.148	0.111
	100	0.013	0.226	0.015	0.003	0.005	0.086	0.017	0.697	0.033	0.006	0.009	0.026
	200	0.009	0.277	0.000	0.000	0.000	0.050	0.007	0.620	0.000	0.000	0.000	0.008
	500	0.002	0.338	0.000	0.000	0.000	0.028	0.003	0.569	0.000	0.000	0.000	0.002
2.5	10	0.081	0.131	1.109	1.140	1.397	0.958	1.105	2.032	3.694	3.919	4.392	2.609
	20	0.065	0.239	0.837	0.869	1.060	0.622	0.485	1.487	2.368	2.560	2.765	1.173
	50	0.058	0.293	0.595	0.695	0.798	0.320	0.186	1.063	1.436	1.723	1.785	0.391
	100	0.027	0.298	0.475	0.602	0.672	0.195	0.091	0.816	1.082	1.373	1.428	0.165
	200	0.013	0.255	0.374	0.525	0.589	0.098	0.045	0.602	0.864	1.139	1.185	0.065
	500	0.007	0.146	0.273	0.484	0.532	0.047	0.018	0.316	0.674	0.979	1.014	0.022
3.0	10	0.019	0.224	1.085	1.103	1.304	0.996	1.574	2.613	3.976	4.160	4.504	3.186
	20	0.005	0.268	0.877	0.924	1.044	0.639	0.751	1.793	2.525	2.701	2.865	1.512
	50	0.003	0.193	0.738	0.844	0.893	0.357	0.313	1.003	1.540	1.728	1.760	0.537
	100	0.004	0.130	0.754	0.868	0.872	0.222	0.163	0.567	1.160	1.352	1.341	0.242
	200	0.001	0.053	0.775	0.887	0.896	0.141	0.084	0.238	0.920	1.112	1.106	0.113
	500	0.001	0.006	0.811	0.929	0.940	0.070	0.035	0.066	0.769	0.965	0.980	0.042
4.0	10	0.049	0.209	1.105	1.023	1.181	1.100	2.450	3.694	5.254	5.458	5.552	4.850
	20	0.043	0.169	0.896	0.891	0.985	0.748	1.339	2.065	3.000	3.234	3.315	2.269
	50	0.031	0.049	0.835	0.847	0.899	0.400	0.585	0.763	1.608	1.739	1.789	0.815
	100	0.008	0.014	0.849	0.878	0.891	0.255	0.291	0.328	1.161	1.240	1.285	0.377
	200	0.007	0.016	0.838	0.877	0.882	0.156	0.147	0.159	0.911	1.005	1.010	0.174
	500	0.004	0.002	0.831	0.877	0.878	0.066	0.058	0.064	0.773	0.861	0.862	0.062

| Show Table

DownLoad: CSV

Table 3. Simulation results of the DRL distribution for different parameter values.

Para.	n	Bias						mse
$\lambda$		mle	mom	ade	cvme	lse	mpse	mle	mom	ade	cvme	lse	mpse
5.0	10	0.106	0.077	1.058	1.013	1.242	1.164	3.881	4.545	6.333	6.748	7.579	6.421
	20	0.061	0.036	0.955	0.915	1.007	0.824	2.026	2.335	3.527	3.991	4.014	3.056
	50	0.046	0.015	0.864	0.869	0.904	0.447	0.821	0.839	1.793	1.906	1.970	1.070
	100	0.013	0.002	0.857	0.866	0.857	0.270	0.399	0.411	1.244	1.325	1.303	0.513
	200	0.013	0.006	0.833	0.843	0.849	0.161	0.198	0.214	0.946	1.006	1.009	0.238
	500	0.001	0.003	0.832	0.846	0.843	0.082	0.079	0.078	0.792	0.829	0.822	0.086
6.0	10	0.167	0.115	1.192	1.073	1.169	1.260	5.339	5.663	8.785	8.977	9.398	8.321
	20	0.049	0.027	0.946	0.952	0.975	0.890	2.794	2.926	4.236	4.861	4.686	3.860
	50	0.008	0.019	0.873	0.889	0.881	0.499	1.081	1.079	2.081	2.225	2.292	1.417
	100	0.007	0.016	0.859	0.843	0.870	0.294	0.523	0.517	1.392	1.421	1.522	0.631
	200	0.007	0.006	0.845	0.839	0.846	0.181	0.261	0.269	1.030	1.061	1.065	0.308
	500	0.001	0.001	0.832	0.832	0.835	0.085	0.107	0.108	0.823	0.841	0.847	0.112
8.0	10	0.141	0.089	1.294	1.077	1.320	1.622	8.569	8.656	12.77	12.61	13.63	13.41
	20	0.074	0.009	1.050	0.991	1.111	1.110	4.265	4.178	6.399	6.736	7.355	6.250
	50	0.029	0.003	0.883	0.860	0.918	0.577	1.685	1.630	2.726	3.019	3.118	2.151
	100	0.015	0.014	0.867	0.861	0.859	0.372	0.834	0.814	1.787	1.862	1.859	1.029
	200	0.008	0.019	0.837	0.808	0.842	0.219	0.416	0.422	1.202	1.204	1.285	0.493
	500	0.006	0.001	0.826	0.814	0.820	0.097	0.169	0.175	0.877	0.879	0.902	0.176
10.0	10	0.106	0.086	1.365	1.333	1.446	1.868	12.89	12.67	17.66	19.87	19.51	20.26
	20	0.001	0.063	1.143	1.079	1.067	1.308	6.232	6.279	8.767	9.553	9.777	9.056
	50	0.025	0.036	0.899	0.892	0.933	0.686	2.483	2.373	3.774	4.076	4.248	3.246
	100	0.018	0.016	0.868	0.886	0.872	0.457	1.247	1.246	2.255	2.421	2.464	1.489
	200	0.001	0.009	0.852	0.827	0.829	0.267	0.636	0.599	1.467	1.501	1.501	0.712
	500	0.002	0.007	0.823	0.801	0.818	0.119	0.246	0.240	0.967	0.971	1.005	0.255

| Show Table

DownLoad: CSV

Based on the simulation criteria, it is observed that all estimation approaches work quite well in estimating the parameter λ of the DRL distribution.

6. Data analysis

In this section, the importance of the proposed distribution is discussed by using data sets from different areas. We shall compare the fits of the DRL distribution with different competitive distributions such as Poisson (Poi), discrete Pareto (DPr), discrete Rayleigh (DR), discrete inverse Rayleigh (DIR), discrete Burr-Hatke (DBH), discrete Bilal (DBi), discrete Lindley (DL), new discrete Lindley (NDL), and discrete Burr-XII (DBXII) distributions. The fitted probability distributions are compared using some criteria, namely, the negative log-likelihood ( $-L$ ), Akaike information criterion (aic), and Kolmogorov-Smirnov (ks) test with its p-value.

6.1. Data set I: Covid-19 data in Pakistan

The first data set represents the number of deaths due to coronavirus in Pakistan during the period March 18, 2020, to April 30, 2020, which were obtained from the public reports of the National Institute of Health (NIH), Islamabad, Pakistan (https://covid.gov.pk/stats/pakistan). The mean, variance, and di of data set I are 9.4773,102.39, and 10.804, respectively. The mle(s) along with standard error(s) "se(s)" and goodness-of-fit measures for this data are presented in Table 4.

Table 4. The mle(s), se(s), and goodness-of-fit measures for data set I.

Model	$\lambda$		$\delta$		Goodness-of-fit measures
Model	mle	se	mle	se	$-L$	aic	ks	p-value
DRL	8.8686	1.5033	-	-	145.22	292.43	0.156	0.2300
Poi	9.4773	0.4641	-	-	283.94	569.89	0.391	< 0.0001
DPr	0.5021	0.0757	-	-	162.19	326.38	0.401	< 0.0001
DR	9.9883	0.7535	-	-	168.85	339.70	0.339	< 0.0001
DIR	7.4291	1.2625	-	-	166.31	334.61	0.382	< 0.0001
DBH	0.9950	0.0115	-	-	175.37	352.74	0.647	< 0.0001
DBi	11.838	1.2932	-	-	151.29	304.59	0.213	0.0370
DL	0.8313	0.0165	-	-	149.17	300.33	0.184	0.1000
NDL	0.1640	0.0161	-	-	148.44	298.89	0.237	0.0140
DBXII	0.9536	0.0434	11.907	11.305	150.70	305.40	0.302	0.0007

| Show Table

DownLoad: CSV

The results in Table 4 show that the DRL distribution provides a better fit over other competing discrete models since it has the minimum aic, and ks values with the highest p-value. Figure 4 shows the probability-probability (pp) plots for all tested models which prove the empirical results listed in Table 4.

Figure 4. The pp plots for all fitted distributions for data set I.

DownLoad: Full-Size Img PowerPoint

6.2. Data set II: Hydrology data

The second data set was reported in [25], which represents the exceedance of flood peaks in m³/s of the Wheaton River near Carcross in Yukon Territory, Canada based on the discretization concept. The mean, variance, and di of this data are 11.806,152.38, and 12.908, respectively. The mle(s), se(s), and goodness-of-fit measures for data set II are reported in Table 5.

Table 5. The mle(s), se(s), and goodness-of-fit measures for data set II.

Model	$\lambda$		$\delta$		Goodness-of-fit measures
Model	mle	se	mle	se	$-L$	aic	ks	p-value
DRL	11.214	1.4497	-	-	252.71	507.43	0.133	0.1600
Poi	11.805	0.4049	-	-	564.38	1130.8	0.408	< 0.0001
DPr	0.4770	0.0563	-	-	276.82	555.64	0.311	< 0.0001
DR	12.280	0.7239	-	-	300.65	603.29	0.323	< 0.0001
DIR	4.7947	0.6303	-	-	331.46	664.92	0.497	< 0.0001
DBH	0.9966	0.0072	-	-	302.29	606.57	0.572	< 0.0001
DBi	14.621	1.2479	-	-	272.50	546.99	0.257	0.0002
DL	0.8592	0.0109	-	-	264.30	530.59	0.232	0.0009
NDL	0.1373	0.0107	-	-	262.09	526.17	0.271	< 0.0001
DBXII	0.8205	0.0591	2.6112	0.9287	270.50	544.99	0.228	0.0011

| Show Table

DownLoad: CSV

It is observed that the DRL model is the best among all competitive distributions. Figure 5 illustrates the pp plots for all tested distributions which prove the empirical results reported in Table 5.

Figure 5. The pp plots for all fitted distributions for data set II.

DownLoad: Full-Size Img PowerPoint

6.3. Data set III: Forest fire in Greece

The third data set was listed in [26] and represents the number of fires in Greece forest districts for the period from 1st July 1998 to 31 August 1998. The mean, variance, and di measures are 5.2, 32.382, and 6.2272, respectively. The mle(s), se(s), and goodness-of-fit measures for data set II are listed in Table 6.

Table 6. The mle(s), se(s), and goodness-of-fit measures for data set III.

Model	$\lambda$		$\delta$		Goodness-of-fit measures
Model	mle	se	mle	se	$-L$	aic	ks	p-value
DRL	4.3673	0.5510	-	-	301.11	604.21	0.1510	0.0140
Poi	5.2000	0.2174	-	-	434.16	870.32	0.282	< 0.0001
DPr	0.6250	0.0597	-	-	339.05	680.10	0.352	< 0.0001
DR	5.6788	0.2714	-	-	352.72	707.45	0.261	< 0.0001
DIR	3.5198	0.3748	-	-	360.90	723.80	0.413	< 0.0001
DBH	0.9833	0.0136	-	-	352.42	706.85	0.532	< 0.0001
DBi	6.7993	0.4693	-	-	310.75	623.49	0.107	< 0.0001
DL	0.7337	0.0156	-	-	303.88	609.75	0.193	0.0100
NDL	0.2567	0.0152	-	-	302.73	607.47	0.169	0.0037
DBXII	0.7486	0.0459	2.4582	0.4938	325.00	654.01	0.287	< 0.0001

| Show Table

DownLoad: CSV

It is found that the new discrete model is the best among all tested distributions. Figure 6 shows the pp plots for all competitive distributions which prove the empirical results listed in Table 6.

Figure 6. The pp plots for all fitted distributions for data set III.

DownLoad: Full-Size Img PowerPoint

6.4. Data set IV: AG-positive leukemia

The fourth data set represents the time to death (in weeks) of AG-positive leukemia patients [27]. The mean, variance, and di values are 62.471, 2954.3, and 47.29, respectively. The estimates and goodness-of-fit measures for all competitive distributions are listed in Table 7.

Table 7. The mle(s), se(s), and goodness-of-fit measures for data set IV.

Model	$\lambda$		$\delta$		Goodness-of-fit measures
Model	mle	se	mle	se	$-L$	aic	ks	p-value
DRL	61.943	15.273	-	-	87.425	176.85	0.152	0.8300
Poi	62.470	1.9169	-	-	475.26	952.52	0.470	0.0011
DPr	0.2838	0.0688	-	-	98.335	198.67	0.324	0.0560
DR	58.076	7.0429	-	-	96.794	195.59	0.309	0.0770
DIR	25.310	6.543	-	-	128.59	259.18	0.681	< 0.0001
DBH	0.9999	0.0029	-	-	119.81	241.62	0.716	< 0.0001
DBi	75.109	13.164	-	-	92.886	187.77	0.219	0.3900
DL	0.9692	0.0052	-	-	91.858	185.72	0.215	0.4100
NDL	0.0306	0.0052	-	-	91.458	184.92	0.218	0.3900
DBXII	0.9975	0.0008	117.30	45.982	96.151	196.30	0.327	0.0530

| Show Table

DownLoad: CSV

It is noted that the DRL is the best for this data. Figure 7 shows the pp plots for all tested distributions which prove the empirical results mentioned in Table 7.

Figure 7. The pp plots for all fitted distributions for data set IV.

DownLoad: Full-Size Img PowerPoint

7. Conclusions

In this article, a new one-parameter discrete model has been proposed entitled a discrete Ramos-Louzada (DRL) distribution. The new model can be used effectively in modeling asymmetric data with overdispersion phenomena. Some of its statistical properties have been derived. It was found that all its properties can be expressed in closed forms, which makes the new model can be utilized in different analysis, especially, in time series and regression. Various estimation techniques including maximum likelihood, moments, least squares, Anderson's-Darling, Cramer von-Mises, and maximum product of spacing estimator, have been investigated to get the best estimator for the real data. The estimation performance of these estimation techniques has been assessed via a comprehensive simulation study. The flexibility of the proposed discrete model has been tested utilizing four distinctive real data sets in various fields. Finally, we hope that the DRL distribution attracts wider sets of applications in different fields.

Conflict of interest

The authors declare that they have no conflict of interest to report regarding the present study.

References

[1]	A. Bekir, E. H. M. Zahran, Bright and dark soliton solutions for the complex Kundu-Eckhaus equation, Optik, 223 (2020), 165233. https://doi.org/10.1016/j.ijleo.2020.165233 doi: 10.1016/j.ijleo.2020.165233
[2]	A. Bekir, E. H. M. Zahran, Three distinct and impressive visions for the soliton solutions to the higher-order nonlinear Schrö dinger equation, Optik, 228 (2021), 166157. https://doi.org/10.1016/j.ijleo.2020.166157 doi: 10.1016/j.ijleo.2020.166157
[3]	A. Bekir, E. H. M. Zahran, New vision for the soliton solutions to the complex Hirota-dynamical model, Phys. Scripta, 96 (2021), 055212. https://doi.org/10.1088/1402-4896/abe889 doi: 10.1088/1402-4896/abe889
[4]	A. Biswas, 1-soliton solution of the K (m, n) equation with generalized evolution, Phys. Lett. A, 372 (2008), 4601–4602. https://doi.org/10.1016/j.physleta.2008.05.002 doi: 10.1016/j.physleta.2008.05.002
[5]	H. Triki, A. M. Wazwaz, Bright and dark soliton solutions for a K (m, n) equation with t-dependent coefficients, Phys. Lett. A, 373 (2009), 2162–2165. https://doi.org/10.1016/j.physleta.2009.04.029 doi: 10.1016/j.physleta.2009.04.029
[6]	H. Triki, A. M. Wazwaz, Bright and dark solitons for a generalized Korteweg-de Vries-modified Korteweg-de Vries equation with high-order nonlinear terms and time-dependent coefficients, Can. J. Phys., 89 (2011), 253–259. https://doi.org/10.1139/P11-015 doi: 10.1139/P11-015
[7]	N. A. Kudryashov, The Painlevé approach for finding solitary wave solutions of nonlinear non-integrable differential equations, Optik, 183 (2019), 642–649. https://doi.org/10.1016/j.ijleo.2019.02.087 doi: 10.1016/j.ijleo.2019.02.087
[8]	A. Bekir, E. H. M. Zahran, Optical soliton solutions of the thin-film ferro-electric materials equation according to the Painlevé approach, Opt. Quantum Electron., 53 (2021), 118. https://doi.org/10.1007/s11082-021-02754-w doi: 10.1007/s11082-021-02754-w
[9]	A. Bekir, E. H. M. Zahran, Painlevé approach and its applications to get new exact solutions of three biological models instead of its numerical solutions, Int. J. Mod. Phys. B, 34 (2020), 2050270. https://doi.org/10.1142/S0217979220502707 doi: 10.1142/S0217979220502707
[10]	A. Bekir, E. H. M. Zahran, New visions of the soliton solutions to the modified nonlinear Schrodinger equation, Optik, 232 (2021), 166539. https://doi.org/10.1016/j.ijleo.2021.166539 doi: 10.1016/j.ijleo.2021.166539
[11]	M. S. M. Shehata, H. Rezazadeh, E. H. M. Zahran, E. Tala-Tebue, A. Bekir, New optical soliton solutions of the perturbed Fokas-Lenells equation, Commun. Theor. Phys., 71 (2019), 1275c1280. https://doi.org/10.1088/0253-6102/71/11/1275 doi: 10.1088/0253-6102/71/11/1275
[12]	A. Bekir, E. H. M. Zahran, New multiple-different impressive perceptions for the solitary solution to the magneto-optic waveguides with anti-cubic nonlinearity, Optik, 240 (2021), 166939. https://doi.org/10.1016/j.ijleo.2021.166939 doi: 10.1016/j.ijleo.2021.166939
[13]	A. Biswas, Y. Yildirim, E. Yasar, M. F. Mahmood, A. S. Alshomrani, Q. Zhou, et al., Optical soliton perturbation for Radhakrishnan-Kundu-Lakshmanan equation with a couple of integration schemes, Optik, 163 (2018), 126–136. https://doi.org/10.1016/j.ijleo.2018.02.109
[14]	N. A. Kudryashov, D. V. Safonova, A. Biswas, Painleve analysis and a solution to the traveling wave reduction of the Radhakrishnan-Kundu-Lakshmanan equation, Regul. Chaotic Dyn., 24 (2019), 607–614. https://doi.org/10.1134/S1560354719060029 doi: 10.1134/S1560354719060029
[15]	H. U. Rehman, M. S. Saleem, A. M. Sultan, M. Iftikhar, Comments on dynamics of optical solitons with Radhakrishnan-Kundu-Lakshmanan model via two reliable integration schemes, Optik, 178 (2019), 557–566. https://doi.org/10.1016/j.ijleo.2018.12.010 doi: 10.1016/j.ijleo.2018.12.010
[16]	T. A. Sulaiman, H. Bulut, G. Yel, S. S. Atas, Optical solitons to the fractional perturbed Radhakrishnan-Kundu-Lakshmanan model, Opt. Quant. Electron., 50 (2018), 372. https://doi.org/10.1007/s11082-018-1641-7 doi: 10.1007/s11082-018-1641-7
[17]	B. Sturdevant, D. A. Lott, A. Biswas, Topological 1-soliton solution of the generalized Radhakrishnan, Kundu, Lakshmanan equation with nonlinear dispersion, Mod. Phys. Lett. B, 24 (2010), 1825–1831. https://doi.org/10.1142/S0217984910024109 doi: 10.1142/S0217984910024109
[18]	S. Arshed, A. Biswas, P. Guggilla, A. S. Alshomrani, Optical solitons for Radhakrishnan-Kundu-Lakshmanan equation with full nonlinearity, Phys. Lett. A, 384 (2020), 126191. https://doi.org/10.1016/j.physleta.2019.126191 doi: 10.1016/j.physleta.2019.126191
[19]	A. Bansal, A. Biswas, M. F. Mahmood, Q. Zhou, M. Mirzazadeh, A. S. Alshomrani, et al., Optical soliton perturbation with Radhakrishnan-Kundu-Lakshmanan equation by Lie group analysis, Optik, 163 (2018), 137–141. https://doi.org/10.1016/j.ijleo.2018.02.104
[20]	A. Biswas, 1-soliton solution of the generalized Radhakrishnan, Kundu, Lakshmanan equation, Phys. Lett. A, 373 (2009), 2546–2548. https://doi.org/10.1016/j.physleta.2009.05.010 doi: 10.1016/j.physleta.2009.05.010
[21]	A. Biswas, Optical soliton perturbation with Radhakrishnan-Kundu-Lakshmanan equation by traveling wave hypothesis, Optik, 171 (2018), 217–220. https://doi.org/10.1016/j.ijleo.2018.06.043 doi: 10.1016/j.ijleo.2018.06.043
[22]	A. Biswas, M. Ekici, A. Sonmezoglu, A. S. Alshomrani, Optical solitons with Radhakrishnan, Kundu, Lakshmanan equation by extended trial function scheme, Optik, 160 (2018), 415–427. https://doi.org/10.1016/j.ijleo.2018.02.017 doi: 10.1016/j.ijleo.2018.02.017
[23]	N. A. Kudryashov, The Radhakrishnan-Kundu-Lakshmanan equation with arbitrary refractive index and its exact solutions, Optik, 238 (2021), 166738. https://doi.org/10.1016/j.ijleo.2021.166738 doi: 10.1016/j.ijleo.2021.166738
[24]	D. D. Ganji, A. Asgari, Z. Z. Ganji, Exp-function based solution of nonlinear Radhakrishnan, Kundu and Laskshmanan (RKL) equation, Acta Appl. Math., 104 (2008), 201–209. https://doi.org/10.1007/s10440-008-9252-0 doi: 10.1007/s10440-008-9252-0
[25]	O. Gonzalez-Gaxiola, A. Biswas, Optical solitons with Radhakrishnan-Kundu-Lakshmanan equation by Laplace-Adomian decomposition method, Optik, 179 (2019), 434–442. https://doi.org/10.1016/j.ijleo.2018.10.173 doi: 10.1016/j.ijleo.2018.10.173
[26]	A. Neirameh, Soliton solutions modeling of generalized Radhakrishnan-Kundu-Lakshmanan equation, J. Appl. Phys., 8 (2018), 71–80. https://doi.org/10.22051/JAP.2019.21375.1099 doi: 10.22051/JAP.2019.21375.1099
[27]	S. S. Singh, Solutions of Kudryashov-Sinelshchikov equation and generalized Radhakrishnan-Kundu-Lakshmanan equation by the first integral method, Int. J. Phys. Res., 4 (2016), 37–42. https://doi.org/10.14419/ijpr.v4i2.6202 doi: 10.14419/ijpr.v4i2.6202
[28]	B. Ghanbari, J. F. Gómez-Aguilar, Optical soliton solutions for the nonlinear Radhakrishnan-Kundu-Lakshmanan equation, Mod. Phys. Lett. B, 33 (2019), 1950402. https://doi.org/10.1142/S0217984919504025 doi: 10.1142/S0217984919504025
[29]	S. Rehman, J. Ahmad, Modulation instability analysis and optical solitons in birefringent fibers to RKL equation without four wave mixing, Alex. Eng. J., 60 (2021), 1339–1354. https://doi.org/10.1016/j.aej.2020.10.055 doi: 10.1016/j.aej.2020.10.055
[30]	Y. Yıldırım, A. Biswas, M. Ekici, H. Triki, O. Gonzalez-Gaxiol, A. K. Alzahrani, et al., Optical solitons in birefringent fibers for Radhakrishnan-Kundu-Lakshmanan equation with five prolific integration norms, Optik, 208 (2020), 164550. https://doi.org/10.1016/j.ijleo.2020.164550 doi: 10.1016/j.ijleo.2020.164550
[31]	Y. Yıldırım, A. Biswas, Q. Zhou, A. S. Alshomrani, M. R. Belic, Optical solitons in birefringentfibers for Radhakrishnan-Kundu-Lakshmanan equation with acouple of strategic integration architectures, Chin. J. Phys., 65 (2020), 341–354. https://doi.org/10.1016/j.cjph.2020.02.029 doi: 10.1016/j.cjph.2020.02.029
[32]	J. H. He, Exp-function method for fractional differential equations, Int. J. Nonlinear Sci. Numer. Simul., 14 (2013), 363–366. https://doi.org/10.1515/ijnsns-2011-0132 doi: 10.1515/ijnsns-2011-0132
[33]	Y. Tian, J. Liu, A modified exp-function method for fractional partial differential equations, Therm. Sci., 25 (2021), 1237–1241. https://doi.org/10.2298/TSCI200428017T doi: 10.2298/TSCI200428017T
[34]	F. Y. Ji, C. H. He, J. J. Zhang, A fractal Boussinesq equation for nonlinear transverse vibration of a nanofiber-reinforced concrete pillar, Appl. Math. Modell., 82 (2020) 437–448. https://doi.org/10.1016/j.apm.2020.01.027 doi: 10.1016/j.apm.2020.01.027
[35]	J. H. He, N. Qie, C. H. He, Solitary waves travelling along an unsmooth boundary, Results Phys., 24 (2021), 104104. https://doi.org/10.1016/j.rinp.2021.104104 doi: 10.1016/j.rinp.2021.104104
[36]	J. H. He, W. F. Hou, C. H. He, T. Saeed, T. Hayat, Variational approach to fractal solitary waves, Fractals, 29 (2021), 2150199. https://doi.org/10.1142/S0218348X21501991 doi: 10.1142/S0218348X21501991
[37]	C. X. Liu, Periodic solution of fractal Phi-4 equation, Therm. Sci., 25 (2021), 1345–1350 https://doi.org/10.2298/TSCI200502032L
[38]	J. H. He, E. D. Yusry, Periodic property of the time-fractional Kundu-Mukherjee-Naskar equation, Results Phys., 19 (2020), 103345. https://doi.org/10.1016/j.rinp.2020.103345 doi: 10.1016/j.rinp.2020.103345
[39]	J. H. He, Variational principle and periodic solution of the Kundu-Mukherjee-Naskar equation, Results Phys., 17 (2020), 103031. https://doi.org/10.1016/j.rinp.2020.103031 doi: 10.1016/j.rinp.2020.103031
[40]	J. H. He, W. F. Hou, N. Qie, K. A. Gepreel, A. H. Shirazi, H. M. Sedighi, Hamiltonian-based frequency-amplitude formulation for nonlinear oscillators, Facta Univ. Ser. Mech. Eng., 19 (2021), 199–208. https://doi.org/10.22190/FUME201205002H doi: 10.22190/FUME201205002H
[41]	Y. Zhao, Y. B. Lei, Y. X. Xu, S. L. Xu, H. Triki, A. Biswas, et al., Vector spatiotemporal solitons and their memory features in cold rydberg gases, Chin. Phys. Lett., 39 (2022), 034202. https://doi.org/10.1088/0256-307X/39/3/034202 doi: 10.1088/0256-307X/39/3/034202
[42]	S. L. Xu, Y. B. Lei, J. T. Du, Y. Zhao, R. Hua, J. H. Zeng, Three-dimensional quantum droplets in spin-orbit-coupled Bose-Einstein condensates, Chaos Solitons Fract., 164 (2022), 112665. https://doi.org/10.1016/j.chaos.2022.112665 doi: 10.1016/j.chaos.2022.112665
[43]	K. Y. Huang, Y. Zhao, S. Q. Wu, S. L. Xu, M. R. Belic, B. A. Malomed, Quantum squeezing of vector slow-light solitons in a coherent atomic system, Chaos Solitons Fract., 163 (2022), 112557. https://doi.org/10.1016/j.chaos.2022.112557 doi: 10.1016/j.chaos.2022.112557
[44]	T. A. Nofal, Simple equation method for nonlinear partial differential equations and its applications, J. Egypt. Math. Soc., 24 (2016), 204–209. https://doi.org/10.1016/j.joems.2015.05.006 doi: 10.1016/j.joems.2015.05.006
[45]	N. A. Kudryashov, V. B. Loguinova, Extended simplest equation method for nonlinear differential equations, Appl. Math. Comput., 205 (2008), 396–402. https://doi.org/10.1016/j.amc.2008.08.019 doi: 10.1016/j.amc.2008.08.019
[46]	Y. L. Ma, C. B. Li, Q. Wang, A series of abundant exact travelling wave solutions for a modified generalized Vakhnenko equation using auxiliary equation method, Appl. Math. Comput., 211 (2009), 102–107. https://doi.org/10.1016/j.amc.2009.01.036 doi: 10.1016/j.amc.2009.01.036
[47]	G. B. Whitham, Comments on periodic waves and solitons, IMA J. Appl. Math., 32 (1984), 353–366. https://doi.org/10.1093/imamat/32.1-3.353 doi: 10.1093/imamat/32.1-3.353

This article has been cited by:

1.	Muhammad Ahsan-ul-Haq, On Poisson Moment Exponential Distribution with Applications, 2022, 2198-5804, 10.1007/s40745-022-00400-0
2.	Mohamed S. Eliwa, Muhammad Ahsan-ul-Haq, Amani Almohaimeed, Afrah Al-Bossly, Mahmoud El-Morshedy, Barbara Martinucci, Discrete Extension of Poisson Distribution for Overdispersed Count Data: Theory and Applications, 2023, 2023, 2314-4785, 1, 10.1155/2023/2779120
3.	Muhammad Ahsan-ul-Haq, Afrah Al-Bossly, Mahmoud El-Morshedy, Mohamed S. Eliwa, Maciej Lawrynczuk, Poisson XLindley Distribution for Count Data: Statistical and Reliability Properties with Estimation Techniques and Inference, 2022, 2022, 1687-5273, 1, 10.1155/2022/6503670
4.	Abdulaziz S. Alghamdi, Muhammad Ahsan-ul-Haq, Ayesha Babar, Hassan M. Aljohani, Ahmed Z. Afify, The discrete power-Ailamujia distribution: properties, inference, and applications, 2022, 7, 2473-6988, 8344, 10.3934/math.2022465
5.	Elebe Emmanuel Nwezza, Uchenna Ugwunnaya Uwadi, Christian Osagie, Modeling the number of component failures: A Poison-geometric distribution, 2022, 16, 24682276, e01206, 10.1016/j.sciaf.2022.e01206
6.	Muhammed Irshad, Christophe Chesneau, Veena D’cruz, Radhakumari Maya, Discrete Pseudo Lindley Distribution: Properties, Estimation and Application on INAR(1) Process, 2021, 26, 2297-8747, 76, 10.3390/mca26040076
7.	Muhammad Ahsan-ul-Haq, Javeria Zafar, A new one-parameter discrete probability distribution with its neutrosophic extension: mathematical properties and applications, 2023, 2364-415X, 10.1007/s41060-023-00382-z
8.	Ahmed Sedky Eldeeb, Muhammad Ahsan-ul-Haq, Ayesha Babar, A new discrete XLindley distribution: theory, actuarial measures, inference, and applications, 2024, 17, 2364-415X, 323, 10.1007/s41060-023-00395-8
9.	Osama Abdulaziz Alamri, Classical and Bayesian estimation of discrete poisson Agu-Eghwerido distribution with applications, 2024, 109, 11100168, 768, 10.1016/j.aej.2024.09.063
10.	Hassan M. Aljohani, Muhammad Ahsan-ul-Haq, Javeria Zafar, Ehab M. Almetwally, Abdulaziz S. Alghamdi, Eslam Hussam, Abdisalam Hassan Muse, Analysis of Covid-19 data using discrete Marshall–Olkinin Length Biased Exponential: Bayesian and frequentist approach, 2023, 13, 2045-2322, 10.1038/s41598-023-39183-6
11.	Safar M. Alghamdi, Muhammad Ahsan-ul-Haq, Olayan Albalawi, Majdah Mohammed Badr, Eslam Hussam, H.E. Semary, M.A. Abdelkawy, Binomial Poisson Ailamujia model with statistical properties and application, 2024, 17, 16878507, 101096, 10.1016/j.jrras.2024.101096
12.	Amani Alrumayh, Hazar A. Khogeer, A New Two-Parameter Discrete Distribution for Overdispersed and Asymmetric Data: Its Properties, Estimation, Regression Model, and Applications, 2023, 15, 2073-8994, 1289, 10.3390/sym15061289
13.	John Kwadey Okutu, Nana K. Frempong, Simon K. Appiah, Atinuke O. Adebanji, Pritpal Singh, A New Generated Family of Distributions: Statistical Properties and Applications with Real-Life Data, 2023, 2023, 2577-7408, 1, 10.1155/2023/9325679
14.	M. R. Irshad, Muhammed Ahammed, R. Maya, Christophe Chesneau, INAR(1) process with Poisson-transmuted record type exponential innovations, 2024, 19, 15741699, 145, 10.3233/MAS-231458
15.	Khaled M. Alqahtani, Mahmoud El-Morshedy, Hend S. Shahen, Mohamed S. Eliwa, A discrete extension of the Burr-Hatke distribution: Generalized hypergeometric functions, different inference techniques, simulation ranking with modeling and analysis of sustainable count data, 2024, 9, 2473-6988, 9394, 10.3934/math.2024458
16.	Mohamed Eliwa, Mahmoud El-Morshedy, Hend Shahen, Modelling dispersed count data under various shapes of failure rates: A discrete probability analogue of odd Lomax generator, 2023, 37, 0354-5180, 6177, 10.2298/FIL2318177E
17.	M. Ahsan-ul-Haq, M. R. Irshad, E. S. Muhammed Ahammed, R. Maya, New Discrete Bilal Distribution and Associated INAR(1) Process, 2023, 44, 1995-0802, 3647, 10.1134/S1995080223090020
18.	Yingying Qi, Dan Ding, Yusra A. Tashkandy, M.E. Bakr, M.M. Abd El-Raouf, Anoop Kumar, A novel probabilistic model with properties: Its implementation to the vocal music and reliability products, 2024, 107, 11100168, 254, 10.1016/j.aej.2024.07.035
19.	Khlood Al-Harbi, Aisha Fayomi, Hanan Baaqeel, Amany Alsuraihi, A Novel Discrete Linear-Exponential Distribution for Modeling Physical and Medical Data, 2024, 16, 2073-8994, 1123, 10.3390/sym16091123
20.	Thanasate Akkanphudit, The Discrete Gompertz–Weibull–Fréchet Distribution: Properties and Applications, 2023, 44, 1995-0802, 3663, 10.1134/S1995080223090032
21.	Mohamed S. Algolam, Mohamed S. Eliwa, Mohamed El-Dawoody, Mahmoud El-Morshedy, A discrete extension of the Xgamma random variable: mathematical framework, estimation methods, simulation ranking, and applications to radiation biology and industrial engineering data, 2025, 10, 2473-6988, 6069, 10.3934/math.2025277
22.	Howaida Elsayed, Mohamed Hussein, A New Discrete Analogue of the Continuous Muth Distribution for Over-Dispersed Data: Properties, Estimation Techniques, and Application, 2025, 27, 1099-4300, 409, 10.3390/e27040409
23.	Amir Mushtaq, Mohamed Kayid, Ghadah Alomani, A new discrete generalized class of distribution with application to radiation and COVID-19 data, 2025, 18, 16878507, 101485, 10.1016/j.jrras.2025.101485

Reader Comments

Your name:*

Email:*
© 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)