Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

A weighted hybrid discrete probability model: Mathematical framework, statistical analysis, estimation techniques, simulation-based ranking, and goodness-of-fit evaluation for over-dispersed data

  • 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

    Related Papers:

    [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 σ, ρ:TT are defined by

    σ()=inf{sTs>} and ρ()=sup{sTs<},

    (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:TR is said to be rd-continuous and is written hCrd(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:TR is differentiable at T whenever

    hΔ:=limsh()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Δ:=limsh(σ())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 hk of a continuously differentiable function h:RR and a (delta) differentiable function k:TR results in

    (hk)Δ={10h(k+sμkΔ)ds}gΔ. (1.3)

    A function h:TR 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:TR is denoted by

    Crd=Crd(T)=Crd(T,R).

    The set of functions h:TR that are differentiable and whose derivative is rd-continuous is denoted by

    C1rd=C1rd(T)=C1rd(T,R).

    Finally, if h:TR is a function, then we define the function hσ:TR 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 xx0 such that xC1([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,pCrd(0,)T, fC(R,R) R is continuous and satisfies uf(u)>0 for u0. Additionally, for each k>0, there exists M=Mk>0 such that f(u)/uM, |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 z0 or

    z()>0,  a()zΔ()>0,  r()(a()zΔ())Δ>0, (r()(a()zΔ())Δ)Δ<0,

    and for this characteristic, we indicate that z2.

    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 limx()A()=0.

    Proof. Let x() be a non-oscillatory solution of Eq (1.4), where x()>0, and x(δ())>0 for 1 for some 10. According to Theorem 2.1, the corresponding function z()=x()A() is a positive solution of (2.2), implying that either z0 or z2 for 1.

    Let us examine the case where z2. In this case, we observe that

    a()zΔ()1r1(s)r(s)(a(s)zΔ(s))ΔΔsr()(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)Δsa()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 z0. Then limz()=k0, 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))Δskγ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: limz()=limx()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 limx()A()=0.

    Proof. Let x() be a non-oscillatory solution of Eq (1.4), where x()>0, and x(δ())>0 for 1 for some 10. According to Theorem 2.1, the corresponding function z()=x()A() is a positive solution of (2.2), implying that either z0 or z2 for 1.

    First, let us assume that z2. 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, z2.

    Subsequently, let us assume that z0. Proceeding similarly to the proof of Theorem 2.7, it becomes evident that condition (2.8) once more leads to the conclusion that limx()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 10, 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 10. According to Theorem 2.1, the corresponding function z()=x()A() is a positive solution of (2.2), implying that either z0 or z2 for 1. Assuming that z()2, we have

    a()zΔ()1r1(s)r(s)(a(s)zΔ(s))ΔΔsr()(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)Δsr()(a()zΔ())Δ2ϕ(s)a(s)Δs=r()(a()zΔ())Δψ(). (2.23)

    There exists 32 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 z0. Integrating (2.2) from to ξ(), we obtain

    r()(a()zΔ())Δξ()P(s)zγ(δ(s))Δszγ(δ(ξ()))ξ()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Δszγ(δ(ξ(ξ())))ξ()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 sup0(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 limx()A()=0.

    Proof. Let x() be a non-oscillatory solution of Eq (1.4), where x()>0, and x(δ())>0 for 1 for some 10. According to Theorem 2.1, the corresponding function z()=x()A() is a positive solution of (2.2), implying that either z0 or z2 for 1.

    Firstly, let us consider z2; then we have r()(a()zΔ())Δ is decreasing, and moreover,

    r()(a()zΔ())ΔP(s)zγ(δ(s))Δszγ((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 z0. Proceeding similarly to the proof of Theorem 2.7, it becomes evident that condition (2.8) once more leads to the conclusion that limx()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 sup0(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 10. According to Theorem 2.1, the corresponding function z()=x()A() is a positive solution of (2.2), implying that either z0 or z2 for 1. Assume z()0. Integrating (2.2) from ν to yields

    r(ν)(a(ν)zΔ(ν))ΔνP(s)z(δ(s))Δsz(δ())ν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())))+p0x(α)=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()=p03/2, the corresponding canonical equation is

    z()+p03/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α(vup0s3/2dsdu)dv=lim(α2+(12α)α+4)3/2>1

    and by choosing ρ()=

    lim sup0(P(s)ρ(s)ψ(δ(s))ϕ(s)ρΔ(s)r(s)4ρ(s))Δs=lim sup0(p0s3/2α2s3s14s)ds=.

    It follows that (2.42) is oscillatory. Also, by Theorems 2.2 and 2.4, Eq (2.42) is oscillatory or limx()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()=p01/2, the corresponding canonical equation is

    Δ(Δ(Δ(z())))+p01/2x(2)=0. (2.45)

    It is clear that (2.1) and (2.8) hold. Further, (2.17) becomes

    lim sup{1321(s+1)p0s1/2s2+2sp02s(s2)   +(2)p0s2s}=.

    Hence, by Theorem 2.3, every solution is oscillatory or limx()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
  • Reader Comments
  • © 2025 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)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(457) PDF downloads(41) Cited by(0)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog