Research article Special Issues

Measuring digital economy in China

  • Received: 17 July 2022 Revised: 21 August 2022 Accepted: 23 August 2022 Published: 31 August 2022
  • JEL Codes: G20, G10, O16

  • The COVID-19 pandemic has highlighted the importance of the digital economy in restoring economic and social development, creating more jobs and improving people's well-being. To inform policy makers about changes to digital strategies, measuring the digital economy is a prerequisite. This study aimed to compile an index of digital economy at the provincial (municipalities, autonomous regions, collectively referred to as "provinces") level to present an accurate and in-depth depiction of how it has developed in China. Our sample covers 31 provinces in China, over the period 2010–2020. This paper firstly constructs the digital economy index system from the four dimensions of digital users, digital platforms, digital industries and digital innovation, and then adopts a combination of entropy weighting method and grey target theory to measure the digital economy index. This paper study revealed that China's digital economy has been on an upward trend from 2010 to 2019 and has a decline in 2020, and the digital innovation is an important driving force for the growth of the digital economy index. The convergence of China's digital economy is decreasing, indicating that the gap in digital economy development between provinces is increasing. The proposed index in this study can be used as a screening tool, decision making tool, benchmarking tool and guidance of high-quality digital economy development.

    Citation: Yanting Xu, Tinghui Li. Measuring digital economy in China[J]. National Accounting Review, 2022, 4(3): 251-272. doi: 10.3934/NAR.2022015

    Related Papers:

    [1] Kawkab Al Amri, Qamar J. A Khan, David Greenhalgh . Combined impact of fear and Allee effect in predator-prey interaction models on their growth. Mathematical Biosciences and Engineering, 2024, 21(10): 7211-7252. doi: 10.3934/mbe.2024319
    [2] Saheb Pal, Nikhil Pal, Sudip Samanta, Joydev Chattopadhyay . Fear effect in prey and hunting cooperation among predators in a Leslie-Gower model. Mathematical Biosciences and Engineering, 2019, 16(5): 5146-5179. doi: 10.3934/mbe.2019258
    [3] Dirk Stiefs, Ezio Venturino, Ulrike Feudel . Evidence of chaos in eco-epidemic models. Mathematical Biosciences and Engineering, 2009, 6(4): 855-871. doi: 10.3934/mbe.2009.6.855
    [4] Yuhong Huo, Gourav Mandal, Lakshmi Narayan Guin, Santabrata Chakravarty, Renji Han . Allee effect-driven complexity in a spatiotemporal predator-prey system with fear factor. Mathematical Biosciences and Engineering, 2023, 20(10): 18820-18860. doi: 10.3934/mbe.2023834
    [5] Yuanfu Shao . Bifurcations of a delayed predator-prey system with fear, refuge for prey and additional food for predator. Mathematical Biosciences and Engineering, 2023, 20(4): 7429-7452. doi: 10.3934/mbe.2023322
    [6] Hongqiuxue Wu, Zhong Li, Mengxin He . Dynamic analysis of a Leslie-Gower predator-prey model with the fear effect and nonlinear harvesting. Mathematical Biosciences and Engineering, 2023, 20(10): 18592-18629. doi: 10.3934/mbe.2023825
    [7] Ranjit Kumar Upadhyay, Swati Mishra . Population dynamic consequences of fearful prey in a spatiotemporal predator-prey system. Mathematical Biosciences and Engineering, 2019, 16(1): 338-372. doi: 10.3934/mbe.2019017
    [8] Shunyi Li . Hopf bifurcation, stability switches and chaos in a prey-predator system with three stage structure and two time delays. Mathematical Biosciences and Engineering, 2019, 16(6): 6934-6961. doi: 10.3934/mbe.2019348
    [9] Rongjie Yu, Hengguo Yu, Chuanjun Dai, Zengling Ma, Qi Wang, Min Zhao . Bifurcation analysis of Leslie-Gower predator-prey system with harvesting and fear effect. Mathematical Biosciences and Engineering, 2023, 20(10): 18267-18300. doi: 10.3934/mbe.2023812
    [10] Wanxiao Xu, Ping Jiang, Hongying Shu, Shanshan Tong . Modeling the fear effect in the predator-prey dynamics with an age structure in the predators. Mathematical Biosciences and Engineering, 2023, 20(7): 12625-12648. doi: 10.3934/mbe.2023562
  • The COVID-19 pandemic has highlighted the importance of the digital economy in restoring economic and social development, creating more jobs and improving people's well-being. To inform policy makers about changes to digital strategies, measuring the digital economy is a prerequisite. This study aimed to compile an index of digital economy at the provincial (municipalities, autonomous regions, collectively referred to as "provinces") level to present an accurate and in-depth depiction of how it has developed in China. Our sample covers 31 provinces in China, over the period 2010–2020. This paper firstly constructs the digital economy index system from the four dimensions of digital users, digital platforms, digital industries and digital innovation, and then adopts a combination of entropy weighting method and grey target theory to measure the digital economy index. This paper study revealed that China's digital economy has been on an upward trend from 2010 to 2019 and has a decline in 2020, and the digital innovation is an important driving force for the growth of the digital economy index. The convergence of China's digital economy is decreasing, indicating that the gap in digital economy development between provinces is increasing. The proposed index in this study can be used as a screening tool, decision making tool, benchmarking tool and guidance of high-quality digital economy development.



    The effect of disease on eco-epidemiology system is a significant topic from both mathematical and ecological perspectives. The disease factor usually leads to a more complex and diverting dynamics than those in the disease-free system [1,2]. Within the interactions between predator and prey, the disease could only spread in prey or predator population, also could spread between prey and predator [3,4,5]. Birds (particularly pelicans) infect vibrio and die by preying on vibrio-infected fish (particularly tilapia) at the Salton Sea in the desert of Southern California [3], which is an example of disease spreads amongst the prey. For the disease in predator, taking fox rabies as an example, foxes (Vulpis) infect rabies and transmit to other foxes or their prey rabbits by biting in Europe and North America [6]. More relevant examples could be found in [7]. From the mathematical epidemiology point of view, one needs much more attention in the dynamics of infected predator to observe whether the presence of the prey allows the survival of a part of the predator population [8].

    A variety of diseased predator models have been proposed to study the complex interaction between prey and predator with infected diseases [2,9,10] and the reference therein. Most common epidemic model applied in predator-prey interactions is the SI-type, i.e., the predator population Y(t) is divided into two sub-classes, namely susceptible predator S(t) and infected predator I(t), respectively [10,11,12]. The infection term could be mass-action term (bilinear form) βSI or saturation form βSIS+I [4]. The infected predators usually behave differently with susceptible ones, and suffer an additional death rate. In a epidemic model, the global dynamics are usually determined by the basic reproduction number R0, i.e., the disease will dies out in the population when R01, and the disease will persist in the population when R0>1. However, the basic reproduction number is no longer a threshold parameter determining the global dynamics in diseased predator models, on the contrary, the dynamics are relatively comprehensive and unexpected.

    Predation is the key force in a prey-predator interaction, which could affect the size of prey population by direct hunting [9,13,14,15], and elicit a variety of anti-predator responses [16,17,18]. Consequently, prey tends to alter behaviors in a certain extent, such as change of habitat, foraging activity, vigilance, physiological changes. This anti-predator behaviors accelerate the extinction, evolution and development of prey population in the long run. Under the risk of predation, prey may reduce its foraging activity in order to stay alert, leading to starvation which impacts on population growth [19,20]. Therefore, an immediately result of anti-predator behaviors is the reduction of prey growth rate, which is the cost for prey in prey defense [19,21,22,23,24,25,26].

    Consider a simple birth-death process of the prey X(t) with the cost of anti-predator behaviors [27]:

    dXdt=[F(k,Y)a]XdX,

    where X,Y represent the density of the prey and predator, respectively. a is the birth rate of prey, d is the natural death rate of prey. F(k,Y) accounts for the cost of anti-predator defence due to fear, the parameter k reflects the level of fear which drives anti-predator behaviors of prey. The fear factor F(k,Y) has some specific assumptions under the ecological motions, for details see [20,27].

    To derive a simple diseased predator model incorporating the anti-predator defence due to fear, we adopted the following fear effect term F(k,Y):

    F(k,Y)=11+kY=11+k(S+I).

    Based on the results in [4,9,11], we can obtain the eco-epidemiological system with cost of anti-predator behaviors as following system of nonlinear differential equations:

    {dXdt=rX1+k(S+I)rX2KaXS1+bX,dSdt=eaXS1+bXd1SβSI,dIdt=βSId2I, (1.1)

    where X,S,I represent the density of prey, susceptible predator and infected predator at time t, respectively. r is the intrinsic growth rate of prey, K is the carrying capacity of the prey, a is the predation coefficient, b is the predators handling time of a prey, e is the biomass conversion constant, β is the transmissibility coefficient. d1 and d2 are the mortality rates of the susceptible predator and infected predator, and naturally d1<d2.

    This paper consists of six sections. In the next section, we prove the positivity and boundedness of the solution of system (1.1). In Section 3, we provide the existence conditions of the equilibria of the model. We analyze the stability of equilibria and show the occurrence of Hopf bifurcation in Section 4. In Section 5, the correctness of the theoretical proof is illustrated by numerical simulation. Finally, we summarize our results with ecological interpretations in Section 6.

    In view of the ecological significance, we only consider the solutions (X(t),S(t),I(t)) of system (1.1) on

    R3+={(X(t),S(t),I(t))R3+:X(t)0,S(t)0,I(t)0}.

    Theorem 2.1. Each solution of system (1.1) with initial value (X(0),S(0),I(0))R3+ is positive and ultimately bounded.

    Proof. Since the right-hand side of system (1.1) is completely continuous and locally Lipschitzian on R3+, the solution (X(t),S(t),I(t)) with initial condition (X(0),S(0),I(0))R3+ exists and is unique on R3+.

    By integrating, it follows from system (1.1) that

    X(t)=X(0)exp{t0(r1+k(S(τ)+I(τ))rX(τ)KaS(τ)1+bX(τ))dτ}0,S(t)=S(0)exp{t0(eaX(τ)1+bX(τ)d1βI(τ))dτ}0,I(t)=I(0)exp{t0(βS(τ)d2)dτ}0.

    Hence, the solution (X(t),S(t),I(t)) of system (1.1) with the initial condition (X(0),S(0),I(0))R3+ remains positive.

    From the first equation of (1.1), we can obtain

    dXdt=rX1+k(S+I)rX2KaXS1+bXrXrX2K=rX(1XK),

    then

    lim suptX(t)K.

    Let N(t)=eX(t)+S(t)+I(t), we can get

    dNdt=erX1+k(S+I)erX2Kd1Sd2IerXerX2Kd1Sd2IerX(1XK)+ed1Xd1NeK(r+d1)24rd1N,

    then

    lim suptN(t)eK(r+d1)24rd1.

    This ends the proof.

    Remark 2.2. From Theorem 2.1, we know that all positive solutions of system (1.1) with initial conditions (X(0),S(0),I(0))R3+ are defined in the following positive bounded invariant:

    Γ:={(X(t),S(t),I(t))R3+:0X(t)K,0eX(t)+S(t)+I(t)eK(r+d1)24rd1}.

    System (1.1) possesses at most three boundary equilibria:

    (i) Trivial equilibrium: E0=(0,0,0);

    (ii) Axial equilibrium: E1=(K,0,0);

    (iii) Planar equilibrium: E2=(X2,S2,0) exists if eabd1>0 and K>d1eabd1, where

    X2=d1eabd1,S2=[K(eabd1)2+rd1ke]+[K(eabd1)2rd1ke]2+4K2kre(eabd1)32Kk(eabd1)2. (3.1)

    For epidemic models, the most critical problem is the threshold property for the extinction and persistence of the disease, which is generally governed by the basic reproduction number R0. The basic reproduction number can be interpreted as the expected number of secondary cases produced, in a completely susceptible population, by a typical infected individual during its entire period of infectiousness. Following [28], we define the basic reproduction number for the predator population in the system (1.1) by

    R0:=βS2d2,

    where S2 is given by (3.1).

    Next, we mainly focus on the existence of positive equilibrium E3=(X3,S3,I3) of system (1.1). The coordinates X3,S3,I3 are positive solutions to the following system of equilibrium equations:

    {r1+k(S3+I3)rX3KaS31+bX3=0,eaX31+bX3d1βI3=0,βS3d2=0.

    Thus,

    S3=d2β,I3=X3(eabd1)d1(bX3+1)β,

    and X3 is the positive root of (3.2) in (X2,+):

    Q(X)=m3X3+m2X2+m1X+m0=0, (3.2)

    where

    m3:=bβr(k(eabd1)+b(kd2+β)),m2:=βr(Kb2βk(eabd1)2b(kd2+β)+bkd1),m1:=(kd2a(eabd1)b(akd22+aβd22β2r))Kβr(kd1+kd2+β),m0:=K(ad2(kd1kd2β)β2r).

    If eabd1>0 and r>ad2(β+k(d2d1))β2, we have

    m3<0,m0>0.

    By Descartes' rule of signs, system (1.1) has at least one positive equilibrium E3.

    Hence, we have the following results on the existence of the positive equilibrium. It is worthy to note that the positive equilibrium is not unique due to the impact of fear effect k.

    Theorem 3.1. If eabd1>0 and r>ad2(β+k(d2d1))β2, then system (1.1) has at least one positive equilibrium E3=(X3,S3,I3), where S3=S2R0, I3=X3(eabd1)d1(bX3+1)β and X3 is the positive root of (3.2) in (X2,+).

    Regarding the local stability of trivial equilibrium E0 and axial equilibrium E1, we have the following results. The proof is standard, so we omit it here.

    Theorem 4.1. For system (1.1),

    (i) The trivial equilibrium E0=(0,0,0) is unstable;

    (ii) If one of the following inequalities holds:

    (ii-1) eabd1<0;

    (ii-2) eabd1>0 and K<d1eabd1,

    then the axial equilibrium E1=(K,0,0) is stable; while E1=(K,0,0) is unstable if eabd1>0 and K>d1eabd1.

    Secondly, we will show the local stability of the planar equilibrium E2 of system (1.1). For convenience, set

    r1:=d2(eabd1)βe,r2:=d2(eabd1)(ea+bd1)aβe2,K1:=ea+bd1b(eabd1),K2:=βerd1(eabd1)(βerd2(eabd1)),k1:=Kb(eabd1)2(Kb(eabd1)(ea+bd1))ae2r(ea+bd1),k2:=β((eabd1)(d2(eabd1)βer)K+βerd1)(Kd2(eabd1)2+βerd1)d2. (4.1)

    Theorem 4.2. For system (1.1), assume that eabd1>0. If one of the following inequalities holds:

    (Ⅰ) rr1 and one of the following inequalities holds:

    (Ⅰ-1) d1eabd1<KK1;

    (Ⅰ-2) K>K1 and k>k1;

    (Ⅱ) r1<r<r2 and one of the following inequalities holds:

    (Ⅱ-1) d1eabd1<KK1;

    (Ⅱ-2) K1<K and k>max{k1,k2};

    (Ⅲ) r>r2 and one of the following inequalities holds:

    (Ⅲ-1) d1eabd1<KK2;

    (Ⅲ-2) K2<K and k>max{k1,k2},

    then equilibrium E2 is stable; otherwise, it is unstable.

    Proof. The Jacobian matrix of system (1.1) at E2 is given by

    J2=(a11a12a13a210βS200βS2d2),

    where

    a11:=X2(rK+abS2(1+bX2)2),a12:=krX2(1+kS2)2aX21+bX2,a13:=krX2(1+kS2)2,a21:=eaS2(1+bX2)2.

    Hence, the characteristic equation of J2 is given as

    f(λ)(λβS2+d2)=0, (4.2)

    where

    f(λ):=λ2a11λa12a21.

    Clearly, one can see that J2 has three eigenvalues λ1, λ2 and λ3=βS2d2. Since a12<0,a21>0, then a12a21>0.

    From (3.1), we can obtain

    a11=X2(rK+abS2(1+bX2)2)=X2Φ2Kka2e2,

    where

    Φ:=ab(K(eabd1)2rd1ke)2+4K2kre(eabd1)3ab(K(eabd1)2+rd1ke)2rka2e2.

    Note that the sign of Φ depends on

    ˜Φ:=a2b2((K(eabd1)2rd1ke)2+4K2kre(eabd1)3)(ab(K(aebd1)2+rd1ke)+2rka2e2)2=4a2kerP(k),

    where

    P(k):=ae2r(ea+bd1)k+Kb(eabd1)2(Kb(eabd1)(ea+bd1)).

    One can obtain that P(k) is decreasing with respect to k. If Kea+bd1b(eabd1) holds, we have P(0)0, which means that P(k)<0 for all k>0; if K>ea+bd1b(eabd1) and k>k1 hold, we can get P(k)<0. Therefore, when one of the following inequalities holds:

    (i) Kea+bd1b(eabd1);

    (ii) K>ea+bd1b(eabd1) and k>k1,

    we can obtain a11<0, which implies that the real parts of λ1 and λ2 are all negative.

    It follows from system (3.1) that

    βS2d2=β[K(eabd1)2+rd1ke]+β[K(eabd1)2rd1ke]2+4K2kre(eabd1)32Kk(eabd1)2d2=Θ2Kk(eabd1)2,

    where

    Θ:=K(eabd1)2(2kd2+β)βekrd1+β(K(eabd1)2rd1ke)2+4K2kre(eabd1)3.

    Note that the sign of Θ depends on

    ˜Θ:=β2(K(eabd1)2rd1ke)2+4β2K2kre(eabd1)3(K(eabd1)2(2kd2+β)+βekrd1)2=4Kk(aebd1)2[(Kd2(eabd1)2+βerd1)d2k+β((eabd1)(d2(eabd1)βer)K+βerd1)].

    Then if one of the following inequalities holds:

    (Ⅰ) eabd1>0 and rd2(eabd1)βe;

    (Ⅱ) eabd1>0, r>d2(eabd1)βe and one of the following inequalities:

    (Ⅱ-1) Kβerd1(eabd1)(βerd2(eabd1));

    (Ⅱ-2) K>βerd1(eabd1)(βerd2(eabd1)) and k>k2:=β((eabd1)(d2(eabd1)βer)K+βerd1)(Kd2(eabd1)2+βerd1)d2,

    we have λ3=βS2d2<0.

    Thus, we can arrive at the conclusion.

    It should be pointed out that another way to state Theorem 4.2 is as follows.

    Remark 4.3. For system (1.1), assume that eabd1>0 and R0<1. If one of the following inequalities:

    (Ⅰ) d1eabd1<KK1;

    (Ⅱ) K>K1 and k>k1

    holds, then the planar equilibrium E2 is stable; otherwise, it is unstable.

    Next, we will show the local stability of the positive equilibrium E3 of system (1.1).

    The Jacobian matrix of system (1.1) at E3 is given by

    J3=(b11b12b13b210d20βI30),

    where

    b11=X3(rK+abS3(1+bX3)2),b12=krX3(1+k(S3+I3))2aX31+bX3<0,b13=krX3(1+k(S3+I3))2<0,b21=eaS3(1+bX3)2>0. (4.3)

    The characteristic equation of J3 is given as

    λ3+A1λ2+A2λ+A3=0, (4.4)

    where

    A1=b11,A2=βd2I3b12b21,A3=b11βd2I3b13b21βI3. (4.5)

    Note that if A1>0 holds, then b11<0, which means that A3>0. According to Routh-Hurwitz criterion, the positive equilibrium E3 is locally asymptotically stable when A1>0 and A1A2A3>0.

    Therefore, we can establish the following statement.

    Theorem 4.4. Assume that eabd1>0 and r>ad2(β+k(d2d1))β2 hold. The positive equilibrium E3 of system (1.1) is locally asymptotically stable if A1>0 and A1A2A3>0, where Ai,i=1,2,3 is defined as in (4.5). Otherwise, it is unstable.

    Remark 4.5. Theorem 4.4 gives a sufficient condition about the stability of the positive equilibrium E3 for system (1.1). However, the complexity of model (1.1) leads to the failure to theoretically demonstrate how the fear factor affects the stability of the positive equilibrium. This will be discussed later through numerical simulations.

    In this subsection, we take k as the bifurcation parameter. The characteristic equation of system (1.1) at E3 is (4.4), and Ai(k),i=1,2,3 are defined as (4.5).

    Theorem 4.6. Hopf bifurcation near the positive equilibrium E3 for system (1.1) occurs whenever the critical parameter k attains the value k=kh in the domain:

    Ω={khR+:Δ(kh):=[A1(k)A2(k)A3(k)]|k=kh=0withA2(kh)>0,[dΔ(k)dk]|k=kh0}.

    Proof. If k=kh, the characteristic Eq (4.4) equals

    λ3+A1(kh)λ2+A2(kh)λ+A3(kh)=0, (4.6)

    then (4.6) can be factorized as

    (λ2+A2(kh))(λ+A1(kh))=0. (4.7)

    Clearly, (4.7) has three roots: λ1=iA2(kh), λ2=iA2(kh) and λ3=A1(kh). The roots are of the form λ1=p1(k)+ip2(k), λ2=p1(k)ip2(k) and λ3=p3(k), where pi(k)(i=1,2,3) are real numbers.

    From the characteristic Eq (4.4), we can get

    dλdk=λ2A1+λA2+A33λ2+2A1λ+A2, (4.8)

    where =ddk. Substituting λ=iA2 into (4.8), we obtain that

    A3A2A1+iA2A22(A2iA1A2)=dΔ(k)dk2(A21+A2)+i[A2A22A2A1A2dΔ(k)dk2A2(A21+A2)],

    which implies that

    [dRe(λ)dk]|k=kh=dΔ(k)dk2(A21+A2)|k=kh.

    By using monotonicity condition in the real part of the complex root dRe(λ)dk|k=kh0, the transversality condition dΔ(k)dk|k=kh0 can be obtained to ensure the existence of Hopf bifurcation.

    Results from numerical simulations are provided in this section to demonstrate our theoretical results. As we will show, the observations shed lights on the impact of fear factor. We choose the parameters of system (1.1) as follows:

    r=0.8,a=0.2,b=0.1,e=0.9,d1=0.05,β=0.1,d2=0.053. (5.1)

    Then we have

    eabd1=0.175>0,d1eabd1=0.286,r1=d2(eabd1)βe=0.103,r2=d2(eabd1)(ea+bd1)aβe2=0.106,K1=ea+bd1b(eabd1)=10.571,K2=βerd1(eabd1)(βerd2(eabd1))=0.328.

    Example 5.1 (The stability of E1).

    We adopt K=0.2,k=0.01, then system (1.1) has trivial equilibrium E0=(0,0,0) and axial equilibrium E1=(0.2,0,0). In this case, one can know that the conditions of Theorem 4.1 are satisfied, which means that E1 is locally asymptotically stable. The numerical results are shown in Figure 1.

    Figure 1.  Population dynamics of X(t), S(t) and I(t) of system (1.1) with K=0.2,k=0.01.

    Example 5.2 (The impacts of K and k on the stability of E2).

    In this example, we will choose three values of carrying capacity K for numerical experiments. We conclude that the carrying capacity and fear effect are other key factors related to the extinction of infected predators, in addition to the basic reproduction number R0.

    Firstly, we take K=0.3<K2, then we have k=0.1 which yields that R0=0.263<1. In this case, system (1.1) has trivial equilibrium E0=(0,0,0), axial equilibrium E1=(0.3,0,0), and planar equilibrium E2=(0.286,0.139,0). By Theorem 4.5, E2 is locally asymptotically stable, see Figure 2(a). Thus, when the carrying capacity of the prey K is small, no matter what the level of fear k is, the small size of prey population will lead to the extinction of infected predators.

    Figure 2.  The impacts of K and k on the stability of E2.

    Secondly, for comparison, we take K2<K=15, then

    k1=Kb(eabd1)2(Kb(eabd1)(ea+bd1))ae2r(ea+bd1)=0.149,k2=β((eabd1)(d2(eabd1)βer)K+βerd1)(Kd2(eabd1)2+βerd1)d2=10.873.

    Choosing k=0.1<max{k1,k2} which yields R0=5.792>1, then we have

    A1=0.25791>0,A1A2A3=0.00523>0.

    In this case, system (1.1) has trivial equilibrium E0=(0,0,0), axial equilibrium E1=(15,0,0), planar equilibrium E2=(0.286,3.070,0), and positive equilibrium E3=(7.351,0.530,7.125). By Theorem 4.5, E2=(0.286,3.070,0) is unstable. On the contrary, E3=(7.351,0.530,7.125) is locally asymptotically stable. The numerical simulation is shown in Figure 2(b).

    Finally, we take K2<K=60, then we have

    k1=Kb(eabd1)2(Kb(eabd1)(ea+bd1))ae2r(ea+bd1)=6.629,k2=β((eabd1)(d2(eabd1)βer)K+βerd1)(Kd2(eabd1)2+βerd1)d2=12.238.

    Choosing k=30>max{k1,k2} which yields that R0=0.649<1, system (1.1) has trivial equilibrium E0=(0,0,0), axial equilibrium E1=(60,0,0), and planar equilibrium E2=(0.286,0.344,0). By Theorem 4.5, E2 is locally asymptotically stable, see Figure 2(c). Thus, when the carrying capacity of the prey K is relatively large, a high level of fear k will lead to the extinction of infected predators.

    Example 5.3 (The impact of k on the stability of E3). We adopt K=60, then we have kh=0.26. In the next, we will choose three values of k, corresponding to the local stability of E3, Hopf bifurcation, and instability of E3, to illustrate the impact of fear factor on the population dynamics.

    Firstly, we take k=0.1<kh which yields that R0=5.881>1, then system (1.1) has trivial equilibrium E0=(0,0,0), axial equilibrium E1=(60,0,0), planar equilibrium E2=(0.286,3.117,0) and a unique positive equilibrium E3=(24.047,0.530,12.213). In this case, we obtain that

    A1=0.29863>0,A1A2A3=0.00065>0,

    which means that E3 is local asymptotically stable. The numerical results are shown in Figure 3.

    Figure 3.  Population dynamics of X(t), S(t) and I(t) of system (1.1) with K=60,k=0.1<kh.

    Secondly, we take k=0.26=kh which yields that R0=4.682>1, then system (1.1) has trivial equilibrium E0=(0,0,0), axial equilibrium E1=(60,0,0), planar equilibrium E2=(0.286,2.481,0) and a unique positive equilibrium E3=(12.975,0.530,9.665). In this case, we obtain that

    A1=0.14694>0,A1A2A3=0,

    which means that system (1.1) undergoes a Hopf bifurcation and there is a limit cycle around E3. The numerical results and the bifurcation diagrams of system (1.1) with respect to the parameter k are shown in Figures 4 and 5, respectively. Comparing Figures 3 and 4(a), one can see that there are two different implications induced by the fear factor k: the first is that the stability of E3 converts from stable into unstable, and the second is the decrease of values of X3 and I3 of E3.

    Figure 4.  Population dynamics of X(t), S(t) and I(t) of system (1.1) with K=60,k=0.26=kh. (a) Time-series plots; (b) Phase portraits in 3-dimensional space.
    Figure 5.  Bifurcation diagram of the system (1.1) with respect to the parameter k. Here K=60, other parameters are taken as in (5.1).

    Finally, we take k=0.5>kh which yields that R0=3.821>1, then system (1.1) has trivial equilibrium E0=(0,0,0), axial equilibrium E1=(60,0,0), planar equilibrium E2=(0.286,2.025,0) and a unique positive equilibrium E3=(7.685,0.530,7.322). In this case, we can obtain that

    A1=0.07642>0,A1A2A3=0.00051<0,

    which means that E3 is unstable. The numerical results are shown in Figure 6. One can find that the difference between Figures 4 and 6 is the decrease of values of E3 from (12.975,0.530,9.665) to (7.685,0.530,7.322), which is induced by the impact of the feat factor.

    Figure 6.  Population dynamics of X(t), S(t) and I(t) of system (1.1) with K=60,k=0.5>kh. (a) Time-series plots; (b) Phase portraits in 3-dimensional space.

    In this paper, we explored a predator-prey model that incorporates infectious disease in predator population and the cost of anti-predator behaviors. The cost of anti-predator behaviors is measured by a fear effect k leading to an reduction of prey's birth rate. We fulfill a complete stability analysis of equilibria for system (1.1) and show that the system (1.1) exhibits the Hopf bifurcation. Biologically, we focus on the impact of fear effect on the population dynamics. As we will see later, the cost of a high level of fear effect is disastrously. The main findings are summarized in the following.

    1) Small size of prey population leads to the extinction of infected predators.

    If the carrying capacity K is relatively small, the planar equilibrium E2 is stable, see Figure 2(a). Thus, no matter what the level of fear effect k is, a small size of prey population will lead to the extinction of infected predators.

    2) Low level of the fear effect doesn't impact on the population dynamics.

    If the level of fear effect k<kh, the positive (coexistence) equilibrium E3 is stable, see Figure 3. Hence, we conclude that a small fear effect k is not the key disturbance and does not change the coexistence dynamics of system (1.1). However, the densities of the prey and infected predator gradually decrease as k increasing.

    3) Certain medium level of the fear effect lead to periodic oscillation.

    If k=kh, the fear effect can destabilize the stability of E3 and will benefit the occurrence of periodic oscillation. In other words, system (1.1)undergoes a limit cycle, see Figures 4 and 5.

    4) High level of the fear effect leads to complex dynamics and the infected predator can go to extinction.

    If k>kh, E3 is unstable, see Figure 6. Therefore, a large fear effect k persistently and dramatically influence the population dynamics of prey and predator. Furthermore, if the level of the fear factor k is extremely high, the planar equilibrium E2 is stable, see Figure 2(c). The prey will respond to perceived predation risk and show a variety of anti-predator responses, dramatically decreasing the recruitment of susceptible predator, which will lead to an extinction of infected predator.

    The authors would like to thank the editor and the referees for their helpful comments. This research was supported by the National Natural Science Foundation of China (Grant No. 12171192, 12031020 and 12071173), the Science and Technology Research Projects of the Education Office of Jilin Province, China (JJKH20211033KJ), the Technology Development Program of Jilin Province, China (20210508024RQ) and Huaian Key Laboratory for Infectious Diseases Control and Prevention, China (HAP201704).

    The authors declare that there are no conflicts of interest regarding the publication of this paper.



    [1] Abrell T, Pihlajamaa M, Kanto L (2016) The role of users and customers in digital innovation: Insights from B2B manufacturing firms. Inf Manag 53: 324–335. https://doi.org/10.1016/j.im.2015.12.005 doi: 10.1016/j.im.2015.12.005
    [2] Ali Research Institute, KPMG (2018) 2018 Global Digital Economy Development Index. Available from: http://i.aliresearch.com/img/20180918/20180918153226.pdf.
    [3] Atkinson R, Foote C (2021) The 2020 State New Economy Index. Available from: https://www2.itif.org/2020-state-new-economy-index.pdf.
    [4] Baboo S, Nunkoo R, Kock F (2021) Social media attachment: Conceptualization and formative index construction. J Bus Res 139: 437–447. https://doi.org/10.1016/j.jbusres.2021.09.064 doi: 10.1016/j.jbusres.2021.09.064
    [5] Baker S, Bloom N, Davis S (2016) Measuring Economic Policy Uncertainty. Q J Econ 131: 1593–1636. https://doi.org/10.1093/qje/qjw024 doi: 10.1093/qje/qjw024
    [6] BEA (2018) Defining and Measuring the Digital Economy. Available from: https://www.bea.gov/digital Economy.
    [7] BEA (2022) Toward a Digital Economy Satellite Account. Available from: https://www.bea.gov/data/special-topics/digital-economy.
    [8] Beaunoyer E, Dupere S, Guitton MJ (2020) COVID-19 and digital inequalities: Reciprocal impacts and mitigation strategies. Comput Hum Behav 111: 106424. https://doi.org/10.1016/j.chb.2020.106424 doi: 10.1016/j.chb.2020.106424
    [9] CAICT (2017) White Paper on China's Digital Economy Development. Available from: http://www.caict.ac.cn/english/research/whitepapers/.
    [10] CAICT (2021) White Paper on China's Digital Economy Development. Available from: http://www.caict.ac.cn/english/research/whitepapers/.
    [11] CCID Consulting (2020) Digital Economic Development Index. Available from: http://www.ccidconsulting.com/.
    [12] Cornell University (2018) Global Innovation Index 2018: Energizing the World with Innovation. Available from: https://www.wipo.int/publications/en/details.jsp?id=4330.
    [13] European Commission (2021) Digital Economy and Society Index. Available from: https://digital-strategy.ec.europa.eu/en/news/digital-economy-and-society-index-2021.
    [14] Golinelli D, Boetto E, Carullo G, et al. (2020) Adoption of Digital Technologies in Health Care During the COVID-19 Pandemic: Systematic Review of Early Scientific Literature. J Med Internet Res 22: e22280. https://doi.org/10.2196/22280 doi: 10.2196/22280
    [15] Government of Australia (2021) Digital economy strategy 2030. Available from: https://apo.org.au/node/312247.
    [16] Hajkowicz S (2006) Multi-attributed environmental index construction. Ecol Econ 57: 122–139. https://doi.org/10.1016/j.ecolecon.2005.03.023 doi: 10.1016/j.ecolecon.2005.03.023
    [17] International Telecommunication Union (2017) Measuring the Information Society Report 2017. Available from: https://www.itu.int/en/ITU-D/Statistics/Pages/publications/mis2017.aspx.
    [18] Lee C, Zhong J (2016) Construction of a responsible investment composite index for renewable energy industry. Renew Sust Energ Rev 51: 288–303. https://doi.org/10.1016/j.rser.2015.05.071 doi: 10.1016/j.rser.2015.05.071
    [19] Li J, Yu F, Deng G, et al. (2017) Industrial Internet: A Survey on the Enabling Technologies, Applications, and Challenges. IEEE Commun Surv Tutor 19: 1504–1526. https://doi.org/10.1109/COMST.2017.2691349 doi: 10.1109/COMST.2017.2691349
    [20] Li T, Ma J (2021) Does digital finance benefit the income of rural residents? A case study on China. Quant Financ Econ 4: 664–688. https://doi.org/10.3934/QFE.2021030 doi: 10.3934/QFE.2021030
    [21] Liao G, Li Z, Wang M, et al. (2022) Measuring China's urban digital finance. Quant Financ Econ 6: 385–404. https://doi.org/10.3934/QFE.2022017 doi: 10.3934/QFE.2022017
    [22] OECD (2014) Measuring the Digital Economy: A New Perspective. http://dx.doi.org/10.1787/9789264221796-en
    [23] Selahattin K, Aykut E, Havvanur F (2021) The effect of COVID-19 pandemic on residential real estate prices: Turkish case. Quant Financ Econ 5: 623–639. https://doi.org/10.3934/QFE.2021028 doi: 10.3934/QFE.2021028
    [24] The 117th Congress (2021) S.1260-United States Innovation and Competition Act of 2021. Available from: https://www.congress.gov/bill/117th-congress/senate-bill/1260.
    [25] The World Economic Forum (2021) The Network Readiness Index 2021. Available from: https://networkreadinessindex.org/.
    [26] USAID (2021) digital Strategy 2020–2024. Available from: https://www.usaid.gov/usaid-digital-strategy.
    [27] Wang B, Tian J, Cheng L, et al. (2018) Spatial Differentiation of Digital Economy and Its Influencing Factors in China. Sci Geogr Sin 38: 859–868. https://doi.org/10.13249/j.cnki.sgs.2018.06.004 doi: 10.13249/j.cnki.sgs.2018.06.004
    [28] Whitmore A (2012) A statistical analysis of the construction of the United Nations E-Government Development Index. Gov Inf Q 29: 68–75. https://doi.org/10.1016/j.giq.2011.06.003 doi: 10.1016/j.giq.2011.06.003
    [29] Woetzel J, Seong J, Wang KW, (2017) China's digital economy: A leading global force. Available from: https://www.mckinsey.com/featured-insights/china/chinas-digital-economy-a-leading-global-force.
    [30] Xiang S, Wu W (2019) Research on the Design of China's Digital Economy Satellite Account Framework. Stat Res 36: 3–16. https://doi.org/10.19343/j.cnki.11-1302/c.2019.10.001 doi: 10.19343/j.cnki.11-1302/c.2019.10.001
    [31] Xu L, Li M, D Q (2017) Dynamic Evaluation Method Based on Grey Target. J Syst Sci Math Sci 37: 112–124. https://doi.org/10.12341/jssms13047 doi: 10.12341/jssms13047
    [32] Xu X, Zhang M (2020) Research on the Scale Measurement of China's Digital Economy: Based on the Perspective of internet Comparison. China Ind Econ 5: 23–41. https://doi.org/10.19581/j.cnki.ciejournal.2020.05.013 doi: 10.19581/j.cnki.ciejournal.2020.05.013
    [33] Zhu J, Hipel K (2012) Multiple stages grey target decision making method with incomplete weight based on multi-granularity linguistic label. Inf Sci 212: 15–32. https://doi.org/10.1016/j.ins.2012.05.011 doi: 10.1016/j.ins.2012.05.011
    [34] Zhu Y, Tian D, Yan F (2020) Effectiveness of Entropy Weight Method in Decision-Making. Math Probl Eng 2020: 3564835. https://doi.org/10.1155/2020/3564835 doi: 10.1155/2020/3564835
  • This article has been cited by:

    1. Chunmei Zhang, The effect of the fear factor on the dynamics of an eco-epidemiological system with standard incidence rate, 2024, 9, 24680427, 128, 10.1016/j.idm.2023.12.002
    2. Hongqiuxue Wu, Zhong Li, Mengxin He, Dynamic analysis of a Leslie-Gower predator-prey model with the fear effect and nonlinear harvesting, 2023, 20, 1551-0018, 18592, 10.3934/mbe.2023825
    3. Zhuoying Zhao, Xinhong Zhang, Unraveling the transmission mechanism of animal disease: Insight from a stochastic eco-epidemiological model driven by Lévy jumps, 2025, 191, 09600779, 115859, 10.1016/j.chaos.2024.115859
    4. Huazhou Mo, Yuanfu Shao, Stability and bifurcation analysis of a delayed stage-structured predator–prey model with fear, additional food, and cooperative behavior in both species, 2025, 2025, 2731-4235, 10.1186/s13662-025-03879-y
  • Reader Comments
  • © 2022 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(3047) PDF downloads(237) Cited by(32)

Figures and Tables

Figures(1)  /  Tables(17)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog