Analysis of the COVID-19 epidemic in french overseas department Mayotte based on a modified deterministic and stochastic SEIR model =================================================================================================================================== * Solym Mawaki Manou-Abi * Julien Balicchi ## Abstract In order to anticipate a future trends in the development of the novel coronavirus COVID-19 epidemic starting early at march 13, in the french overseas department Mayotte, we study in this paper a deterministic and stochastic epidemic model through several possible states: susceptible, exposed, infected including symptomatics and asymptomatics cases, simple or mild removed, severe removed (including hospitalized) and death cases. We first fit the constant transmission rate parameter to the epidemic data in Mayotte during an early exponential growth phase using an algorithm with a package *R*0 and based on the ML method. Thus, after the intervention and control measures imposed by the gouvernement, we introduce a simple temporally varying and exponentially decreasing transmission rate with a control parameter. This allows us to understand wheter these interventions are having the desired impact of controlling the epidemic by pointing out some percentages to maintain control which is critical in Mayotte given the fragility of its health infrastructure and the significant fraction of the population without access to water. Key words and phrases * SEIR compartmental model * Coronavirus * Mathematical model * Reproductive number * Epidemiology ## 1. Introduction A novel coronavirus *Covid* −19 spread from the capital of the Hubei province in China to the rest of the country, then to most of the world and in particular the french overseas department Mayotte. Following the emergence of this novel coronavirus and its spread outside of China, Europe is now experiencing large epidemics. In response, France and many European countries have implemented unprecedented non-pharmaceutical interventions including case isolation of symptomatic individuals and their contacts, the closure of schools and universities, banning of mass gatherings and some public events, and most recently, widescale social distancing including local and national lockdowns of populations with all but essential internal travel banned. Also in response to the rising numbers of cases and deaths, and to maintain the capacity of health systems to treat as many severe cases as possible, France like European countries and other continents, have implemented some process measures to control the epidemic including his overseas department outside the metropolitan France. The epidemic began slightly later in France, from February. To prevent further spread of *Covid*− 19, all travel into Mayotte for tourism and family visits was prohibited from Friday, March 20. Only essential travel for selected individuals, such as residents returning home, professionals from essential services, and patients with health conditions was allowed to entry. Local authorities have advised all travelers arriving from the mainland France or areas affected by *Covid* −19 to self-isolate at home for 14 days and avoid contact with others. In addition, the government has urged the public to observe good personal hygiene, such as regular hand washing. Individuals who develop respiratory symptoms, such as cough, fever, and difficulty breathing, within 14 days of arrival into Mayotte are advised not to visit the doctor or emergency room directly, but to contact the SAMU (French Emergency Medical Services) immediately by dialing 15. Further international spread of the virus was expected until March 29; but from March 30 the government has decided to cancel all flights both departure and arrival in Mayotte. Many mathematical models of the COVID-19 coronavirus epidemic in China, USA, Italy and France have been developed, and some of these are listed in th following papers []. Mathematical models can be defined as a method of emulating real life situations with mathematical equations to expect their future behavior. In epidemiology, mathematical models play role as a tool in analyzing the spread and control of infectious diseases. It is well known that one of the most useful parameters concerning infectious diseases is called basic reproduction number. It can be specific to each strain of an epidemic model. To anticipate future trends in the development of the epidemic in Mayotte starting early at March 13, we explore in this paper a deterministic and stochastic epidemic model to describe the number of cases. To this end, we describe the evolution of an individual through several possible states: susceptible, exposed, infected including symptomatics and asymptomatics cases, removed (isolated or cured), hospitalized and death cases. Estimating the basic reproduction numbers for SARS-CoV-2 presents challenges due to the high proportion of infections not detected by health systems and regular changes in testing policies, resulting in different proportions of infections being detected over time and between countries. Most countries so far only have the capacity to test a small proportion of suspected cases and tests are reserved for severely ill patients or for high-risk groups (e.g. contacts of cases appeared in Mulhouse). Looking at case data, therefore, gives a systematically biased view of trends. However, in France and other contries in Europe the basic reproduction numbers for SARS-CoV-2 is betwen 2 and 3 at the begining of the epidemic with an average around 2.5. In our case, firstly we fit the constant transmission rate parameter of the model at the begining of the epidemic in Mayotte to observed data. After that we focused on the effects of the France government imposed public policies designed to contain this epidemic in his overseas department Mayotte from March 20, wich correspond to just a few days after that effective date of control measures which really started in practice. We assume that from March 20, the transmission rate parameter is not only time-depending and exponentially decreasing but also does not affect the recovery rate unless the disease is curable, which is not the case for Covid-19 actually. Our model can thus be used to project the time-line of the model forward in time. We then propose various scenarios on this control to show the impact of strict respect of the government measure can lead to have a control on the capacity of the hospital staff in Mayotte. The novelty of our paper is based on the new model (1) This paper is organized as follows. Section 2 deals with some preliminaries and basis aspects of the deterministic SEIR model intended to clarify the compuations and simulayions. ## 2. Preliminaries and basis aspects of the deterministic SEIR model A deterministic susceptible-exposed-infectious-recovered (SEIR) model for infectious diseases is developed here with the aim to make experimental simulations in next sections. Because the classical model assumes that the infected person’s incubation period is not infectious, this assumption is quite different from the infection characteristics of the new coronavirus infection. Therefore we will use a revised SEIR model to analyse and predict the trend of the outbreak. The model meets the following assumption. (H1) The mortality rate *µ* induced by the disease is considered. (H2) The transmission method is person-to-person and after a short infectious period individual becomes permanently immune. (H3) Once exposed, individuals go through a latent period ![Graphic][1] of approximately 6 days after which they become infectious for a period ![Graphic][2] that is estimated to be 8 days for slight or mild infectious cases and for a period ![Graphic][3] for severe infectious cases that is estimated to be 21 days. (H4) We assume that there exists a proportion *p* of mild infections and a mortality rate *µ*. The model divides the total population into six categories: * Susceptibles **S** not yet infected but a risk of being infected * Exposed **E** infected but not yet be infectious, in a latent period; Infectious * **I** confirmed to be infectious (with infectious capacity and not yet be quarantined) and contain reported and unreported cases. * Removed with mild infections **R***M*, Removed **R***S* with severe infections and **Death** death cases. Define *S*(*t*) *E*(*t*), *I*(*t*), *R**M* (*t*), *R**S*(*t*) and *D*(*t*) as the number of susceptible, exposed, infectious, removed with mild infections, severe infections with hospitalization and death cases in the population at time *t*, respectively. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/17/2020.04.15.20062752/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/F1) Figure 1. Compartments of the SEIR model Consider a time interval [*t, t* + *dt*] where the small change *dt* represents the length between the two points at which measurements are taken. Let *β* denote the disease transmission coefficient. This can be though as the rate at which each infectious individual makes potentially infections contacts with each other individual, where a potential infectious contacts will transmit the disease if the contact is made by an infectious individual with a susceptible individual. It is the product of the contact rate and the probability of infection given a contact per unit time. The number of new infections in the time interval [*t, t* + *dt*] is ![Formula][4] where *N* is the total population size at the beginning date *t* of the epidemic. Given initial conditions *S*(0) = *N*, initial average of exposed individuals *E* > 0 and *I*(0) = *I* > 0, the model consists of the following system of ordinary differential equations: ![Formula][5] The total population size is denoted by ![Formula][6] so that the population size remains constant. Remark 2.1. *It is worth to mentionned that the initial average of exposed individuals in front of an infectious individuals E* *is not E*(0). One can also look, by this model, to keep track of the cumulative number of Covid-19 cases from the time of onset of symptoms for *C*(*t*) which is not a compartent and defined as follows: ![Formula][7] With contact tracing, a proportion of individuals exposed to the virus is quarantined. The quarantined individuals can either move to the compartment *E* or S q, depending on whether they are effectively infected or not., while the other proportion, 1 q, consists of individuals exposed to are effectively infected or not [20], while the other proportion, 1 q, consists of individuals exposed the virus who are missed from the contact tracing and move to the exposed compartment, E, once to the virus who are missed from the contact tracing and Remark 2.2. *One can assume that some international infectious travelers individuals y*(*t*) *negligible in front of the total population size S*(*t*) (*S*(*t*) + *y*(*t*) ∼ *S*(*t*)) *but not negligible in front of Infectious individuals I*(*t*) *subsequently travel at time t and are eventually detected in their destination to Mayotte*. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/17/2020.04.15.20062752/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/F2) Figure 2. Compartments of the SEIR model including exposed travellers ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/17/2020.04.15.20062752/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/F3) Figure 3. A map of Mayotte Hence this leads to the new equations of the model: ![Formula][8] The dynamic of new exposed travellers *y*(*t*) can be thought as follows: ![Formula][9] where *α* is the average number of imported exposed individuals. The following table lists the parameters of the model we’ll use for numerical simulation and estimates. These estimates are often wide due to wide variations from one infection to another. Only the parameter *β* will be estimated. View this table: [Table1](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/T1) Remark 2.3. *In the above model, we only consider constant-valued parameters. However, many of these quantities actually vary with time. Take for example the transmission rate coeffcient β. One can then thing naturally about a possible time varying transmission coefficient β*(*t*) *that allows the effect of control interventions and by the way we will explore how the basic reproduction number R* *change*. According to this remark, the previous system of ordinary differential equations becomes: ![Formula][10] This is a realistic model that can be used by epidemiologists to study the dynamics of the diseases and the effects of control interventions. The basic reproduction number *R* is defined as the average number of secondary cases generated by a primary case over his/her infectious period when introduced into a large population of susceptibles individuals [19]. Note that the constant *R* thus measures the initial growth rate of the epidemic and for the model above it can be shown (see the next section) that ![Formula][11] where *γ*1 = *pγ**s* and *γ*2 = (1 *p*)*γ**h*. Furthermore Chowell et al. [2] define the time-dependent effective reproductive number ![Formula][12] of secondary cases per infections cases at time *t*. One can see that if *S*(*t*) ∼ *N* it follows that ![Graphic][13] is a function proportional to the time varying transmission coefficient *β*(*t*). The time point at which *R*(*t*) assumes values smaller that 1 indicates when control measures have becomes effective in controlling the epidemic. The intervention strategies to control the spread of Covid include surveillance, placement of suspected cases in quarantine for 2 weeks (the maximum estimated length of the incubation period), education of hospital personnel and community members on the use of strict barrier nursing techniques (i.e. protective clothing and equipment, patient management), and the rapid burial or cremation of patients who die from the disease. In order to account for the control intervention we assume that the transmission rate function *β*(*t*) is constant up to the time point *t*∗ when the control measures are introduced and after that decays exponentially. Assumes that a fraction *c* ∈ [0, 1) of new infections is reduced during a period *T* = *t* −*t*∗ where *t*∗ is the time point at which control measures are introduced. Then we can write: ![Formula][14] ![Formula][15] ![Formula][16] where *q* = − ln(1 − *c*) ≥ 0 is a kind of quarantine control paramter. Thus in the rest of this paper we will assume that there exists *q* ∈ [0, 1) such that ![Formula][17] where the date *t*∗ denote the starting time where the control intervention start and *n* equal to the period where the control parameter relied for example one can assume that the percentage of a control measure is maintened about a 7-day period. This will be the approach for simulations in the paper. Note that such kind of parametrization of *β*(*t*) was already consider in [**?**] and [**?**]. Let us point out that the geometry of the likelihood function can not permit identification of this quarantine parameter *q* as its estimator is correlated. Our way for the transmission rate is a more parsimonious parametrization and it does not affect *γ* the recovery rate unless the disease is curable, which is not the case for Covid-19 actually. Now, let us discussed about the importance of *c* and *p*. * If *c* = 0 (*q* = 0) one can see that *β*(*t*) = *β* which means that 0% of new infections are reduced so that there is no effort to contain new infections which mean that the control intervention failed. * Note that in the time interval [0, *t* −*t*∗] if we consider *q* = 0.2 (*c* ∼ 0.18) and *β* = 0.422 this is equivalent to a higher level of control measures since this mean that after time *T* the control measure started, approximately 80% of new interventions are reduced. One can then thing on 1− *p* as the percentage of succefull intervention measures during the major public policies process introduced to avoid new infections. * If a fraction *c* of new infections is reduced during every period *T* = *t* − *t*∗ then the epidemic will go to extinction after *nT* time where *n* is the interger part of ![Formula][18] To close this paragraph, let us say some words above a possibility to have an almost periodicity or asymptotic periodicity in the data. Some disease have seasonal patterns such that measles and flu, which may be well represented by an almost periodic or asymptotic periodictime varying transmission rate resulting from, for example, changes in contact rates corresponding to when control measures or reported numbers of infectious fluctuate on some day the day of the week. In the spirit of the above parametrization, we will consider according to observation in real data where we think that there is a slight increasing in a sequence of change times *t*1, …*t**K* and followed by a slowly varying exponential decrease, a piecewise linear continuous transimission coefficient *β*(*t*) in the form: ![Formula][19] where *β* = *β*(0) at the begining and ![Formula][20] and *n* = 1, 2 denoting the change of one week to another. ## 3. About the french overseas department Mayotte and epidemic Data Mayotte, is an island and overseas department of France located in the northern Mozambique Channel in the Indian Ocean off the coast of Southeast Africa, between northwestern Madagascar and northeastern Mozambique. It is also located between the two southeasternmost islands of the Comoros archipelago. Mayotte is made up of one main island (also called Grande Terre), one smaller island called Petite Terre which lies about 1.5 miles east of the main island and islets. The capital, Mamoudzou, is located on the eastern coast of the main island and in the neighbouring smaller island of Petite-Terre, the DzaoudziPamandzi International Airport is located. The territory is also known as Maore, the native name of its main island, especially by advocates of its inclusion in the Union of the Comoros. Mayotte has a total land area of about 144 square miles or 374 square kilometres (sq. km), and with its almost 285000 people according to official estimates, is very densely populated at 747 inhabitants per sq.km. Mayotte is the poorest department of France. However, many immigrants from other Mozambique Channel nations enter Mayotte illegally, as it is one of the most prosperous regions in the area. From 1976 into the 21st century, Mayotte had a special status with France as a territorial collectivity, conceived as being midway between an overseas territory and an overseas department. Its status was changed to a departmental collectivity in 2001 and then to overseas department in 2011. Mayotte is represented in the French National Assembly by a deputy and in the French Senate by two senators. It is administered by a French-appointed prefect and an elected Departmental Council. The judiciary is modeled on the French system. Although, as a department, Mayotte is now an integral part of France, and despite being in France, the vast majority of Mayottes population do not speak French as a first language, but a majority of the people 14 years and older report in the census that they can speak French (with varying levels of fluency). The majority language is Shimaore, a language related closely to those found in the neighbouring Comoros islands. Kibushi, a Malagasy language, is the second most widely spoken. Both have been influenced by Shimaore. The vast majority of the population is Muslim. As concerned the land, note that a volcanic mountain range forms a north-south chain on Mayotte island, with summits from about 1, 600 to 2, 000 feet (500 to 600 metres) in elevation. Protected waters for shipping and fishing are created by surrounding coral reefs some distance from the shore. Mayotte age structure is as follows: 46.6% of the population is under age 15, 51% of the population between 15 and 64 years old and 2.4% of the population is older than 65. As we can see the Mayotte population is young. This type of age structure is common for developing countries with high birth and death rates. Relatively short life expectancy, as well as low level of education and poor health care. ### 3.1 Epidemic Data process in Mayotte Mayottes Regional Health Agency and local authorities confirms the first case of new coronavirus infection, an imported case, on march 14 and reported on march 13. He was an individual who traveled from the metropolitan France and back to Mayotte three days ago. Sub-sequently, five other imported case was confirmed on march 20. After that 8 others imported on March 26. Thus a total of 14 imported cases were confirmed. Three others imported cases were confirmed from march 27 to 29. One can thus think that almost 1 imported case is introduced per day. In this case and according to Remark 8, we can assume that among susceptibles individuals a negligible fraction of international infectious individuals subsequently travel and are detected in their destination to Mayotte. Thus, according to the model (8) one may take *y*(*t*) = *t* a variable number *y*(*t*) of exposed international travellers to Mayotte at time *t*. Clearly, this is reasonable according to the epidemic history data in Mayotte at the begining. But since March 30, the authorities of Mayotte department and France have decided to close and cancel all flights to and from Mayotte; which therefore limits the introduction of new imported cases but cannot prevent the spread of the virus. Hence for *t* ∈ [*t*, *t**c*], this leads to the new equations of the model: ![Formula][21] where *t* = *march*13 and *t**c* = *march*29. But for *t* ≥ *t**c* +1 we assume that *y*(*t*) = 0. We obtain information and data from the regional french agency of health called ARS of Mayotte where the second author of this paper belongs. View this table: [Table2](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/T2) According to ARS Mayotte and Santé Publique France in Mayotte, the probable contact number with a confirmed case is around 15 − 16 so we take for our model *E*(0) = 15 at the begining of the epidemic. More informations can be found at the following website of the Mayotte local authorities: [http://www.mayotte.gouv.fr](http://www.mayotte.gouv.fr) and [www.mayotte.ars.sante.fr](http://www.mayotte.ars.sante.fr) in Coronavirus COVID-19: Point de situation View this table: [Table3](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/T3) So for simulation purposes we will have that *I*(0) = 3 from march 13 to 20, *I*(0) = 14 from march 20 to 26 and *I*(0) = 17 from march 27 to nowadays. The following table presents the data observed in Mayotte. Data information includes the cumulative number of reported (confirmed), imported, hospitalized cases including critical cases. For more informations about the collection of the data, we refer the reader to the following website of the Mayotte local authorities: [http://www.mayotte.gouv.fr](http://www.mayotte.gouv.fr) or [www.mayotte.ars.sante.fr](http://www.mayotte.ars.sante.fr) in Coronavirus COVID-19 Point de situation. View this table: [Table 1.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/T4) Table 1. Covid cases in Mayotte from march 13 to 19 ## 4. Illustration and Numerical Simulations of the epidemic COVID-19 in Mayotte The following picture illustrate both the number of cumulative confirmed, cured, hospitalized and death cases from march 13 to april 14. ### 4.1 Estimation of the basic reproduction number *R* in the begining of the epidemic The Reproduction number *R* may be estimated at different times during an epidemic. Some importants review found that the estimated mean *R* for COVID-19 is in [2, 3] with a median of2.5 which is considerably higher than the SARS estimate or the influenza. These estimates of *R* depend on the estimation method used as well as the validity of the underlying assumptions. Due to insufficient data and short onset time, current estimates of *R* for *COV ID* − 19 are possibly biased. However, as more data are accumulated, estimation error can be expected to decrease. View this table: [Table 2.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/T5) Table 2. Covid cases in Mayotte from march 20 to 26 View this table: [Table 3.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/T6) Table 3. Covid cases in Mayotte from march 27 to april 3 View this table: [Table 4.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/T7) Table 4. Covid cases in Mayotte from april 4 to 10 Methods for estimating the basic reproduction number at the beginning of the out-break and the time-dependent reproduction number of secondary cases at any time during the outbreak, are usually bases on the following methods, but all required some hypotheses in order to apply the methods. View this table: [Table 5.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/T8) Table 5. Covid cases in Mayotte from april 11 to 17 For instance, the Maximum Likewood (ML) that we are interested in this paper as an estimator according to [9] and [27] propose extensions and options implemented in the software *R*. Others methods like Exponential growth rate (EG) and Sequential Bayesian method also exists. Note that the ML methods require the user to select the time period over which growth is exponential and to determine the serial intervall meaning the duration between symptoms onset of the primary cases and symptoms onset of the secondary cases in a transmission chain. In the case of *Covid*− 19 in Mayotte, the epidemic curve can be analyzed the first week from the case on march 15 to march 21 where over consecutive time units and a generation time interval with mean 2 days and standard deviation 0.5 day. To estimate the early dynamics of *Covid* − 19 transmission in Mayotte, we fitted the value of the serial interval with a gamma distribution with mean 2 and standard deviation 0.5. By using the *R* Package, we estimate that the basic reproduction number for the epidemic in Mayotte ![Formula][22] We presented the code using the *R* Package in the following lines. Now, we are able to apply the preceding compartmental model for the prediction of the *Covid* − 19 in Mayotte. Note that, according to our assumptions we have ![Formula][23] Here is a Typical session code. >**library(R0)** # **Loading package** ># **epidemic curve can be input as a list of dates** >**epid** = **c(“2020-03-15”, “2020-03-16”, “2020-03-17”**,**”2020-03-18”, “2020**- **03-19”**,**”2020-03-20”**,**”2020-03-21”)** ># **create generation time: gamma distribution, with mean** 2 **time units and standard deviation** 0.5 **time units** >**mGT** = **generation.time (“gamma”, c(2, 0.5))** ># **outbreak during the 2020 Coronavirus in Mayotte** >**Mayotte.2020**=c(0,1,1,1,2,3,2) ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/17/2020.04.15.20062752/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/F4) Figure 4. Covid19 in Mayotte: cumulative confirmed, cured, hospitalized and death cases from march 13 to april 14. >**data(Mayotte.2020)** >**est.R0.ML(Mayotte.2020, mGT, begin=1, end=7, range=c(0**,**3))** ># **Reproduction number estimate using Maximum Likelihood method**. >**plot(res.R)** # **diplays results** >**res=est.R0.ML(Mayotte.2019, mGT, begin=1, end=7, range=c(0**,**3))** >**plot(res)** ### 4.2 Prediction of the Covid-19 in Mayotte Assume that there was any intervention control, of course this is not true beacause from march 16, the french government imposed public policies designed to contain this epidemic in all departments including overseas department, for instance in Mayotte. The following table illustrative cumulative confirmed compared to predicted data with our mathematical model. View this table: [Table 6.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/T9) Table 6. Reported Covid and predicted cases in Mayotte from march 13 to april 10 Now Remark 4.1. *Note that during in this period the model predict that at most two hospitalized individual will died. This was the case at april 10 in Mayotte*. ## Discussion Comparing model predictions with observed confirmed cases reported in Mayotte, we found that, from april 3 to 10 the model predicted almost two times higher cases than those were reported. Thus, we deduce that the control measures and public policies really started around April 1 instead of March 16 the official date. ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/17/2020.04.15.20062752/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/F5) Figure 5. Predictions cases of Covid19 in Mayotte from march 26 to april 10 without any public policies. The following picture illustrate the number of cumulative predictions cases in Mayotte as weel as the pic and the end of the epidemic without any control measure. ## Discussion Without any control measure, we can see that the pic will be attained around June 25 an the extinction of the epidemic can be esperate at the end of august. Remark 4.2. *Note also that the cancelation of travel restrictions on the peak time and peak value has an effect on the peak time and peak value*. ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/17/2020.04.15.20062752/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/F6) Figure 6. Predictions cases of Covid19 in Mayotte without any public control. ## 5. Prediction when varying the quarantine control parameter*q* Now let us assume that from march 20, the transmission rate parameter is time-depending, does not affect the recovery rate unless the disease is curable, which is not the case for Covid-19 actually and exponentially decreasing as described in 7 with a control quarantine paramter *q*. This should be actually the case since the control measure started at this time. ![Formula][24] with *n* = 1 means march 20 and *n* = 2 means march 27 and so one. The following table show the evolution of the epidemic when the control paramter *q* = 0.2 under the following time varying function ![Figure 7.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/17/2020.04.15.20062752/F7.medium.gif) [Figure 7.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/F7) Figure 7. Predictions cases of Covid19 in Mayotte with a control paramter *q* = 0.2. ## Discussion Under control measure with a control paramter *q* = 0.2, one observed that the pic is attained later around July 17 with a maximum of infectious people less that the one with no control measure. ![Figure 8.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/17/2020.04.15.20062752/F8.medium.gif) [Figure 8.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/F8) Figure 8. Predictions cases of Covid19 in Mayotte with various control parameter ## Discussion Now, let us remark that at this stage with the exponential decreasing parameter *β* if *q* = 0.5 every 7 days one can observe a rapide extinction of the epidemic but with *q* = 0.05 or *q* = 0.07 every 7 days the extinction come later. At this stage, we hink that a different parametrization need to be introduced. So one would like to imagine a slowly varying parameter in the form ![Formula][25] But since we already observed that the control measures and public policies really started around April 1 instead of March 16 the official date, we will assume that there was an increasing of the transmission parameter *β*(*t*) during the second week (from March 20 to 26 and from March 26 to April 3) and after that the decrease phase started Now we assume in the rest of the prediction scenarios that there was a slowly varying in transmission on the form ![Formula][26] where *β* = *β*(0) at the begining and ![Formula][27] and *n* = 1, 2 …. denoting the change of one week to another. Let us see what becomes the simulation with a specific choosen quantine parameter *q* = 0.1 compared with reported Covid-19 cases in Mayotte. ![Figure 9.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/17/2020.04.15.20062752/F9.medium.gif) [Figure 9.](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/F9) Figure 9. Predictions cases of Covid19 in Mayotte with the new control parameter To examine the possible impact of enhanced interventions on disease infections, we plotted the number of cumulative infected individuals *I*(*t*) and the predicted cumulative number of reported cases with varying the quarantine parameter *q*. This analysis shows that incresing the parameter persistently decreases the peak value but may either delay or bring forward the peak, as shown in the above pictures. Remark 5.1. *Note that, our analysis shows that increasing quarantine rate will bring forward the peak and lead to a reduction of the peak value in terms of the number of infected individuals. This indicates that enhancing quarantine and isolation following contact tracing and reducing the contact rate can signicantly lower the peak and reduce the cumulative number of predicted reported* Remark 5.2. *Considering the spreading of the virus and in order to examine the impact the impact of the travel restriction We find that with travel restriction* (*no imported exposed individuals*), *the number of infected individuals compared with the scenario of no travel restriction will decrease. This means that the effect of a travel restriction to Mayotte on the* 2019 −*nCoV infection is almost equivalent to increasing quarantine parameter*. ### 5.1 Discussion about the relevance about stochastic models At this stage, let us point out the fact that a stochastic pertubation is necessary since in general, susceptible individual does not have the same probability of being infected and depends on many factors depending on the behavior of individuals which changes over the course of the epidemic for instance the confinement will changes the habit. Therefore suppose that some stochastic environmental factor acts simultaneously on each individual in the population so that *β* changes to a a stochastic process (*β*(*t*))*t*≥0 such that ![Formula][28] where *B* is a standard brownian motion. Hence the number of potential infectious contacts that a single infected individual makes with another individual during a period *dt* is normally distributed with mean *βdt* and variance *σ*2*dt*. Hence ![Formula][29] As *V ar*(*β*(*t*)) → 0 when *dt* → 0, this is reasonable model for biologist and it is a well-established way of introducing stochastic environmental noise into biologically realistic population dynamic models. To motivate the above assumption on *β*, let us we argue as follows. Suppose that thenumbers of potentially infectious contacts between an infectious individual and another individual in successive time intervals [*t, t*+*T*),[*t*+*T, t*+2*T*), …, [*t*+(*n*1)*T, t*+*nT*) are independent identically distributed random variables and that *n* is very large. By the Central limit Theorem, the total number of potentially infectious contacts made in [*t, t* + *nT*) has approximately a normal distribution with mean *nx* and ![Graphic][30], where *x* and ![Graphic][31] are, respectively, the mean and variance of the underlying distribution in each of the separate time intervals of length *T*. Hence it is reasonable to assume that the total number of potentially infectious contacts has a normal distribution whose mean and variance scale as the total length of the time interval as in our assumptions. Therefore we replace *β*(*t*) in (7) by ![Formula][32] where *B* is a standard brownian motion and *σ* is the diffusion paramter that allows the effect of random control measures. Imitating the arguments of proof in [], we can show that under the above assumption on *β*, that the underline SDE rewrittes as follows admits a positive solution. ![Formula][33] ## 6. Towards a generalized continuous time Markov SEIR model The analogous stochastic model (continuous timeMarkov chain) is constructed by considering the following events: exposure *E*, infection *I* and removal (mild or severe)*R*. The transition rates are defined as View this table: [Table10](http://medrxiv.org/content/early/2020/04/17/2020.04.15.20062752/T10) When an infection occurs, the number of suscptibles decreases by 1 and the number of exposed increases by 1; when the latency period *γ* = (*pγ**s* + (1− *p*)*γ**h*) ends, the number of number of exposed decreases by 1 and the number of infectives increases by 1; and finally when there is a recovery, the number of infectives decreases by 1, and finally when there is a recovery, the number of infectives decreases by 1 and the number of recovered increases by 1. increments distributed exponentially: Now, to understand the following probabilistic modelisation, assumes first that the latency period *λ* and the duration of infection *γ* satisfies an exponential distribution *E*∼ *E* (*λ*), *I* ∼ (*γ*) whe the mean retired rate *γ* = *pγ**s* + (1 *p*)*γ**h*. In that case, the above three types of events happen as follows: * Infection of a susceptible (such an event decreases *S*(*t*) by one, and increases *I*(*t*) by one, so the event happen at rate: ![Formula][34] * Exposed (contact) of a susceptible with an infective (such an event increases *E*(*t*) by one, and decreases *S*(*t*) by one, so the event happen at rate: ![Formula][35] * Removed of an infective with mild infection *R**M* (such an event decreases *I*(*t*) by one and increases *R**M* (*t*) by 1 so the event) happen at rate ![Formula][36] * Retired *R**S* with severe infection (such an event decreases *I*(*t*) by one and increases *R**S*(*t*) by 1 so the event) happen at rate ![Formula][37] Hence we have the following equations with *P**s,e*, *P**e,i*, *P**i,rs*, *P**i,rh* *P**rh,d* standard mutually independent poisson processes: ![Formula][38] Note that one can clearly forget about the last equation since If we set ![Formula][39] Now define ![Formula][40] We have ![Formula][41] The above model equations assumes that *β, γ**s*, *γ**h* and *λ* are constant. But as explained above a temporally varying parameter is relevant in order to explained stochastic or deterministic fluctuations in *β* so we consider *β*(*t*) and the model equations can be read as follows. ![Formula][42] Note that here again the population is closed: ![Formula][43] and the above stochastic model can reads as follows: ![Formula][44] Where ![Formula][45] and ![Formula][46] The *P**j* are standard mutually independent poisson processes defined above. Note that we forgot about the last equation in *D*(*t*) since it is not relevant in the study of the previous model equations since it can be deduced easily. In order to estabilished some relevant results let us recall the following Large Numbers for Poisson processes. Proposition 6.1. *Let* (*P* (*t*))*t*≥0 *be Poisson process with rate λ. Then* ![Formula][47] Imitating the arguments in Theorem 2.2.7 in [22], we can established the following Theorem. Theorem 6.1. *Law of Large Numbers Assume that* ![Formula][48] *is a locally Lipschitz function in x, locally uniformly in time t such* ![Formula][49] *that and the unique solution of the ODE* ![Formula][50] *does not explose in finite time. Let* ![Graphic][51] *denote the solution of the above Stochastic Differential Equation. Then* ![Graphic][52] *converges almost surely and locally uniformly in t to* ![Formula][53] *which is exactly our deterministic model normalized with N with s* = *S*/*N, i* = *I*/*N, e* = *E*/*N r**m* = *R**M* /*N, r**s* = *R**S*/*N and z* = (*s, e, i, r**m*, *r**h*). Now, since the above stochastic model converges to our deterministic diffential system model, it is natural to look at fluctuations of the difference between the stochastic epidemic process and its deterministic limit. In this way, we introduce the following rescaled difference between ![Graphic][54] and *z*(*t*): ![Formula][55] Define the continuous time martingales ![Formula][56] We have ![Formula][57] Now, consider the following 5-dimensional process whose *j*-th component is defined as ![Formula][58] We wish to show that ![Graphic][59] converges in law to a Gaussian process. It is clear that ![Formula][60] Now let us established the main result of this section. The proof is contained in Theorem 2.3.2 in [22]. Theorem 6.2. *Central Limit Theorem Under the assumption of the previous theorem, assume furthermore that the function b is locally of class C*1 *in x, locally uniformly in t. Then as N*→ + ∞ *the above rescaled process* ![Graphic][61] *converges in law to a Gaussian process V defined by the following stochastic differential equation*. ![Formula][62] *where* (*B*1, *B*2, *B*3, *B*4, *B*5) *are mutually independent Brownian motion*. ## 7. Methods for computing*R*and*R*(*t*) The basic reproduction number *R* is defined as the average number of secondary cases generated by a primary case over his/her infectious period when introduced into a large population of susceptibles individuals [19]. The constant *R* thus measures the initial growth rate of the epidemic. A typical first step in analysing a system of differential equations is finding the equilibra. At equilibrium we do not mean that individuals are fixed within particular compartments, but rather the rate of individuals entering and leaving a compartment is excatly balanced. This is equivalent to setting the right hand size of our equations, i.e. the rate of change of individuals into and out of each compartment equal to zero. In epidemiology, models generally have two important equilibria: * the disease-free equilibrium (DFE) and the endemic equilibrium. The DFE requires that there are no infected individuals in the population, or in the case of SIR model *I*∗ = 0, when the star designates an equilibrium solution. In contrast the endemic equilibrium corresponds to the state in which infected individuals persist indefinitely such that *I*∗ > 0. Finding equilibria is only the first step to understanding long-term behavior in a system. We must also determine which of the behavior are typically realized. This involves determining the stability of the equilibria and will show us whether we will approach the equilibria or more away from it, assuming we have started near by. Now coming back to the method for computing *R*, there is a rich mathematical theory that describe how this quantity can be computed for a long range of SEIR-type models. Let us introduce a commonly used method for finding *R*: the next generation matrix method. of [**?**]. A nice feature of the next generation matrix method is that it only requires use of the DFE which is often easy to compute. In the SEIR model the DFE is given by ![Formula][63] Set ![Formula][64] After determining the DFE, we must create a sub-model that only considers the “disease” compartments, a subset of the equations in the modified SEIR model. The disease compartment are those that include individuals that are in any stage of infection which, for the SEIR model, includes both the exposed and infections individuals namely equations: ![Formula][65] Set ![Formula][66] and we can write the disease compartment in the form ![Formula][67] ![Formula][68] Here ![Graphic][69] contains any terms that directly lead to new infections entering each compartment *j* (here *j* = 2 because disease compartments are *E* and *I*) and ![Graphic][70] a vector of the *j* disease compartments. Notice that the second element of ℱ is zero because no new infections enter the compartment *I* rather their transition from to the *E* into *I*. ![Formula][71] where ![Graphic][72] contains all of the outputs and ![Graphic][73] contains all the inputs from each disease class. This includes term such as mortality or transition between class. ![Formula][74] The Jacobian matrix of our equivalent sub-model equations evaluated at the DFE is given by ![Formula][75] As a result, we can factor out the vector ![Graphic][76] on the right hand size, leaving us with ![Formula][77] so that ![Formula][78] From linear stability analysis we know that the DFE is equivalent to the real part of all eigen values of *F* − *V* being less than zero. ![Formula][79] means that the DFE is stable where *ρ* is the spectral radius of *FV* −1. ![Formula][80] so that ![Formula][81] ## Data Availability https://www.mayotte.ars.sante.fr/coronavirus-covid-19-point-de-situation-et-conduite-tenir [https://www.mayotte.ars.sante.fr/coronavirus-covid-19-point-de-situation-et-conduite-tenir](https://www.mayotte.ars.sante.fr/coronavirus-covid-19-point-de-situation-et-conduite-tenir) ## Conflict of Interests The authors are gratefull to the staff of ARS-Mayotte for usefull informations and to Benoîte De Saporta and Etienne Pardoux for many relevant discussions during the prepatory of this paper.The authors declare that there is no conflict of interest regarding the publication of this paper. * Received April 15, 2020. * Revision received April 15, 2020. * Accepted April 17, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. [1].Britton and E. Pardoux, Stochastic Epidemic Models with Inference. Mathematical Biosciences Subseries, Springer, Doi: 10.1007/978-3-030-30900-8; (2020) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-3-030-30900-8&link_type=DOI) 2. [2]. G. Chowell et al. The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and Uganda. Journal of Theoretical Biology; Volume 229 (1), 119–126 (2004) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2004.03.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15178190&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F17%2F2020.04.15.20062752.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222192000011&link_type=ISI) 3. [3].Obadia et al. The R0 package: a toolbox to estimate reproduction numbers for epidemic outbreaks. BMC Medical Informatics and Decision Making 2012, 12:147. 4. [4]. J. Balicchi and Anne Barbail Structural and Predictive Analysis of the Birth Curve in Mayotte from 2011 to 2017. In Wiley Book, February 2020: [https://doi.org/10.1002/9781119706922.ch112020](https://doi.org/10.1002/9781119706922.ch112020) 5. [5].Komatsu Estimation of the reproductive number of novel coronavirus Covid − 19 and the probable outbreak size on the Diamond Princess cruise ship: A data-driven analysis International Journal of Infectious Diseases (2020) 6. [6]. Joseph T Wu et al. Nowcasting and forecasting the potential domestic and international spread of the 2019 −nCoV outbreak originating in Wuhan, China: a modelling study. Lancet Infect Dis 2020; 395: 68997 7. [7]. A. J. Kucharskir et al. Early dynamics of transmission and control of Covid − 19: a mathematical modelling study. Lancet Infect Dis, [https://doi.org/10.1016/](https://doi.org/10.1016/) S1473-3099(20)30144-4; (2020) 8. [8]. Ying Liu et al. The reproductive number of COVID-19 is higher compared to SARS coronavirus. Journal of Travel Medicine, (2020), 14 doi: 10.1093/jtm/taaa021. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jtm/taaa021&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32052846&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F17%2F2020.04.15.20062752.atom) 9. [9].Li et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus Infected Pneumonia N Engl J Med (2020); 382:1199–1207 DOI: 10.1056/NEJMoa2001316 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001316&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31995857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F17%2F2020.04.15.20062752.atom) 10. [10]. Toshikazu Kuniya. Prediction of the Epidemic Peak of Coronavirus Disease in Japan, 2020.J. Clin. Med. (2020), 9, 789; doi:10.3390/jcm9030789 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/jcm9030789&link_type=DOI) 11. [11].Bulletin d’informations Covid-19 Agence Rgionale de Sant de Mayotte (2020). 12. [12].COV ID − 19: point pidmiologique du 10 mars 2020 Sant Publique France (2020). 13. [13]. Seth Flaxman et al. Estimating the number of infections and the impact of nonpharmaceutical interventions on COVID-19 in 11 European countries Imperial College COVID-19 Response Team (2020). 14. [14].Covid−19: Point pidmiologique - Situation au 15 mars 2020. Sant Publique France (2020). 15. [15]. P. Lekone et al. Statistical Inference in a Stochastic Epidemic SEIR Model with Control Intervention: Ebola as a Case Study Biometrics 62, 11701177 (2005). 16. [16]. Biao Tang et al. Estimation of the Transmission Risk of the 2019−nCoV and Its Implication for Public Health Interventions. J. Clin. Med. (2020), 9, 462; doi:10.3390/jcm9020462 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/jcm9020462&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32046137&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F17%2F2020.04.15.20062752.atom) 17. [17]. Biao Tang et al. An updated estimation of the risk of transmission of the novel coronavirus (2019-nCov) Infectious Disease Modelling 5 (2020) 248e255 18. [18]. Jonathan M. Read et al. Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions. medRxiv preprint doi: [https://doi.org/10.1101/2020.01.23.20018549](https://doi.org/10.1101/2020.01.23.20018549), (2020). 19. [19]. O. Dieckmann. On the definition and the computation of the basic reproduction ratio R0in models for infectious diseases in heterogeneous populations. J. Math. Biol. (1990) 28: 365-382, (1990). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF00178324&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2117040&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F17%2F2020.04.15.20062752.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DG35000001&link_type=ISI) 20. [20]. P. Magal and G. Webb. Predicting the number of reported and unreported cases for the COVID-19 epidemic in South Korea, Italy, France and Germany Prepublication (2020). 21. [21]. Z. Liu, P. Magal, O. Seydi and G. Webb. Understanding unreported cases in the 2019-nCov epidemic outbreak in Wuhan, China, and the importance of major public health interventions. MPDI Biology, 2020, 9(3), 50. 22. [22].Britton and E. Pardoux, Stochastic Epidemic Models with Inference. Mathematical Biosciences Subseries, Springer, Doi: 10.1007/978-3-030-30900-8; (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-3-030-30900-8&link_type=DOI) 23. [23]. Z. Liu, P. Magal, O. Seydi, and G. Webb, Predicting the cumulative number of cases for theCOV ID − 19 epidemic in China from early data medRxiv, 2020. 24. [24]. Solym Manou-Abi et al. Asymptotically periodic solution for a stochastic fractional integro-differential equation In Mathematical Modeling of Random and Deterministic Phenomena [https://doi.org/10.1002/9781119706922.ch6](https://doi.org/10.1002/9781119706922.ch6), Wiley ISTE London (2020) 25. [25]. Solym Manou-Abi et al. Asymptotically periodic solution of a stochastic differential equation Bulletin of the Malaysian Mathematical Sciences Society (2018) 26. [26]. Solym Manou-Abi et al. Asymptotically periodic functions in the stepanov sense and its application for an advanced differential equation with piecewise constant argument in a Banach space Mediterr. J. Math. (2018) 15: 25 27. [27].Groupe de modlisation de l’quipe ETE Laboratoire MIVEGEC, CNRS, IRD, Universit de Montpellier Rapport 1, (2020) 28. [28].Zhou et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study Lancet 2020; DOI:[https://doi.org/10.1016/S0140-6736(20)30566-3](https://doi.org/10.1016/S0140-6736(20)30566-3) [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif [4]: /embed/graphic-2.gif [5]: /embed/graphic-3.gif [6]: /embed/graphic-4.gif [7]: /embed/graphic-5.gif [8]: /embed/graphic-8.gif [9]: /embed/graphic-9.gif [10]: /embed/graphic-11.gif [11]: /embed/graphic-12.gif [12]: /embed/graphic-13.gif [13]: /embed/inline-graphic-4.gif [14]: /embed/graphic-14.gif [15]: /embed/graphic-15.gif [16]: /embed/graphic-16.gif [17]: /embed/graphic-17.gif [18]: /embed/graphic-18.gif [19]: /embed/graphic-19.gif [20]: /embed/graphic-20.gif [21]: /embed/graphic-21.gif [22]: /embed/graphic-29.gif [23]: /embed/graphic-30.gif [24]: /embed/graphic-35.gif [25]: /embed/graphic-38.gif [26]: /embed/graphic-39.gif [27]: /embed/graphic-40.gif [28]: /embed/graphic-42.gif [29]: /embed/graphic-43.gif [30]: /embed/inline-graphic-5.gif [31]: /embed/inline-graphic-6.gif [32]: /embed/graphic-44.gif [33]: /embed/graphic-45.gif [34]: /embed/graphic-47.gif [35]: /embed/graphic-48.gif [36]: /embed/graphic-49.gif [37]: /embed/graphic-50.gif [38]: /embed/graphic-51.gif [39]: /embed/graphic-52.gif [40]: /embed/graphic-53.gif [41]: /embed/graphic-54.gif [42]: /embed/graphic-55.gif [43]: /embed/graphic-56.gif [44]: /embed/graphic-57.gif [45]: /embed/graphic-58.gif [46]: /embed/graphic-59.gif [47]: /embed/graphic-60.gif [48]: /embed/graphic-61.gif [49]: /embed/graphic-62.gif [50]: /embed/graphic-63.gif [51]: /embed/inline-graphic-7.gif [52]: /embed/inline-graphic-8.gif [53]: /embed/graphic-64.gif [54]: /embed/inline-graphic-9.gif [55]: /embed/graphic-65.gif [56]: /embed/graphic-66.gif [57]: /embed/graphic-67.gif [58]: /embed/graphic-68.gif [59]: /embed/inline-graphic-10.gif [60]: /embed/graphic-69.gif [61]: /embed/inline-graphic-11.gif [62]: /embed/graphic-70.gif [63]: /embed/graphic-71.gif [64]: /embed/graphic-72.gif [65]: /embed/graphic-73.gif [66]: /embed/graphic-74.gif [67]: /embed/graphic-75.gif [68]: /embed/graphic-76.gif [69]: /embed/inline-graphic-12.gif [70]: /embed/inline-graphic-13.gif [71]: /embed/graphic-77.gif [72]: /embed/inline-graphic-14.gif [73]: /embed/inline-graphic-15.gif [74]: /embed/graphic-78.gif [75]: /embed/graphic-79.gif [76]: /embed/inline-graphic-16.gif [77]: /embed/graphic-80.gif [78]: /embed/graphic-81.gif [79]: /embed/graphic-82.gif [80]: /embed/graphic-83.gif [81]: /embed/graphic-84.gif