Analysis and Parametrical Estimation with Real COVID-19 Data of a New Extended SEIR Epidemic Model with Quarantined Individuals

During the recent COVID-19 pandemic, quarantine and testing policies have been of vital importance since the causative agent has been a novel virus and no vaccine was developed at the time. In this work, a new epidemiological deterministic model is proposed, analyzed, and discussed. Such a model includes quarantine periods of people with symptoms that have been tested positive, and it will trigger a trace of their close contact, who will be also tested and put in quarantine if the result is positive. Moreover, how the model parameters affect its stability is analyzed with the basic reproduction number R0. Since the COVID-19 outbreak in Spain (approximately 13/03/2020) until 25/04/2021, different restrictions have been applied. For discussion of a real case study, data have been gathered and used from the Spanish Autonomous Community of Cantabria to estimate the parameters and to see how the restrictions have affected their values. In the parameter estimation process, it has been assumed that the constructed model follows the structure of an ARX model. Finally, by considering that the gathered data are subject to certain errors, the paper discusses how to adequate the model usefulness for its use in fitting and processing data through an estimation mechanism involving the provided daily total and positive performed tests.


Introduction
In late December of 2019, China reported to the World Health Organization (WHO) several cases of pneumonia apparently linked to the Hunan Wholesale Seafood Market in Wuhan City, and they pointed out that this disease might be caused by a new virus. Chinese scientists analyzed the virus from a hospitalized person, and the 11th of January, China confirmed that the abnormal arising of pneumonia cases was due to a novel coronavirus, similar to SARS and MERS, so they named it SARS-CoV-2 (afterwards it will be renamed to . It was also found out that fever and sore cough were the common symptoms, and that pneumonia was a less common consequence of the disease causing virus [1]. An abundant background literature on many medical aspects and transmission mechanisms of COVID-19 is available, some of them being linked with general description of epidemic models. See, for instance,  and references there in and also [24][25][26][27][28] for some general results on epidemic model parameterizations of interest. Some specific related details are now described. By late January of 2020, some cities of China locked down so as to slow down the increase in hospitalized people, and some airlines decided to suspend flights to (or from) China. On the 30th of January, the WHO (World Health Organization) reported that there were 7,818 confirmed cases all over the world, from which 7,736 cases were in China, 1,370 were severe and 170 deaths [21]. Moreover, the Centers for Disease Control and Prevention (CDC) wrote a press release which outlined the first case of person-to-person spread in U.S, and consequently world countries geared up [8]. e applied prevention measures were not enough to avoid the virus spreading and on March 11 determined that COVID-19 was a global pandemic [19] since there were 132,758 confirmed cases in the world and 4,955 deaths [18]. On March 17, Europe and Russia closed the borders, and by the time, many countries announced lock downs (for instance, Spain, France, and so on) [1].
When Spain imposed a quarantine period on the 13th of March of 2020, its epidemiological situation was critical: there were 4,209 confirmed cases from which 575 were new ones, and 42% of the total cases were hospitalized [5]. Even if the rules were very restrictive (only on necessary situations it was allowed going outside and wearing the face mask was mandatory), the cases continued raising until 25/03/2020 with a pick value of 720 new cases and 825 deaths [6], and it was not until 25/05/2020, when the number of new cases was approximately 2,000 and 68 the number of deaths [6], that more permission was given to the inhabitants (i.e., it was possible working out outside at certain time ranges). On the 21th of June 2020, the Spanish Government announced the so-called "new normality," the mask face and social distancing (1.5 m) were still mandatory although at this point each autonomous community could impose more constrains but their capability was limited [7]. During summer vacation, the situation maintained stable; however, when people started working (approximately after 09/09/2020), the epidemiological situation worsen. To tackle this problem, each Autonomous Community (AC) combined and applied differently restrictions (as, for instance, curfew or bars close down), and considering the later results, inhabitants cast doubt on whether the imposed rules were worth or not.
erefore, a mathematical model which could represent the COVID-19 behavior over certain conditions could be very useful.
e aim of epidemiological models is to understand and predict the dynamics of the propagation of a disease, and it must exhibit an equilibrium between three elements [11], namely, accuracy, transparency, and flexibility. Accuracy is the capacity to reproduce the real life behavior and a necessary feature to establish good control actions, perhaps the complexity of the model is proportional to the mentioned accuracy. On the other hand, transparency considers how the model components interact in the dynamics, and it has to be noticed that it is a property that is opposed to complexity and consequently to accuracy. Finally, flexibility is the ability to adapt the model to other circumstances.
If the model fits well with the real disease propagation dynamics, it can be used to infer optimal control actions, so it makes possible implementing more efficient measures and optimizing the use of resources. In 2001, the epidemic footand-mouth outbreak in the UK and the result of three different models applied to this case predicted a large-scale spreading and a better control of the epidemic if the infected and exposed animals were culled [26]. In many research studies, it has been considered a spreading of a very high death rate virus such as smallpox [11] to verify whether a massive vaccination, which could generate health problems, is advantageous or not. However, the obtained results differ from each other since there are many epidemiological uncertainties. As far as COVID-19 is a pandemic, there are many data sources that provide information about the epidemiological situation, so it is feasible to create a parametric model from real data and reduce epidemiological uncertainties. us, clear control strategies could be built up and implemented.
In this work, a deterministic mathematical model is built up, which spans the well-known SEIR epidemic model by adding people in quarantine and those being asymptomatic. To reduce the model dimension, the infected people with symptoms and asymptomatic will be gathered, so a SEIQR epidemic model is obtained which is integrated by the following subpopulations: Susceptible (S), Exposed (E), Infectious (I), Quarantined (Q), and Recovered or Immune (R). In this model, it is assumed that only exposed and infectious people will cause new infections since people in Q, which represent the individuals who stay in hospital and at home during some period due to the disease, are concerned about their illness, and they try not being in contact with others. On the other hand, it will be considered that some percentage of the individuals in the infectious subpopulation I, more precisely the symptomatic ones, will be reported and carried to Q, and it will trigger a tracing method as it is shown in [3,23]. us, we can say that there are the subpopulations including infected individuals which are the exposed, the infectious, and the quarantined subpopulations. In particular, the asymptomatic infected subpopulation is included in the exposed one.
After constructing the SEIQR model schematic, defining how individuals interact, and defining the necessary parameters, the respective nonlinear system of differential equations will be obtained and analyzed. e theoretical expression for the basic reproduction number, which is a feature dependent on the system's parameters and a threshold that determines whether the disease is out of control or not, will be obtained through the next generation matrix [27]. For stability analysis purposes, the equilibrium points, that is, the disease-free one and the endemic equilibrium one, which determine the states in which the disease is extinguished from the society and when the disease prevail (i.e., the usual flu disease), respectively, are calculated. e system is linearized around the equilibrium points, and the eigenvalue problem is solved to evaluate when their real parts become positive/negative (condition instability/stability of the system). Some of the model parameters (for instance, the virus incubation period and the probability to get infected by symptomatic/asymptomatic people) are a disease characteristic, so their values have been collected via other studies [2,3,12,13,27,29]. On the other hand, in [30], a model with susceptible, exposed, infectious, quarantined, and removed is discussed with parametrical values fitted to Saudi Arabia data which test the model and analyze the effects of hospital beds availability, quarantines, and media effects on the foreseen dynamics. In [31], a sliding-mode vaccination controller is designed to fight against COVID-19 and testing simulations are performed on data of Iran and Russia. In [32], machine learning designs are used in a susceptible-infected-recovered model to discuss the influence of vaccination in the COVID-19 pandemic evolution in Brazil as well as the impacts of the immunization speed and vaccines efficacy.
However, parameters that depend on the society behavior (i.e., the average number of contacts per day) must be determined with the respective real data. Since the autonomous community of Cantabria gathers wide data about the epidemiological situation, it has been used to estimate the unknown parameters, such as the average contact rate and the tracing effectiveness. It has been assumed that the discrete differential expressions of variables such as Q and I follow an ARX model structure.
us, the ARX model parameters have been estimated with the least square method and their values have been used to infer the unknown epidemiological values of the parameters. Finally, the differential equations have been solved with the obtained parameters, and the results have been compared with real data. e paper is organized as follows. Section 2 describes the new proposed epidemic SEQIR (Susceptible-Exposed-Quarantined-Infectious-Recovered) model. Such a model has several specific novel characteristics, whose effects are integrated and analyzed together in a coordinated way, namely, (a) e proportions of asymptomatic and symptomatic individuals in the quarantined group and in the infectious group have, in general, different average transmission rates of contagion to susceptible individuals as it is also potentially distinct from the above mentioned ones the average transmission rate of the exposed individuals to the susceptible ones. (b) Various average transmission rates are parameterized by the infective contacts between individuals. In particular, the values of the average transmission rates are directly proportional to the average potentially infectious contacts. ose contacts occur among the various infected types of individuals (that is, those being currently allocated in the exposed subpopulation, including the asymptomatic infective, in the symptomatic infectious subpopulation, and in the quarantined subpopulation) with susceptible individuals. (c) Various average transmission rates are also parameterized in a direct proportion by the corresponding probabilities of infecting susceptible individuals. (d) e model is also parameterized by the average rate of attendance to the doctor for disease positivity testing or evaluation of symptoms checking and by the effectiveness in tracing the positive cases. (e) e increase in the attendance rate of symptomatic suspect individuals contributes positively to the incorporation to the quarantined subpopulation of those being positively tested. (f ) e total day-to-day quarantined subpopulation incorporates the individuals who have been tested positive. is includes the hospitalized individuals including those under intensive care and also those quarantined at their homes having no serious symptoms.
e positivity of the model, in terms of the nonnegativity of all the components of the state trajectory solution under any given finite nonnegative initial conditions, is also proved and discussed in that section. Section 2 also includes a brief discussion on other approaches given for the pandemic study including the potential vaccination and treatment control relevance in the damping of the infection as well as the influences of the intervention measures in the pandemic evolution. On the other hand, Section 3 deals with the disease-free and the endemic equilibrium points which are proved to be unique. e disease-free equilibrium point is proved to be locally asymptotically stable if the basic reproduction number is less than unity, while the endemic one is locally asymptotically stable if such a value exceeds unity. In addition, the endemic equilibrium point is proved to be nonreachable in the first case since its infective components are negative so that the eventual convergence of the state trajectory solution to it would contradict the previously proved positivity of the model performed in the former section. Furthermore, it is proved that only one of the two equilibrium points is a global asymptotic attractor depending again on the value of the basic reproduction number. Section 4 gives some simulations with the proposed model being evaluated for the COVID-19 pandemic with demographic data from Spain. Section 5 develops in more detail a case study against COVID-19 which is performed with official data acquisition taken from the Cantabria Autonomous Community (a region in North Spain). In particular, the estimations of the average periods lasted by people to get to a Health Center or to a doctor for eventual positivity inspection, the efficiency of the positivity tracing, and the average number of contagion contacts are estimated from real data along the initial confinement period and also for a subsequent period. In this subsequent period, the formerly applied strong public intervention measures, of almost total general confinement of the whole population, were left while still public intervention rules were maintained, for instance, a duty of general use of masks, social distance keeping, partial restrictions of mobility, closing of nocturnal amusement places or very significant limitation rules on the number of attendees to spectacles and restaurants, and so on. Finally, conclusions end the paper. A subsection with a glossary of the main symbols is listed below, and the proofs of the stated mathematical results are given in appendixes.

Nomenclature
S, E, Q, I, R: subpopulations of susceptible, exposed, quarantined, infectious, and recovered individuals. It is assumed that the exposed subpopulation E includes the asymptomatic infective individuals, while the symptomatic infective (or infectious) integrated in the subpopulation I.
Discrete Dynamics in Nature and Society 3 N: total population Λ: average recruitment parameter β e , β a , β s : exposed to susceptible, asymptomatic to susceptible, and symptomatic to susceptible average transmission rates σ: average incubation rate c: average recovery rate ρ: average immunity loss rate δ: average quarantine rate μ: average natural death rate α: average disease mortality rate p 1 : average probability of evolution from exposed to symptomatic p 2 : average probability of evolution from exposed to asymptomatic λ: average rate of symptomatic visits to doctor for testing c e � β e /p e , c a � β a /p a , c s � β s /p s : average numbers of close contacts per day of a member from a group of (either exposed, asymptomatic or symptomatic infectious) to susceptible; and p e , p a , and p s are the respective probabilities to infect a susceptible individual t , p: daily total and positive tests r: average number of recovered quarantined individuals who were confined because positive testing

Model Setting and Its Description and Positivity
Epidemic models can operate in either deterministic or in stochastic frameworks depending on how the population of one block transforms into the others; that is, deterministic models consider that the rates (namely, the chance to be infected) maintain constant with respect to time, and therefore there is a unique solution for each initial condition. In contrast, stochastic models add random or probabilistic variables into the rates, so the model will give a set of probable solutions. Even if human behavior and infections have stochastic components, when the real system (stochastic) is composed by a large group of people, deterministic models could be used to represent the system. en, the ratios are constant and could be determined by statistical results [11]. e graphical representation of the system is made by a group of blocks (or compartments), each one represents a certain part of the population, and the arrows make the connection between two blocks (A ⟶ B) and specify the way in which part of the population from A transforms into the inhabitants of the block B. To build up the model flow chart, it is important to know the features of the tackling disease which are the stages of the disease, incubation period, immunity period, and so on.
Regarding COVID-19, the WHO reported that there are two main ways in which the disease is transmitted: presymptomatic and asymptomatic transmissions [20]. Presymptomatic transmission, which is possible during the incubation period (time between exposure and the symptoms onset) which is 5-6 days [17,20], is more likely to happen 1-3 days before symptoms appear [26]. e preliminary way of transmission is via symptomatic cases, but even if the percentage of asymptomatic cases is 16%-17%, [2,13], it has been seen that they show similar viral loads [10,13,26,27,29]. Some research studies have spotlighted the importance of presymptomatic and asymptomatic cases since they contribute to the virus high spread [27,29]. erefore, many countries were forced to include quarantine periods to reduce the spreading. e flow chart depicted in Figure 1 shows a deterministic model of COVID-19 which includes the characteristics mentioned before. Note that Figure 1 have been in contact with infected people and therefore are infected too, but they do not present symptoms because they are still incubating the virus. ey move freely; thus, they will maintain an average number of close contacts per day c free . (iii) Symptomatic (S y ): this group consists of people who got the disease and present symptoms. Soon after the incubation time, an individual of S y will probably move freely (the average number of close contacts per day will be equal to c free ) since the host manifests presymptoms, which commonly are confused with tiredness or symptoms of a common cold/flu. After the list of symptoms will become larger and their intensity will increase, therefore, he/ she will go to the doctor and put in quarantine. It is assumed that the average number of close contacts of people from S y is lower than c free .
(iv) Quarantined (Q): this group is integrated by subgroups coming from S and S y . People from S y who presented clear symptoms will have tested positive and put in quarantine. eir positive results will trigger a tracing through their close contacts, which is assumed that are individuals from S. Traced contacts will be tested, and those with positive result will be put in quarantine. It is assumed that the average number of contacts per day of people in Q is nearly zero because they are aware of their disease and they will interact less with others. Note that Q covers people who have tested positive and are staying at home or hospital; therefore, only people from Q might die due to the disease.  Discrete Dynamics in Nature and Society e newborns are introduced in the class S with a recruitment Λ, and from the classes, it removed a portion with a rate μ (daily natural death rate). It is assumed that the population is homogeneously mixed, so it is possible to determine the transmission rates β e , β a , and β s , which indicate the probability to infect susceptible people when they are in contact with exposed, asymptomatic, and symptomatic individuals, respectively. e corresponding expressions are β e � c e p e , β a � c a p a , and β s � c s p s where c i is the average number of close contacts per day of a member from exposed, asymptomatic, and symptomatic, respectively, to susceptible, while p e , p a , and p s are the respective probabilities to infect a susceptible individual. As it was mentioned before, it will be assumed that c e � c a � c free , c e > c s > 0, and p a ≈ p s > p e . ese are necessary parameters to specify how susceptible people get infected; one individual from S has a probability E/N, I/N, and A/N of contacting a person from E, I, and A, respectively, so the chance to get infected is (β e E + β s S y + β a A)/N. erefore, all susceptible people who will get disease per day is (1) Exposed people after the incubation time (σ − 1 days) will evolve to asymptomatic people with a probability p 2 . ose ones evolving to symptomatic individuals may not consider, at the early beginning, that they are still sick. However, after λ − 1 days (λ indicates the rate at which people from group S y visit the doctor, and it is supposed that its value will be larger as the test stock, sanitary resources, and the disease impact on society are increased), their symptoms will increase and they will go to the doctor, test positive, and put in quarantine. is will start a tracing, which effectiveness is determined by p 1 , of the positive tested symptomatic people with close contacts c s λp 1 S y S/N, and those who test positive β s λp 1 S y S/N will go to Q. erefore, the number of new exposed per day will be equal to the susceptible people who get the disease per day (see expression (1) minus traced people who have tested positive): People in Q include those who have been tested positive, which means hospitalized ones, people in Intensive Care Unit (ICU), and those who are in their homes. Some people from Q will remain there during the quarantine period (δ − 1 days), and others will die from COVID-19 with a rate α. After c − 1 days, both asymptomatic and symptomatic individuals will recover from the disease, and they will be part of the recovered population R during a period ρ − 1 . Now, six blocks are shown in the whole scheme from Figure 1(a), so that six differential equations are needed to represent the system. Managing with six equations might be a difficult task; therefore, a new compartment I involving together A and S y is defined to reduce the number of equations from six to five, that is, and in Figure 1(b), the block I (infectious) substitutes the blocks A and S y . In the following, the dependence on time is omitted in the equations for the sake of exposition and simplicity when no confusion is expected. e differential system representing the deterministic model is given by the following set of first-order differential equations: and R(0) � R 0 . e above model is being considered in the sequel for analysis. All the parameters are nonnegative and, furthermore, p 1 , p 2 ∈ [0 , 1]. e total population and its time derivative become e following result holds on nonnegativity of the trajectory solution for any given finite nonnegative initial conditions. e proof is given in Appendix A.

1.
Assume that en, the following properties hold: for all t ∈ R + so that the five subpopulations are nonnegative for all time (ii) e total population is nonnegative and bounded for all time and then the subpopulations are also bounded for all time.
Close techniques for evaluating the nonnegativity of the solution have been used in some background literatures for other epidemic models. See, for instance, [27][28][29][30][31][32]. Also, it can be pointed out that other related alternative mathematical models for infective disease have been proposed in the recent literature with applications to parameterizations related to COVID-19. For instance, in [33], mixed contagions of susceptible individuals from asymptomatic and symptomatic infectious ones are studied based on provided official data taken from some European countries. In [34], SARS-CoV-2 (COVID-19) transmission data from Spain and Italy were analyzed and processed to estimate the transmission rate levels through SIR models with undelayed and delayed propagation infection related to the periods of different public intervention measures. In [35], confinements and quarantines have been considered as impulsive actions which decrease at certain time instants, and along certain duration time periods, the amount of individuals is able to produce contagions. In [36], an analysis of the contagion propagation levels from small amounts of exposed and infected outsiders coming in a certain habitat was studied under eventual delayed resusceptibility. In [37], a model for COVID-19 including a potential herd immunity has been developed for Austria, Luxembourg, and Sweden. Also, the effect of delaying the vaccination while assuming unlimited vaccine units supply is studied in [38]. A very complete recent work on epidemic modeling is given in [39]. Such a well-organized work discusses different deterministic, stochastic, and optimization models and includes also a discussion on the pandemic evolution in France. On the other hand, some mathematical techniques of interest for its application in mathematical epidemic models related to positivity and stability as well as numerical methods, estimation methods, or probability and statistics tools can be found in [39][40][41][42][43][44][45][46][47].

Equilibrium Points and Stability Analysis
Now, the disease-free and the endemic equilibrium points of the above model are seen to be unique, and their components are characterized. It is also seen that the endemic equilibrium is not reachable (or, roughly speaking, it does not exist in view of eorem 1 (i)) for small values of the transmission rate since the infective components are negative.

Equilibrium Points.
e equilibrium points are found by zeroing the time derivatives of the subpopulations in (4)- (8).
e disease-free equilibrium point satisfies E dfe � Q dfe � I dfde � 0. Looking at (8), it is seen that those constraints imply that R dfe � 0 and, from the first one, S dfe � N dfe � Λ/μ so that x dfe � (Λ/μ, 0, 0, 0, 0) T is indeed unique. e interpretation is that, in the absence of disease, the whole population is susceptible. Note that, in the absence of recruitment, i.e., if Λ � 0, then the disease-free equilibrium point implies the population extinction. 6 Discrete Dynamics in Nature and Society e existence, uniqueness, and characterization of the endemic equilibrium point are discussed in the two following results which are, respectively, proved in Appendixes B and C.

Lemma 1. Assume that
en, the following properties hold: (i) We have (ii) We have Let us define the relative transmission rates related to some transmission rate reference value β leading to β e � β er β, β a � β ar β, and β s � β sr β. en, a necessary condition for the endemic equilibrium point to exist is that β > β c . (iii) e quarantined value at the endemic equilibrium point Q ee , if it exists (i.e., if Q ee > 0 ), satisfies the subsequent constraint: and (v) We have where Note that the reference transmission rate β is fixed for each given triple of transmission rates β e , β a , and β s with their relative values being β er � β e /β, β ar � β a /β, and β sr � β s /β. e following brief discussion concludes that a unique feasible endemic equilibrium point can also exist irrespective of the parameter λ which describes the attendance rate to the doctor or Health Service for checking or treatment.
Discrete Dynamics in Nature and Society Lemma 1 leads directly to the next result whose proof is given in Appendix C.

Theorem 2.
Assume that λ ∈ [0, λ c ] and β > β c . en, the endemic equilibrium point is unique subject to an endemic equilibrium population given by And, as a result, there is a unique reachable endemic equilibrium point given by (11), (12), (14), (15), and (17). If β � β c , then the endemic equilibrium point coincides with the disease-free one.

Remark 1.
eorem 2 establishes that the endemic equilibrium point does not exist for sufficiently small transmission rates so that, in this case, the disease-free equilibrium point is the only feasible equilibrium point. For instance, note that if the three transmission rates are zero, then the susceptible individuals at the endemic equilibrium would be infinity what contradicts eorem 1(ii) so that, for zero transmission rates, the endemic equilibrium does not exist, as expected.
Remark 2. Note that it becomes attractive to refer the three defined transmission rates to a prescribed reference one β by using the relative values of the defined one with respect to such a reference transmission rate. In particular, it is easy in that way to characterize the existence of the endemic equilibrium point referring to such a reference transmission rate exceeding a threshold which depends on the remaining parameters including the three relative transmission rates. It will be seen later on in Section 3.2 that the necessary condition β > β c for the reachability of the endemic equilibrium point (equivalent to its existence from its necessary positivity for that concern) is also equivalent to the basic reproduction number to fulfill R 0 > 1, which implies in turn the instability of the disease-free equilibrium point. us, one will conclude that if the disease-free equilibrium point is locally asymptotically stable, then the endemic equilibrium one is not reachable. Simulations corroborate that if R 0 � 1, then the disease-free equilibrium point is locally asymptotically stable. As a result, the disease-free equilibrium point is a global attractor if R 0 ≤ 1 and equivalently if β ≤ β c .
Remark 3. Note that the sufficient conditions λ ∈ [0, λ c ] and β > β c of eorem 2 for the uniqueness of the endemic equilibrium point are independent of the equilibrium subpopulations amount since it is based on the fact that I ee < N ee .
For the particular case, λ � 0 which reflects that people, in average, do not go to the doctor under symptoms since the delayed time to go to the doctor is λ − 1 � +∞. For such a case, it can be deduced as well that the endemic equilibrium is unique. However, the reasoning has to be specific for the discussion of the particular model (4)-(8) obtained by zeroing the terms affected by a factor λ in the differential equations. Since λ � 0 yields to division by zero in some relevant equations in the proof of Lemma 1, the next related result is stated and proved separately from Lemma 1 and eorem 2. Its proof is given in Appendix D.
Theorem 3. If λ � 0, then there are a unique disease-free equilibrium point and a unique endemic equilibrium. e first one is x dfe � (Λ/μ, 0, 0, 0, 0) T ; similarly for λ ≠ 0, and the endemic one satisfies the subsequent relations: where Furthermore, the endemic equilibrium point is reachable if the reference transmission rate exceeds a minimum threshold according to

Local Asymptotic Stability Analysis.
e next technical result, proved in Appendix E, relies on the evaluation of the Jacobian matrix with respect to any equilibrium point. It will be then used to characterize the local asymptotic stability of both equilibrium points.
with entries as follows: 8 Discrete Dynamics in Nature and Society e above Jacobian matrix can be decomposed as J * � T * − V with the transmission and transition matrices, respectively, − V (being independent of the transmission rates and of the considered equilibrium point) and by using the relative transmission rates β er � β e /β, β sr � β s /β, and β ar � β a /β with respect to the reference transmission rate β.
e above result is useful for the local stability analysis around the disease-free and endemic equilibrium points, respectively, x df and x ee , via the respective Jacobian ma- In the calculation of the entries of the Jacobian matrix around x * , it has been taken into account that the total population is not constant with respect to time due to the mortality associated to the disease. For instance, generates one of the additive terms in J * 11 and so on for the remaining entries of the Jacobain matrix around the considered equilibrium point. e following remark gives some considerations of conceptual interest related to the Jacobain matrix around the equilibrium points and its partition into the sum of the transition and transmission matrices.
Remark 4. Note that the transmission matrix (− V) is a Metzler stability matrix, i.e., its diagonal entries are negative and its off-diagonal entries are nonnegative so that all its eigenvalues are real negative and its associated fundamental (or transition) matrix e − Vt is nonsingular and has nonnegative entries, for all time. Note also that the transmission matrix T * is zero if the transmission rates are zero (i.e., it is zero in the absence of disease). e fact that e − Vt has nonnegative entries for all time implies that any solution trajectory is nonnegative for all time under any given nonnegative initial conditions, in the absence of disease which is a necessary condition for the nonnegativity of the whole model to have a nonnegative solution trajectory under nonnegative initial conditions ( eorem 1).
On the other hand, is stable if and only if the spectral radius of V − 1 T * , identical to that of T * V − 1 is small enough. According to equation (23), this fact happens for the disease-free equilibrium point if β is small enough. It can be proved that this is equivalent also to β < β c and it will be seen later on that this is also further equivalent to the basic reproduction number to be less than one. Finally, according to (10), this also implies that the endemic equilibrium point is not reachable.
Local asymptotic stability and instability results around the equilibrium points are stated in the subsequent result, whose proof is given in Appendix F.

Theorem 4.
e following properties hold: (i) e Jacobian matrix around the disease-free equilibrium point is stable if and only if the basic reproduction number with the average critical reference transmission rate β c being defined in (10) and then the diseasefree equilibrium point is locally asymptotically stable. Moreover, the spectral radius of V − 1 T df equalizes the basic reproduction number where, such that for the disease-free equilibrium point, the transmission and transition matrices are as follows:

Remark 5.
It is now discussed how the calculation of the basic reproduction number can be simplified by manipulating the infective variables, that is, the exposed, infectious, and quarantine subpopulations, only. Note from equations (27)-(28) that the Jacobian matrix at the disease-free equilibrium point is as follows: which makes obvious by its direct inspection that the linearized subsystem of the infective components E, I, and Q around the disease-free equilibrium points is not coupled to the linearized subsystem of noninfective components S and R. Note that the Jacobian matrix of the infective linearized subsystem around the disease-free equilibrium point is given by where the corresponding transition and transition matrices are as follows: which concludes that the local stability around the diseasefree equilibrium point holds if R 0 � ρ(T i df V i df ) < 1 which guarantees that J i df is a stability matrix since − V i df is a stability matrix of eigenvalues − (μ + σ), − (μ + α + δ), and − (λ(1 − p 2 ) + μ + c). erefore, the local stability around the disease-free equilibrium point can be analyzed more easily from the Jacobian matrix of the infective linearized subsystem since if E(t), I(t), Q(t) ⟶ 0 as t ⟶ ∞ for any initial conditions sufficiently close to the disease-free equilibrium point, then the disease-free equilibrium point is locally asymptotically stable and conversely this happens if the eigenvalues of J i df have negative real parts.
Discrete Dynamics in Nature and Society

Global Asymptotic Stability Analysis.
e subsequent result extends eorem 4 by proving that the given asymptotic stability results are also global and that only one of the two equilibrium points is a global asymptotically stable attractor depending on the range of values of the basic reproduction number, equivalently, on that of the reference transmission rate for any given fixed triple of relative transmission rates. Its proof is given in Appendix G.

Computer Simulations on the
Epidemic Evolution e subsequent parameterization data to test the epidemic model for the COVID-19 pandemic correspond to Spain. It is assumed that the contact rates c e , c a , and c s which conform the respective transmission rates β e , β a , and β s are supposed identical and equal to c r . e transmission probability of the exposed p e is 10 − 3 times lower than that of the symptomatic p s which equalizes that of the asymptomatic individuals p a . e average time before an infectious individual goes to the doctor is λ − 1 � 2 days. e tracing effectiveness p 1 is 0.5 while the average time of staying in quarantine is equal to δ − 1 � σ − 1 + c − 1 . e remaining model parameters and the initial conditions are displayed in Tables 1 and 2. e obtained simulation results are depicted in Figure 2(a); it shows the dynamics of the variables S, R, Q, E, and I and the cumulative number of deaths with respect to time. e variables show characteristics of an underdamped dynamical system; its response is a sum of a transient response (which is damped sinusoidal) and the steady state. In this case, the disease persists, which is an expected result since R 0 � 2.36. It is important to note that, approximately, 300 days after the outbreak, the total number of deaths reaches approximately 1 million, and at the end of the simulation (600 days after the outbreak) is doubled. is high mortality can be attributed to the high number of contagious contacts. e disease spreading can be decreased by reducing the value of c r � c e � c s � c a ; therefore, in an alternative simulation, c r was reduced to 4 (R 0 � 0.89), and as it can be seen in Figure 2(b), the disease tends to disappear; the variables Q, E, and I show an overdamped response, and their values continuously decrease with respect to time. As long as their values and the total number of deaths are small compared with S long over the simulation, the variable S seems to be unchanging (see the top graph of Figure 2(b)). e results of both above simulations correspond to the expected behavior in the cases of which the basic reproduction number is greater and lower than unity. Besides, considering the cumulative number of deaths of both simulations, it is clear that reducing the value of the parameter of average contact rate c r could prevent a substantial increase in the total deaths.
On the other hand, Figure 3 shows the disease evolution as the basic reproduction number is unity of close to unity. It is seen that if the reproduction number is unity so that the disease-free and the endemic equilibrium points coincide, the disease-free is asymptotically stable as in the case when it takes values below unity as expected from the theoretical discussion performed in the former section. It is also seen that, for values of the basic reproduction number which are slightly higher than unity, the endemic equilibrium is a globally asymptotically stable attractor.

Data Acquisition.
is section evaluates the proposed model with real discrete data from Cantabria Community. e tuple of daily discrete data of the k − th day is defined y(k) � (t(k), p(k), d(k), r(k)) where t(k) is the number of performed tests, p(k) is the number of positive tested results, d(k) is the number of deaths, and r(k) is the number of recovered individuals. Over the above array, the subindices d and c are used for daily and cumulative data. From cumulative data, one can obtain the daily ones by using the derivative of the accumulated ones by Euler forward approximation which results in y d (k + h) � (y c (k + h) − y c (h))/h, and typically h will be taken as a period of one day. Some relevant data from Cantabria, from February 29, 2020, to February 14, 2021, are displayed in Figure 4.
It can be pointed out that it is of interest to evaluate from the data the number of individuals in quarantine, which also includes the hospitalized ones since, approximately, the deaths are proportional to the quarantined individuals, i.e., d(k) � αQ(k) and the fraction of recovered individuals is proportional to the quarantined individuals, i.e., r(k) � δQ(k). us, by discretizing (8) for one day period, one gets In the first propagation wave of the disease, since the number of infectious I(k) individuals compared with recovered individuals R(k) is close to zero and ρ + μ ≪ 1, it can be stated that r(k) is approximately the daily increment of the recovered since from the above relation, r(k) � δQ(k) � R(k + 1) − R(k). Note that r(k) is a part of the recovered individuals, that is, those recovered individuals from quarantine, since most of the recovered individuals from the slight or the asymptomatic infection are not accounted for the data processing. erefore, we keep in the following lowercase notation r(k) for this partial account of the estimated recovered daily increment R(k + 1) − R(k). e period from February 29, 2020, to February 14, 2021, gave a cumulative of recovered and death people which was 25,500 which is only a 5% of the population of Cantabria. erefore, S/N has been taken unity along that period. Regarding the number of positive tests, death individuals, and recovered tested daily individuals p(k), d(k), and r(k), some average delays τ 1 , τ 2 , and τ 3 exist from the contributions of the daily infectious and quarantined individuals caused by the transmission process and the delays in data processing. e used evolution equations are as follows: where the superscript hat stands for the various estimates and w p (k)      moving-average filter which takes 14 samples has been applied to (33)- (35). By combining the above equations, one gets for the discretization of the estimate of the quarantined individuals through (6) for h � 1 via the Euler forward approximation that

Estimation of the Delays.
Note that as far as the coexistence with COVID-19 increases, it is expected faster responses in the performance of positivity tests and public supply of data. erefore, the potential estimation of delays can be adapted to various disease evolution periods. Now, define τ 4 � τ 2 − τ 3 and τ 5 � τ 1 − τ 2 . Algorithm 1 has been performed to estimate them from the registered data.
Step 2. For each subset, seek for delays between r(.) and d(.) (in particular, the delays between the corresponding additive contributions of d(k + τ 2 ) and r(k + τ 3 ) to (45)) with MATLAB finddelay() (for delays between signals) to get τ (i) Step 5. Calculate (36) under the equivalent form as follows: Step 6. Apply the linear regression method to the data , d(k) and compute the coefficients of determination R (2r) ℓi and Step 8. While ℓ ≤ 2n do ℓ←ℓ + 1 and go to Step 4.

Remark 6. Both coefficients of determination R (2r)
ℓi and R (2d) ℓi , i � 1, 2, are the proportions of variability in the set of given data due to prediction [35,36], and they range from zero to unity. e data become better fitted as the corresponding coefficient of determination becomes close to unity. In contrast, values close to zero lead to a high variability and a bad prediction as a result.

Estimation of Average Mortality Rate and Stay in Quarantine Periods.
e estimations of α i and δ i , i � 1, 2, are displayed in Table 3. Note that α 1 > α 2 as expected from a higher hospitalization percentage along the first studied period. Regarding the recovery rate, after the first wave, tracing techniques were implemented when someone with symptoms went to hospital and his/her close contacts were tested to detect new infections, and those tested positive were quarantined which can explain that δ 2 < δ 1 . Figures 6 and 7 display the linear regressions along the first and second evaluation time periods to estimate the mortality rates and stay in quarantine period.

5.4.
Estimations of λ , c r , and p 1 . By approximating S/N to unity, we take equations (5) and (7), the infective subsystem of the epidemic model, and we rewrite the set of third-order equations of first-order as a differential equation of thirdorder with the quarantined subpopulation as variable. is leads to Q ... where where p e � β e /c r and it has been considered that p a � p s � β s /c r depending on the average transmission rates β e � β a and β s and that the average contact rates satisfy c r � c e � c s � c a . It has been also assumed that λ is the average rate of visits to the doctor or to hospital for the particular disease COVID-19. Now, apply the Euler forward approximation to (38) for discretization which leads to and equalize the above expression by equating corresponding coefficients to the prediction values of the quarantined individuals given by   which leads to the system of algebraic equations J 1 x 1 � v 1 with x 1 � (λ, c r , λc r ) T which are in expanded form: where the parameters which are not unknowns are taken from Table 1. us, since the third row is not independent of the two former ones, one takes the minimum square error solution (λ , c r ) of (42) given by λ, c r � Arg min with stopping computational value error less than 10 − 4 . Now, proceed in the same way to get a second-order differential equation of I(t), by using the combination of (5) and (7), resulting in with By applying the Euler forward discretization method, one gets In this case, one gets proceeding in the same way as done with (40) to get (42) the algebraic system J 2 p 1 � v 2 with the equivalent explicit expanded scalar relationships as follows:

18
Discrete Dynamics in Nature and Society And the estimation p 1 of the effectiveness tracing p 1 is calculated by replacing the necessary values of Table 1 and the solution (λ , c r ) of (43) results to be p 1 � Arg min with stopping computational error less than 10 − 4 . e solution is found through the Newton-Raphson method [42,43].

Relationship between Positive and Total Tests
Interpreted via the Estimations of c r and p 1 . e total and positive tests are described, respectively, by As a consequence of testing symptomatic people with COVID-19 (1 − p 2 )λI(k), we traced contacts c r λp 1 (k)(1 − p 2 )(IS/N) and other people that present similar symptoms to COVID-19 but are not infected ε(k). By taking S(k)/N(k) � 1, one gets where e estimations of the values of λ, c r , and p 1 for each one of the mentioned time intervals are displayed in Table 4. e smaller values of the estimation of λ were found in the period from 19/12/2020 to 06/01/2021 which matches with winter holidays when a significant number of people moved to their relative houses so that people under infection symptoms have found more difficulty to go to the outpatient. Regarding the average number of contacts per day c r , as tighter restrictions were applied, smaller values of c r were observed. e smallest value c r (3) � 2.07 is found from 19/03/2020 to 28/05/2020 approximately as the quarantine was imposed while the next closer smaller one c r (12) � 3.14 was reached from 25/11/2020 to 13/12/2020 as the state of alertness was still valid. us, these values are closely related to the levels of imposed restrictions.
Note that the parameter λ − 1 of time of meeting doctor or health service for testing varies approximately from 1.58 days to 1.07 days. Assuming that, at first, COVID-19 was underestimated and its symptoms were mixed up with cold/ flu ones, infected people with symptoms might typically have to visit the doctors at lower rates. On the other hand, holiday seasons could also be associated with lower values of λ since people might be outside their usual homes so that health assistance might be reduced. e tracing effectiveness p 1 was never before implemented in Spain on big scale cases such as in this pandemic; therefore, there has been notably uncertainty regarding how to implement it. In that way, as the tested positivity was high, people were tested randomly. e variances of the parameters λ(k), c r (k), and p 1 (k) are calculated through the values of Table 4. On the other hand, Figures 8-10 display the quarantined, death, and recovered individuals in both periods compared with their estimations taking account of the variances of the involved parameters, for instance, and concerning Figure 8, the corresponding variances, or its associated typical deviations, of Q(k) are obtained via (36) as follows. For each k − th sample, corresponding to a particular day, one takes data from two weeks (that is, 14 days) neighboring days, namely, the former seven ones, the current ones, and the next six ones so there is one delay of six days to provide the evaluation of the mean value of the involved data over each day. One proceeds in this way for the estimations on the mentioned groups of two weeks for the estimations of d(k), r(k), and p(k). With those data, one calculates the k-th variance, or its associate typical deviation, of the estimation of the quarantined subpopulation Q(k). Figure 11 shows a linear regression of total tests t(k) versus positive tests p(k) in two periods of time from 01/03/ 2020 to16/07/2020 and from 17/07/2020 to 02/02/2021. e adjustment of the corresponding linear regression equation t(k) � r a (k)p(k) + ε(k) via least-squares gives the pairs (r a (k) , ε(k)). Note from (54) that r a (k) ∈ (0, 1] and that t(k) � r a (k)p(k) + ε(k) ≥ p(k) so that ε(k) ≥ (1 − r a (k))p a (k) ≥ 0 so that, for potentially got negative values of r a (k) and ε(k), the corresponding data are not collected.  # death people dp (k -τ 2 ) dp (t -τ 2 ) dp (t -τ 2 ) ± σ dp (t -τ 2 ) --dp (k -τ 2 ) dp (t -τ 2 ) dp (t -τ 2 ) ± σ dp (t -τ 2 ) -- Figure 9: Comparison between measured death people per day (black dots) with respect to the estimation (blue solid line) carried out with the estimated parameters. e dashed red lines show the standard error of the simulation.

Conclusions
is paper has proposed a SEQIR model which incorporates a quarantined subpopulation to those of the standard ones in the common SEIR-type models. e quarantined subpopulation contains the individuals who have been tested positive.
is includes the hospitalized individuals including those under intensive care and also those quarantined at their homes having no serious symptoms. e model parameterizes the average transmission depending on the infective contacts as the average rate of attendance to the doctor for disease positivity testing or symptom evaluation checking and the average effectiveness tracing on positive cases. e proposed model considers that the proportions of asymptomatic and symptomatic individuals in the quarantined and the infectious subpopulations can potentially    have different average transmission rates of contagion to the susceptible group. Also, various average transmission rates are assumed to be parameterized by the infective contacts between individuals which allows to evaluate the infection transmission reduction under public intervention measures like confinements or mobility restrictions, social distance maintenance, or use of masks. A main novelty of the proposed model is that it is parameterized by the average rate of attendance to the doctor for disease positivity testing or evaluation of symptoms checking and by the effectiveness in tracing the positive cases. It is proved that the increase in the attendance rate for testing of symptomatic suspect individuals contributes positively to the incorporation to the quarantined subpopulation of those being positively tested so that the infection transmission might be controlled more efficiently.
e model is also proved to be positive, an obvious need for the coherency of any epidemic model, in the sense that the state trajectory solution is seen to be nonnegative for all time for any given set of nonnegative initial conditions while it is proved to have unique disease-free and endemic equilibrium points. e disease-free equilibrium point is proved to be locally asymptotically stable if the basic reproduction number is less than unity, and under these conditions, the endemic one is not reachable, while it is locally asymptotically stable if such a value exceeds unity. Both points are globally asymptotically stable under the abovementioned conditions, and they are locally asymptotically stable. e mathematical condition for the basic reproduction number to be less than or larger than unity is proved to be equivalent to the transmission rates to be less than or larger than a certain critical transmission rate under the assumption that the remaining model parameters are kept unmodified. is helps to interpret the eventual transmission attenuation in terms of reduction of the disease transmission rates, implying as a result a decrease in the infective contagion contacts, as being dependent on the potential intervention measures. e model is tested by computer simulations for COVID-19 with data from Spain. Subsequently, a detailed case study was performed against COVID-19 through official data taken from the Spanish Cantabria Autonomous Community. In particular, the estimations of the average periods lasted by people to get to a health center or to a doctor for eventual positivity inspection, the positivity efficiency of tests, and the average number of contagion contacts are estimated from real data under the confinement period and also for a subsequent period. In this second period, the first strong public intervention measures of almost total general confinement were left, while still the population was under severe rules like general use of masks, social distance keeping, closing of nocturnal amusement places, or very important limitation rules on the number of attendees to spectacles and restaurants.
In the proposed model, there is a good trade-off between accuracy, transparency, and flexibility. Even if the quarantined subpopulation gathers people who were tested positive and were staying either at home or in hospital, reducing so the model complexity, a quantitative and qualitative fit of data, with respect to real data, has been checked. Regarding the transparency, the calculated basic reproduction number shows a simple relation with respect to parameters such as the average number of close contacts per day or the tracing effectiveness. e proposed model has been then evaluated with real data acquisition from the Spanish Autonomous Region of Cantabria. From the data source, the number of positive test result per day and the number of death people due to the disease and recovered people were gathered.
A noteworthy delay between the quarantined with respect to death and recovered can be interpreted in terms of delays in providing data and in the eventual attendance for checking eventual positivity testing. e analysis of the imposed restriction with respect to time gives an overview of how the limitations might affect the value of contagion contact rates.

A. Proofs of the Mathematical Results (Proof of Theorem 1)
It can be directly seen from (4) Now, proceed by contradiction. Assume that there is some t 1 ≥ 0 for which at least one component of x(t 1 ) is negative and x: [0, t 1 ) ⟶ R + 5 with x(0) � x 0 ∈ R + 5 . In view of (4)-(8) and since the trajectory solution is everywhere continuous on R because it is everywhere time-differentiable on R, one concludes that for t − 1 � lim ε⟶0 and all ε ∈ (0, t 1 ) since those components are negative at t 1 and nonnegative for t ∈ [0, t 1 ). It also follows that this event can only happen for any nontrivial trajectory solution if t 1 > 0 since, otherwise, the claimed ε ∈ (0, t 1 ) such that x i (t 1 − ε) ≥ 0 and x i (t − 1 ) � 0; ∀i ∈ 5 t 1 cannot exist for t 1 � 0. However, in view of (A.1), _ x i (t − 1 ) ≥ 0 and ∀i ∈ 5 t 1 so that x i (t 1 ) ≥ 0 and ∀i ∈ 5 t 1 and then x i (t 1 ) ≥ 0 and ∀i ∈ 5. So, there is no t 1 ≥ 0 (and then 5 t 1 � ∅) for which some component of x(t 1 ) is negative if x 0 ∈ R + 5 . Property (i) has been proved. Now, one has from (9) and Property (i) since Q(t) ≥ 0; ∀t ∈ R + that _ N(t) ≤ Λ − μN(t) and Furthermore, again for Property (i), one has that N(t) � S(t) + E(t) + I(t) + Q(t) + R(t) � 5 i�1 x i (t) ≥ 0 and x i (t) ≥ 0; ∀t ∈ R + , ∀i ∈ 5.

B. Proof of Lemma 1 and Associated Remark
Property (i) follows directly by zeroing _ I(t) in (7). Property (ii) follows by zeroing _ E(t) in (5) that S ee � (σ + μ)N ee E ee β e E ee + p 2 β a + 1 − p 2 1 − p 1 λ β s I ee , (B.1) which leads to (12) from (11) since E ee ≠ 0. On the other hand, note that the total population cannot asymptotically extinguish if Λ > 0 either for the disease-free equilibrium point (from (9) sinceQ dfe � 0) or for the endemic one since the endemic equilibrium infected subpopulations should be positive. en, for E ee > 0, S ee < N ee since the sum of all infected subpopulations at the equilibrium is necessarily positive so that one gets from (10) and (12) that g S � S ee N ee � (σ + μ) λ 1 − p 2 + μ + c λ 1 − p 2 + μ + c β e + p 2 β a + 1 − p 2 1 − p 1 λ β s σ < 1.

G. Proof of Theorem 5
Note the following facts under the assumption λ ∈ [0, λ c ]: Fact a. If R 0 < 1, then the endemic equilibrium point is not reachable (or it does not exist in the first orthant of the state space) since it is not positive while the diseasefree one is locally asymptotically stable. So, for R 0 < 1, there is only a unique locally stable attractor and any potential limit cycle surrounding it at a phase plane of any two components, if any, would be unstable (so, it would vanish under any perturbation which would not be detectable) because of the stability/instability alternation between equilibrium points and limit cycles from Poincaré's theory of qualitative differential equations concerning the presence of mixed equilibrium point and limit cycles. Also, no limit oscillation could affect to the disease-free equilibrium values of the exposed, infectious, or quarantined individuals since it would reach negative values contradicting eorem 1. As a result, the disease-free equilibrium point is globally asymptotically stable if R 0 (β) < 1(β < β c ). Fact b. If R 0 > 1, then the endemic equilibrium point is unique, reachable, and locally asymptotically stable while the disease-free one is unstable. Potential limit cycles at any plane of any two components, if any, should surround it but not both equilibrium points since their joint Poincaré index would be "+2" (if the disease-free equilibrium point is not a saddle point) or "0" (provided that the disease-free equilibrium point is a saddle point). Such a limit cycle should also be unstable in the case of existence since the endemic equilibrium point is locally asymptotically stable so that it could only be potentially surrounded by an unstable limit cycle in any phase plane of two components. us, the endemic equilibrium point is a globally asymptotically stable attractor if R 0 > 1(β > β c ). Fact c. e equilibrium point is still unique for R 0 � 1 because both equilibrium points coincide.
As a result of the above Facts a-c, there is only a globally asymptotically stable attractor which is the disease-free equilibrium point if R 0 (β) ≤ 1(β ≤ β c ) and the endemic one if R 0 > 1(β > β c ).

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.