The influence of an appropriate reporting time and publicity intensity on the spread of infectious diseases

: We present a stochastic time-delay susceptible-exposed-asymptomatic-symptom-vaccinated-recovered (SEAQVR) model with media publicity e ff ect in this study. The model takes into account the impacts of noise, time delay and public sensitivity on infectious illness propagation. The stochastic dynamics of the system are analyzed at the Hopf bifurcation, using time delay and noise intensity as bifurcation parameters, and the theoretical conclusions are validated using numerical simulation. Increasing the time delay and sensitivity coe ffi cient can e ff ectively delay the occurrence of the peak number of infected individuals and mitigate the extent of infection. Additionally, time delay and noise intensity are shown to have specific thresholds, beyond which periodic infections occur. Notably, heightened public sensitivity reduces the threshold for time delay, and media publicity directly a ff ects public sensitivity. The numerical simulation reveals that increasing media publicity intensity does not always yield better results, and that the sensitivity of the public at present is an important reference index for setting an appropriate publicity intensity.


Introduction
In recent years, the dynamics of infectious diseases have changed fundamentally as a result of population increase, faster urbanization, climate change and economic globalization. Particularly, people have become more prone to infectious illness epidemics, and particular microbes have gained medication resistance. We have also discovered substantial changes in vector ecology, such as increased Aedes activity, which leads to increased disease dissemination. Furthermore, because of modern society's high degree of connectedness and urbanization, diseases are more prone to spread among people [1]. Dengue fever broke out globally in 2016, resulting in 100 million illnesses and 38,000 deaths [2]. When an infectious disease spreads, it not only harms people's bodily and mental health, but it also necessitates the use of human, material and financial resources to manage it. In some developing countries, serious infectious diseases will generate civil discontent and impede social development. As a result, analyzing the transmission trend of infectious diseases and controlling measures is critical.
In 1972, Kermack and McKendrick constructed the SIR cabin model [3] by using the dynamic system modeling method, which divided the population into three types and studied the disease's spread law and epidemic mechanism in detail, laying the theoretical foundation for the dynamic model of infectious diseases. In recent years, Korolev developed an SEIRD model which outlines how to use supplementary information from random tests to calibrate the model's initial parameters and restrict the range of probable forecasts for future deaths [4]. The SEAIR model was developed by Basnarkov, who discovered that the centrality of feature vectors roughly determines the chance of infection [5]. He et al. proposed a fractional discrete-time SIR model with vaccination, which proved and quantified the complex dynamics of the system [6]. Meng et al. developed the SEIRV model and an evolutionary game model to investigate the differences between forced and voluntary inoculation methods in heterogeneous networks [7]. Goel developed an SIR model with a Holling type-II treatment rate and used the Lyapunov approach to investigate the stability of the equilibrium point [8].
The media is crucial in the prevention and control of infectious diseases [9]. When an epidemic spreads in a country, the government employs the media to educate citizens on how to appropriately respond. The execution of certain policies and the dissemination of information will have an impact on human behavior [10,11]. Public preventative measures that are implemented on time and effectively can significantly reduce the infection rate [12]. The intensities of publicity, cultural level and social duty may influence people's sensitivity to information published in the media [13,14]. Individualism, collectivism and ethnic diversity will each have an impact on how people respond to health emergencies [15].
Examining the impact of time delay on system stability constitutes a crucial aspect of investigating system dynamics. Studying the epidemic model with time delay can better depict the transmission mechanism of illnesses because most biological processes involve time delay. Misra et al. analyzed the stability and direction of the Hopf bifurcation, as well as the time delay of carrying out an awareness plan [16]. According to Gutierrez et al., a delay in reporting the death toll resulted in a significant increase in illness severity [17]. The public's behavior and psychological state will be impacted by the media's delayed reporting of the pandemic scenario, which is also influenced by technology, capability, resources and other variables.
Random factors strongly impact the propagation of infectious illnesses in nature. Factors such as immunity, temperature and humidity, for example, will affect the infection rate. People will take different countermeasures based on the heterogeneity of information, which will interfere with the infection rate. The deterministic equation has been idealized. In comparison to the deterministic equation, the stochastic equation can more accurately represent the actual scenario [18][19][20][21][22][23]. According to the Din et al. stochastic model, noise assured the extinction of the hepatitis B virus [24]. Krause et al. extended the random SIS epidemic model to the spatial network, thus obtaining the random epidemic ensemble population model; they discovered that Gaussian white noise can be used to offset the cure rate [25]. These random factors can be approximately simulated by Gaussian white noise and introduced into the deterministic epidemic model to study the influence of noise on the dynamic behavior of the system [26]. Based on the model presented in [27], this paper examines the influence of reporting time delay and Gaussian white noise on infection rates. It introduces the sensitivity coefficient, denoted as "k," which represents the public's responsiveness to the number of newly infected individuals displaying symptoms. A stochastic differential equation model incorporating time delay is developed. Excessive promotion or non-dissemination of information about infectious diseases is generally not recommended. By conducting dynamic analyses of time delay, noise intensity and the sensitivity coefficient, this research provides valuable insights for media organizations to devise appropriate publicity strategies. Numerical simulations were performed to investigate the impact of different levels of public sensitivity on the selection of effective publicity intensity.
This paper is structured as follows: In the second section, we provide a detailed explanation of the SEAVQR model, the calculation of the basic reproductive number and the necessary conditions for the existence of the endemic equilibrium point. We further investigate the Hopf bifurcation conditions and present the stability analysis. The third section focuses on deriving the stochastic Itô equation and examining the occurrence of stochastic bifurcation. In the fourth section, we discuss the numerical simulation to analyze the respective influences of reporting time delay, the sensitivity coefficient and noise intensity on the spread of infectious diseases. Finally, in the fifth part, we summarize the findings from the preceding sections.

The analysis of the model
Building upon the model proposed in [27], we consider the effects of noise, reporting time delay and the sensitivity coefficient on the infection rate β, and introduce the SEAQVR model. Figure 1 illustrates the schematic diagram of the model, while the model itself is described as follows: .
(2.1) Population has been divided into six categories: susceptible (S ), exposed (E), asymptomatic infected (A), symptomatic infected (Q), vaccinated (V) and recovered (R). Ω denotes the supplementary population, µ represents the natural mortality rate, δ signifies the mortality rate, φ denotes the rate of transition from exposed individuals to infected individuals and γ refers to the recovery rate of infected individuals, both symptomatic and asymptomatic. η denotes the vaccination rate, while Γ represents the rate of vaccine failure. α represents the proportion of symptomatic individuals among the infected population, whereas β and β g denote the infection rates post-vaccination and without vaccination, respectively. We posit that β is influenced by the collective awareness of prevention among the general public. As the number of infections escalates, individuals often proactively adopt measures to safeguard themselves and impede the transmission chain, while also bolstering their own immunity. k is the coefficient of sensitivity of the public to the number of new symptomatic infections. Naturally, the level of protection corresponds to the sensitivity of the public. The greater the sensitivity, the more active and vigilant are the protective measures. However, in real-world situations, the public's sensitivity to the number of newly infected individuals varies. For instance, if the public lacks accurate information, society becomes permeated with feelings of fear and confusion, consequently disrupting the public's sensitivity. With a range of − π 2 , π 2 , arc tan x is a monotonically increasing odd function with a flat change and global boundedness. The model becomes more logical because it can adjust the infection rate β within (0, 1). ξ (t) is Gaussian white noise with a power spectral density, and ε is a sufficiently small-scale parameter.
The basic reproductive number R 0 represents the average number of infections in a susceptible population for an infected person in a disturbance-free system. There is always a disease-free equilibrium point P 1 (S 1 , The basic reproductive number of the undisturbed model, as calculated via the next-generation matrix method [28], is Next, the dynamic behavior near the endemic equilibrium point is considered. Theorem 2.1. If R 0 > 1, system (2.1) has an endemic equilibrium. If R 0 < 1, system (2.1) has no endemic equilibrium.
Proof. Assuming that the right-hand term of the first six equations in (2.1) are equal to 0, we can determine the endemic equilibrium point P * (S * , E * , A * , Q * , V * , R * ), where A * satisfies the following equation: Clearly, considering the condition m 1 < 0, let us assume that A 1 and A 2 represent two roots of (2.2). If the condition R 0 > 1 is met, we have that m 3 > 0. Then, A 1 A 2 < 0 is obtained, indicating the existence of a positive root among the two roots, and, consequently, an endemic equilibrium point. □ In the case of R 0 < 1, m 3 < 0. So, we have that A 1 A 2 > 0, indicating that both roots can be either positive or negative. In this scenario, the following formula holds.
Consequently, A 1 + A 2 < 0, indicating that both roots are negative and there is no existence of an endemic equilibrium point.
Considering that the first five equations of (2.1) do not account for the recovered individuals (R), the following only analyzes the first five equations and linearizes the system at the endemic equilibrium point. Let The equilibrium point of (2.1) is shifted to the origin, and the equation is expressed in vector form as follows: Without Gaussian white noise, (2.1) can be linearized as follows.
The characteristic equation corresponding to (2.3) is and f 2 (0) = 0 , whether or not 0 is the root of the characteristic equation relies on whether f 1 (0) = 0. Since f 1 (0) > 0, it follows that 0 cannot be the root of the characteristic equation given by (2.4).
In order to determine if endemic diseases are stable at their equilibrium point when τ = 0, the characteristic equation given by (2.4) is then transformed into the following.
All roots of the characteristic equation (2.4) have negative real parts, and the equilibrium point of endemic diseases is locally asymptotically stable.
Theorem 2.2. If R 0 > 1 and H1 and H2: 0 are satisfied, the system will generate Hopf bifurcation at the endemic equilibrium point.
Proof. Differentiating the two sides of (2.5) with respect to τ, the Hopf transversality condition is as follows: dλ dτ Re dλ dτ and Re dλ dτ −1 have the same signs. Therefore, if condition (H2) holds, according to Hopf bifurcation theory, Hopf bifurcation occurs at τ = τ 0 , and the theorem is proved. □
where −τ 0 ⩽ θ ⩽ 0. We use the extended equation of z t (θ) to represent (2.1) and calculate the center manifold as follows: where L [z t (θ)] and F [z t (θ)] are the linear and nonlinear parts of (2.1), respectively. After combining (3.1) with (3.2), we get Considering <Ψ (s) , h (θ, y (t)) >= 0, here are the calculated stochastic ordinary differential equations:ẏ N , Carry out the polar coordinate transformation by using the stochastic averaging method [29]: where R (t) and φ (t) are the amplitude and phase of the solution, respectively. We can obtain the stochastic ODEs with R (t) and φ (t): Since R (t) and φ (t) are both slow varying processes, the random averaging method is used to average time on the pseudo-period 2π ω 0 . The amplitude process after smoothing R (t) is a Markov diffusion process, and the Itô equation is obtained as shown below.
R = 0 is the equilibrium point of the system, whose Lyapunov exponent is When λ < 0, that is, µ 1 < 1 2 µ 3 , the trivial solution to the equation is locally asymptotically stable. When λ > 0, that is, µ 1 > 1 2 µ 3 , the trivial solution to the equation is unstable. Next, the global dynamic properties are obtained via boundary classification and the three-exponential method. When µ 2 < 0, the right boundary R → +∞ is an entry boundary, and this condition is the premise of the following discussion. Because m (0) = σ (0) = 0, the left boundary belongs to the first kind of singular boundary. Its diffusion coefficient α = 2, drift coefficient β = 1 and characteristic value c = 2µ 1 µ 3 . When 2µ 1 µ 3 > 1, the nontrivial stationary probability density exists in the system.
Because α − β = 1, the stationary probability density can be simplified as p (R) = O (R c−α ). Dynamic bifurcation, also known as D-bifurcation, refers to a change in integrability. Phenomenological bifurcation, P-bifurcation, refers to a change of probability density image shape. Through analysis, D-bifurcation occurs at µ 1 = 1 2 µ 3 , and P-bifurcation occurs at µ 1 = µ 3 , as shown in Figure 2.

Numerical simulation
In this section, in order to reveal the dynamic behavior of the system, we utilize a sample size of 100 and employ the Euler-Maruyama method to carry out the simulations. Based on the theoretical study in the preceding sections, time delay plays a vital role in studying the spread of infectious diseases as a bifurcation parameter. The stability of the system changes as the time delay rises, and the system has a periodic solution.
The basic reproductive number of the system in an undisturbed state R 0 = 2.02741 > 1 is determined using the numbers in Table 1, so the equilibrium point of endemic diseases exists in system (2.1) and is recorded as P * (S * , E * , A * , V * , Q * , R * ), in which S * = 552030.54, E * = 1405470.79, A * = 87494.34, V * = 4741153.57, Q * = 71666.036, R * = 3142180. L (0) = −10 −12 and m = 25.1664, so both H1 and H2 are satisfied. Through the analysis in the second section, it is calculated that ω 0 = 0.302386, τ 0 = 12.5623 and Hopf bifurcation occurs at τ = τ 0 . The noise intensity was set to K = 0 and we drew the time series diagrams of S , E, A, V, Q, R with τ = 10 and τ = 15, as illustrated in Figures 3 and 4. When τ = 10 < τ 0 , the endemic equilibrium point is stable, and when τ = 15 > τ 0 , an infection periodically develops at the endemic equilibrium point. It should be mentioned that there is a time interval when the media reports the progress of infectious diseases. The system may oscillate if the delay of reporting the number of newly infected individuals with symptoms exceeds the threshold.  Suppose that k = 0.1; then, τ 0 = 1.75 can be determined; we applied τ = 10, 15, 20 in system (2.1) correspondingly. Figure 5 depicts the oscillation of the Q/N time-series curve. Observing Figure 5, the number of symptomatic infections increases as the reporting time delay grows. By maintaining a reporting time delay below the threshold, it ensures that individuals can comprehend the epidemic situation and implement timely prevention and control measures. However, when the reporting time is delayed beyond the threshold, the public lacks access to relevant information, hampering their ability to make scientifically informed decisions and implement effective prevention measures. As a consequence, the number of infected individuals increases. As a result, the threshold τ 0 is extremely useful for the study, control and eradication of infectious illnesses. We have chosen various sensitivity coefficients, k, and computed τ 0 accordingly. Based on these calculations, we have generated a line chart, as shown in Figure 6, illustrating the relationship between τ 0 and the sensitivity coefficient, k. It is discovered that τ 0 decreases as the sensitivity coefficient k increases. As the sensitivity coefficient increases, there is a corresponding rise in the level of attention people pay to the number of new symptomatic infections. Consequently, this increased sensitivity generates a higher demand for real-time reporting.  Postponing the occurrence of the peak period can be advantageous for the government, as it enables the allocation of resources, thereby bolstering medical capabilities and providing additional time to realize the control of infectious diseases. By reducing the peak, the burden on medical resources can be mitigated and infected individuals can be brought to a manageable stage, ensuring adherence to medical standards during an outbreak. Hence, it is crucial for the media to exercise reasonable control over the timing of reporting and publicity efforts. When the sensitivity coefficient, denoted as k, remains constant, a reporting delay close to τ 0 can yield significant results with minimal effort. Additionally, it should be noted that the value of the sensitivity coefficient, k, can be subject to variation.
The random dynamic behavior close to τ 0 will now be then examined. Select noise intensity as the bifurcation parameter, set τ = 12.562 < τ 0 , ε = 0.00000001 and refer to Table 1 for additional parameters. The third section's theoretical study reveals that the system will randomly bifurcate when µ 1 = µ 3 , that is, when K=1.613. For Figure 9, we have selected a specific time period ranging from 3000d to 5000d to conduct simulations and generate the joint probability density diagram and marginal probability density diagram for susceptible persons (S ) and symptomatic infected persons (Q) with noise intensity increasing from 1 to 200. The marginal probability map transforms from a single peak to a multi-peak structure as the noise intensity rises, and a random P-bifurcation takes place. A time-series plot of Q/N between K = 1 and K = 200 is shown in Figure 10(a) and 10(b). When the noise intensity is minimal, the system usually remains stable. When the noise intensity is high, the system oscillates, which implies that, as the noise intensity rises, so does the population of infected people. Let τ = 10, ε = 0.01 and the noise intensity K traverse 1, 10, 15 and 20 in Figure 11. The variance of the random process considerably rises with the noise level over a period ranging from 4000d to 5000d. As can be seen from Figure 12(a) and 12(b), the noise intensity increases and the system loses stability over a period ranging from 2000d to 5000d. Social networking platforms like Twitter, Weibo and Facebook are often used, which has facilitated the spread of rumors that might potentially frighten the population [30,31]. Furthermore, the dissemination of false information will impede the reporting of information and may encourage people to select inefficient preventative and control measures that cannot successfully control infectious diseases [32]. Therefore, media should convey the information to the public in a stable, objective and scientific manner.   We normally believe that, when an epidemic strikes, the greater the media coverage, the better. However, as shown in Figure 6, τ 0 drops as the sensitivity coefficient increases, and the media's reporting time delay is restricted by human resources, technical equipment and so on. In some circumstances, the reporting delay will exceed the threshold τ 0 . As the amplitude of the periodic oscillation increases, control methods must be changed over time, posing significant obstacles to epidemic prevention and control. Furthermore, before the health department can supply the information, vast volumes of noise information in the data must be collected and pre-processed. A minor reporting time delay may reduce the accuracy of the reporting information and influence the system's stability. Governments should review whether publicity intensity matches correlative sectors' abilities to respond to emergencies and gather and analyze data. It will be self-defeating if the correlative sectors' capacity does not fulfill the threshold criteria. It is assumed that the correlative sectors' capacity can meet the τ 0 requirement, and we choose τ 0 as the reporting time.
With increased publicity intensity, the peak might occur considerably sooner, as indicated in Figure 13(a). In Figure 13(a), an increase in the value of k results in the peak arrival time being shifted forward by 20 days. Increasing the threshold by reducing the intensity of publicity can reduce work pressure in the health industry, but the peak size may increase, as illustrated in Figure 13(b). According to Figure 13(b), reducing the value of k leads to an increase in the peak number of symptomatic infections by 600 individuals. It should be emphasized that cultural factors influence the sensitivity coefficient k. It cannot be changed randomly with a change in publicity intensity, but it can change within an interval. In a group with a sense of duty, k is frequently greater. Figure 13(c) depicts a time-series plot mimicking k at various publicity intensities ranging from 0.25 to 0.5. We discovered that, assuming that the health department's capability can match the publicity intensity, the larger the publicity intensity, the smaller the peak, and the later the peak occurs. Of course, the sensitivity coefficient's growth has an upper bound. Following the previous research, we suspect that there is a k 0 . When k > k 0 , we should raise the degree of publicity. When k < k 0 , the arrival time and peak value cannot be simultaneously optimized, and trade-offs must be made based on the actual scenario. The magnitude and timing of spikes have a substantial societal and economic impact, with overburdened healthcare systems increasing infectious mortality, direct morbidity and medical consequences [33]. As a result, the government should make appropriate publicity intensity based on current circumstances, such as medical resources, social stability and public welfare. The current sensitivity of society which can roughly reflect the change interval of the sensitivity coefficient under the action of media publicity, is also an important reference index for the government.

Conclusions
In this paper, a SEAVQR model with Gaussian white noise is established. Due to the process of data collection and preprocessing, there exists a temporal gap in reporting infectious diseases. We discovered that an acceptable reporting time delay aids in the control of infectious illness spread. Hopf bifurcation occurs at τ = τ 0 . Setting τ 0 as the reporting time-delay threshold, and under the condition that the reporting time delay is less than the threshold, the public can gain valuable insight into the epidemic situation and take prompt preventative and control actions. When the reporting time delay exceeds the threshold, a lack of public information on responses can lead to an increase in the number of infected people. If the reporting time delay exceeds the threshold, the larger the time delay, the greater the obstacle to the control of infectious diseases. The system takes into account the effect of noise on the transmission of infectious diseases. The dissemination of inaccurate information can hinder the public's ability to make informed judgments, and the proliferation of rumors can exacerbate the "fear effect," leading to an unstable system. When the reporting time delay approaches the threshold, the system may oscillate periodically if the noise intensity becomes larger. Therefore, the accuracy of reporting information and the timeliness of reporting should be considered by the government in the early control of the epidemic.
Furthermore, we discovered that increasing the reporting delay could postpone the arrival of the peak of symptomatic infected patients and minimize the peak of infection, giving us more precious time to control the epidemic. The same result can be obtained by increasing the sensitivity coefficient. In a stable state, adopting a longer reporting delay and the government enhancing public sensitivity through media awareness will help enhance medical treatment quality and living standards. As a result, while the sensitivity coefficient is fixed, it is better to set the reporting time delay around the threshold. However, we observe that increasing the sensitivity coefficient results in a drop in the stated delay threshold. The media's publicity intensity for the epidemic is not that the higher the intensity, the better. When the correlative departments' work abilities cannot meet the stronger publicity intensity, the public may be unable to adopt timely preventative and control measures due to low accuracy or big reporting delays, potentially leading to the spread of infectious illnesses becoming out of control. We also discovered the existence of k 0 . When k < k 0 , government increases in publicity intensity may cause the peak to arrive earlier. If the government decreases publicity intensity to raise the threshold at the expense of lowering the sensitivity coefficient, it will result in a peak increase and, to some extent, undermine public welfare. When k > k 0 , however, the increase in publicity intensity could be helpful in postponing and lowering the peak.
To sum up, during the pandemic, the government should correctly understand the public's attitude toward this infectious disease. If the sensitivity k exceeds k 0 , it should try its best to improve the sensitivity. If k is less than k 0 , it is necessary to formulate appropriate publicity intensity for media publicity under the guidance of social and economic conditions and the ability of correlative departments in order to maximize benefits. In addition, the appropriate reporting time should be taken near the threshold. The media should grasp timeliness, accuracy and effectiveness in reporting. Code is available at https://github.com/Wangqiubao/HC num sim.

Use of AI tools declaration
The authors declare that they have not used artificial intelligence tools in the creation of this article.