
Piper guineense and honey contain antioxidative, anti-inflammatory, and antimicrobial properties that can help restore neuronal and other cell damage. To investigate the neuromodulatory roles of p. guineense and honey against lead toxicity on the hippocampus and cerebellum, impairing social behaviors and motor activities.
Thirty Wistar rats were separated into six groups of five rats each, marked with dye. Group A served as control; B was untreated lead; C was a medium dose of the extract (50 mg/kg) and honey (1000 mg/kg); D was a high dose of the extract (80 mg/kg) and honey (1500 mg/kg); E received extract (80 mg/kg), and F received honey (1500 mg/kg). All groups received 110 mg/kg of lead orally, except the control. Social interaction, antidepressant effects, and motor activities were studied using a sociability chamber (SC), Forced Swim Test (FST), and String methods. A blood sample was used to evaluate glutathione peroxidase (GPx) and glutathione oxide transaminase (GOT), while the lipid level was estimated using cerebellar homogenate. Neuronal damage, vacuolation, necrosis, cell degeneration, and alterations in both hippocampus and cerebellum marked untreated group, with decreased GPx and GOT activities followed by impaired motor activities, social behavior, memory, and motivation. Using SCT, group B spent significantly lesser time (47.60 ± 47.60) with stranger 1 compared to A (138.20 ± 34.05), while group C spent considerably more time with stranger 1 (86.80 ± 30.32) than group B at P ≥ 0.05. The treatment increased the enzyme level and restored histoarchitecture (
Citation: UCHEWA O. Obinna, EMECHETA S. Shallom, EGWU A. Ogugua, EDE C. Joy, IBEGBU O. Augustine. Neuromodulatory roles of PIPER GUINEENSE and honey against Lead-Induced neurotoxicity in social interactive behaviors and motor activities in rat models[J]. AIMS Neuroscience, 2022, 9(4): 460-478. doi: 10.3934/Neuroscience.2022026
[1] | Jun Pan, Haijun Wang, Feiyu Hu . Revealing asymmetric homoclinic and heteroclinic orbits. Electronic Research Archive, 2025, 33(3): 1337-1350. doi: 10.3934/era.2025061 |
[2] | Bin Long, Shanshan Xu . Persistence of the heteroclinic loop under periodic perturbation. Electronic Research Archive, 2023, 31(2): 1089-1105. doi: 10.3934/era.2023054 |
[3] | Guowei Yu . Heteroclinic orbits between static classes of time periodic Tonelli Lagrangian systems. Electronic Research Archive, 2022, 30(6): 2283-2302. doi: 10.3934/era.2022116 |
[4] | Xianjun Wang, Huaguang Gu, Bo Lu . Big homoclinic orbit bifurcation underlying post-inhibitory rebound spike and a novel threshold curve of a neuron. Electronic Research Archive, 2021, 29(5): 2987-3015. doi: 10.3934/era.2021023 |
[5] | Minzhi Wei . Existence of traveling waves in a delayed convecting shallow water fluid model. Electronic Research Archive, 2023, 31(11): 6803-6819. doi: 10.3934/era.2023343 |
[6] | Anna Gołȩbiewska, Marta Kowalczyk, Sławomir Rybicki, Piotr Stefaniak . Periodic solutions to symmetric Newtonian systems in neighborhoods of orbits of equilibria. Electronic Research Archive, 2022, 30(5): 1691-1707. doi: 10.3934/era.2022085 |
[7] | Xianjun Wang, Huaguang Gu . Post inhibitory rebound spike related to nearly vertical nullcline for small homoclinic and saddle-node bifurcations. Electronic Research Archive, 2022, 30(2): 459-480. doi: 10.3934/era.2022024 |
[8] | Xiaochen Mao, Weijie Ding, Xiangyu Zhou, Song Wang, Xingyong Li . Complexity in time-delay networks of multiple interacting neural groups. Electronic Research Archive, 2021, 29(5): 2973-2985. doi: 10.3934/era.2021022 |
[9] | Peng Mei, Zhan Zhou, Yuming Chen . Homoclinic solutions of discrete p-Laplacian equations containing both advance and retardation. Electronic Research Archive, 2022, 30(6): 2205-2219. doi: 10.3934/era.2022112 |
[10] | Caiwen Chen, Tianxiu Lu, Ping Gao . Chaotic performance and circuitry implement of piecewise logistic-like mapping. Electronic Research Archive, 2025, 33(1): 102-120. doi: 10.3934/era.2025006 |
Piper guineense and honey contain antioxidative, anti-inflammatory, and antimicrobial properties that can help restore neuronal and other cell damage. To investigate the neuromodulatory roles of p. guineense and honey against lead toxicity on the hippocampus and cerebellum, impairing social behaviors and motor activities.
Thirty Wistar rats were separated into six groups of five rats each, marked with dye. Group A served as control; B was untreated lead; C was a medium dose of the extract (50 mg/kg) and honey (1000 mg/kg); D was a high dose of the extract (80 mg/kg) and honey (1500 mg/kg); E received extract (80 mg/kg), and F received honey (1500 mg/kg). All groups received 110 mg/kg of lead orally, except the control. Social interaction, antidepressant effects, and motor activities were studied using a sociability chamber (SC), Forced Swim Test (FST), and String methods. A blood sample was used to evaluate glutathione peroxidase (GPx) and glutathione oxide transaminase (GOT), while the lipid level was estimated using cerebellar homogenate. Neuronal damage, vacuolation, necrosis, cell degeneration, and alterations in both hippocampus and cerebellum marked untreated group, with decreased GPx and GOT activities followed by impaired motor activities, social behavior, memory, and motivation. Using SCT, group B spent significantly lesser time (47.60 ± 47.60) with stranger 1 compared to A (138.20 ± 34.05), while group C spent considerably more time with stranger 1 (86.80 ± 30.32) than group B at P ≥ 0.05. The treatment increased the enzyme level and restored histoarchitecture (
Many countries already recognize measles as a highly contagious and dangerous disease for infants and children, which can cause death in the pre-immunization period [1,2]. World Health Organization (WHO) has estimated that there were around 535,000 human deaths due to measles in 2000, and most of them occurred in tropical and sub-tropical such as Indonesia, China, Malaysia, Colombia, and Myanmar. Since then, several measles immunization programs have been launched in high-risk countries which have an impact on reducing global measles deaths by 78% in 2012 [3]. Dramatic decreases in the number of cases have also been reported in China since 1986 and can reach around 0.5 per million in the 1995–2005 period. It is an impact of the implementation of the two-dose measles vaccination program, which has been carried out during that period [4,5].
Meanwhile, a different condition for measles cases was reported in Australia. The outbreaks occur even with a reasonably large reproduction number (R0) of around 15 to 20, whereas efforts to eliminate measles have been launched [6]. Consequently, the coverage of the immunization level was increased to 92–95% in order to reduce the number of measles cases [7]. Measles outbreaks have also been reported in Italy, with 4347 confirmed cases and 5404 notified cases in 2017 [8].
Measles are caused by enveloped RNA virus with the genus Morbillivirus in the family paramyxovirus. Transmission usually occurs through direct contact with sufferers or via air such as mucus contact, coughing, and sneezing [9]. Sharing tableware with an infected person may also have a high potential for measles. Other names for this disease are morbilli or rubeola. Measles transmission starts from infected individuals coughing or sneezing who can make the virus spread through the air and can survive up to several hours.
The virus infects the respiratory tract and then spreads throughout the body. The period of measles transmission begins four days before and four days after the rash appears, the peak of transmission at the time of the initial symptoms (prodromal phase), namely in the first 1–3 days of illness [11]. Complications that often occur are diarrhea and bronchopneumonia [12]. Specific symptoms of this disease are reddish skin rashes that appear within the incubation period (7–18 days) and can last for 5–6 days [10].
Symptoms of measles are characterized by fever and are accompanied by one or more symptoms of cough, runny nose, red eyes, or watery eyes and then followed by reddish spots or rash that usually appears starting from behind the ear. Furthermore, symptoms appear on the body in the form of maculopapular shape and found grayish-white patches on the inner cheeks. In children, this disease can cause deadly severe complications if not handled properly. Therefore, early treatment, such as consultation with a doctor or the nearest health care provider, is needed to avoid complications.
The disease can be prevented by immunization using two doses of MMR (measles, mumps, and rubella) vaccine, which is highly effective and safe [13]. However, the implementation of immunization strategy has failed to cover rural populations in developing countries [14]. One in five children in the world (around 24 million children) did not receive the first dose of measles vaccination in a routine immunization program in 2008 [2]. Efforts to overcome the spread of measles globally by providing safe and cost-effective vaccines have been improving rapidly to reduce the mortality of children from measles. The impact was a dramatic decrease in the number of measles cases [10]. World Health Assembly endorsed Global Vaccine Action Plan (GVAP) for eliminating measles, congenital rubella syndrome (CRS), and rubella in 2012 [15], followed by all WHO regions who agreed to eradicate measles [10]. In Indonesia, there were around 23,164 cases of measles and 30,463 cases of rubella from 2010 to 2015. The number of these cases was allegedly lower than the reality, given the number of cases that were not reported [16]. In 2016, there were around 3112 measles cases reported in Jakarta [17].
Modeling of measles with vaccination provides a better understanding and interpreting collected data and dynamics of measles to produce effective and efficient control strategies. A spatial pattern of misaligned disease data was investigated to estimate and map the risk of measles by utilizing ecological Bayesian regression models in Namibia [18]. A discrete-time seasonal model of measles transmission was proposed by Zhenguo Bai and Liu in 2015. They described the seasonal fluctuation of measles cases in China [19]. Muhammad Farman et al. have presented analysis and numerical simulations of an SEIR model represented the measles transmission with non-integer time-fractional derivatives in 2017 and considered the measles control [20]. An SEIR epidemic model was performed to investigate the impact of exposed individuals in the dynamics of measles transmission [21]. The optimal policy for controlling measles was examined through an SIRV model using supplemental immunization activities with separation of vaccinated compartment [22]. In 2018, O.J. Peter et al. proposed an SVEIR model of measles transmission by considering natural recovery and recovery due to vaccination [23].
The progress in measles transmission modeling has given challenges for researchers in shaping up the health policy decisions. This paper aims to construct two measles transmission models involving vaccination and to investigate the most influential parameters for proposing a control strategy. The proposed model introduces a hospitalized compartment wherein infected individuals entering hospitals are considered isolated and cannot transmit the infection to the susceptible group. First, we presented the weekly data and the measles progress in the data section. We then extracted measles incidence data to be the number of incidences per total population. Second, two mathematical models were constructed, and then afterward, the equilibrium points, as well as their stability conditions, were analyzed. Further, there exist oscillatory solutions if the parameters are satisfying a certain condition. Third, the unobserved parameters and initial conditions were estimated from the data along their feasible regions. Fourth, we did a sensitivity analysis to analyze the effect of changing parameter values and the most influential parameter towards the dynamics. We then conducted several simulations of control strategies to minimize the number of infected by taking different vaccinated proportions in the last section of this paper.
Following the policy of the Ministry of Health in achieving the elimination target for measles in 2020. the Jakarta Provincial Health Office has targeted for no endemic area of measles for more than 12 months with the implementation of adequate measles surveillance. The indicator is to conduct measles surveillance based on individual cases or CBMS (Case-Based Measles Surveillance), which has been carried out gradually since 2014. In this surveillance, for every case found with symptoms of fever, rash on the body, accompanied by one or more symptoms of cough, runny nose, and red-eye, blood specimens (serum) from the individual are examined at a national reference laboratory to confirm the diagnosis of Measles. Another indicator of the surveillance is discarded ratio (misclassification of measles) obtained from laboratory tests [24].
Based on a leaflet from Jakarta Provincial Health Office, hospitals are instructed to identify people groups who are vulnerable to measles by considering the mapping of cases and immunization status at the rural level. It is implemented as a global target for measles elimination and Rubella's control in 2020. Besides, several attempts were made to follow up on each reported case in a health service facility. They are examining, treating, and verifying alerts by carrying out epidemiological investigations in less than 24 hours. Then they are taking and sending specimens for laboratory examinations, reporting each case of measles to Jakarta Health Office surveillance website by submitting individual data (form C1), and monitoring measles cases in the local area or various sources of infection, including cases reported at the local hospital.
Achievement of the measles target is carried out by laboratory sampling in 2017 amounted to 77%, 2018 amounted to 75%, and until July 2019 amounted to 70% based on the evaluation results of measles program holder in Jakarta. Even though the proportion of measles until July 2019 had only reached 70%, the achievement of discarded measles in Jakarta had exceeded 2 per 100,000 population, equal to 4.17 per 100,000 population [25].
In this study, we validated our model with actual data of measles in Jakarta. The data recorded here is the number of hospitalized individuals. We chose Jakarta as our observed area due to the most available of completed data collection compared to other cities in Indonesia. Jakarta Health Office has collected measles data from 44 sub-district health centers, 283 village health centers, and 152 hospitals in Jakarta, which has verified and passed the validation process. Here, the data used is a weekly basis of incidence data started from the first to the last week in 2017, see Figure 1. The left figure is a plot of time series measles data in all districts of Jakarta. The right one is the weekly total incidence of measles in Jakarta, produced by adding up the data in the left data plot.
We developed a standard SIR model (Kermack–McKendrick epidemic model[26]) for measles transmission involving vaccination and the addition of a hospitalized compartment to analyze the dynamics of measles transmission. Hereafter a stochastic model was constructed from a theory of stochastic differential equations containing a Wiener process.
We formulated a deterministic model (SIHR) to describe the measles transmission dynamic. The total population is assumed to be constant in this model. To illustrate the model equations, we divide the total population into four compartments. Let ˉS,ˉI,ˉH, and ˉR are the number of susceptible, infected, hospitalized, and recovered individuals including individuals who were successfully vaccinated, respectively.
The susceptible compartment ˉS increased by birth at a rate of Λ. The individuals probably got their infection by contact with infected individuals at rate β. Vaccination process is given against the susceptible population ˉS at rate rp, where r and p are a probability of success in preventing measles and the proportion of vaccinated individuals, respectively. Consequently, the susceptible individuals might be moved to the infected compartment ˉI at rate β(1−p)+β(1−r)p=β(1−rp). Then, the infected individuals tend to go to the hospital at a rate of δ for getting treatment. The infected and hospitalized individuals recovered naturally at a rate of γ1 and due to treatment at a rate of γ2, respectively. Logically, it should be γ2>γ1. All compartments decreased by a natural death at a rate of μ.
The transmission dynamics of measles can be expressed by a system of differential equations as follows.
dˉS(t)dt=Λ−βN(1−rp)ˉS(t)ˉI(t)−(rp+μ)ˉS(t),dˉI(t)dt=βN(1−rp)ˉS(t)ˉI(t)−(δ+γ1+μ)ˉI(t),dˉH(t)dt=δˉI(t)−(γ2+μ)ˉH(t),dˉR(t)dt=γ1ˉI(t)+γ2ˉH(t)+rpˉS(t)−μˉR(t), | (3.1) |
where dˉN(t)dt=Λ−μN(t). We assumed that the population is constant then Λ=μN,N∈N. We then normalized the variables by substituting S(t)=ˉS(t)N,I(t)=ˉI(t)N,H(t)=ˉH(t)N,R(t)=ˉR(t)N. Consequently, it reduced the system (3.1) into three-dimensional space in (3.2) by ignoring the recovered compartment due to the fact that recovered and successfully vaccinated individuals did not contribute to the transmission.
dS(t)dt=μ−β(1−rp)S(t)I(t)−(rp+μ)S(t),dI(t)dt=β(1−rp)S(t)I(t)−(δ+γ1+μ)I(t),dH(t)dt=δI(t)−(γ2+μ)H(t). | (3.2) |
Since, the parameters and state variables are non negative ie t≥0, we investigated model (3.2) in the following feasible region:
Ω={(S(t),I(t),H(t))∈R3+|0≤S(t),I(t),H(t)≤1,∀t} |
All of the parameters used here are positive and described in Table 1. The non-negative octant Ω is positively invariant of model (3.2). This follows immediately evaluating normal vector at each plane S=0,I=0, and H=0, respectively. The model has two equilibria; we called it as a disease-free equilibrium (DFE) E0=(S0,I0,H0)=(μrp+μ,0,0) and a endemic equilibrium E1 that can be expressed as follows.
S∗=μR0(rp+μ),I∗=(R0−1)μ(δ+γ1+μ)R0,H∗=δ(R0−1)μ(δ+γ1+μ)R0(γ2+μ), | (3.3) |
No. | Parameters | Description | Value | Reff. |
1. | μ | Death rate and birth rate | 152∗70 | [27] |
2. | β | Infection rate | Estimated | - |
3. | r | Probability of succesed vaccination | 90% | [28] |
4. | p | Proportion of vaccinated individuals | Varied | - |
5. | δ | Hospitalized rate | Estimated | - |
6. | γ1 | Recovery rate of infected | 1/3 | [10] |
7. | γ2 | Recovery rate of hospitalized | 1 | - |
where R0=β(1−rp)μ(μ+rp)(δ+γ1+μ) is a basic reproduction number for the model. It can be defined as the spectral radius of matrix ρ(FV−1) which is
F=[β(1−rp)μrp+μ000] and V=[δ+γ1+μ0−δγ2+μ]. |
Further, we give the proof of stability of the equilibria both locally and globally, in Theorem 3.1.1 and 3.1.2. The existence of oscillatory solutions is also presented in Proposition 3.1.1.. Although it is trivial, we stated it explicitly to support simulations and to provide easy understanding.
Theorem 3.1.1 Local stability of the equilibria.
1. If R0<1, then the disease-free equilibrium E0 is locally asymptotically stable in Ω.
2. If R0>1, then the endemic equilibrium E1 is locally asymptotically stable in Ω.
Proof.
1. The Jacobian matrix of model (3.2) at E0 is
J(E0)=[−(rp+μ)−β(1−rp)μpr+μ00β(1−rp)μpr+μ−(δ+γ1+μ)00δ−(γ2+μ)]. |
The eigen values of the matrix are −(μ+γ2),−(rp+μ), and (δ+γ1+μ)(R0−1). If R0<1 then all eigen values are real negative, implies the disease-free equilibrium E0 is locally asymptotically stable in Ω.
2. The Jacobian matrix of model (3.2) at E0 is
J(E1)=[−β(1−rp)(R0−1)μ(δ+γ1+μ)R0−(rp+μ)−β(1−rp)μ(rp+μ)R00β(1−rp)(R0−1)μ(δ+γ1+μ)R0β(1−rp)μ(pr+μ)R0−(δ+γ1+μ)00δ−(γ2+μ)]. |
The characteristic polynomial of J(E1) is given by f(λ)=λ3+Aλ2+Bλ+C where A=(pr+μ)R0+(μ+γ2),B=(R0−1)(pr+μ)(δ+γ1+μ)+R0(pr+μ)(μ+γ2), and C=(γ2+μ)(R0−1)(pr+μ)(δ+γ1+μ). If R0>1 then coefficients of the polynomial are positive and AB>C. It indicates that Routh-Hurwitz conditions for the polynomial are satisfied. Therefore the endemic equilibrium E1 is locally asymptotically stable in Ω.
Theorem 3.1.2 Global stability of the equilibria.
1. If R0≤1, then the disease-free equilibrium E0 is globally asymptotically stable in Ω.
2. If R0>1, then the endemic equilibrium E1 is globally asymptotically stable in the sub-region Ω2⊆Ω that satisfy H≥max{H∗,δIˆγ} or H≤min{H∗,δIˆγ}.
Proof.
1. We first considered the Lyapunov function V(S,I,H)=S0(δ+γ1+μ)I+1δH and got
˙V=S0δ+γ1+μ˙I+1δ˙H=S0(δ+γ1+μ)[β(1−rp)SI−(δ+γ1+μ)I]+1δ[δI−(γ2+μ)H]=[R0S−1]I−S0I−γ2+μδH, |
which is less than zero in Ω if R0≤1. A dot represents differentiation with respect to time. Furthermore, ˙V=0 if I=I0 and H=H0. Hence, the largest compact invariant set contained in the set Ψ={(S,I,H)∈Ω,˙V(S,I,H)=0} when R0<1 is the singleton E0. Therefore, the DFE is globally asymptotically stable in Ω based on the LaSalle invariance principle theorem in [29].
2. For simplifying this proof, we assumed that ˆβ=β(1−rp),ˆα=rp+μ,ˆδ=δ+γ1+μ and ˆγ=γ2+μ. Consequently, the system (3.2) can rewrite as follows.
˙S=μ−ˆβSI−ˆαS˙I=ˆβSI−ˆδI˙H=δI−ˆγH, | (3.4) |
where R0=ˆβμˆαˆδ,S∗=μR0ˆα,I∗=(R0−1)μR0ˆδ, and H∗=δ(R0−1)μˆγR0ˆδ. We then considered a Lyapunov function
U(S,I,H)=S−S∗−S∗ln(SS∗)+I−I∗−I∗ln(II∗)+H−H∗−H∗ln(HH∗). |
Clearly U:Ω→R and continuously differentiable. In addition, U(S∗,I∗,H∗)=0 while U(S,I,H)>0,∀(S,I,H)≠(S∗,I∗,H∗). The derivative with respect to the time can be derived as follows.
˙U=(1−S∗S)˙S+(1−I∗I)˙I+(1−H∗H)˙H=(1−S∗S)(μ−ˆβSI−ˆαS)+(1−I∗I)(ˆβSI−ˆδI)−(1−H∗H)(δI−ˆγH)=(μ−ˆαS∗SS∗−μS∗S+ˆβS∗I+ˆαS∗−ˆδI−ˆβI∗S+ˆδI∗)−1H(H−H∗)(δI−ˆγH)=(μ−μR0SS∗−μS∗S+ˆβμR0ˆαI+μR0−ˆδI−ˆβμ2(R0−1)R20ˆαˆδSS∗+μ(R0−1)R0)−1H(ˆγH2−δIH−ˆγH∗H+δH∗I)=(μ−ˆαˆδˆβSS∗−μS∗S+ˆδI+ˆαˆδˆβ−ˆδI−μSS∗+ˆαˆδˆβSS∗+μ−ˆαˆδˆβ)−ˆγH(H2−(δIˆγ+H∗)H+δH∗Iˆγ)=(μ−μS∗S−μSS∗+μ)+(−ˆαˆδˆβSS∗+ˆαˆδˆβ+ˆαˆδˆβSS∗−ˆαˆδˆβ)−ˆγH(H(H−δIˆγ−H∗)+δH∗Iˆγ)=μ(2−S∗S−SS∗)+ˆαˆδˆβ(−SS∗+1+SS∗−1)−ˆγH(H(H−δIˆγ−H∗)+δH∗Iˆγ)=−μS∗S((SS∗)2−2SS∗+1)−ˆγH(H(H−δIˆγ−H∗)+δH∗Iˆγ)=−μS∗S(1−SS∗)2−ˆγH(H−H∗)(H−δIˆγ). |
If R0>1 and H≥max{H∗,δIˆγ} or H≤min{H∗,δIˆγ} are satisfied then ˙U≤0. Therefore, E1 is a globally asymptotically stable within Ω2.
Proposition 3.1.1 There exists oscillatory solutions of system (3.2) to the endemic equilibrium if R0>1 and I∗>μ(rp+μ)R04(δ+γ1+μ)2.
Proof. The characteristic polynomial of the Jacobian matrix at E1 can be used as a sign of the existence of oscillatory solutions around the equilibrium of E1. Remember that we can check it through its discriminant. The eigenvalue of the matrix has non-zero imaginary parts if the discriminant Δ=A3C27−A2B2108−ABC6+B327+C24 is bigger than zero [30,31]. Therefore, we have Δ=−1108(pr+μ)Q21Q2 where Q1=R0(pr+μ)(δ+γ1−γ2)+μ(2γ2−pr−δ−γ1)−pr(δ+γ1)+γ22≥0 and Q2=R20(pr+μ)+(1−R0)4(μ+δ+γ1).
If R0>1 and I∗>μ(pr+μ)R04(δ+γ1+μ)2 then Q2<0, implies Δ>0.
In this part, we constructed a system of Itˆo SDEs (stochastic differential equations) with continuous-valued random variables by utilizing the probability transition rates in Table 2 and assuming △t (time step) approach zero [32,33]. A diffusion approximation is a designation of the process resulting from the SDE model by Kurtz in 1978 [34]. Let y(t) denote the number of compartments in the population at generation t, x(t) denote the proportion of compartments in a population, and x(t)=y(t)/N. The random vectors x(t)=(s(t)i(t)h(t))T and y(t)=(ˉs(t)ˉi(t)ˉh(t))T are vectors of continuous random variables with state space s(t),i(t),h(t)∈[0,1] and ˉs(t),ˉi(t),ˉh(t)∈[0,N],t∈[0,∞).
i | Change, (△y)i | Probability, Pi |
1 | [1,0,0]T | μN△t |
2 | [−1,1,0]T | βN(1−rp)ˉs(t)ˉi(t)△t |
3 | [−1,0,0]T | rpˉs(t)△t |
4 | [−1,0,0]T | μˉs(t)△t |
5 | [0,−1,1]T | δˉi(t)△t |
6 | [0,−1,0]T | γ1ˉi(t)△t |
7 | [0,−1,0]T | μˉi(t)△t |
8 | [0,0,−1]T | γ2ˉh(t)△t |
9 | [0,0,−1]T | μˉh(t)△t |
Based on the transitions and rates in Table 2, the corresponding SDE model can be written as the form
dx=ˉf(x(t))dt+G(x(t))dˉW(t), | (3.5) |
where G is a diffusion matrix and ˉW(t)=(ˉW1(t),ˉW2(t),⋯,ˉW9(t))T is a 9×1 matrix (vector) of independent Wiener processes. To order △t,
ˉf(x(t))△t=E(△x(t)|x(t))=1NE(△y(t)|x(t)), |
where ˉf is drift vector and
GGT△t=M△t=E(△x(t)(△x(t))T|x(t))=1N2E(△y(t)(△y(t))T|x(t)). |
Next, we computed the expectation and the covariance matrix based on Table 2. For a sufficiently small △t, the expectation
E(△x)=[μ−β(1−rp)s(t)i(t)−(rp+μ)s(t)β(1−rp)s(t)i(t)−(δ+γ1+μ)i(t)δi(t)−(γ2+μ)h(t)]△t. | (3.6) |
The covariance matrix can be written as
M=1N[m11m120m12m22m230m23m33]△t, | (3.7) |
where m11=μ+β(1−rp)s(t)i(t)+(rp+μ)s(t), m22=β(1−rp)s(t)i(t)+(δ+γ1+μ)i(t), m33=δi(t)+(γ2+μ)h(t), m12=−β(1−rp)s(t)i(t) and m23=−δi(t). The diffusion matrix G(x(t)) is given by
1√N[√μ−√−m12−√rps(t)−√μs(t)000000√−m1200−√δi(t)−√γ1i(t)−√μi(t)000000√δi(t)00−√γ2h(t)−√μh(t)] |
such that GGT=V.
We defined an optimization problem for parameter estimation in model (3.2). For simplification of the optimization problem, the model can be rewritten
˙X(t)=f(X(t),t,ϑ),X(t0)=X0 | (3.8) |
where X(t)=(S(t)I(t)H(t))T and ϑ={β,δ,I(0),H(0)} is sets of estimated parameters and initial conditions. In this case, we have two unobserved parameters, β and δ, and two unknown initial conditions, I(0) and H(0). p is vaccination related parameter, which is not included into consideration for parameter estimation. This parameter is then used to compare the impact of vaccination against the measles epidemic. Hence, the optimization problem can be written as follows.
minimizeθ:Ψ(x(t),ϑ)=√n∑t=1(H(t)−Data(t))2subjectto:system(3.8),0≤I(0)≤1,0≤H(0)≤1,β≥0,0≤δ≤1. | (3.9) |
Ψ is squared residual norm (resnorm). The main idea of parameter estimation is how to find the best combination of ϑ such that minimizing the resnorm. In this case, we used nonlinear least squares method to look for the best combination of ϑ.
In this part, we discuss a local sensitivity analysis. It gives us information on which parameters are the most influential to the dynamics of the system, with the initial value of the given parameter. We investigate the sensitivity of parameters that are meaningful in the biological interpretation of measles transmission. The idea of this sensitivity analysis is to determine the effect of changing parameter values on the solution of the system ˙X=f(t,X,θ) in Eq (3.8) where X(t,θ)∈R3 depend on time t and θ is the set of parameters, θ={β,p,δ,γ1,γ2}∈R5. We first define a sensitivity function of X(t,θ) as follows.
ζ(t,θ)=(ζik(t,θ))=ddθiXk(t,θ). | (3.10) |
where i=1,...5,k=1,2,3, and ζ is time dependent t [35]. By differentiating ζ with respect to time t, then the first order solution sensitivity is a solution of differential equation,
˙ζ=ddtdXdθ=ddθdXdt=ddθf(t,X,θ)=∂X∂θ⋅df(t,X,θ)dX+∂f(t,X,θ)∂θ=ζ⋅J(t,X,θ)+Q(t,X,θ) | (3.11) |
with ζ(t=0)=0. J(t,X,θ) and Q(t,X,θ) are 3-by-3 Jacobian matrix and 5-by-3 matrix, respectively. The entries of matrix Q are derivatives of the right side of Eq (3.10) with respect to the parameters.
In the past decade, some countries in the European Union (EU) are burdened with cases of high measles due to low immunization coverage [36]. Several factors potentially reduced coverage of the measles vaccine, such as low family income [37], parent education, lack of access to health care networks, migration [37,38,39,40], and beliefs about the risk of vaccination [41]. Therefore, high vaccine coverage is needed to prevent outbreaks in the general population as well as in the poorly pooled population [42]. The effectiveness of immunization also depends on reducing disparities associated with socioeconomic levels and access to health care [43,44].
Measles vaccine is available and circulating in society. Jakarta provincial office also supports measles vaccination programs to minimize the cases as much as possible by implementing several government programs that have been launched. There are three stages of vaccine implementation in Jakarta. The stages of vaccination are given to susceptible individuals aged 9–11 months, 18–24 months, and 7–15 years[45]. However, the efficiency of this intervention program has not been fully assessed. Hence it is necessary to analyze it to get the best strategy that is efficient and adjusted for the budget. Here, we discussed the optimal way to control measles transmission via model (3.2) by a term of vaccination p. The parameter values obtained from the estimation parameter are used in this simulation by varying the value of p to observe the impact of vaccination.
We confirmed the constructed model based on the above analysis results and the actual data in this part. The numerical method used in this simulation is Runge-Kutta method. First, the oscillatory solution was performed by taking parameter values from Table 1, β=2,p=0, and δ=0.5. The other values are taken according to either references or assumptions, see Table 1. Second, the solution oscillates and tends to the equilibrium point, see Figure 2. It is a numerical simulation of Proposition 3.1.1 about the existence of oscillatory solutions around E1.Third, we estimated parameters in the model satisfying the minimum error between data and hospitalized solution and used the parameters for stochastic simulation. There are two unobserved parameters in model (3.2), they are β and δ. For the initial condition, we assumed N=9.5⋅106 as a total population in Jakarta, S(0)=1−(I(0)+H(0)),R(0)=0. The unknown initial conditions, H(0) and I(0) are estimated together with the unobserved parameters. Fourth, we presented a study of sensitivity analysis and then performed control a strategy which is simulated with several levels of immunization coverage.
The result of parameter estimation shows that ϑ converges to β=0.67,δ=0.3583,I(0)=2.0157⋅10−5,H(0)=2⋅10−6 with a resnorm of 1.0486⋅10−10, see Figure 4(a). The residual plot and fitting the distribution of the residual are shown in Figure 4(b) and Figure 4(c), respectively. The results show that the absolute errors are mostly below 2⋅10−6, indicated in the yellow area, which can be interpreted that the average estimated error is under two people per million population. Further, the residual follows a normal distribution with mean equal to zero and 1.4242⋅10−6 variance by using a one-sample t−test at significance level 99 %. Afterward, we used the estimated parameter for stochastic simulation with △t=10−4. The results are shown in Figure 5. The simulation results show that the data are within the 95 % confidence interval of the stochastic model and the average of the stochastic solutions across a thousand runs is relatively close to the deterministic solution.
In sensitivity analysis, we are interested in observing the parameter sensitivity under endemic conditions (R0>1) such that the sensitivity analysis was conducted locally in neighborhood of θ⋆ where,
θ⋆={β⋆=2,p⋆=2×10−7,δ⋆=0.5,γ⋆1=1/3,γ⋆2=1}. |
The results of this sensitivity analysis are shown in Figure 6 by taking fixed biological parameters μ=1/(52⋅70) and r=0.9. We then compared which parameters are dominantly contributing to the change of the solutions. The parameters discussed here for sensitivity analysis are β,p,δ,γ1, and γ2, respectively. Recalling the Eq (3.10), if ζ is large enough, it means that a small change in parameter θi makes a big contribution to Xk. Therefore, significant parameters have large values of ζ. A positive or negative value of ζ indicates the effect which is directly or inversely proportional, respectively.
In Figure 6(a)–(c), the most sensitive parameter to the susceptible population is p. The second and third are δ and γ1. They have a positive effect on the susceptible population even though changes to these parameters do not greatly affect the susceptible population. It indicates that when the parameters increase, then the susceptible population also increases. Unlike the β parameter, this parameter is the fourth sensitive parameter after δ and γ1. But the effect of changing the parameter β and p have an opposite effect on the population. It means that by increasing infection rates or the proportion of vaccinated individuals, it decreases the susceptible population. Meanwhile, the parameter γ2, the effect of the changes is not significant on the susceptible population.
For the sensitivity of the infected population with respect to the parameters, δ and γ1 are the most sensitive parameter and have a opposite effect on the infected population. The parameter p is more sensitive than the parameter β. But the parameter p has a opposite impact on the infected population, see Figure 6(d)–(f). It can be interpreted that when the rate of treatment to the hospital increases, the infected population decreases. The government (health department) can issue policies such as the health insurance program to increase the rate of treatment at the hospitals to reduce the number of infected populations. Natural healing is more sensitive than vaccination parameters. It indicates that when an outbreak is occurring, providing treatment accesses easier for infected individuals is better than vaccinating.
The most sensitive parameter in the hospitalized population is the recovery rate caused by treatments (γ2). The second one is the hospitalized rate (δ). The recovery rates, both natural healing, and healing caused by treatments have an opposite effect on the hospitalized population. It can be seen that the increase in the recovery rate will further reduce the number of hospitalized population. The infection rate (β) has a positive effect on the hospitalized population even though the sensitivity level is lower than p, see Figure 6(g)–(i). Vaccination parameter (p) has an opposite impact, which means that giving vaccines with a large proportion will reduce the number of hospitalized population.
Jakarta has about 2.5 million children under 15 years old or around 26% of the total population. They are the most dominant measles vaccination target. Therefore, the variation of vaccination coverage is less than 30% in this simulation. Variation of the level of vaccine coverage is represented by the parameter p and the results are presented in Figure 7. The far-left graphics represent the dynamics infected with respect to the vaccine coverage 0.2% to 1% and 2% to 6% in the right image. Both provide information that increasing vaccine coverage, the number of cases is getting smaller. Nevertheless, this level of vaccine coverage must be adjusted and synchronized with the intended budget to overcome this disease.
We have constructed deterministic and stochastic models for a measles transmission with vaccination. Two equilibria are obtained with their stability, both locally and globally. Oscillatory solutions are also found around the endemic equilibrium when certain conditions are satisfied, are shown in Figures 2 and 3. In the case of parameter estimation, the residuals are mostly less than two people per one million population and follow a normal distribution with a mean of 0 at the significance level of 99%. The stochastic model is shown to capture the randomness of the transmission, represented by independent Wiener process variables and diffusion matrices. It is caused by the influence of other factors in which not all of them can be accommodated into the model. However, the simulation results show that the actual data are within the 95% confidence interval of the stochastic model, which means that the model can reach the random transmission of this disease contained in the data (see Figure 5). The average of the stochastic solutions across a thousand runs is relatively close to the deterministic solution. Moreover, the sensitivity analysis results show that the most influential parameters in the infected population are hospitalized rate and natural recovery rate, then followed by the vaccinated proportion. It can be a valuable information for decision-makers that improving access to health care networks is one of the dominant factors in reducing the number of measles cases. Vaccine coverage is also very influential in reducing the spread of the virus as a preventive measure. We hence proposed a control strategy which is providing treatment accesses easier is better than vaccinating when an outbreak occurs.
There are several limitations to this study. The first problem is the unique solution for parameter estimation. It brings us to the global optimization problem in high dimensional space, which is still a challenge for researchers in the optimization field. The assumption that the total population is almost constant over time and some uniformity assumptions on parameter values taken are also limitations of this study. In the case of sensitivity analysis, it is a local sensitivity in the neighborhood of θ∗, which is time-dependent dynamics, and its sensitivities have adhered to the initial parameters.
MF and HF gratefully acknowledge the financial support provided by Indonesia Ministry of Research and Technology through PMDSU Program [grant No. 1511/E4.4/2015]. Parts of this work of ES and NN are funded by ITB Research Grant 2019 and ITB Research, Community Service, and Innovation Program (P3MI-ITB) 2019, respectively [grant No. 1000M/I1.C01/PL/2019].
MF performed the mathematical analysis, stochastic model, and constructed the manuscript. DS did the mathematical analysis and sensitivity analysis. S collected and interpreted data. HF did the sensitivity analysis. NN reviewed the manuscript, control strategy, and explained the results. ES reviewed the manuscript and analyzed the mathematical method.
The authors declare there is no conflicts of interest.
[1] |
Anand KS, Vikas D (2012) Hippocampus in health and disease: An overview. Ann Indian Acad Neurol 15: 239-246. https://doi.org/10.4103/0972-2327.104323 ![]() |
[2] |
Sweatt JD (2016) Neural plasticity and behavior–sixty years of conceptual advances. J Neurochem 139: 179-199. https://doi.org/10.1111/jnc.13580 ![]() |
[3] |
Arakawa H, Arakawa K, Blanchard DC, et al. (2008) A new test paradigm for social recognition evidenced by urinary scent marking behavior in C57BL/6J mice. Behav Brain Res 190: 97-104. https://doi.org/10.1016/j.bbr.2008.02.009 ![]() |
[4] |
Silverman JL, Yang M, Lord C, et al. (2010) Behavioural phenotyping assays for mouse models of autism. Nat Rev Neurosci 11: 490-502. https://doi.org/10.1038/nrn2851 ![]() |
[5] |
Pintana H (2012) Effects of metformin on learning and memory behaviors and brain mitochondrial functions in high fat diet induced insulin resistant rats. Life Sci 91: 409-414. https://doi.org/10.1016/j.lfs.2012.08.017 ![]() |
[6] |
Martinos MM, Yoong M, Patil S, et al. (2012) Recognition memory is impaired in children after prolonged febrile seizures. Brain 135: 3153-3164. https://doi.org/10.1093/brain/aws213 ![]() |
[7] |
Blake RV, Wroe SJ, Breen EK, et al. (2000) Accelerated forgetting in patients with epilepsy: evidence for impairment in memory consolidation. Brain 123: 472-483. https://doi.org/10.1093/brain/123.3.472 ![]() |
[8] |
Scola C, Bourjade M, Jover M (2015) Social interaction is associated with changes in infants' motor activity. Socioaffect Neurosci Psychol 5: 1. https://doi.org/10.3402/snp.v5.28256 ![]() |
[9] |
Couture SM, Penn DL, Losh M, et al. (2010) Comparison of social cognitive functioning in schizophrenia and high functioning autism: more convergence than divergence. Psychol Med 40: 569-579. https://doi.org/10.1017/S003329170999078X ![]() |
[10] | Rubin R, Strayer DS (2008) Environmental and Nutritional pathology, Rubins pathology; Clinicopathologic Foundations of Medicine. Lippincot Williams and Wilkins 57-69. |
[11] | Owalabi J, Williams F, Fabiyi O (2014) Evaluation of moringa's effects against lead induced disruption of the hippocampus in animal models. World J Life Sci Med Res 3: 39-45. |
[12] |
Kaidanovich-Beilin O, Lipina TV, Takao K, et al. (2009) Abnormalities in brain structure and behavior in GSK-3alpha mutant mice. Mo Brain 2: 35. https://doi.org/10.1186/1756-6606-2-35 ![]() |
[13] | Iliyasu MO, Ibegbu AO, Sambo JS, et al. (2015) Histopathological changes on the hippocampus of adult wistar rats exposed to lead acetate and aqueous extract of psidium guajava leaves. Afr J Cell Pathol 5: 26-31. |
[14] | Rachael DR, Patrick DW, Melissa CD, et al. (2014) The role of hippocampus in flexible cognition and social behavior. Front Hum Neurosci 8: 742. https://doi.org/10.3389/fnhum.2014.00742 |
[15] |
Wani AL, Ara A, Usmani JA (2015) Lead toxicity: a review. Interdiscipl Toxicol 8: 55-64. https://doi.org/10.1515/intox-2015-0009 ![]() |
[16] |
Zhang H, Liu Y, Zhang R, et al. (2014) Binding Mode Investigations on the Interaction of Lead (II) Acetate with Human Chorionic Gonadotropin. J Phys Chem B 118: 9644-9650. https://doi.org/10.1021/jp505565s ![]() |
[17] | Wasswa JN, Omorodion TN, Avwioro GO, et al. (2017) Histological effect of piper guineense (uziza) leaves on the liver of wistar rats. Int J Res Rev 4: 36-41. |
[18] | Kummer K, Klement S, Eggart V, et al. (2011) Conditioned place preference for social interaction in rats: contribution of sensory components. Front Behav Neurosci 5: 80. |
[19] |
Yang B, Ren Q, Ma M, et al. (2016) Antidepressant effects of (C)-MK-801 and MK-801 in the social defeat stress model. Int J Neuropsychopharmacol 19: pyw080. https://doi.org/10.1093/ijnp/pyw080 ![]() |
[20] | Kaidanovich-Beilin O, Lipina T, Vukobradovic I, et al. (2011) Assessment of Social Interaction Behaviors. J Vis Exp 48. https://doi.org/10.3791/2473 |
[21] |
Tariq M, Khan HA, Deeb SAl, et al. (2005) Neuroprotective Effect of Nicotine Against 3-Nitropropionic Acid (3-Np)-Induced Experimental Huntington's Disease in Rats. Brain Res Bull 67: 161-168. https://doi.org/10.1016/j.brainresbull.2005.06.024 ![]() |
[22] | Cole TG, Klotzsch SG, Namara MC (1997) Measurement of triglyceride concentration. Handbook of lipoprotein testing . Washington: AACC Press 115-26. |
[23] |
Okada M, Matsui H, Ito Y, et al. (1998) Low-density lipoprotein cholesterol can be chemically measured: a new superior method. J Lab Clin Med 132: 195-201. https://doi.org/10.1016/S0022-2143(98)90168-8 ![]() |
[24] |
Ochei JO, Enitan SS, Effedua HI, et al. (2017) Libido enhancement potential of piper guineense in male Wistar rats. Asian J Bio 4: 1-9. https://doi.org/10.9734/AJOB/2017/36542 ![]() |
[25] |
Ercal N, Gurer-Orhan H, Aykin-Burns N (2001) Toxic metals and oxidative stress. Part1. Mechanisms involved in metal induced oxidative damage. Curr Top Med Chem 1: 529-539. https://doi.org/10.2174/1568026013394831 ![]() |
[26] |
Lucke-Wold BP, Naser ZJ, Logsdon AF, et al. (2015) Amelioration of Nicotinamide Adenine Dinucleotide Phosphate-Oxidase Mediated Stress Reduces Cell Death after Blast-Induced Traumatic Brain Injury. Transl Res 166: 509-528.e1. https://doi.org/10.1016/j.trsl.2015.08.005 ![]() |
[27] |
Sharma P, Chambial S, Shukla KK (2016) Lead and neurotoxicity. Indian J Clin Biochem 30: 1-2. https://doi.org/10.1007/s12291-015-0480-6 ![]() |
[28] | Samarghandian ST, Farkhondeh M, Azimi-Nezhad AM, et al. (2019) Antidotal or protective effects of honey and chrysin, its major polyphenols, against natural and chemical toxicities. Acta Biomed 90(4): 533-550. |
[29] | Amadi G Effects of methanolic leaf extract of piper guineense on some biochemical parameters of diabetic female albino rats (2018). Available at https://futospace.futo.edu.ng and accessed on 202 April 11th |
[30] |
Blake GJ, Otvos JD, Rifai N, et al. (2002) Low-density lipoprotein particle concentration and size as determined by nuclear magnetic resonance spectroscopy as predictors of cardiovascular disease in women. Circulation 106: 1930-1937. https://doi.org/10.1161/01.CIR.0000033222.75187.B9 ![]() |
[31] | Yankelevitch-Yahav R, Franko M, Huly A, et al. (2015) The Forced Swim Test as a Model of Depressive-Like Behavior. J Vis Exp : 52587. https://doi.org/10.3791/52587 |
[32] | Can A, Dao DT, Arad M, et al. (2012) The Mouse Forced Swim Test. J Vis Exp : e3638. https://doi.org/10.3791/3638 |
[33] |
Weilch KD, Pfister JA, Lima FG (2013) Effect Of Α7 Nicotinic Acetylcholine Receptor Agonists and Antagonists on Motor Function in Mice. Toxicol Appl Pharm 266: 366-374. https://doi.org/10.1016/j.taap.2012.11.024 ![]() |
[34] | Ijomone OM, Olaibe OK, Biose IJ, et al. (2014) Performance of motor associated behavioral tests following chronic nicotine administration. Ann Neurosci 21: 42-46. https://doi.org/10.5214/ans.0972.7531.210203 |
[35] |
Kropff E, Yang SM, Schinder AF (2015) Dynamic role of adult born dentate granule cells in memory processing. Curr Opin Neurobiol 35: 21-26. https://doi.org/10.1016/j.conb.2015.06.002 ![]() |
[36] |
Moy SS, Nadler JJ, Perez A, et al. (2004) Sociability and preference for social novelty in five inbred strains: an approach to assess autistic-like behavior in mice. Genes Brai Behav 3: 287-302. https://doi.org/10.1111/j.1601-1848.2004.00076.x ![]() |
[37] |
Ahmad F, Salahuddin M, Alsamman K, et al. (2018) Developmental lead (Pb)-induced deficits in hippocampal protein translation at the synapses are ameliorated by Ascorbate supplementation. Neuropsych Dis Treat 14: 3289-3298. https://doi.org/10.2147/NDT.S174083 ![]() |
[38] |
Musa SA, Omoniye MI, Hamman WO, et al. (2012) Preventive activity of Ascorbic acid on lead acetate induced cerebellar damage in adult wistar rats. Med Health Sci J 13: 99-104. https://doi.org/10.15208/mhsj.2012.68 ![]() |
[39] | Hsiang J, Díaz E (2011) Lead and developmental neurotoxicity of the central nervous system. Curr Neurobiol 2: 35-42. |
[40] |
Lucke-Wold BP, Turner RC, Logsdon AF, et al. (2016) Endoplasmic reticulum stress implicated in chronic traumatic encephalopathy. J Neurosurg 124: 687-702. https://doi.org/10.3171/2015.3.JNS141802 ![]() |
1. | Feng Ding, Ling Xu, Xiao Zhang, Yihong Zhou, Xiaoli Luan, Recursive identification methods for general stochastic systems with colored noises by using the hierarchical identification principle and the filtering identification idea, 2024, 57, 13675788, 100942, 10.1016/j.arcontrol.2024.100942 | |
2. | Guiyao Ke, Jun Pan, Feiyu Hu, Haijun Wang, Dynamics of a New Four-Thirds-Degree Sub-Quadratic Lorenz-like System, 2024, 13, 2075-1680, 625, 10.3390/axioms13090625 | |
3. | Ling Xu, Huan Xu, Feng Ding, Adaptive Multi-Innovation Gradient Identification Algorithms for a Controlled Autoregressive Autoregressive Moving Average Model, 2024, 43, 0278-081X, 3718, 10.1007/s00034-024-02627-z | |
4. | Siyu Liu, Yanjiao Wang, Feng Ding, Ahmed Alsaedi, Tasawar Hayat, Joint iterative state and parameter estimation for bilinear systems with autoregressive noises via the data filtering, 2024, 147, 00190578, 337, 10.1016/j.isatra.2024.01.035 | |
5. | Haijun Wang, Jun Pan, Guiyao Ke, Feiyu Hu, A pair of centro-symmetric heteroclinic orbits coined, 2024, 2024, 2731-4235, 10.1186/s13662-024-03809-4 | |
6. | Lijuan Liu, Decomposition‐based maximum likelihood gradient iterative algorithm for multivariate systems with colored noise, 2024, 34, 1049-8923, 7265, 10.1002/rnc.7344 | |
7. | Jiayun Zheng, Feng Ding, A filtering‐based recursive extended least squares algorithm and its convergence for finite impulse response moving average systems, 2024, 34, 1049-8923, 6063, 10.1002/rnc.7307 | |
8. | Huanqi Sun, Weili Xiong, Feng Ding, Erfu Yang, Hierarchical estimation methods based on the penalty term for controlled autoregressive systems with colored noises, 2024, 34, 1049-8923, 6804, 10.1002/rnc.7323 | |
9. | Haoming Xing, Feng Ding, Xiao Zhang, Xiaoli Luan, Erfu Yang, Highly-efficient filtered hierarchical identification algorithms for multiple-input multiple-output systems with colored noises, 2024, 186, 01676911, 105762, 10.1016/j.sysconle.2024.105762 | |
10. | Yanjiao Wang, Yiting Liu, Jiehao Chen, Shihua Tang, Muqing Deng, Online identification of Hammerstein systems with B‐spline networks, 2024, 38, 0890-6327, 2074, 10.1002/acs.3792 | |
11. | Guiyao Ke, Creation of Three‐Scroll Hidden Conservative Lorenz‐Like Chaotic Flows, 2024, 7, 2513-0390, 10.1002/adts.202400247 | |
12. | Haijun Wang, Jun Pan, Guiyao Ke, Revealing More Hidden Attractors from a New Sub-Quadratic Lorenz-Like System of Degree 6 5, 2024, 34, 0218-1274, 10.1142/S0218127424500718 | |
13. | Ling Xu, Feng Ding, Xiao Zhang, Quanmin Zhu, Novel parameter estimation method for the systems with colored noises by using the filtering identification idea, 2024, 186, 01676911, 105774, 10.1016/j.sysconle.2024.105774 | |
14. | Saida Bedoui, Kamel Abderrahim, Feng Ding, Iterative parameter identification for Hammerstein systems with ARMA noises by using the filtering identification idea, 2024, 38, 0890-6327, 3134, 10.1002/acs.3865 | |
15. | Jun Pan, Haijun Wang, Feiyu Hu, Revealing asymmetric homoclinic and heteroclinic orbits, 2025, 33, 2688-1594, 1337, 10.3934/era.2025061 | |
16. | Haijun Wang, Jun Pan, Feiyu Hu, Guiyao Ke, Asymmetric Singularly Degenerate Heteroclinic Cycles, 2025, 35, 0218-1274, 10.1142/S0218127425500725 | |
17. | Jun Pan, Haijun Wang, Guiyao Ke, Feiyu Hu, Creation of single-wing Lorenz-like attractors via a ten-ninths-degree term, 2025, 23, 2391-5471, 10.1515/phys-2025-0165 |
No. | Parameters | Description | Value | Reff. |
1. | μ | Death rate and birth rate | 152∗70 | [27] |
2. | β | Infection rate | Estimated | - |
3. | r | Probability of succesed vaccination | 90% | [28] |
4. | p | Proportion of vaccinated individuals | Varied | - |
5. | δ | Hospitalized rate | Estimated | - |
6. | γ1 | Recovery rate of infected | 1/3 | [10] |
7. | γ2 | Recovery rate of hospitalized | 1 | - |
i | Change, (△y)i | Probability, Pi |
1 | [1,0,0]T | μN△t |
2 | [−1,1,0]T | βN(1−rp)ˉs(t)ˉi(t)△t |
3 | [−1,0,0]T | rpˉs(t)△t |
4 | [−1,0,0]T | μˉs(t)△t |
5 | [0,−1,1]T | δˉi(t)△t |
6 | [0,−1,0]T | γ1ˉi(t)△t |
7 | [0,−1,0]T | μˉi(t)△t |
8 | [0,0,−1]T | γ2ˉh(t)△t |
9 | [0,0,−1]T | μˉh(t)△t |
No. | Parameters | Description | Value | Reff. |
1. | μ | Death rate and birth rate | 152∗70 | [27] |
2. | β | Infection rate | Estimated | - |
3. | r | Probability of succesed vaccination | 90% | [28] |
4. | p | Proportion of vaccinated individuals | Varied | - |
5. | δ | Hospitalized rate | Estimated | - |
6. | γ1 | Recovery rate of infected | 1/3 | [10] |
7. | γ2 | Recovery rate of hospitalized | 1 | - |
i | Change, (△y)i | Probability, Pi |
1 | [1,0,0]T | μN△t |
2 | [−1,1,0]T | βN(1−rp)ˉs(t)ˉi(t)△t |
3 | [−1,0,0]T | rpˉs(t)△t |
4 | [−1,0,0]T | μˉs(t)△t |
5 | [0,−1,1]T | δˉi(t)△t |
6 | [0,−1,0]T | γ1ˉi(t)△t |
7 | [0,−1,0]T | μˉi(t)△t |
8 | [0,0,−1]T | γ2ˉh(t)△t |
9 | [0,0,−1]T | μˉh(t)△t |