Data modeling played a crucial role in a variety of research domains due to its widespread practical applications, especially when handling complex datasets. This study explored a specific discrete distribution, characterized by a single parameter, developed using the weighted combining discretization method. The statistical properties of this distribution were rigorously derived and expressed mathematically, covering essential aspects such as moments, skewness, kurtosis, covariance, index of dispersion, order statistics, entropies, mean residual life, residual coefficient of variation function, stress-strength models, and premium principles. These properties highlighted the model's suitability for analyzing right-skewed data with heavy tails, making it a powerful tool for probabilistic modeling in situations where data exhibited overdispersion and increasing failure rates. The research introduced a range of estimation techniques, including maximum product of spacings, method of moments, Anderson-Darling, right-tail Anderson-Darling, maximum likelihood, least squares, weighted least squares, Cramer-Von-Mises, and percentile, each explained in detail. A ranking simulation study was performed to assess the performance of these estimators, with ranking techniques used to determine the most effective estimator across various sample sizes. The study further applied the proposed model to real-world datasets, demonstrating its ability to address complex data scenarios and showcasing its superior performance in comparison to traditional models such as the geometric, Poisson, and negative binomial distributions. Overall, the results emphasized the proposed model's potential as a versatile and effective tool for modeling over-dispersed and skewed data, with promising implications for future research in diverse fields.
Citation: Mahmoud El-Morshedy, Mohamed S. Eliwa, Mohamed El-Dawoody, Hend S. Shahen. A weighted hybrid discrete probability model: Mathematical framework, statistical analysis, estimation techniques, simulation-based ranking, and goodness-of-fit evaluation for over-dispersed data[J]. Electronic Research Archive, 2025, 33(4): 2061-2091. doi: 10.3934/era.2025091
[1] | Amparo Pascual-Ahuir . Annual Report 2021. AIMS Bioengineering, 2022, 9(1): 64-67. doi: 10.3934/bioeng.2022006 |
[2] | Xu Guo . Annual Report 2023. AIMS Bioengineering, 2024, 11(1): 18-23. doi: 10.3934/bioeng.2024002 |
[3] | Xu Guo . Annual Report 2022. AIMS Bioengineering, 2023, 10(1): 62-66. doi: 10.3934/bioeng.2023006 |
[4] | Sven Sölmann, Anke Rattenholl, Hannah Blattner, Guido Ehrmann, Frank Gudermann, Dirk Lütkemeyer, Andrea Ehrmann . Mammalian cell adhesion on different 3D printed polymers with varying sterilization methods and acidic treatment. AIMS Bioengineering, 2021, 8(1): 25-35. doi: 10.3934/bioeng.2021004 |
[5] | Fabrizio Belleggia . Hard and soft tissue augmentation of vertical ridge defects with the “hard top double membrane technique”: introduction of a new technique and a case report. AIMS Bioengineering, 2022, 9(1): 26-43. doi: 10.3934/bioeng.2022003 |
[6] | Maria Júlia Bento Martins Parreira, Bruna Trazzi Pagani, Matheus Bento Medeiros Moscatel, Daniela Vieira Buchaim, Carlos Henrique Bertoni Reis, Beatriz Flávia de Moraes Trazzi, Acácio Fuziy, Rogerio Leone Buchaim . Effects of systemic administration of the retinoid Isotretinoin on bone tissue: A narrative literature review. AIMS Bioengineering, 2024, 11(2): 212-240. doi: 10.3934/bioeng.2024012 |
[7] | Ayub Ahmed, Bashdar Salam, Mahmud Mohammad, Ali Akgül, Sarbaz H. A. Khoshnaw . Analysis coronavirus disease (COVID-19) model using numerical approaches and logistic model. AIMS Bioengineering, 2020, 7(3): 130-146. doi: 10.3934/bioeng.2020013 |
[8] | Daria Wehlage, Hannah Blattner, Al Mamun, Ines Kutzli, Elise Diestelhorst, Anke Rattenholl, Frank Gudermann, Dirk Lütkemeyer, Andrea Ehrmann . Cell growth on electrospun nanofiber mats from polyacrylonitrile (PAN) blends. AIMS Bioengineering, 2020, 7(1): 43-54. doi: 10.3934/bioeng.2020004 |
[9] | Mohammad Z. Ansari, Anil K. Nirala . Assessment of Fevicol (adhesive) Drying Process through Dynamic Speckle Techniques. AIMS Bioengineering, 2015, 2(2): 49-59. doi: 10.3934/bioeng.2015.2.49 |
[10] | Qing Song, Yu Mao, Mark Wilkins, Fernando Segato, Rolf Prade . Cellulase immobilization on superparamagnetic nanoparticles for reuse in cellulosic biomass conversion. AIMS Bioengineering, 2016, 3(3): 264-276. doi: 10.3934/bioeng.2016.3.264 |
Data modeling played a crucial role in a variety of research domains due to its widespread practical applications, especially when handling complex datasets. This study explored a specific discrete distribution, characterized by a single parameter, developed using the weighted combining discretization method. The statistical properties of this distribution were rigorously derived and expressed mathematically, covering essential aspects such as moments, skewness, kurtosis, covariance, index of dispersion, order statistics, entropies, mean residual life, residual coefficient of variation function, stress-strength models, and premium principles. These properties highlighted the model's suitability for analyzing right-skewed data with heavy tails, making it a powerful tool for probabilistic modeling in situations where data exhibited overdispersion and increasing failure rates. The research introduced a range of estimation techniques, including maximum product of spacings, method of moments, Anderson-Darling, right-tail Anderson-Darling, maximum likelihood, least squares, weighted least squares, Cramer-Von-Mises, and percentile, each explained in detail. A ranking simulation study was performed to assess the performance of these estimators, with ranking techniques used to determine the most effective estimator across various sample sizes. The study further applied the proposed model to real-world datasets, demonstrating its ability to address complex data scenarios and showcasing its superior performance in comparison to traditional models such as the geometric, Poisson, and negative binomial distributions. Overall, the results emphasized the proposed model's potential as a versatile and effective tool for modeling over-dispersed and skewed data, with promising implications for future research in diverse fields.
A recently developed theory gaining significant interest is time scale T, introduced by Stefan Hilger to bridge the gap between continuous and discrete analysis [1]. In simpler terms, it aims to unify the study of differential equations (governing continuous change) and difference equations (modeling discrete jumps) [2].
The core concept of time scale T involves defining a time domain as any non-empty, closed set of real numbers. The familiar differential and difference equations emerge as special cases when the time scale is the set of all real numbers or integers, respectively.
To gain a comprehensive understanding, it is necessary to review some basic concepts of time scale theory. The forward and backward jump operators σ, ρ:T→T are defined by
σ(ℓ)=inf{s∈T∣s>ℓ} and ρ(ℓ)=sup{s∈T∣s<ℓ}, |
(supplemented by inf∅=supT and sup∅=infT). A point ℓ∈T is called right-scattered, right-dense, left-scattered or left-dense if σ(ℓ)>ℓ, σ(ℓ)=ℓ, ρ(ℓ)<ℓ, ρ(ℓ)=ℓ holds, respectively. The set Tκ is defined to be T if T does not have a left-scattered maximum; otherwise, it is T without this left-scattered maximum. The graininess function μ:T→[0,∞) is defined by μ(ℓ)=σ(ℓ)−ℓ. Hence, the graininess function is constant 0 if T=R, while it is constant ℓ for T=Z. However, a time scale T could have nonconstant graininess. A function h:T→R is said to be rd-continuous and is written h∈Crd(T,R), provided that h is continuous at right dense points and at left dense points in T, left hand limits exist, and are finite. We say that h:T→R is differentiable at ℓ∈T whenever
hΔ:=lims→ℓh(ℓ)−h(s)ℓ−s |
exists when σ(ℓ)=ℓ (here by s→ℓ it is understood that s approaches ℓ in the time scale), and when h is continuous at ℓ and σ(ℓ)>ℓ it is
hΔ:=lims→ℓh(σ(ℓ))−h(ℓ)μ(ℓ). |
The product and quotient rules [3, Theorem 1.20] for the derivative of the product hk and the quotient h/k of two differentiable functions h and k are as follows:
(hk)Δ(ℓ)=hΔ(ℓ)k(ℓ)+h(σ(ℓ))kΔ=h(ℓ)kΔ(ℓ)+hΔ(ℓ)k(σ(ℓ)), | (1.1) |
(hk)Δ(ℓ)=hΔ(ℓ)k(ℓ)−h(ℓ)kΔ(ℓ)k(ℓ)k(σ(ℓ)). | (1.2) |
The chain rule [3, Theorem 1.90] for the derivative of the composite function h∘k of a continuously differentiable function h:R→R and a (delta) differentiable function k:T→R results in
(h∘k)Δ={∫10h′(k+sμkΔ)ds}gΔ. | (1.3) |
A function h:T→R is called rd-continuous provided it is continuous at right-dense points in T and its left-sided limits exist (finite) at left-dense points in T. The set of rd-continuous functions f:T→R is denoted by
Crd=Crd(T)=Crd(T,R). |
The set of functions h:T→R that are differentiable and whose derivative is rd-continuous is denoted by
C1rd=C1rd(T)=C1rd(T,R). |
Finally, if h:T→R is a function, then we define the function hσ:T→R by hσ(ℓ)=h(σ(ℓ)) for all ℓ∈T.
Bohner and Peterson's book [3] provides a comprehensive overview and organization of this new calculus. Beyond these basic cases, numerous other time scales can be defined, leading to a wealth of applications. One such application is the study of population dynamic models, as explored in [4]. To delve deeper into the theory, readers can consult the referenced papers [5,6] and monographs [3,7].
Recent years have seen a surge in research on the oscillation and non-oscillation of solutions to dynamic equations on time scales. For further exploration, readers can refer to the references provided [8,9,13,14,15,16,17,18,19].
The present paper investigates the asymptotic behavior of solutions to the semi-canonical third dynamic equation
(a2(ℓ)(a1(ℓ)xΔ(ℓ))Δ)Δ+p(ℓ)xγ(δ(ℓ))=0, ℓ∈[ℓ0,∞)T, | (1.4) |
where γ is the ratio of positive odd integers.
In this paper, we consider the following conditions:
(ⅰ) a1(ℓ)∈C2rd([ℓ0,∞)T,(0,∞)), a2(ℓ)∈C1rd([ℓ0,∞)T,(0,∞)), p(ℓ)∈C([ℓ0,∞)T,(0,∞)) and Eq (1.4) is in semi-canonical form, i.e.,
∫∞ℓ0Δsa2(s)=∞ and ∫∞ℓ0Δsa1(s)<∞; | (1.5) |
(ⅱ) δ∈C1rd([ℓ0,∞)T), δΔ(ℓ)≥0, and limℓ→∞δ(ℓ)=∞.
Let us recall that a solution of Eq (1.4) is a nontrivial real-valued function x satisfying the equation for ℓ≥ℓx for some ℓx≥ℓx0 such that x∈C1([ℓx,∞)T,R), a1xΔ∈C1([ℓx,∞)T,R), and a2(a1(xΔ))Δ∈C1([ℓx,1)T,R). We exclude solutions that vanish identically in some neighborhood of infinity, assuming that such solutions exist for Eq (1.4). A solution x(ℓ) of Eq (1.4) is termed oscillatory if it exhibits arbitrarily large zeros on [ℓx,∞)T; otherwise, it is classified as non-oscillatory.
The study of oscillatory behavior in Eq (1.4) often hinges on its form. Equation (1.4) is in canonical form if
∫∞ℓ0Δsa1(s)=∫∞ℓ0Δsa2(s)=∞, |
and it is in non-canonical form if
∫∞ℓ0Δsa1(s)<∞ and ∫∞ℓ0Δsa2(s)<∞. |
If either
∫∞ℓ0Δsa1(s)<∞ and ∫∞ℓ0Δsa2(s)=∞, | (S1) |
or
∫∞ℓ0Δsa1(s)=∞ and ∫∞ℓ0Δsa2(s)<∞, | (S2) |
then we will say that (1.4) is in semicanonical form.
The groundwork for studying third-order dynamic equations on general time scales was laid by Erbe et al. [13], who focused on equations of the form
(a2(ℓ)((a1(ℓ)xΔ(ℓ))Δ))Δ+p(ℓ)f(x(ℓ))=0, ℓ∈[ℓ0,∞)T, | (1.6) |
where a1,a2,p∈Crd(ℓ0,∞)T, f∈C(R,R) R is continuous and satisfies uf(u)>0 for u≠0. Additionally, for each k>0, there exists M=Mk>0 such that f(u)/u≥M, |u|≥k. Using the Riccati transformation technique, they established sufficient conditions that guarantee every solution to this equation either oscillates or converges to zero.
Building on Erbe et al.'s work [13], Hassan [14] investigated a more general form of the third-order equation
(a2(ℓ)((a1(ℓ)xΔ(ℓ))Δ)α)Δ+f(ℓ,x(δ(ℓ)))=0, ℓ∈[ℓ0,∞)T, | (1.7) |
where α≥1 and δ(ℓ)≤ℓ, in the canonical form.
In the particular case of T=R and γ=1, Chatzarakis et al. [11] established new oscillation criteria for the differential equation
(a2(ℓ)(a1(ℓ)x′(ℓ))′)′+p(ℓ)x(δ(ℓ))=0, |
in the canonical form. Recently, techniques have been developed to study the oscillatory behavior of solutions to third-order equations. Moaaz et al. [21,22] extended the improved methods used in studying second-order equations [23,24]. The development of oscillation criteria for delay differential equations of odd orders can also be observed through the works [25,26].
Our literature review indicates a scarcity of research on the oscillatory behavior of solutions to Eq (1.4) when it takes the semi-canonical form (S1). This paper tackles Eq (1.4) in its less-studied semi-canonical form. We begin by transforming it into the more common canonical form. This transformation allows us to then establish new criteria for determining when solutions to Eq (1.4) oscillate.
To enhance readability, we'll use the following symbols:
A(ℓ):=∫∞ℓΔsa1(s), a(ℓ):=a1(ℓ)A(ℓ)Aσ(ℓ), r(ℓ):=a2(ℓ)Aσ(ℓ), |
P(ℓ):=p(ℓ)Aγ(ℓ), ϕ(ℓ):=∫ℓℓ1Δsr(s), ψ(ℓ):=∫ℓℓ1ϕ(s)Δsa(s), and z(ℓ):=x(ℓ)A(ℓ). |
Lemma 2.1. [27] Assume that x is an eventually positive solution of (1.4) satisfying (1.5). Then there exists ℓ1∈[ℓ0,∞)T such that x satisfies one of the following three cases:
(I) xΔ>0, (a1(ℓ)(xΔ(ℓ)))Δ>0, (a2(ℓ)((a1(ℓ)(xΔ(ℓ)))Δ)Δ<0;
(II) xΔ<0, (a1(ℓ)(xΔ(ℓ)))Δ>0, (a2(ℓ)((a1(ℓ)(xΔ(ℓ)))Δ)Δ<0;
(III) xΔ<0, (a1(ℓ)(xΔ(ℓ)))Δ<0, (a2(ℓ)((a1(ℓ)(xΔ(ℓ)))Δ)Δ<0.
Theorem 2.1. Assume that
∫∞ℓΔsr(s)=∞. | (2.1) |
Then the semi-canonical dynamic Eq (1.4) has a solution x(ℓ) if and only if the corresponding canonical equation
(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ+P(ℓ)zγ(δ(ℓ))=0, | (2.2) |
admits the solution z(ℓ)=x(ℓ)A(ℓ).
Proof. Referring back to σ(ℓ) as the forward jump operator and performing differentiation yields
a2(ℓ)Aσ(ℓ)(a1(ℓ)A(ℓ)Aσ(ℓ)(x(ℓ)A(ℓ))Δ)Δ=a2(ℓ)Aσ(ℓ){a1(ℓ)A(ℓ)Aσ(ℓ)[xΔ(ℓ)A(t)−x(t)AΔ(t)A(t)Aσ(ℓ)]}Δ=a2(ℓ)Aσ(ℓ){a1(ℓ)xΔ(ℓ)A(ℓ)−a1(ℓ)x(ℓ)AΔ(ℓ)}Δ=a2(ℓ)Aσ(ℓ){(a1(ℓ)xΔ(ℓ))ΔAσ(ℓ)+a1(ℓ)xΔ(ℓ)AΔ(ℓ)+xΔ(ℓ)}=a2(ℓ)(a1(ℓ)xΔ(ℓ))Δ. | (2.3) |
From (2.1), we have
∫∞ℓAσ(s)a2(s)Δs=∞, | (2.4) |
and
∫∞ℓΔsa1(s)A(s)Aσ(s)=∫∞ℓ(1A(s))ΔΔs=limℓ→∞(1A(ℓ)−1A(ℓ0))=∞. | (2.5) |
Combining (2.3) with (1.4), we obtain
(a2(ℓ)(a1(ℓ)xΔ(ℓ))Δ)Δ+p(ℓ)xγ(δ(ℓ))=0(a2(ℓ)Aσ(ℓ)(a1(ℓ)A(ℓ)Aσ(ℓ)(x(ℓ)A(ℓ))Δ)Δ)Δ+p(ℓ)Aγ(ℓ)xγ(δ(ℓ))Aγ(ℓ)=0(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ+P(ℓ)zγ(δ(ℓ))=0. | (2.6) |
It is clear that x(ℓ)A(ℓ) is a solution of (2.6). Moreover, considering (2.4) and (2.5), it is apparent that Eq (2.6) is in canonical form and from [28] this canonical from is unique.
Theorem 2.1 significantly streamlines the analysis of Eq (1.4) by reducing it to the scope of (2.2), thereby directing our focus towards only two classifications of solutions that ultimately exhibit positivity, i.e., either
z(ℓ)>0, a(ℓ)zΔ(ℓ)<0, r(ℓ)(a(ℓ)zΔ(ℓ)>0, (r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ<0, |
and in this case, we denote z∈ℵ0 or
z(ℓ)>0, a(ℓ)zΔ(ℓ)>0, r(ℓ)(a(ℓ)zΔ(ℓ))Δ>0, (r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ<0, |
and for this characteristic, we indicate that z∈ℵ2.
Theorem 2.2. Let γ≥1 and (2.1) hold. Suppose that
lim supℓ→∞{1ϕγ(δ(ℓ))∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)Δs+1ϕγ−1(δ(ℓ))∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))Δs +ϕ(δ(ℓ))∫∞ℓP(s)ψγ(s)Δs}={∞,γ>1,1,γ=1, | (2.7) |
and
∫∞ℓ01a(u)∫∞u1r(v)∫∞vP(s)ΔsΔvΔu=∞. | (2.8) |
Then every non-oscillatory solution z(ℓ) of (1.4) satisfies limℓ→∞x(ℓ)A(ℓ)=0.
Proof. Let x(ℓ) be a non-oscillatory solution of Eq (1.4), where x(ℓ)>0, and x(δ(ℓ))>0 for ℓ≥ℓ1 for some ℓ1≥ℓ0. According to Theorem 2.1, the corresponding function z(ℓ)=x(ℓ)A(ℓ) is a positive solution of (2.2), implying that either z∈ℵ0 or z∈ℵ2 for ℓ≥ℓ1.
Let us examine the case where z∈ℵ2. In this case, we observe that
a(ℓ)zΔ(ℓ)≥∫ℓℓ1r−1(s)r(s)(a(s)zΔ(s))ΔΔs≥r(ℓ)(a(ℓ)zΔ(ℓ))Δ∫ℓℓ1Δsr(s)≥r(ℓ)(a(ℓ)zΔ(ℓ))Δϕ(ℓ). |
Hence,
(a(ℓ)zΔ(ℓ)ϕ(ℓ))Δ=ϕ(ℓ)(a(ℓ)zΔ(ℓ))Δ−(a(ℓ)zΔ(ℓ))ϕΔ(ℓ)ϕ(ℓ)ϕσ(ℓ)=ϕ(ℓ)r(ℓ)(a(ℓ)zΔ(ℓ))Δ−(a(ℓ)zΔ(ℓ))r(ℓ)ϕ(ℓ)ϕσ(ℓ)≤0. | (2.9) |
Consequently, it can be inferred from (2.9) that
z(ℓ)≥∫ℓℓ1zΔ(s)Δs=∫ℓℓ1a(s)zΔ(s)ϕ(s)ϕ(s)a(s)Δs≥a(ℓ)zΔ(ℓ)ϕ(ℓ)ψ(ℓ). | (2.10) |
Combining (2.10) with (2.2), we see that a(ℓ)zΔ(ℓ)ϕ(ℓ)ψ(ℓ) is a positive solution to the dynamic inequality
(r(ℓ)χΔ(ℓ))Δ+P(ℓ)ψγ(ℓ)ϕγ(δ(ℓ))χγ(δ(ℓ))≤0, | (2.11) |
where χ(ℓ):=a(ℓ)zΔ(ℓ). Integration (2.11) from ℓ to ∞ and considering the nonincreasing nature of χ(ℓ)/ϕ(ℓ), we obtain
χΔ(ℓ)≥1r(ℓ)∫∞ℓP(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs. |
Therefore,
χ(ℓ)≥∫ℓℓ11r(s)∫∞sP(u)ψγ(u)ϕγ(δ(u))χγ(δ(u))ΔuΔs=∫ℓℓ11r(s)∫ℓsP(u)ψγ(u)ϕγ(δ(u))χγ(δ(u))ΔuΔs+∫ℓℓ11r(s)∫∞ℓP(u)ψγ(u)ϕγ(δ(u))χγ(δ(u))ΔuΔs. | (2.12) |
Integrating by parts, we obtain
χ(ℓ)≥∫ℓℓ1ϕ(σ(s))P(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs+ϕ(ℓ)∫∞ℓP(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs. | (2.13) |
It follows that
χ(δ(ℓ))≥∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs+ϕ(δ(ℓ))∫∞δ(ℓ)P(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs=∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs+ϕ(δ(ℓ))∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs +ϕ(δ(ℓ))∫∞ℓP(s)ψγ(s)ϕγ(δ(s))χγ(δ(s))Δs. | (2.14) |
Utilizing the monotonicity characteristics of χ(ℓ) and χ(ℓ)/ϕ(ℓ), we have χ(δ(ℓ))≤χ(δ(s)) and χ(δ(s))ϕ(δ(s))≥χ(δ(ℓ))ϕ(δ(ℓ)) for s≥ℓ, hence (2.13) takes the form
χ(δ(ℓ))≥χγ(δ(ℓ))ϕγ(δ(ℓ))∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)Δs+χγ(δ(ℓ))ϕγ−1(δ(ℓ))∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))Δs +ϕ(δ(ℓ))χγ(δ(ℓ))∫∞ℓP(s)ψγ(s)Δs, | (2.15) |
χ1−γ(δ(ℓ))≥1ϕγ(δ(ℓ))∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)Δs+1ϕγ−1(δ(ℓ))∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))Δs +ϕ(δ(ℓ))∫∞ℓP(s)ψγ(s)Δs. | (2.16) |
This contradicts (2.7). Subsequently, let us assume that z∈ℵ0. Then limℓ→∞z(ℓ)=k≥0, and we propose that k=0. If not, it would imply z(ℓ)≥k>0. Integrating (2.2) from ℓ to ∞ yields
r(ℓ)(a(ℓ)zΔ(ℓ))Δ≥∫∞ℓP(s)zγ(δ(s))Δs≥kγ∫∞ℓP(s)Δs. |
Therefore,
−a(ℓ)zΔ(ℓ)≥kγ∫∞ℓ1r(u)∫∞uP(s)ΔsΔu, |
and
z(ℓ1)≥kγ∫∞ℓ11a(u)∫∞u1r(v)∫∞vP(s)ΔsΔvΔu. |
This leads to a contradiction to (2.8). Thus, we conclude: limℓ→∞z(ℓ)=limℓ→∞x(ℓ)A(ℓ)=0, and, the proof of the theorem is complete.
Theorem 2.3. Let 0<γ<1 and (2.1) hold. If (2.8) and
lim supℓ→∞{1ϕ(δ(ℓ))∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)Δs+∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))Δs +ϕγ(δ(ℓ))∫∞ℓP(s)ψγ(s)Δs}=∞ | (2.17) |
hold, then every non-oscillatory solutionz(ℓ) of (1.4) satisfies limℓ→∞x(ℓ)A(ℓ)=0.
Proof. Let x(ℓ) be a non-oscillatory solution of Eq (1.4), where x(ℓ)>0, and x(δ(ℓ))>0 for ℓ≥ℓ1 for some ℓ1≥ℓ0. According to Theorem 2.1, the corresponding function z(ℓ)=x(ℓ)A(ℓ) is a positive solution of (2.2), implying that either z∈ℵ0 or z∈ℵ2 for ℓ≥ℓ1.
First, let us assume that z∈ℵ2. Proceeding similarly to the proof of Theorem 2.7, we arrive at (2.15). Dividing (2.16) by ϕ1−γ(δ(ℓ), we obtain
(χ(δ(ℓ))ϕ(δ(ℓ)))1−γ≥1ϕ(δ(ℓ))∫δ(ℓ)ℓ1ϕ(σ(s))P(s)ψγ(s)Δs+∫ℓδ(ℓ)P(s)ψγ(s)ϕγ(δ(s))Δs +ϕγ(δ(ℓ))∫∞ℓP(s)ψγ(s)Δs. | (2.18) |
In view of the decreasing nature of χ(δ(ℓ))/ϕ(δ(ℓ)) and the fact that 0<γ<1, there exists a constant C>0 such that
(χ(δ(ℓ))ϕ(δ(ℓ)))1−γ≤C. |
Taking the lim sup as ℓ→∞, we establish a contradiction to (2.18), and consequently, z∉ℵ2.
Subsequently, let us assume that z∈ℵ0. Proceeding similarly to the proof of Theorem 2.7, it becomes evident that condition (2.8) once more leads to the conclusion that limℓ→∞x(ℓ)A(ℓ)=0. This completes the proof.
Theorem 2.4. Suppose that conditions (i), (ii), and δΔ(ℓ)>0 are satisfied on [ℓ0,∞)T, γ≤1, and there exists a function ξ(ℓ) such that
ξΔ(ℓ)≥0, ξ(ℓ)>ℓ ,and θ(ℓ)=δ(ξ(ξ(ℓ)))<ℓ. | (2.19) |
If
lim infℓ→∞∫ℓδ(ℓ)P(s)ψγ(δ(s))Δs{=∞,γ<1,>1/e,γ=1, | (2.20) |
and
lim infℓ→∞∫ℓθ(ℓ)(1a(s)∫ξ(ℓ)s1r(u)∫ξ(u)uP(v)ΔvΔu)Δs{=∞,γ<1,>1/e,γ=1, | (2.21) |
for all ℓ1≥ℓ0, then Eq (1.4) is oscillatory.
Proof. Let x(ℓ) be a non-oscillatory solution of Eq (1.4), where x(ℓ)>0, and x(δ(ℓ))>0 for ℓ≥ℓ1 for some ℓ1≥ℓ0. According to Theorem 2.1, the corresponding function z(ℓ)=x(ℓ)A(ℓ) is a positive solution of (2.2), implying that either z∈ℵ0 or z∈ℵ2 for ℓ≥ℓ1. Assuming that z(ℓ)∈ℵ2, we have
a(ℓ)zΔ(ℓ)≥∫ℓℓ1r−1(s)r(s)(a(s)zΔ(s))ΔΔs≥r(ℓ)(a(ℓ)zΔ(ℓ))Δϕ(ℓ). |
It follows that
zΔ(ℓ)≥r(ℓ)(a(ℓ)zΔ(ℓ))Δϕ(ℓ)a(ℓ). | (2.22) |
Integrating the above inequality from ℓ2 to ℓ, we obtain
z(ℓ)≥∫ℓℓ2r(s)(a(s)zΔ(s))Δϕ(s)a(s)Δs≥r(ℓ)(a(ℓ)zΔ(ℓ))Δ∫ℓℓ2ϕ(s)a(s)Δs=r(ℓ)(a(ℓ)zΔ(ℓ))Δψ(ℓ). | (2.23) |
There exists ℓ3≥ℓ2 such that δ(ℓ)≥ℓ2 for all ℓ≥ℓ3. Then, we have
z(δ(ℓ))≥r(δ(ℓ))(a(δ(ℓ))zΔ(δ(ℓ)))Δψ(δ(ℓ)), for all ℓ≥ℓ3. |
Combining this with (2.2) yields
YΔ(ℓ)+P(ℓ)ψγ(δ(ℓ))Yγ(δ(ℓ))≤0, for ℓ≥ℓ3, | (2.24) |
where Y(ℓ):=r(ℓ)(a(ℓ)zΔ(ℓ))Δ. Integrating (2.24) from δ(ℓ) to ℓ, we have
Y(δ(ℓ))≥Y(δ(ℓ))−Y(ℓ)≥Yγ(δ(ℓ))∫ℓδ(ℓ)P(s)ψγ(δ(s))Δs. | (2.25) |
Hence,
Y1−γ(δ(ℓ))≥∫ℓδ(ℓ)P(s)ψγ(δ(s))Δs for ℓ≥ℓ3. |
According to [29, Theorem 1], we reach the intended contradiction.
Now, consider z∈ℵ0. Integrating (2.2) from ℓ to ξ(ℓ), we obtain
r(ℓ)(a(ℓ)zΔ(ℓ))Δ≥∫ξ(ℓ)ℓP(s)zγ(δ(s))Δs≥zγ(δ(ξ(ℓ)))∫ξ(ℓ)ℓP(s)Δs, |
where θ(ℓ):=δ(ξ(ξ(ℓ))). Consequently,
(a(ℓ)zΔ(ℓ))Δ≥zγ(δ(ξ(ℓ)))r(ℓ)∫ξ(ℓ)ℓP(s)Δs. | (2.26) |
Integrating (2.26) from ℓ to ξ(ℓ), we have
−a(ℓ)zΔ(ℓ)∫ξ(ℓ)ℓzγ(δ(ξ(s)))r(s)∫ξ(s)sP(u)ΔuΔs≥zγ(δ(ξ(ξ(ℓ))))∫ξ(ℓ)ℓ1r(s)∫ξ(s)sP(u)ΔuΔs=zγ(θ(ℓ))∫ξ(ℓ)ℓ1r(s)∫ξ(s)sP(u)ΔuΔs. | (2.27) |
It follows that
zΔ(ℓ)+(1a(ℓ)∫ξ(ℓ)ℓ1r(s)∫ξ(s)sP(u)ΔuΔs)zγ(θ(ℓ))≤0. | (2.28) |
The remainder of the proof follows a similar pattern to the one described above and is therefore omitted.
Theorem 2.5. Let (2.1) hold. Assume that there exists a function ρ(ℓ)∈C1rd(T,R+), such that
lim supℓ→∞∫ℓℓ0(P(s)ρ(s)ψ(δ(s))ϕ(s)λγ−1−ρΔ(s)r(s)4ρ(s))Δs=∞, | (2.29) |
and (2.8) hold. Then every solution z(ℓ) of (1.4) is oscillatory or satisfies limℓ→∞x(ℓ)A(ℓ)=0.
Proof. Let x(ℓ) be a non-oscillatory solution of Eq (1.4), where x(ℓ)>0, and x(δ(ℓ))>0 for ℓ≥ℓ1 for some ℓ1≥ℓ0. According to Theorem 2.1, the corresponding function z(ℓ)=x(ℓ)A(ℓ) is a positive solution of (2.2), implying that either z∈ℵ0 or z∈ℵ2 for ℓ≥ℓ1.
Firstly, let us consider z∈ℵ2; then we have r(ℓ)(a(ℓ)zΔ(ℓ))Δ is decreasing, and moreover,
r(ℓ)(a(ℓ)zΔ(ℓ))Δ≥∫∞ℓP(s)zγ(δ(s))Δs≥zγ(ℓ(s))∫∞ℓP(s)Δs. | (2.30) |
Let us define the generalized Riccati substitution
ω(ℓ)=ρ(ℓ)r(ℓ)(a(ℓ)zΔ(ℓ))Δa(ℓ)zΔ(ℓ). | (2.31) |
Applying both the product rule and the quotient rule, we obtain
ωΔ(ℓ)=(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ(ρ(ℓ)a(ℓ)zΔ(ℓ))+(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)σ(ρ(ℓ)a(ℓ)zΔ(ℓ))Δ=(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)Δ(ρ(ℓ)a(ℓ)zΔ(ℓ)) +(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)σ((a(ℓ)zΔ(ℓ))ρΔ(ℓ)−ρ(ℓ)(a(ℓ)zΔ(ℓ))Δ(a(ℓ)zΔ(ℓ))(a(ℓ)zΔ(ℓ))σ)≤−P(ℓ)ρ(ℓ)(zγ(δ(ℓ))a(ℓ)zΔ(ℓ))+ρΔ+(ℓ)ρ(σ(ℓ))ω(σ(ℓ)) −ρ(ℓ)(r(ℓ)(a(ℓ)zΔ(ℓ))Δ)σ(a(ℓ)zΔ(ℓ))Δ(a(ℓ)zΔ(ℓ))(a(ℓ)zΔ(ℓ))σ. | (2.32) |
Using the monotonicity of r(ℓ)(a(ℓ)rΔ(ℓ))Δ and (a(ℓ)rΔ(ℓ)), we have
(a(ℓ)zΔ(ℓ))Δ≥r(σ(ℓ))r(ℓ)(a(σ(ℓ))zΔ(σ(ℓ)))Δ, | (2.33) |
and
1a(ℓ)zΔ(ℓ)≥1a(σ(ℓ))zΔ(σ(ℓ)). | (2.34) |
Combining (2.33) and (2.34) with (2.32), we obtain
ωΔ(ℓ)≤−P(ℓ)ρ(ℓ)(zγ(δ(ℓ))a(ℓ)zΔ(ℓ))+ρΔ+(ℓ)ρ(σ(ℓ))ω(σ(ℓ))−ρ(ℓ)r(ℓ)ρ2(σ((ℓ))ω2(σ(ℓ)). | (2.35) |
From (2.9), (2.10), and the fact that δ(ℓ)≤ℓ, we have
z(δ(ℓ))a(ℓ)zΔ(ℓ)≥ψ(δ(ℓ))ϕ(ℓ), for ℓ≥ℓ3. | (2.36) |
This, together with (2.33), leads to
ωΔ(ℓ)≤−P(ℓ)ρ(ℓ)ψ(δ(ℓ))ϕ(ℓ)zγ−1(δ(ℓ))+ρΔ+(ℓ)ρ(σ(ℓ))ω(σ(ℓ))−ρ(ℓ)r(ℓ)ρ2(σ((ℓ))ω2(σ(ℓ)). | (2.37) |
Since zΔ(ℓ)>0, then there exists a constant λ>0 such that z(ℓ)≥λ for ℓ≥ℓ3. Consequently, (2.37) can be expressed as
ωΔ(ℓ)≤−P(ℓ)ρ(ℓ)ψ(δ(ℓ))ϕ(ℓ)λγ−1+ρΔ+(ℓ)ρ(σ(ℓ))ω(σ(ℓ))−ρ(ℓ)r(ℓ)ρ2(σ((ℓ))ω2(σ(ℓ))≤−P(ℓ)ρ(ℓ)ψ(δ(ℓ))ϕ(ℓ)λγ−1+ρΔ(ℓ)r(ℓ)4ρ(ℓ). | (2.38) |
Integrating both sides of (2.38) from ℓ4>ℓ3 to ℓ, we obtain
∫ℓℓ4(P(s)ρ(s)ψ(δ(s))ϕ(s)λγ−1−ρΔ(s)r(s)4ρ(s))Δs≤ω(ℓ4), | (2.39) |
which contradicts (2.29). Now, assume that z∈ℵ0. Proceeding similarly to the proof of Theorem 2.7, it becomes evident that condition (2.8) once more leads to the conclusion that limℓ→∞x(ℓ)A(ℓ)=0. This completes the proof.
Theorem 2.6. Let γ=1,
lim supℓ→∞∫ℓδ(ℓ)(1a(v)∫ℓv1r(u)∫ℓuP(s)ΔsΔu)Δv>1, | (2.40) |
and assume that there exists a function ρ(ℓ)∈C1rd(T,R+), such that
lim supℓ→∞∫ℓℓ0(P(s)ρ(s)ψ(δ(s))ϕ(s)−ρΔ(s)r(s)4ρ(s))Δs=∞. | (2.41) |
Then every solution of (1.4) is oscillatory.
Proof. Let x(ℓ) be a non-oscillatory solution of Eq (1.4), where x(ℓ)>0, and x(δ(ℓ))>0 for ℓ≥ℓ1 for some ℓ1≥ℓ0. According to Theorem 2.1, the corresponding function z(ℓ)=x(ℓ)A(ℓ) is a positive solution of (2.2), implying that either z∈ℵ0 or z∈ℵ2 for ℓ≥ℓ1. Assume z(ℓ)∈ℵ0. Integrating (2.2) from ν to ℓ yields
r(ν)(a(ν)zΔ(ν))Δ≥∫ℓνP(s)z(δ(s))Δs≥z(δ(ℓ))∫ℓνP(s)Δs. |
Integrating again twice from ν to ℓ, we obtain
z(ν)≥z(δ(ℓ))∫ℓν(1a(v)∫ℓv1r(u)∫ℓuP(s)ΔsΔu)Δv. |
Replacing ν with δ(ℓ) leads to contradiction to (2.40). Hence, every positive solution z(ℓ) does not satisfy ℵ0. Therefore, if (2.40) holds, then z(ℓ)∈ℵ2. Proceeding as in Theorem (2.29) with γ=1, completes the proof.
Example 2.1. Consider the third order linear differential equation
(1ℓ(ℓ2(x′(ℓ)))′)′+p0√ℓx(αℓ)=0, ℓ≥1, | (2.42) |
where p0 is a constant and α∈(0,1). Here a2(ℓ)=1ℓ, a1(ℓ)=ℓ2, p(ℓ)=p0√ℓ and δ(ℓ)=αℓ. It is clear that (2.42) is semi-canonical. Since A(ℓ)=1ℓ, a(ℓ)=r(ℓ)=1, and P(ℓ)=p0ℓ3/2, the corresponding canonical equation is
z′′′(ℓ)+p0ℓ3/2x(αℓ)=0. | (2.43) |
It is clear that (2.1) holds. Applying Theorem 2.6, we have
lim supℓ→∞∫ℓδ(ℓ)(1a(v)∫ℓv1r(u)∫ℓuP(s)ΔsΔu)Δv=lim supℓ→∞∫ℓαℓ(∫ℓv∫ℓup0s3/2dsdu)dv=limℓ→∞(−α2+(1−2α)√α+4)ℓ3/2>1 |
and by choosing ρ(ℓ)=ℓ
lim supℓ→∞∫ℓℓ0(P(s)ρ(s)ψ(δ(s))ϕ(s)−ρΔ(s)r(s)4ρ(s))Δs=lim supℓ→∞∫ℓℓ0(p0s3/2α2s3s−14s)ds=∞. |
It follows that (2.42) is oscillatory. Also, by Theorems 2.2 and 2.4, Eq (2.42) is oscillatory or limℓ→∞x(ℓ)A(ℓ)=0.
Remark 2.1. It is worth noting that the existing results in [10,30,31,32] cannot be directly applied to Eq (2.42) due to the fact that a1(ℓ)≠1.
Example 2.2. Consider the second order difference equation
Δ(1ℓ+1Δ(ℓ(ℓ+1))Δx(ℓ))+p0x1/2(ℓ−2)=0, ℓ≥1, | (2.44) |
where p0 is a constant. Here a2(ℓ)=1ℓ+1, a1(ℓ)=ℓ(ℓ+1), p(ℓ)=p0 and δ(ℓ)=ℓ−2. It is clear that (2.44) is semi-canonical. Since A(ℓ)=1ℓ, a(ℓ)=r(ℓ)=1, and P(ℓ)=p0ℓ1/2, the corresponding canonical equation is
Δ(Δ(Δ(z(ℓ))))+p0ℓ1/2x(ℓ−2)=0. | (2.45) |
It is clear that (2.1) and (2.8) hold. Further, (2.17) becomes
lim supℓ→∞{1ℓ−3ℓ−2∑1(s+1)p0s1/2s2+ℓ∑ℓ−2sp0√2s(s−2) +√(ℓ−2)∞∑ℓp0s√2s}=∞. |
Hence, by Theorem 2.3, every solution is oscillatory or limℓ→∞x(ℓ)A(ℓ)=0.
The results of this study are presented in a novel and generalizable framework, highlighting their broad applicability. Our approach involves a unique transformation that converts the equation from the semi-canonical form to the more tractable canonical form. This transformation facilitates the derivation of new oscillation criteria with fewer restrictions compared to the existing literature. Theorems 2.19 and 2.40 illustrate our criteria, ensuring that all solutions oscillate. The results obtained are consistent with the results in [11,13,14] and can be extended to non linear difference equations. Our approach has the potential to be extended to both non-canonical and semi-canonical forms (as defined in (S2)), potentially leading to new oscillation conditions.
Ahmed M. Hassan: Writing-original draft, Writing-review and editing, Making major revisions; Clemente Cesarano: Supervision, Writing-review and editing; Sameh S. Askar: Formal analysis, Writing-original draft; Ahmad M. Alshamrani: Writing-original draft, Making major revisions. All authors have read and approved the final version of the manuscript for publication.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors present their appreciation to King Saud University for funding this research through Researchers Supporting Project number (RSPD2024R533), King Saud University, Riyadh, Saudi Arabia.
The authors declare that there is no conflict of interest regarding the publication of this paper.
[1] |
T. Nakagawa, S. Osaki, The discrete Weibull distribution, IEEE Trans. Reliab., 24 (1975), 300–301. https://doi.org/10.1109/TR.1975.5214915 doi: 10.1109/TR.1975.5214915
![]() |
[2] |
K. B. Kulasekera, D. W. Tonkyn, A new discrete distribution, with applications to survival, dispersal and dispersion, Commun. Stat.- Simul. Comput., 21 (1992), 499–518. https://doi.org/10.1080/03610919208813032 doi: 10.1080/03610919208813032
![]() |
[3] |
H. Sato, M. Ikota, A. Sugimoto, H. Masuda, A new defect distribution metrology with a consistent discrete exponential formula and its applications, IEEE Trans. Semicond. Manuf., 12 (1999), 409–418. https://doi.org/10.1109/66.806118 doi: 10.1109/66.806118
![]() |
[4] | J. D. Smith, A review of Finn, Fischer, and Handler (Eds.), collaborative/therapeutic assessment: A casebook and guide, JPA, 95 (2012), 234–235. https://doi.org/10.1080/00223891.2012.730086 |
[5] |
M. Roederer, A. Treister, W. Moore, L. A. Herzenberg, Probability binning comparison: A metric for quantitating univariate distribution differences, Cytometry, 45 (2001), 37–46. https://doi.org/10.1002/1097-0320(20010901)45:1<37::AID-CYTO1142>3.0.CO;2-E doi: 10.1002/1097-0320(20010901)45:1<37::AID-CYTO1142>3.0.CO;2-E
![]() |
[6] |
A. Barbiero, A. Hitaj, Discrete approximations of continuous probability distributions obtained by minimizing Cramer-von Mises-type distances, Stat. Papers, 64 (2023), 1669–1697. https://doi.org/10.1007/s00362-022-01356-2 doi: 10.1007/s00362-022-01356-2
![]() |
[7] | T. Ghosh, D. Roy, N. K. Chandra, Reliability approximation through the discretization of random variables using reversed hazard rate function, Int. J. Math. Comput. Stat. Nat. Phys. Eng., 7 (2013), 96–100. |
[8] |
S. Chakraborty, Generating discrete analogues of continuous probability distributions-A survey of methods and constructions, J. Stat. Distrib. Appl., 2 (2015), 6. https://doi.org/10.1186/s40488-015-0028-6 doi: 10.1186/s40488-015-0028-6
![]() |
[9] | S. Kotsiantis, D. Kanellopoulos, Discretization techniques: A recent survey, GESTS Int. Trans. Comput. Sci. Eng., 32 (2006), 47–58. |
[10] | G. Casella, R. L. Berger, Statistical Inference Vol. 70, Duxbury Press, 1990. Available from: https://philpapers.org/rec/CASSIV. |
[11] | C. M. Bishop, Pattern Recognition and Machine Learning, Springer, 2006. Available from: https://link.springer.com/book/9780387310732#bibliographic-information. |
[12] | A. Gelman, J. B. Carlin, H. S. Stern, D. B. Dunson, A. Vehtari, D. B. Rubin, Bayesian Data Analysis, CRC Press, 2013. https://doi.org/10.1201/b16018 |
[13] | R. E. Barlow, Statistical theory of reliability and life testing, 1975. https://cir.nii.ac.jp/crid/1571980074720917504 |
[14] | C. D. Lai, M. Xie, Stochastic Ageing and Dependence for Reliability, Springer, 2006. https://doi.org/10.1007/0-387-34232-X |
[15] |
B. Singh, R. P. Singh, A. S. Nayal, A. Tyagi, Discrete inverted Nadarajah-Haghighi distribution: Properties and classical estimation with application to complete and censored data, Stat. Optim. Inf. Comput., 10 (2022), 1293–1313. https://doi.org/10.19139/soic-2310-5070-1365 doi: 10.19139/soic-2310-5070-1365
![]() |
[16] |
D. Roy, Discrete rayleigh distribution, IEEE Trans. Reliab., 53 (2004), 255–260. https://doi.org/10.1109/TR.2004.829161 doi: 10.1109/TR.2004.829161
![]() |
[17] | T. Hussain, M. Ahmad, Discrete inverse Rayleigh distribution, Pak. J. Stat., 30 (2014). |
[18] | S. D. Poisson, Probabilité des Jugements en Matière Criminelle et en Matière Civile, Précédées des Règles Générales du Calcul des Probabilités, Paris, France: Bachelier, 1837. |
[19] |
M. El-Morshedy, M. S. Eliwa, E. Altun, Discrete Burr-Hatke distribution with properties, estimation methods and regression model, IEEE Access, 8 (2020), 74359–74370. https://doi.org/10.1109/ACCESS.2020.2988431 doi: 10.1109/ACCESS.2020.2988431
![]() |
[20] |
H. Krishna, P. S. Pundir, Discrete Burr and discrete Pareto distributions, Stat. Methodol., 6 (2009), 177–188. https://doi.org/10.1016/j.stamet.2008.07.001 doi: 10.1016/j.stamet.2008.07.001
![]() |
[21] | D. J. Hand, F. Daly, K. J. McConway, A. D. Lunn, E. O. Ostrowski, A Hand Book of Small Data Sets, Chapman and Hall/CRC, 1993. https://doi.org/10.1201/9780429246579 |
[22] | J. F. Lawless, Statistical Models and Methods for Lifetime Data, John Wiley & Sons, 2011. |
[23] |
P. Damien, S. Walker, A Bayesian non-parametric comparison of two treatments, Scand. J. Stat., 29 (2002), 51–56. https://doi.org/10.1111/1467-9469.00891 doi: 10.1111/1467-9469.00891
![]() |