key: cord-0893514-potnj7ls authors: Momani, Shaher; Ibrahim, Rabha W.; Hadid, Samir B. title: Susceptible-Infected-Susceptible Epidemic Discrete Dynamic System Based on Tsallis Entropy date: 2020-07-14 journal: Entropy (Basel) DOI: 10.3390/e22070769 sha: 6d86e6949b78d2aa1df0951dc72818a1074da084 doc_id: 893514 cord_uid: potnj7ls This investigation deals with a discrete dynamic system of susceptible-infected-susceptible epidemic (SISE) using the Tsallis entropy. We investigate the positive and maximal solutions of the system. Stability and equilibrium are studied. Moreover, based on the Tsallis entropy, we shall formulate a new design for the basic reproductive ratio. Finally, we apply the results on live data regarding COVID-19. Discrete dynamic systems of SISE were extensively discussed for a long historical period, that successfully described the procedure in disease diffusion (see [1] ). A decay ago, in the year 1927, the traditional SISE was offered [2] . After that, there established an enormous number of periodicals on SISE [3] [4] [5] . In overall, SISEs are considered to be homogeneously combined, which indicates that susceptible persons are infected with the same information. Nevertheless, there are various systems of populations in human culture [6] , and the joining between persons is not identical. The stability and convergence of the systems are studied by using the basic reproductive ratio. This ratio is given in different formula based on the system and the situation of the solution. In our discussion, we suggest new formal of this ratio based on the entropy concept. For COVID-19, the researchers established a suitable ratio called the case fatality rate (CFR). In [7] , the researchers studied SISE at level-liberated networks; it designates that under the suitable parameters, there is probably a threshold at which the disease will persevere. In view of [7] , some inoculation approaches are investigated, which additional develop the mechanisms of SISE on networks [8] [9] [10] . Obviously, the classes of difference equations have numerous practices in SISE system [11] [12] [13] . In reality, a positive interval, and discrete simulations often give information about disease [14] . On the other hand, a difference equation is the discretion of the continuous model [15, 16] , which indicates it practical to respond to the approximation method. Particularly, discrete simulations show a more complex dynamical conduct than the conforming continuous representations [17] [18] [19] [20] . Under these compensations in attention, the state of the discrete SISE system of networks is about excessive investigation care. From the above-mentioned details, we shall deal with a discrete-time SISE system involving Tsallis entropy, which will be important work. We apply the results to live data regarding COVID-19. In traditional statistical methods, the entropy function formerly presented by Rudolf Clausius is construed as statistical entropy utilizing probability theory. The statistical entropy view was introduced in 19th century with the work of physicist Ludwig Boltzmann. This entropy was generalized by Tsallis as follows [21] : Consider a discrete set of probabilities {ρ  } satisfying the condition ∑  ρ  = 1, and α any real number, the Tsallis entropy is formulated by the terms where α is a real parameter which is known as the entropy-index. One of the most important property of Tsallis entropy is that it has a maximum value determining when each micro-state is equiprobable (ρ  = 1/Φ for all ) and then we get [22] ). Machado [23, 24] presented novel formulas for entropy inspired by using the behavior of fractional calculus. The results of the generalized fractional entropy are examined both in usual probability distributions and data series. Moreover, by using the quantum deformed calculus, Hasan et al. [25] introduced a generalized q-entropy. Numerous issues rule the transmissibility of the infection from the affected to the unaffected. In addition, disease dynamical systems can be investigated at altered rules: the single distinct, small collections of people, and among whole people. Different representations are selected given by the complexity of available data. In their contemporary avatar, computers that generate the numbers and distribution designs of infections simulate systems (see [26] [27] [28] [29] ). The SISE system is formulated with N patrons and all the patrons are separated into n groups by their joints (junctions) ( = 1, 2, ..., n). Consequently, it has N = ∑ n =1 N  , where N  represents the total number of the patron with position . It is considered that every patron has two positions, the first position is infected (I) and the second position is the susceptible (S). The susceptible patron may be infected with transmission ratio τ, and the infected patron may improve to a susceptible patron with repossession ratio . Hence, we obtain the equation N  (t) = S  (t) + I  (t) and the discrete system is the Tsallis entropy introduced by the probability that any given connect points to an infected node and T indicates the time-step measure. It is a value indicating out that scheme (1) is recognized by employing the forward Euler pattern to the continuous SISE system and the equilibrium points (in discrete system they are equal to the fixed points) of structure (1) are similar as for the continuous equivalent. By letting Approximate (2) to entropy system, we have where The definition of the function Λ is more general description of the interaction of susceptible and infected individual using entropy. Entropy is a powerful implement for analysis telling the probability distributions of the potential formal of a system, and hence the information encoded in it. Nevertheless, significant information may also be organized in the time-based dynamics, a feature that is not typically taken into account. The notion of scheming entropy based on non-linear designs is utilized to discover spatial structures and processes. Usually, spatial processes have been supposed to be linear, characteristically in relations of a linear auto regressive or heartrending average process. Nevertheless, additional spatial dynamics are possible to display nonlinear types in a technique that is similar to time-based systems. The capacities of nonlinear systems are progressively documented in science, as the restrictions of equilibrium representations in clarifying real-world phenomena convert more seeming. As a consequence, interest is growing in the growth studies in science. We proceed to conclude the existence of solution of (3). Consider the entropy discrete system of SISE (3) . Then it has bounded non-negative solutions if the following hypotheses are achieved Proof. By the maximum value of the Tsallis entropy, System (3) implies that By letting α → ∞, we have lim α→∞ τT (α−1) 2 = 0 for all fixed parameters τ and T . Thus, I(1) is bounded by N. Moreover, since I  (0) = N  with T 1 then this yields that for α → ∞, the initial solution becomes I(0) = N and consequently the step one of solution becomes I(1) = N(1 − T ) ≥ 0, which leads to the non-negative solution I. Hence, by induction, one can prove that 0 ≤ I(t) ≤ N for In this section, we aim to study the stability of SISE (1). By substituting I  (t) = N  (t) − S  (t) in the first equation of System (1), we have which is equivalent to the following system The disease free equilibrium of SISE (6) can be computed by the following construction Ξ 0 (S 1 (0), ..., S n (0), I 1 (0), ..., I n (0)) = (N 1 , ..., N n , 0, ..., 0). By employing the linearization matrix method [19] on the system (6) at the point Ξ 0 , we obtain where ψ (the vector of new infections) and σ ( the vector of all other transitions including disease-connected deaths) are non-negative such that ψ + σ is irreducible and Ψ indicates the Jacobi matrix at Ξ 0 (we assume that this point is unique). Note that Hence, System (6) approximates to the form • We used the maximum value of entropy in our system because our suggested system is formulated only for the infected and susceptible persons. We did not include the removed cases R(t) (death and recovery). This variable may be defined by using the maximum entropy R(t) = Υ max α (I(t)) (see [30] ). Note that entropy index is strongly connected to the number of individuals N and the number of groups n, (1 ≤ n ≤ N), so that when n = N, one would expect the SIS model consequence with non-linear incidence. Very recently, Tsallis and Tirnakli [31] proposed a q-statistical functional arrangement that acts to describe acceptably the existing information for all countries. Reliably, calculations of the dates and altitudes of those peaks in rigorously affected countries become likely unless well-organized actions or vaccines, or functional modifications of the accepted epidemiological approaches, arise. The basic reproductive ratio (λ 0 ) can be explained as the predictable ratio of cases openly produced by one case in a resident where all persons are subject to infection. Mathematically, it is known as the spectral radius of the matrix ψ(I d − σ) −1 (the largest absolute number of the eigenvalues). There are other different definitions and formulas can describe the situation properly. This ratio plays an important role to achieve the stability. It has been shown in many studies if λ 0 > 1 then we indicate an unstable situation and if λ 0 < 1 then the situation is asymptotically stable, while the case λ 0 = 1 indicates the stability, but not being asymptotic [19] . Recently, for COVID-19, researchers suggested the case fatality rate (CFR, the aim is to reduce this ratio) [32] where D indicates the number of dying people. For example, if the number D = 10 and I = 500, then the ratio is 2%. Simultaneously, if it is recorded that there are 500 susceptible persons then The aim is to reduce the rate λ 0 or CFR by isolated position and keep cleaning the environment of the person. In our discussion, we suggest to involve the entropy evaluation for this rate. In [33] the authors formulated λ 0 by using the probability of the survival function P as follows (for discrete data): From the above example I = 500, D = 10 and N = 1000, we have The idea of the probability of the survival function is not suitable for COVID-19. Therefore, based on our SISE system, we suggest to use the maximum entropy Υ max α (I) as follows: , t = 0, 1, 2, ... α = 1. Based on the above data, we indicate the following results (λ 0 ) 2 = 0.5 ≈ λ 0 , that is lim α→2 (λ 0 ) α = λ 0 . Note that , (λ 0 ) 3 = 0.25, (λ 0 ) 4 = 0.16, .... We conclude that the ratio decreases whenever α increases. Hence, the SISE system is stable, while for α < 1, the system is unstable. For example, when α = 0.5 and the probability P(I) = 0.1 this implies that Φ = 1 P(I) = 10, we get Υ max α (I(t) = 1−10 0.5 −0.5 = 4.324; which leads (λ 0 ) 0.5 = 4.324 × 0.5 = 2.162 > 1. From above, we conclude that Theorem 1 can be extended to include the stability as follows: Consider the entropy discreet system of SISE (3). If the conditions hold. Then every solution of (3) is bounded non-negative and stable satisfying the basic reproductive ratio , t = 0, 1, 2, ... . The survival function (or it called reliability function) is a function that offers the probability that a patient, scheme, or other thing of concern will survive further than any indicated time and it is one of the techniques to define and show survival data. It states as the probability that a subject survives longer than time t. The distribution of survival times may be approximated well by a function such as the exponential distribution. Numerous distributions are usually utilized in survival analysis, containing the exponential, Tsallis entropy, gamma, normal and log-logistic. These distributions are formulated by parameters. The entropy optimization principle (includes the maximum entropy) converts it from a measure of information into an implement of statistics conclusively [34] . Since the higher maximum entropy goes to Tsallis entropy (see [35] ), then it is a confidence to employ this fact to define CFR. In this section, we use live data to examine our theoretical results, especially the stability of the SISE system by using (λ) α . Table 1 shows the data from the first infected countries until the end of May. The rate of death is given by using CFR N (t). The basic reproduction ratio is evaluated by using (λ 0 ) α for some α > 1. We use the conditions of Theorem 2 to get non-negative bounded and stable solution. One can consider the following system where a := τT , b := T . Let a = b = 0.4, initial condition (S 0 , I 0 ) = (0.1, 0.1) and α = 2, we have a stable limit cycle for the system of period one (see Figure 1 ). The red line shows the values of each considered case. Example 2. Consider the following system where a := τT , b = T and α = 2 with (0.2, 0.6) as initial point. In all figures, the red line indicates the values of each considered case as follows: • For a = 0.9, b = 0.3, the system has a limit cycle with period 4, while for a = 0.5, b = 0.3 the system has a limit cycle with period 2. Note that the simplest situation is the case where there is no recuperation rate. This leads to an SI-like model, so that the pathogen infects all individuals on the long run. The simple continuous SI model has the logistic function as a solution and its discretized version is the logistic map, which presents the traditional bifurcation diagram as stable solutions. The correct balance between short-and long-term data loading in the world of big data has the following strategies: • An applied perception based on the main usages for the data; • The data constructions wanted for analysis; • The relevancy of the data over time; • Development with little organization is essential. We prepared a new formula of the basic reproductive ratio (λ 0 ) α , which is defined by the Tsallis entropy. The formula is useful for the stability of SISE system involving the Tsallis entropy. It is related to long time data ( by taking in account the above stratifies, it may modify for short time data ). We applied the suggested system by using live data regarding COVID-19. A Short History of Mathematical Population Dynamics A contribution to the mathematical theory of epidemics A new way of investigating the asymptotic behaviour of a stochastic sis system with multiplicative noise Dynamics of an SIS reaction-diffusion epidemic model for disease transmission Dynamics of the impact of twitter with time delay on the spread of infectious diseases The structure and function of complex networks Epidemic spreading in scale-free networks Immunization of complex networks Behaviors of susceptible-infected epidemics on scale-free networks with identical infectivity Immunization and epidemic threshold of an SIS model in complex networks Some discrete-time SI, SIR, and sis epidemic models Effect of discretization on dynamical behavior of SIR and sir models with nonlinear incidence Stability and bifurcation analysis in a discrete sir epidemic model An Introduction to Difference Equations Discretizations of nonlinear differential equations using explicit nonstandard methods Difference equations from discretization of a continuous epidemic model with immigration of infectives Bifurcation and chaos in an epidemic model with nonlinear incidence rates Discrete-time SIS models with simple and complex population dynamics The basic reproduction number in some discrete-time epidemic models Dynamical analysis of a discrete-time SIS epidemic model on complex networks Possible generalization of Boltzmann-Gibbs statistics Domínguez-Jiménez, I. Determining the entropic index q of Tsallis entropy in images through redundancy Fractional order generalized information Fractional Renyi entropy Classification of Covid-19 Coronavirus, Pneumonia and Healthy Lungs in CT Scans Using Q-Deformed Entropy and Deep Learning Features The fractional differential polynomial neural network for approximation of functions Utility function for intelligent access web selection using the normalized fuzzy fractional entropy New Texture Descriptor Based on Modified Fractional Entropy for Digital Image Splicing Forgery Detection Analytic study of complex fractional Tsallis' entropy with applications in CNNs Maximum entropy method for estimating the reproduction number: An investigation for COVID-19 in China Predicting COVID-19 peaks around the world. Front. Phys. 2020, 8, 217 COVID19 in Numbers-R 0 , the Case Fatality Rate and Why We Need to Flatten the curve.webm Perspectives on the basic reproduction ratio A new entropy optimization model for graduation of data in survival analysis Tsallis entropy theory for modeling in water engineering: A review The authors would like express their sincere appreciation to the reviewers for their very careful review of our paper and rich it in information. The authors declare no conflict of interest. The work here is supported by the University Ajman grant: 2020-COVID-19-08.