
Coronary microvascular dysfunction (CMD) is emerging as an important cause of myocardial ischemia, but there is a lack of a non-invasive method for reliable early detection of CMD.
To develop an electrocardiogram (ECG)-based machine learning algorithm for CMD detection that will lay the groundwork for patient-specific non-invasive early detection of CMD.
Vectorcardiography (VCG) was calculated from each 10-second ECG of CMD patients and healthy controls. Sample entropy (SampEn), approximate entropy (ApEn), and complexity index (CI) derived from multiscale entropy were extracted from ST-T segments of each lead in ECGs and VCGs. The most effective entropy subset was determined using the sequential backward selection algorithm under the intra-patient and inter-patient schemes, separately. Then, the corresponding optimal model was selected from eight machine learning models for each entropy feature based on five-fold cross-validations. Finally, the classification performance of SampEn-based, ApEn-based, and CI-based models was comprehensively evaluated and tested on a testing dataset to investigate the best one under each scheme.
ApEn-based SVM model was validated as the optimal one under the intra-patient scheme, with all testing evaluation metrics over 0.8. Similarly, ApEn-based SVM model was selected as the best one under the intra-patient scheme, with major evaluation metrics over 0.8.
Entropies derived from ECGs and VCGs can effectively detect CMD under both intra-patient and inter-patient schemes. Our proposed models may provide the possibility of an ECG-based tool for non-invasive detection of CMD.
Citation: Xiaoye Zhao, Yinlan Gong, Lihua Xu, Ling Xia, Jucheng Zhang, Dingchang Zheng, Zongbi Yao, Xinjie Zhang, Haicheng Wei, Jun Jiang, Haipeng Liu, Jiandong Mao. Entropy-based reliable non-invasive detection of coronary microvascular dysfunction using machine learning algorithm[J]. Mathematical Biosciences and Engineering, 2023, 20(7): 13061-13085. doi: 10.3934/mbe.2023582
[1] | Jing-An Cui, Fangyuan Chen . Effects of isolation and slaughter strategies in different species on emerging zoonoses. Mathematical Biosciences and Engineering, 2017, 14(5&6): 1119-1140. doi: 10.3934/mbe.2017058 |
[2] | Rodolfo Acuňa-Soto, Luis Castaňeda-Davila, Gerardo Chowell . A perspective on the 2009 A/H1N1 influenza pandemic in Mexico. Mathematical Biosciences and Engineering, 2011, 8(1): 223-238. doi: 10.3934/mbe.2011.8.223 |
[3] | Keguo Ren, Xining Li, Qimin Zhang . Near-optimal control and threshold behavior of an avian influenza model with spatial diffusion on complex networks. Mathematical Biosciences and Engineering, 2021, 18(5): 6452-6483. doi: 10.3934/mbe.2021321 |
[4] | Xiaomeng Wang, Xue Wang, Xinzhu Guan, Yun Xu, Kangwei Xu, Qiang Gao, Rong Cai, Yongli Cai . The impact of ambient air pollution on an influenza model with partial immunity and vaccination. Mathematical Biosciences and Engineering, 2023, 20(6): 10284-10303. doi: 10.3934/mbe.2023451 |
[5] | Julien Arino, Fred Brauer, P. van den Driessche, James Watmough, Jianhong Wu . A final size relation for epidemic models. Mathematical Biosciences and Engineering, 2007, 4(2): 159-175. doi: 10.3934/mbe.2007.4.159 |
[6] | Yiwei He, Qianqian Cui, Zengyun Hu . Modeling and analysis of the transmission dynamics of cystic echinococcosis: Effects of increasing the number of sheep. Mathematical Biosciences and Engineering, 2023, 20(8): 14596-14615. doi: 10.3934/mbe.2023653 |
[7] | Ning Bai, Juan Zhang , Li Li, Zhen Jin . Evaluating the effect of virus mutation on the transmission of avian influenza H7N9 virus in China based on dynamical model. Mathematical Biosciences and Engineering, 2019, 16(5): 3393-3410. doi: 10.3934/mbe.2019170 |
[8] | Muntaser Safan, Bayan Humadi . An SIS sex-structured influenza A model with positive case fatality in an open population with varying size. Mathematical Biosciences and Engineering, 2024, 21(8): 6975-7011. doi: 10.3934/mbe.2024306 |
[9] | Haitao Song, Lirong Guo, Zhen Jin, Shengqiang Liu . Modelling and stability analysis of ASFV with swill and the virus in the environment. Mathematical Biosciences and Engineering, 2022, 19(12): 13028-13049. doi: 10.3934/mbe.2022608 |
[10] | Kasia A. Pawelek, Anne Oeldorf-Hirsch, Libin Rong . Modeling the impact of twitter on influenza epidemics. Mathematical Biosciences and Engineering, 2014, 11(6): 1337-1356. doi: 10.3934/mbe.2014.11.1337 |
Coronary microvascular dysfunction (CMD) is emerging as an important cause of myocardial ischemia, but there is a lack of a non-invasive method for reliable early detection of CMD.
To develop an electrocardiogram (ECG)-based machine learning algorithm for CMD detection that will lay the groundwork for patient-specific non-invasive early detection of CMD.
Vectorcardiography (VCG) was calculated from each 10-second ECG of CMD patients and healthy controls. Sample entropy (SampEn), approximate entropy (ApEn), and complexity index (CI) derived from multiscale entropy were extracted from ST-T segments of each lead in ECGs and VCGs. The most effective entropy subset was determined using the sequential backward selection algorithm under the intra-patient and inter-patient schemes, separately. Then, the corresponding optimal model was selected from eight machine learning models for each entropy feature based on five-fold cross-validations. Finally, the classification performance of SampEn-based, ApEn-based, and CI-based models was comprehensively evaluated and tested on a testing dataset to investigate the best one under each scheme.
ApEn-based SVM model was validated as the optimal one under the intra-patient scheme, with all testing evaluation metrics over 0.8. Similarly, ApEn-based SVM model was selected as the best one under the intra-patient scheme, with major evaluation metrics over 0.8.
Entropies derived from ECGs and VCGs can effectively detect CMD under both intra-patient and inter-patient schemes. Our proposed models may provide the possibility of an ECG-based tool for non-invasive detection of CMD.
Domestic pigs originated from the Eurasian wild boar (Sus scrofa), which first appeared about 9000 years ago [1]. They are essential for the transmission of swine influenza. Human beings raise domestic pigs, and then slaughter them for pork [2]. Domestic pigs grow in the food and environment provided by human beings, while human beings get the necessary nutrients by eating them. Consequently, in the breeding process, the swine flu virus is transmitted to human beings through domestic pig-human contact [2,3]. According to this process, a new breed-slaughter model with swine influenza transmission can be proposed as a model (1.1).
{∂S1(x,t)∂t=D1∂2S1(x,t)∂x2+(B12N2(x,t)−ω0N1(x,t))N1(x,t)−s0S1(x,t) −β11I1(x,t)S1(x,t)+γ1I1(x,t),x∈R,t>0,∂I1(x,t)∂t=D1∂2I1(x,t)∂x2+β11I1(x,t)S1(x,t)−(s0+γ1)I1(x,t),x∈R,t>0,∂S2(x,t)∂t=D2∂2S2(x,t)∂x2+(b2−r2N2(x,t)K2+B21N1(x,t))N2(x,t)−d2S2(x,t) −2∑j=1β2jIj(x,t)S2(x,t)+γ2I2(x,t),x∈R,t>0,∂I2(x,t)∂t=D2∂2I2(x,t)∂x2+2∑j=1β2jIj(x,t)S2(x,t)−[e2+γ2+d2]I2(x,t),x∈R,t>0,Ni(x,t)=Si(x,t)+Ii(x,t),i=1,2,x∈R,t>0. | (1.1) |
Domestic pig population N1(x,t) and human population N2(x,t) are assumed to be divided into 2 epidemiological compartments: susceptibles (Si(x,t)) and infectives (Ii(x,t)) at time t and location x, i=1,2. Susceptibles can become infected by means of intra-species or inter-species transmission and then recover as new susceptibles. The notation B12 represents the human breeding parameter for the population growth of domestic pigs, while B21 represents the nutrients from eating domestic pigs to increase the birth rate of human beings. The notation s0 represents the slaughter rate of domestic pigs. It's noteworthy that domestic pigs cannot survive independently without human beings, but human beings can still survive well without the supply of pork [2]. Restrictions on the development of human population mainly come from intra-species competition.
For humans, the notation r2=b2−d2 is the intrinsic growth rate of humans, where b2 and d2 represents the natural natality rate and mortality rate, respectively. K2 is the environmental carrying capacity of human population without domestic pig supply. e2 is the additional mortality rate of humans caused by swine flu. For domestic pigs, ω0 represents the intraspecific competition. During the spread of swine flu, the parameters βij represent the per capita incidence rate from species j to species i, where i,j=1,2. γi denote the recovery rate for domestic animals and humans, i=1,2. D1 and D2 are the diffusion coefficients for domestic animals and humans. It is noteworthy that all parameters mentioned above is positive.
The main purpose of this paper is to propose a new breed-slaughter model with swine influenza transmission, and study the dynamic behavior of it. And then, we focus on the invasion process of infected domestic animals into the habitat of humans. Under certain conditions as in Theorem 2, we construct a propagating terrace linking human habitat to animal-human coexistent habitat, then to swine flu natural foci, which is divided by spreading speeds. Firstly, we calculate the equilibrium points of the model without spatial heterogeneity as a model (1.2) and analyze the existence of them by the persistence theory. Secondly, we discuss their stability by the basic reproduction number. Thirdly, we use these equilibrium points to construct a propagating terrace linking them by spreading speeds.
{dS1(t)dt=(B12N2(t)−ω0N1(t))N1(t)−s0S1(t)−β11I1(t)S1(t)+γ1I1(t),dI1(t)dt=β11I1(t)S1(t)−(s0+γ1)I1(t),dS2(t)dt=(b2−r2N2(t)K2+B21N1(t))N2(t)−d2S2(t)−2∑j=1β2jIj(t)S2(t)+γ2I2(t),dI2(t)dt=2∑j=1β2jIj(t)S2(t)−[e2+γ2+d2]I2(t),Ni(t)=Si(t)+Ii(t),i=1,2. | (1.2) |
In model (1.2), domestic pig population N1(t) and human population N2(t) are assumed to be divided into 2 epidemiological compartments: susceptibles (Si(t)) and infectives (Ii(t)) at time t, i=1,2. Other parameters are the same with model (1.1).
At first, we focus on the breed-slaughter system without swine flu transmission and spatial heterogeneity.
If I1(0)=I2(0)=0, N1(0)>0 and N2(0)>0, model (1.2) turns to a new breed-slaughter system without swine influenza transmission as model (2.1).
{dN1(t)dt=(B12N2(t)−ω0N1(t))N1(t)−s0N1(t),dN2(t)dt=r2(1−N2(t)K2)N2(t)+B21N1(t)N2(t),N1(0)>0,N2(0)>0, | (2.1) |
Similar to the competition system in [4], breed-slaughter system also has abundant dynamic results. For the positive equilibrium point
E∗=(N∗1,N∗2)=(r2(s0−B12K2)B12B21K2−ω0r2,K2(s0B21−ω0r2)B12B21K2−ω0r2) |
of model (2.1), we have three cases:
(a). If B12B21<ω0r2K2 and s0<min{ω0r2B21,B12K2}, the positive equilibrium point E∗ of model (2.1) is stable (Figure 2(a)).
(b). If B12B21>ω0r2K2 and s0>max{ω0r2B21,B12K2}, the positive equilibrium point E∗ of model (2.1) is unstable (Figure 2(b)).
(c). Other than the condition as (a) or (b), the positive equilibrium point E∗ of model (2.1) does not exist.
In order to reflect the effect of interspecific interaction on swine influenza transmission during breeding process as model (1.2), we suppose that B12B21<ω0r2K2 and s0<min{ω0r2B21,B12K2} to guarantee the existence and the stability of the boundary equilibrium point
E3=(r2(s0−B12K2)B12B21K2−ω0r2,0,K2(s0B21−ω0r2)B12B21K2−ω0r2,0) |
with I1(0)=I2(0)=0 in model (1.2).
After calculation, we summarize that there are at most 6 equilibrium points in R4+ of the system (1.2): E0=(0,0,0,0), E1=(0,0,K2,0), E2=(0,0,¯S2,¯I2), E3=(N∗1,0,N∗2,0), E4=(N′1,0,S′2,I′2), E5=(S∗1,I∗1,S∗2,I∗2), where ¯S2=e2+γ2+d2β22, ¯I2=β22K2−(e2+γ2+d2)β22, N∗1=N′1=r2(s0−B12K2)B12B21K2−ω0r2, N∗2=K2(s0B21−ω0r2)B12B21K2−ω0r2, S′2=e2+γ2+d2β22, I′2=β22K2(1+s0B21−B12B21K2B12B21K2−ω0r2)−(e2+γ2+d2)β22. The exact expression of E5 is unknown. However, under certain conditions as in Theorem 2, we can obtain its existence by persistence theory [5,6,7].
If there is no domestic pigs participation, namely N1(0)=S1(0)=I1(0)=0, The persistence and the stability of boundary equilibrium E2=(0,0,¯S2,¯I2) has been proved in [8]. Similarly, we define R0=β22K1b2+e2+γ2. And then, we can get the following lemma.
Lemma 1. If N1(0)=S1(0)=I1(0)=0 and I2(0)>0, {0}×{0}×R2+ is a invariant set of system (1.2). The trivial equilibrium point E0 in model (1.2) is unstable, and we have following two cases:
(a) If R0≤1, the disease-free equilibrium point E1 of model (1.2) is stable;
(b) If R0>1, model (1.2) has a unique equilibrium point E2 in the interior of {0}×{0}×R2+, which is stable, and E1 is unstable.
Furthermore, we consider the transmission process of human influenza with domestic pigs participating, but not infected from them. Namely I1(0)=0, I2(0)>0 and Ni(0)>0, i=1,2. The persistence and the stability of boundary equilibrium E4=(N′1,0,S′2,I′2) is similar to Lemma 1. Taking E3 as the original point by coordinate translation, we can get the following lemma by the persistence theory [5,9,10], when B12B21<ω0r2K2 and s0<min{ω0r2B21,B12K2}.
According to the definition of basic reproduction number in a single population as [5,11,12], we define R1=β1N∗1s0+γ1 as the basic reproduction number of swine flu transmission in demotic pig population and R2=β22N∗2e2+γ2+d2 as the basic reproduction number of swine flu transmission in human population.
Lemma 2. If N1(0)=S1(0)>0, I1(0)=0 and I2(0)>0, R+×{0}×R2+ is a invariant set of system (1.2). The trivial equilibrium point E0 and the boundary equilibrium point E1, E2 in model (1.2) are unstable when B12B21<ω0r2K2 and s0<min{ω0r2B21,B12K2}, and we have following two cases:
(a) If R2≤1, the disease-free equilibrium point E3 of model (1.2) is stable;
(b) If R2>1, model (1.2) has a unique equilibrium point E4 in the interior of R+×{0}×R2+, which is stable, and E3 is unstable.
Next we focus on the discussion about the existence and the stability of the positive equilibrium point E5=(S∗1,I∗1,S∗2,I∗2). At first, we define Rs=max{R1,R2}. Then, we get the theorem as the following.
Theorem 1. If Ni(0)>0 and Ii(0)>0, i=1,2, R4+ is a invariant set of system (1.2). The trivial equilibrium point E0 and the boundary equilibrium point E1, E2 in model (1.2) are unstable when B12B21<ω0r2K2 and s0<min{ω0r2B21,B12K2}, and we have following three cases:
(a) If Rs≤1, the disease-free equilibrium point E3 of model (1.2) is stable;
(b) If Rs>1, R1<R2 and R1≤1, model (1.2) has a unique equilibrium point E4 except for E0, E1, E2 and E3, which is stable, and E3 is unstable;
(c) If Rs>1 and R1≥R2 (or R2>R1>1) model (1.2) has a unique equilibrium point E5 in the interior of R4+, which is stable, and E3 and E4 are unstable.
Proof. If Rs≤1, E4 and E5 do not exist. Similar to the results of Lemma 2 (a), the disease-free equilibrium point E3 of model (1.2) is stable.
Then we consider the results of system (1.2) when Rs>1 and B12B21<ω0r2K2 and s0<min{ω0r2B21,B12K2}. At first, we define
D={(S1,I1,S2,I2)|0≤Ii≤Si+Ii≤N∗i,i=1,2}, |
D1={(S1,I1,S2,I2)|I1=0 or I2=0 ,0≤Si+Ii≤N∗i,i=1,2}, |
D2=D∖D1,˜D2={(S1,I1,S2,I2)|0<Ii≤Si+Ii≤N∗i,i=1,2}. |
D2 and ˜D2 are forward invariant.
Let Ω∗ consists of equilibria E0, E1, E2, E3 and E4. These equilibria cannot be chained to each other in D1. By analyzing the flow in neighborhood of each equilibrium, it is easy to see that Ω∗ is isolated in D and D1 is a uniform strong repeller for ˜D2.
If x(t)=(S1(t),I1(t),S2(t),I2(t)) stays close to E2, we have two cases: if I1(0)=I2(0)=0, then I1(t)=I2(t)=0; if I1(0)>0 or I2(0)>0, then I2(t)>0. Therefore, E2 is isolated in D. Similarly, we can prove that E0, E1 and E3 are isolated in D.
For E4 and E5, we have two cases: (A). R1<R2 and R1≤1; (B). R1≥R2 or R2>R1>1.
(A). R1<R2 and R1≤1
If R1<R2 and R1≤1, E5 do not exist. Similar to the results of Lemma 2 (b), the boundary equilibrium point E4 of model (1.2) is stable.
(B). R1≥R2 or R2>R1>1
If x(t)=(S1(t),I1(t),S2(t),I2(t)) stays close to E4, we have two cases: if I1(0)=0, then I1(t)=0; if I1(0)>0, then I1(t)>0. Since (S1(t),I1(t),S2(t),I2(t)) satisfying system (1.2) has no invariant subset other than E4 in its neighborhood. E4 is isolated in D and a uniform weak repeller for ˜D2. Therefore, we can prove that E0, E1, E2, E3 and E4 are isolated in D.
Using Proposition 4.3 in [5], we can prove that D1 is a uniform weak repeller for ˜D2; and using Theorem 4.5 in [5], we can prove that D1 is a uniform strong repeller for ˜D2.
Then we get that there exists an ϵ>0 such that
liminft→∞min{I1(t),I2(t)}>ϵ, |
with Ni(0)>0 and Ii(0)>0, i=1,2.
Therefore, if B12B21<ω0r2K2, s0<min{ω0r2B21,B12K2} and R1≥R2 (R2>R1>1), there exists at least one internal equilibrium of system (1.2) [9,10,13].
Next, we use Theorem 2 in [11] to discuss the basic reproduction number of system (1.2).
The Jacobian matrix of (I1,I2) is
J=(β11S1−(s0+γ1)0β21S2β22S2−(e2+γ2+d2)), |
Let J=F−V, F be the rate of appearance of new infections in compartment I, V be the rate of transfer of individuals out of compartment I. Then, we get
F=(β11S10β21S2β22S2), |
V= diag(s0+γ1e2+γ2+d2). |
We call FV−1 be the next generation matrix for the model (1.2) and set Rs=ρ(FV−1|E3), where ρ(A) denotes the spectral radius of a matrix A.
Then we get
Rs=max{β11N∗1s0+γ1,β22N∗2e2+γ2+d2}. |
Finally, using Theorem 2 in [11], we can prove Theorem 1.
The basic reproduction number is an important threshold value in the research of the epidemic mathematical model, which determines the disease to break out or not. However, it is not sufficient to discuss the breed-slaughter model with interspecific interaction. The main purpose of this paper is to investigate invasion process of infected domestic animal into human habitat. And we construct a propagating terrace linking human habitat E1 to animal-human coexistent habitat E3, then to swine flu natural foci E4 (or E5), which is divided by certain spreading speeds. The propagating terrace can describe the spatio-temporal continuous change of the transmission of swine flu.
Based on the heterogeneity of the population structure and the temporal and spatial continuity of the mammal movement, the population's spatial factor is considered in the spread of swine flu. If the swine flu host populations are distributed differently in space, the diffusion term may change their local population structure, thus change the swine flu epidemic. In order to describe the population invasion process, we set the initial value is zero in the area x∈(−∞,−x0)∪(x0,∞). The area of (−x0,x0) is the original habitat of N, and N will invade to the area of x∈(−∞,−x0)∪(x0,∞) at the spreading speed s [14].
The definition of spreading speed of a single population is the positive value s satisfied with the conditions as follows,
limt→+∞{sup|x|>ctN(x,t)}=0,∀c>s |
and
limt→+∞inf{inf|x|<ctN(x,t)}>0,∀c<s, |
in the model [4.1]
{∂N(x,t)∂t=D∂2N(x,t)∂x2+rN(x,t)(1−N(x,t)K),x∈R,t>0,N(x,0)=N0>0,x∈[−x0,x0],N(x,0)=0,x∈(−∞,x0)∪(x0,∞). | (4.1) |
The biological description of spreading speed s has been shown in the third figure of Figure 3. The value of s approximates the inverse of the slope of the color lines. It is easy to see that the co-effect of diffusion and reproduction leads to the population territory expansion, in which the local diffusion rate D guarantees the population spatial invasion to new areas and the reproduction rate r guarantees its development on occupied areas. The spreading speed of a single population in the model [4.1] is expressed by s:=2√Dr by [14]. However, it is not enough to study the swine flu with more than one host species [15,16,17,18,19,20]. We redefine the spreading speeds at the human-animal interface, as shown below.
s1:=2√D1(B12K2−s0), |
s2:=max{2√D1(β11N∗1−s0−γ1),2√D2(β22N∗2−e2−γ2−d2)}. |
Due to the participation of two populations, some notations need to be redefined. The notations s and x0 are replaced by si, xi, with i=1,2, corresponding to the two swine flu host populations.
Then we construct a propagating terrace linking human habitat E1 to animal-human coexistent habitat E3, then to swine flu natural foci E4 (or E5), which is divided by certain spreading speeds. The propagating terrace can describe the spatio-temporal continuous change of the transmission of swine flu, which can be show in Theorem 2.
Theorem 2. For system (1.1), if B12B21<ω0r2K2, B21N∗1<r2 and s0<min{ω0r2B21,B12K2}, the initial conditions satisfy that 0<S1(x,0)<N∗1, x∈[−x1,x1]; S1(x,0)=0, x∈(−∞,x1)∪(x1,∞), for some x1>0; 0<I1(x,0)<N∗1, x∈[−x2,x2]; I1(x,0)=0, x∈(−∞,x2)∪(x2,∞), for some x2>0; S2(x,0)=K2, I2(x,0)=0, x∈R.
We set
s1:=2√D1(B12K2−s0),s2:=max{2√D1(β11N∗1−s0−γ1),2√D2(β22N∗2−e2−γ2−d2)}. |
Suppose that s1>s2, x1>x2, then there are three cases about the invasion process as following:
(a) Rs≤1,
limt→+∞sup|x|>ct{|S1(x,t)|+|I1(x,t)|+|S2(x,t)−K2|+|I2(x,t)|}=0, ∀c>s1, |
limt→+∞sup|x|<ct{|S1(x,t)−N∗1|+|I1(x,t)|+|S2(x,t)−N∗2|+|I2(x,t)|}=0, ∀c<s1. |
The system (1.1) forms a propagating terrace, linking E1 to E3.
(b) If Rs>1, R1<R2 and R1≤1,
limt→+∞sup|x|>ct{|S1(x,t)|+|I1(x,t)|+|S2(x,t)−K2|+|I2(x,t)|}=0, ∀c>s1, |
limt→+∞supc2t<|x|<c1t{|S1(x,t)−N∗1|+|I1(x,t)|+|S2(x,t)−N∗2|+|I2(x,t)|}=0, ∀s2<c2<c1<s1, |
limt→+∞sup|x|<ct{|S1(x,t)−N′1|+|I1(x,t)|+|S2(x,t)−S′2|+|I2(x,t)−I′2|}=0, ∀c<s2. |
The system (1.1) forms a propagating terrace, linking E1 to E3, then to E4.
(c) If Rs>1 and R1≥R2 (or R2>R1>1),
limt→+∞sup|x|>ct{|S1(x,t)|+|I1(x,t)|+|S2(x,t)−K2|+|I2(x,t)|}=0, ∀c>s1, |
limt→+∞supc2t<|x|<c1t{|S1(x,t)−N∗1|+|I1(x,t)|+|S2(x,t)−N∗2|+|I2(x,t)|}=0, ∀s2<c2<c1<s1, |
limt→+∞sup|x|<ct{|S1(x,t)−S∗1|+|I1(x,t)−I∗1|+|S2(x,t)−S∗2|+|I2(x,t)−I∗2|}=0, ∀c<s2. |
The system (1.1) forms a propagating terrace, linking E1 to E3, then to E5.
Proof. The epidemic of swine flu originates in the interaction between humans and domestic animals in the breeding process, so the breaking out of swine flu would lag behind this process. Therefore, we first confirm the propagating terrace linking E1 and E3.
The breed-slaughter system without swine flu transmission can be transferred to model (4.2).
{∂N1(x,t)∂t=D1∂2N1(x,t)∂t2+(B12N2(x,t)−ω0N1(x,t))N1(x,t)−s0N1(x,t),∂N2(x,t)∂t=D2∂2N2(x,t)∂t2+r2(1−N2(x,t)K2)N2(t)+B21N1(x,t)N2(x,t). | (4.2) |
Let (N1,N2) be a solution to system (4.2) with the initial condition 0<N1(x,0)<N∗1, x∈[−x1,x1]; N1(x,0)=0, x∈(−∞,x1)∪(x1,∞), for some x1>0; N2(x,0)=K2, x∈R.
If B12B21<ω0r2K2 and s0<min{ω0r2B21,B12K2}, we claim that (N1(x,t),N2(x,t))∈Σ, ∀x∈R,t∈[0,∞), where
Σ:={(N1,N2)∈[0,N∗1]×[0,N∗2]:B12N2(x,t)−ω0N1(x,t)−s0≥0,r2(1−N2(x,t)K2)+B21N1(x,t)≥0}. |
By the strong maximum principle, N1≥0 for t>0. Then we get
∂N2(x,t)∂t≥D2∂2N2(x,t)∂t2+r2(1−N2(x,t)K2)N2(x,t). |
By a comparison, N2≥X, where X is the solution to
{∂X(x,t)∂t=D2∂2X(x,t)∂t2+r2(1−X(x,t)K2)X(x,t),X(x,0)=N2(x,0). | (4.3) |
Then we get the result
limt→+∞infN2(x,t)≥limt→+∞infX(x,t)=K2. |
Set u:=N1 and v:=N2−K2. Then ∂N2(x,t)∂t can be rewritten as
∂v(x,t)∂t=D2∂2v(x,t)∂t2−r2v(x,t)K2(v(x,t)+K2)+B21u(x,t)(v(x,t)+K2). |
Due to u=N1∈[0,N∗1] and v=N2−K2≥0 then
∂v(x,t)∂t≤D2∂2v(x,t)∂t2−(r2−B21N∗1)v(x,t)+B21K2u(x,t). |
By the strong maximum principle, if follows that v≤Y in R×[0,∞), where Y is the solution to
{∂Y(x,t)∂t=D2∂2Y(x,t)∂t2−(r2−B21N∗1)Y(x,t)+B21K2u(x,t),X(x,0)=0,x∈R. | (4.4) |
Then we have
Y(x,t)=B21K2∫t0{e−(r2−B21N∗1)(t−s)∫Re−(x−y)2/[4(t−s)]u(y,s)dy}ds. |
Given ϵ>0, we choose δ>0 small enough such that 2√D1(B12K2−s0+B12δ)<s1+ϵ.
For this δ, we claim that there is τ≫1 such that Y(x,t)<δ+Mu(x,t), ∀x∈R, t≥τ, for some positive constant M. Then it follows that N1 satisfies
∂N1(x,t)∂t≤D1∂2N1(x,t)∂t2+(B12K1−s0+B12δ−(B12M+ω0)N1(x,t))N1(x,t). |
Therefore, according to the comparison principle and the definition of spreading speed [14,15,16,19,21], for any c∈(2√D1(B12K2−s0+B12δ),s1+ϵ), it follows that limt→+∞sup|x|>ctN1(x,t)=0, and then limt→+∞sup|x|>ctN2(x,t)=K2.
Because of the arbitrariness of ϵ, we get
limt→+∞sup|x|>ct{|N1(x,t)|+|N2(x,t)−K2|}=0,∀c>s1. |
Thus, if the swine flu does not break out, namely Rs≤1, for system (1.1),
limt→+∞sup|x|>ct{|S1(x,t)|+|I1(x,t)|+|S2(x,t)−K2|+|I2(x,t)|}=0, ∀c>s1. |
Then we set U:=N∗1−N1 and V:=N∗2−N2. Similar to the proof before, we can get
limt→+∞sup|x|<ct{|N1(x,t)−N∗1|+|N2(x,t)−N∗2|}=0, ∀c<s1. |
If Rs≤1, for system (1.1),
limt→+∞sup|x|<ct{|S1(x,t)−N∗1|+|I1(x,t)|+|S2(x,t)−N∗2|+|I2(x,t)|}=0, ∀c<s1. |
Next we consider the propagating terrace linking E3 to (E4 or E5). Let (S1,I1,S2,I2) be a solution to system (1.1) with the initial condition S1(x,0)=N∗1, S2(x,0)=N∗2, I2(x,0)=0, x∈R. I1(x,0)>0, x∈[−x2,x2]; I1(x,0)=0, x∈(−∞,x2)∪(x2,∞), for some x2>0.
If B12B21<ω0r2K2 and s0<min{ω0r2B21,B12K2}, we claim that (S1(x,t)+I1(x,t),S2(x,t)+I2(x,t))∈Σ, ∀x∈R,t∈[0,∞).
For the spreading speed when Rs>1, comparison principle and strong maximum principle are no longer applicable due to the complexity of system (1.1). However, we can calculate the minimum wave speed from largest eigenvalue of its linearized system at E3 as [22] to link E3 and E4 (or E5).
For the following eigenvalue problem
1λAληλ=cηλ, |
where
Aλ=diag(Diλ2)+J|E3. |
J is the jacobian matrix,
J=(B12N2−2ω0N1−s0−β11I1−β11S1+γ1B12N10β11I1β11S1−(s0+γ1)00B12N2−β21S2r2(1−2N2K2)+B21N1−(β21I1+β22I2)−β22S2+γ20β21S2β21I1+β22I2β22S2−(e2+γ2+d2)). |
For λ≥0, the eigenvalues of the matrix
Aλ=(D1λ2−ω0N∗1−β11N∗1+γ1B12N∗100D1λ2+β11N∗1−(s0+γ1)00B12N∗2−β21N∗2D2λ2−r2N∗2K2−β22N∗2+γ20β21N∗20D2λ2+β22N∗2−(e2+γ2+d2)). |
are D1λ2+β11N∗1−(s1+γ1), D2λ2+β22N∗2−(e2+γ2+d2) and other two impossible results, which cannot define positive wave speed.
Thus, the minimum wave speed can be defined as follows, which can be divided the propagating terrace, linking E3 to E4 (or E5).
s2=max{infλ>0D1λ2+β11N∗1−(s0+γ1)λ,infλ>0D2λ2+β22N∗2−(e2+γ2+d2)λ}=max{2√D1(β11N∗1−s0−γ1),2√D2(β22N∗2−e2−γ2−d2)}. |
If Rs>1, there are two cases: (A). R1<R2 and R1≤1; (B). R1≥R2 or R2>R1>1.
(A). If R1<R2 and R1≤1, then E5 does not exist. s2=2√D2(β22N∗2−e2−γ2−d2), then we get
limt→+∞sup|x|>c2t+x2{|S1(x,t)−N∗1|+|I1(x,t)|+|S2(x,t)−N∗2|+|I2(x,t)|}=0, ∀c>s2, |
limt→+∞sup|x|<ct+x2{|S1(x,t)−N′1|+|I1(x,t)|+|S2(x,t)−S′2|+|I2(x,t)−I′2|}=0, ∀c<s2. |
Combining the results before, linking E1 to E3, then
limt→+∞sup|x|>ct+x1{|S1(x,t)|+|I1(x,t)|+|S2(x,t)−K2|+|I2(x,t)|}=0, ∀c>s1, |
limt→+∞supc2t+x2<|x|<c1t+x1{|S1(x,t)−N∗1|+|I1(x,t)|+|S2(x,t)−N∗2|+|I2(x,t)|}=0, ∀s2<c2<c1<s1, |
limt→+∞sup|x|<ct+x2{|S1(x,t)−N′1|+|I1(x,t)|+|S2(x,t)−S′2|+|I2(x,t)−I′2|}=0, ∀c<s2. |
The system (1.1) forms a propagating terrace, linking E1 to E3, then to E4.
(B). If R1≥R2 or R2>R1>1, set s2=max{2√D1(β11N∗1−s0−γ1),2√D2(β22N∗2−e2−γ2−d2)}. then we get
limt→+∞sup|x|>ct+x1{|S1(x,t)|+|I1(x,t)|+|S2(x,t)−K2|+|I2(x,t)|}=0, ∀c>s1, |
limt→+∞supc2t+x2<|x|<c1t+x1{|S1(x,t)−N∗1|+|I1(x,t)|+|S2(x,t)−N∗2|+|I2(x,t)|}=0, ∀s2<c2<c1<s1, |
limt→+∞sup|x|<ct+x2{|S1(x,t)−S∗1|+|I1(x,t)−I∗1|+|S2(x,t)−S∗2|+|I2(x,t)−I∗2|}=0, ∀c<s2. |
The system (1.1) forms a propagating terrace, linking E1 to E3, then to E5.
If s1>s2, x1>x2 and Rs>1, R1>R2>1, then in Figure 4, the blue area represents the original habitat area of humans at the population size of E1. After domesticating pigs, the red part will be shared with the two species at E3. While after swine flu transmitting between domestic pigs and humans, the internal red part will be shared again with two populations at E5 with swine flu transmission. It is a biological description of propagating terrace of humans with swine flu transmission, which is the local spacial variation of the population.
If s1>s2,x1>x2 and Rs≤1, then in Figure 5, the blue area represents the original habitat area of humans at the population size of E1. After domesticating pigs, the red part will be shared with the two species at E3. Because Rs≤1, there is no swine flu transmission during the breed and slaughter process. Then the propagating terrace links unstable equilibrium E1 and stable equilibrium E2.
We establish a new swine flu mathematical model to reflect the dynamic process of swine flu transmission with interspecific action between domestic pigs and humans, in which the roles of different species will no longer be at the same level. Domestic pigs cannot survive independently without human beings, but human beings can still survive well without the supply of pork. By our new swine flu model, we find that the human-animal interface has promoted the cross-species transmission of swine flu and resulted in the prevalence of flu in humans. In addition, the threshold values of population development and disease transmission are also discussed in order to provide a scientific basis for future health decision makers in swine flu prevention and control. We propose the zoonotic basic reproduction number Rs, which is more applicable to the study of swine flu transmission. Then, it is analyzed that the spreading speed of different species forming propagating terraces is influenced by the intrinsic growth rate r and diffusion rate D.
In this paper, the equilibrium points of the model are calculated and we analyze the existence of the equilibrium points by the persistence theory. Then we discuss their stability by the basic reproduction number. In addition, after redefining the spreading speed, we divide the propagating terrace with two populations, which is an unprecedented task. We concern with the invasion process of infected domestic animals into the habitat of humans. Under certain conditions as in Theorem 2, we construct a propagating terrace linking human habitat to animal-human coexistent habitat, then to swine flu natural foci, which is divided by spreading speeds.
This work is supported by the National Natural Science Foundation of China, 11771044.
The authors declared that they have no conflicts of interest to this work.
[1] |
P. Severino, A. D'Amato, M. Pucci, F. Infusino, F. Adamo, L. I. Birtolo, et al., Ischemic heart disease pathophysiology paradigms overview: from paque activation to microvascular dysfunction, Int. J. Mol. Sci., 21 (2020), 8118. https://doi.org/10.3390/ijms21218118 doi: 10.3390/ijms21218118
![]() |
[2] |
A. E. Moran, M. H. Forouzanfar, G. A. Roth, G. A. Mensah, M. Ezzati, C. J. L. Murray, et al., Temporal trends in ischemic heart disease mortality in 21 world regions, 1980 to 2010, Circulation, 129 (2014), 1483-1492. https://doi.org/10.1161/CIRCULATIONAHA.113.004042 doi: 10.1161/CIRCULATIONAHA.113.004042
![]() |
[3] | L. Dai, Y. Zang, G. Shou, L. Xia, Simulation of MCG signal in 2D cardiac tissue sheet with ischemic condition, in 2011 Computing in Cardiology, (2011), 21-24. Available from: http://cinc.mit.edu/archives/2011/pdf/0021.pdf. |
[4] |
J. Zhang, H. Li, Z. Pu, H. Liu, T. Huang, H. Cheng, et al., Early diagnosis of coronary microvascular dysfunction by myocardial contrast stress echocardiography, Math. Biosci. Eng., 20 (2023), 7845-7858. https://doi.org/10.3934/mbe.2023339 doi: 10.3934/mbe.2023339
![]() |
[5] |
F. Mangiacapra, M. M. Viscusi, G. Verolino, L. Paolucci, A. Nusca, R. Melfi, et al., Invasive assessment of coronary microvascular function, J. Clin. Med., 11 (2021), 228. https://doi.org/10.3390/jcm11010228 doi: 10.3390/jcm11010228
![]() |
[6] |
V. Kunadian, A. Chieffo, P. G. Camici, C. Berry, J. Escaned, A. H. E. M. Maas, et al., An EAPCI expert consensus document on Ischaemia with non-obstructive coronary arteries in collaboration with European society of cardiology working group on coronary pathophysiology & microcirculation endorsed by coronary vasomotor disorders international study group, Eur. Heart J., 41 (2020), 3504-3520. https://doi.org/10.1093/eurheartj/ehaa503 doi: 10.1093/eurheartj/ehaa503
![]() |
[7] |
S. H. Lee, D. Shin, J. M. Lee, T. P. van de Hoef, D. Hong, K. H. Choi, et al., Clinical relevance of Ischemia with nonobstructive coronary arteries according to coronary microvascular dysfunction, J. Am. Heart Assoc., 11 (2022), e025171. https://doi.org/doi:10.1161/JAHA.121.025171 doi: 10.1161/JAHA.121.025171
![]() |
[8] |
J. C. Kaski, F. Crea, B. J. Gersh, P. G. Camici, Reappraisal of ischemic heart disease, Circulation, 138 (2018), 1463-1480. https://doi.org/10.1161/CIRCULATIONAHA.118.031373 doi: 10.1161/CIRCULATIONAHA.118.031373
![]() |
[9] |
M. G. del Buono, R. A. Montone, M. Camilli, S. Carbone, J. Narula, C. J. Lavie, et al., Coronary microvascular dysfunction across the spectrum of cardiovascular diseases: JACC state-of-the-art review, J. Am. Coll. Cardiol., 78 (2021), 1352-1371. https://doi.org/10.1016/j.jacc.2021.07.042 doi: 10.1016/j.jacc.2021.07.042
![]() |
[10] | J. Zhou, Y. Onuma, S. Garg, N. Kotoku, S. Kageyama, S. Masuda, et al., Angiography derived assessment of the coronary microcirculation: is it ready for prime time? Expert Rev. Cardiovasc. Ther., 20 (2022), 549-566. https://doi.org/10.1080/14779072.2022.2098117 |
[11] |
C. N. B. Merz, C. J. Pepine, M. N. Walsh, J. L. F. P. G. Camici, W. M. Chilian, J. A. Clayton, et al., Ischemia and no obstructive coronary artery disease (INOCA), Circulation, 135 (2017), 1075-1092. https://doi.org/10.1161/CIRCULATIONAHA.116.024534 doi: 10.1161/CIRCULATIONAHA.116.024534
![]() |
[12] |
P. Ong, P. G. Camici, J. F. Beltrame, F. Crea, H. Shimokawa, U. Sechtem, et al., International standardization of diagnostic criteria for microvascular angina, Int. J. Cardiol., 250 (2018), 16-20. https://doi.org/10.1016/j.ijcard.2017.08.068 doi: 10.1016/j.ijcard.2017.08.068
![]() |
[13] |
B. Tjoe, L. Barsky, J. Wei, B. Samuels, B. Azarbal, C. N. B. Merz, et al., Coronary microvascular dysfunction: considerations for diagnosis and treatment, Cleve. Clin. J. Med., 88 (2021), 561-571. https://doi.org/doi:10.3949/ccjm.88a.20140 doi: 10.3949/ccjm.88a.20140
![]() |
[14] |
M. A. Marinescu, A. I. Löffler, M. Ouellette, L. Smith, C. M. Kramer, J. M. Bourque, Coronary microvascular dysfunction, microvascular angina, and treatment strategies, JACC Cardiovasc. Imaging, 8 (2015), 210-220. https://doi.org/10.1016/j.jcmg.2014.12.008 doi: 10.1016/j.jcmg.2014.12.008
![]() |
[15] |
V. M. Pereyra, A. Seitz, H. Mahrholdt, R. Bekeredjian, U. Sechtem, P. Ong, Coronary microvascular dysfunction in patients with mild-to-moderate aortic stenosis - insights from intracoronary acetylcholine testing, Int. J. Cardiol. Heart Vasc., 31 (2020), 100658. https://doi.org/10.1016/j.ijcha.2020.100658 doi: 10.1016/j.ijcha.2020.100658
![]() |
[16] |
E. Aribas, J. E. Roeters van Lennep, S. E. Elias-Smale, J. J. Piek, M. Roos, F. Ahmadizar, et al., Prevalence of microvascular angina among patients with stable symptoms in the absence of obstructive coronary artery disease: a systematic review, Cardiovasc. Res., 118 (2022), 763-771. https://doi.org/10.1093/cvr/cvab061 doi: 10.1093/cvr/cvab061
![]() |
[17] | Branch of Cardiovascular Physicians, Chinese Medical Doctor Association, 2020 expert consensus on the prevention and treatment of heart failure after myocardial infarction, Chin. Circ. J., 35 (2020), 1166-1180. Available from: https://corestarbio.com/uploadfile/1669779085_2.pdf. |
[18] |
H. Liu, S. Ou, P. Liu, Y. Xu, Y. Gong, L. Xia, et al., Effect of microcirculatory resistance on coronary blood flow and instantaneous wave-free ratio: a computational study, Comput. Methods Programs Biomed., 196 (2020), 105632. https://doi.org/10.1016/j.cmpb.2020.105632 doi: 10.1016/j.cmpb.2020.105632
![]() |
[19] |
J. Knuuti, W. Wijns, A. Saraste, D. Capodanno, E. Barbato, C. Funck-Brentano, et al., 2019 ESC Guidelines for the diagnosis and management of chronic coronary syndromes, Eur. Heart J., 41 (2019), 407-477. https://doi.org/10.1093/eurheartj/ehz425 doi: 10.1093/eurheartj/ehz425
![]() |
[20] |
Y. Geng, X. Wu, H. Liu, D. Zheng, L. Xia, Index of microcirculatory resistance: state-of-the-art and potential applications in computational simulation of coronary artery disease, J. Zhejiang Univ. Sci. B, 23 (2022), 123-140. https://doi.org/10.1631/jzus.B2100425 doi: 10.1631/jzus.B2100425
![]() |
[21] |
P. G. Camici, G. d'Amati, O. Rimoldi, Coronary microvascular dysfunction: mechanisms and functional assessment, Nat. Rev. Cardiol., 12 (2015), 48-62. https://doi.org/10.1038/nrcardio.2014.160 doi: 10.1038/nrcardio.2014.160
![]() |
[22] |
J. D. Sara, A. Sugrue, V. Kremen, B. Qiang, Y. Sapir, Z. I. Attia, et al., Electrocardiographic predictors of coronary microvascular dysfunction in patients with non-obstructive coronary artery disease: utility of a novel T wave analysis program, Int. J. Cardiol., 203 (2016), 601-606. https://doi.org/10.1016/j.ijcard.2015.10.228 doi: 10.1016/j.ijcard.2015.10.228
![]() |
[23] |
J. Herrmann, J. C. Kaski, A. Lerman, Coronary microvascular dysfunction in the clinical setting: from mystery to reality, Eur. Heart J., 33 (2012), 2771-2783. https://doi.org/10.1093/eurheartj/ehs246 doi: 10.1093/eurheartj/ehs246
![]() |
[24] |
Z. Zhou, Y. Liu, Z. Wang, Y. Wang, J. Zhang, C. Yang, The value of standard deviation of normal RR Interval in predicting coronary microvascular dysfunction in patients with ischemia with non obstructive coronary arteries (in Chinese), Chin. Circ. J., 37 (2022), 804-809. https://doi.org/10.3969/j.issn.1000-3614.2022.08.008 doi: 10.3969/j.issn.1000-3614.2022.08.008
![]() |
[25] | E. Fabris, A. W. J. van't Hof, 7- Angiography and electrocardiography (ECG) for the assessment of coronary microvascular obstruction, in Coronary Microvascular Obstruction in Acute Myocardial Infarction (eds. G. Niccoli and I. Eitel), New York: Academic Press, (2018), 109-126. |
[26] |
P. Ong, A. Athanasiadis, S. Hill, T. Schäufele, H. Mahrholdt, U. Sechtem, Coronary microvascular dysfunction assessed by intracoronary acetylcholine provocation testing is a frequent cause of ischemia and angina in patients with exercise-induced electrocardiographic changes and unobstructed coronary arteries, Clin. Cardiol., 37 (2014), 462-467. https://doi.org/10.1002/clc.22282 doi: 10.1002/clc.22282
![]() |
[27] |
N. Dose, M. M. Michelsen, N. D. Mygind, A. Pena, C. Ellervik, P. R. Hansen, et al., Ventricular repolarization alterations in women with angina pectoris and suspected coronary microvascular dysfunction, J. Electrocardiol., 51 (2018), 15-20. https://doi.org/10.1016/j.jelectrocard.2017.08.017 doi: 10.1016/j.jelectrocard.2017.08.017
![]() |
[28] |
J. D. Sara, R. J. Lennon, M. J. Ackerman, P. A. Friedman, P. A. Noseworthy, A. Lerman, Coronary microvascular dysfunction is associated with baseline QTc prolongation amongst patients with chest pain and non-obstructive coronary artery disease, J. Electrocardiol., 49 (2016), 87-93. https://doi.org/10.1016/j.jelectrocard.2015.10.006 doi: 10.1016/j.jelectrocard.2015.10.006
![]() |
[29] |
R. Alizadehsani, M. Abdar, M. Roshanzamir, A. Khosravi, P. M. Kebria, F. Khozeimeh, et al., Machine learning-based coronary artery disease diagnosis: a comprehensive review, Comput. Biol. Med., 111 (2019), 103346. https://doi.org/10.1016/j.compbiomed.2019.103346 doi: 10.1016/j.compbiomed.2019.103346
![]() |
[30] |
J. S. Richman, J. R. Moorman, Physiological time-series analysis using approximate entropy and sample entropy, Am. J. Physiol. Heart Circ. Physiol., 278 (2000), H2039-H2049. https://doi.org/10.1152/ajpheart.2000.278.6.H2039 doi: 10.1152/ajpheart.2000.278.6.H2039
![]() |
[31] | R. K. Udhayakumar, C. Karmakar, M. Palaniswami, Entropy profiling to detect ST change in heart rate variability signals, in 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), (2019), 4588-4591. https://doi.org/10.1109/EMBC.2019.8857297 |
[32] | H. Rabbani, M. P. Mahjoob, E. Farahabadi, A. Farahabadi, A. M. Dehnavi, Ischemia detection by electrocardiogram in wavelet domain using entropy measure, J. Res. Med. Sci., 16 (2011), 1473-1482. Available from: https://www.researchgate.net/publication/230843943. |
[33] | J. Liu, C. Zhang, T. Ristaniemi, F. Cong, Detection of myocardial infarction from multi-lead ECG using dual-Q tunable Q-factor wavelet transform, in 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), 2019 (2019), 1496-1499. https://doi.org/10.1109/EMBC.2019.8857775 |
[34] |
C. Chen, Y. Jin, I. L. Lo, H. Zhao, B. Sun, Q. Zhao, et al., Complexity change in cardiovascular disease, Int. J. Biol. Sci., 13 (2017), 1320-1328. https://doi.org/10.7150/ijbs.19462 doi: 10.7150/ijbs.19462
![]() |
[35] |
J. M. Ahn, S. Arora, O. G. Solberg, O. Angerås, K. Rolid, M. Rafique, et al., Prognostic value of comprehensive intracoronary physiology assessment early after heart transplantation, Eur. Heart J., 42 (2021), 4918-4929. https://doi.org/10.1093/eurheartj/ehab568 doi: 10.1093/eurheartj/ehab568
![]() |
[36] |
E. A. P. Alday, A. Gu, A. J. Shah, C. Robichaux, A. K. I. Wong, C. Liu, et al., Classification of 12-lead ECGs: the physioNet/computing in cardiology challenge 2020, Physiol. Meas., 41 (2021), 124003-124003. https://doi.org/10.1088/1361-6579/abc960 doi: 10.1088/1361-6579/abc960
![]() |
[37] | X. Zhao, J. Zhang, Y. Gong, L. Xu, H. Liu, S. Wei, et al., Reliable detection of myocardial ischemia using machine learning based on temporal-spatial characteristics of electrocardiogram and vectorcardiogram, Front. Physiol., 13 (2022). https://doi.org/10.3389/fphys.2022.854191 |
[38] |
C. E. Kossmann, D. Brody, G. Burch, H. Hecht, F. Johnston, C. Kay, et al., Report of committee on electrocardiography, American Heart Association, Recommendations for standardization of leads and of specifications for instruments in electrocardiography and vectorcardiography, Circulation, 35 (1967), 583-602. https://doi.org/10.1161/01.cir.35.3.583 doi: 10.1161/01.cir.35.3.583
![]() |
[39] |
J. A. Kors, G. van Herpen, A. C. Sittig, J. H. van Bemmel, Reconstruction of the Frank vectorcardiogram from standard electrocardiographic leads: diagnostic comparison of different methods, Eur. Heart J., 11 (1990), 1083-1092. https://doi.org/10.1093/oxfordjournals.eurheartj.a059647 doi: 10.1093/oxfordjournals.eurheartj.a059647
![]() |
[40] |
J. Song, H. Yan, Z. Xiao, X. Yang, X. Zhang, A robust and efficientalgorithm for ST–T complex detection in electrocardiograms, J. Mech. Med. Biol., 11 (2012), 1103-1111. https://doi.org/10.1142/s0219519411004198 doi: 10.1142/s0219519411004198
![]() |
[41] |
A. Molina-Picó, D. Cuesta-Frau, M. Aboy, C. Crespo, P. Miró-Martínez, S. Oltra-Crespo, Comparative study of approximate entropy and sample entropy robustness to spikes, Artif. Intell. Med., 53 (2011), 97-106. https://doi.org/10.1016/j.artmed.2011.06.007 doi: 10.1016/j.artmed.2011.06.007
![]() |
[42] |
M. A. Busa, R. E. A. van Emmerik, Multiscale entropy: a tool for understanding the complexity of postural control, J. Sport Health Sci., 5 (2016), 44-51. https://doi.org/10.1016/j.jshs.2016.01.018 doi: 10.1016/j.jshs.2016.01.018
![]() |
[43] |
J. M. Yentes, N. Hunt, K. K. Schmid, J. P. Kaipust, D. McGrath, N. Stergiou, The appropriate use of approximate entropy and sample entropy with short data sets, Ann. Biomed. Eng., 41 (2013), 349-365. https://doi.org/10.1007/s10439-012-0668-3 doi: 10.1007/s10439-012-0668-3
![]() |
[44] | D. U. Uguz, F. Berief, S. Leonhardt, C. H. Antink, Classification of 12-lead ECGs using gradient boosting on features acquired with domain-specific and domain-agnostic methods, in 2020 Computing in Cardiology, (2020), 1-4. https://doi.org/10.22489/CinC.2020.363 |
[45] | A. U. Haq, J. Li, M. H. Memon, M. H. Memon, J. Khan, S. M. Marium, Heart disease prediction system using model of machine learning and sequential backward selection algorithm for features selection, in 2019 IEEE 5th International Conference for Convergence in Technology (I2CT), 2019. https://doi.org/10.1109/I2CT45611.2019.9033683 |
[46] |
M. Alnowami, F. Abolaban, E. Taha, A wrapper-based feature selection approach to investigate potential biomarkers for early detection of breast cancer, J. Radiat. Res. Appl. Sci., 15 (2022), 104-110. https://doi.org/10.1016/j.jrras.2022.01.003 doi: 10.1016/j.jrras.2022.01.003
![]() |
[47] |
L. Li, X. Sui, J. Lian, F. Yu, Y. Zhou, Vehicle interaction behavior prediction with self-attention, Sensors, 22 (2022), 429. https://doi.org/10.3390/s22020429 doi: 10.3390/s22020429
![]() |
[48] |
C. Xie, R. Du, J. W. K. Ho, H. H. Pang, K. W. H. Chiu, E. Y. P. Lee, et al., Effect of machine learning re-sampling techniques for imbalanced datasets in 18F-FDG PET-based radiomics model on prognostication performance in cohorts of head and neck cancer patients, Eur. J. Nucl. Med. Mol. Imaging, 47 (2020), 2826-2835. https://doi.org/10.1007/s00259-020-04756-4 doi: 10.1007/s00259-020-04756-4
![]() |
[49] |
T. Saito, M. Rehmsmeier, The precision-recall plot is more informative than the ROC plot when evaluating binary classifiers on imbalanced datasets, PLoS One, 10 (2015), e0118432. https://doi.org/10.1371/journal.pone.0118432 doi: 10.1371/journal.pone.0118432
![]() |
[50] |
S. Ansari, N. Farzaneh, M. Duda, K. Horan, H. B. Andersson, Z. D. Goldberger, et al., A review of automated methods for detection of myocardial ischemia and infarction using electrocardiogram and electronic health records, IEEE Rev. Biomed. Eng., 10 (2017), 264-298. https://doi.org/10.1109/RBME.2017.2757953 doi: 10.1109/RBME.2017.2757953
![]() |
[51] |
R. Correa, P. D. Arini, L. S. Correa, M. Valentinuzzi, E. Laciar, Novel technique for ST-T interval characterization in patients with acute myocardial ischemia, Comput. Biol. Med., 50 (2014), 49-55. https://doi.org/10.1016/j.compbiomed.2014.04.009 doi: 10.1016/j.compbiomed.2014.04.009
![]() |
[52] |
Y. L. Zang, L. Xia, Cellular mechanism of cardiac alternans: an unresolved chicken or egg problem, J. Zhejiang Univ. Sci. B, 15 (2014), 201-211. https://doi.org/10.1631/jzus.B1300177 doi: 10.1631/jzus.B1300177
![]() |
[53] |
M. J. Janse, A. L. Wit, Electrophysiological mechanisms of ventricular arrhythmias resulting from myocardial ischemia and infarction, Physiol. Rev., 69 (1989), 1049-1169. https://doi.org/10.1152/physrev.1989.69.4.1049 doi: 10.1152/physrev.1989.69.4.1049
![]() |
[54] |
C. C. ter Haar, A. C. Maan, S. G Warren, M. Ringborn, B. M. Horáček, M. J. Schalij, et al., Difference vectors to describe dynamics of the ST segment and the ventricular gradient in acute ischemia, J. Electrocardiol., 46 (2013), 302-311. https://doi.org/10.1016/j.jelectrocard.2013.04.004 doi: 10.1016/j.jelectrocard.2013.04.004
![]() |
[55] |
A. Feeny, L. G. Tereshchenko, Beat-to-beat determinants of the beat-to-beat temporal and spatial variability of repolarization, J. Electrocardiol., 49 (2016), 417-422. https://doi.org/10.1016/j.jelectrocard.2016.01.007 doi: 10.1016/j.jelectrocard.2016.01.007
![]() |
[56] | M. Wei, J. Z. Song, H. Yan, Electrocardiogram ST-segment morphology variability analysis base on correlation coefficient entropy and inverse correlation coefficient entropy, Appl. Mech. Mater., 195-196 (2012), 550-554. https://doi.org/10.4028/www.scientific.net/AMM.195-196.550 |
[57] |
M. Ferrario, M. G. Signorini, G. Magenes, S. Cerutti, Comparison of entropy-based regularity estimators: application to the fetal heart rate signal for the identification of fetal distress, IEEE Trans. Biomed. Eng., 53 (2006), 119-125. https://doi.org/10.1109/TBME.2005.859809 doi: 10.1109/TBME.2005.859809
![]() |
[58] |
S. D. Fihn, J. C. Blankenship, K. P. Alexander, J. A. Bittl, J. G. Byrne, B. J. Fletcher, et al., 2014 ACC/AHA/AATS/PCNA/SCAI/STS Focused Update of the Guideline for the Diagnosis and Management of Patients With Stable Ischemic Heart Disease: A Report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines, and the American Association for Thoracic Surgery, Preventive Cardiovascular Nurses Association, Society for Cardiovascular Angiography and Interventions, and Society of Thoracic Surgeons, J. Am. Coll. Cardiol., 64 (2014), 1929-1949. https://doi.org/10.1016/j.jacc.2014.07.017 doi: 10.1016/j.jacc.2014.07.017
![]() |
[59] |
H. A. Mayala, W. Yan, H. Jing, S. Liu, G. Yi, C. Qin, et al., Clinical characteristics and biomarkers of coronary microvascular dysfunction and obstructive coronary artery disease, J. Int. Med. Res., 47 (2019), 6149-6159. https://doi.org/10.1177/0300060519859134 doi: 10.1177/0300060519859134
![]() |