Research article

Fat tails arise endogenously from supply/demand, with or without jump processes

  • We show that the quotient of Levy processes of jump-diffusion type has a fat-tailed distribution. An application is to price theory in economics, with the result that fat tails arise endogenously from modeling of price change based on an excess demand analysis resulting in a quotient of arbitrarily correlated demand and supply whether or not jump discontinuities are present. The assumption is that supply and demand are described by drift terms, Brownian (i.e., Gaussian) and compound Poisson jump processes. If P1dP/dt (the relative price change in an interval dt) is given by a suitable function of relative excess demand, (DS)/S (where D and S are demand and supply), then the distribution has tail behavior F(x)xζ for a power ζ that depends on the function G in P1dP/dt=G(D/S). For G(x)|x|1/q one has ζ=q. The empirical data for assets typically yields a value, ζ˜=3, or ζ[3,5] for many financial markets.

    Many theoretical explanations have been offered for the disparity between the tail behavior of the standard asset price equation and empirical data. This issue never arises if one models price dynamics using basic economics methodology, i.e., generalized Walrasian adjustment, rather than the usual starting point for classical finance which assumes a normal distribution of price changes. The function G is deterministic, and can be calibrated with a smaller data set. The results establish a simple link between the decay exponent of the density function and the price adjustment function, a feature that can improve methodology for risk assessment.

    The mathematical results can be applied to other problems involving the relative difference or quotient of Levy processes of jump-diffusion type.

    Citation: Gunduz Caginalp. Fat tails arise endogenously from supply/demand, with or without jump processes[J]. AIMS Mathematics, 2021, 6(5): 4811-4846. doi: 10.3934/math.2021283

    Related Papers:

    [1] José Rodrigo da-Silva, Ana Alexandre, Clarisse Brígido, Solange Oliveira . Can stress response genes be used to improve the symbiotic performance of rhizobia?. AIMS Microbiology, 2017, 3(3): 365-382. doi: 10.3934/microbiol.2017.3.365
    [2] Yen Thi Ngoc Phan, Minh Thuy Tang, Tu Thi Minh Tran, Van Huu Nguyen, Trang Hien Nguyen, Takeshi Tsuruta, Naoki Nishino . Diversity of lactic acid bacteria in vegetable-based and meat-based fermented foods produced in the central region of Vietnam. AIMS Microbiology, 2017, 3(1): 61-70. doi: 10.3934/microbiol.2017.1.61
    [3] Xavier Cruz-González, Nereha Laza-Pérez, Pedro F. Mateos, Raúl Rivas . Analysis and effect of the use of biofertilizers on Trifolium rubens L., a preferential attention species in Castile and Leon, Spain, with the aim of increasing the plants conservation status. AIMS Microbiology, 2017, 3(4): 733-746. doi: 10.3934/microbiol.2017.4.733
    [4] Ángeles Hidalgo, Francisco-Javier López-Baena, José-Enrique Ruiz-Sainz, José-María Vinardell . Studies of rhizobial competitiveness for nodulation in soybean using a non-destructive split-root system. AIMS Microbiology, 2017, 3(2): 323-334. doi: 10.3934/microbiol.2017.2.323
    [5] Jasmine S. Ritschard, Lea Amato, Yadhu Kumar, Britta Müller, Leo Meile, Markus Schuppler . The role of the surface smear microbiome in the development of defective smear on surface-ripened red-smear cheese. AIMS Microbiology, 2018, 4(4): 622-641. doi: 10.3934/microbiol.2018.4.622
    [6] Caroline L. Marden, Ryan McDonald, Harold J. Schreier, Joy E.M. Watts . Investigation into the fungal diversity within different regions of the gastrointestinal tract of Panaque nigrolineatus, a wood-eating fish. AIMS Microbiology, 2017, 3(4): 749-761. doi: 10.3934/microbiol.2017.4.749
    [7] Neda Askari, Hassan Momtaz, Elahe Tajbakhsh . Acinetobacter baumannii in sheep, goat, and camel raw meat: virulence and antibiotic resistance pattern. AIMS Microbiology, 2019, 5(3): 272-284. doi: 10.3934/microbiol.2019.3.272
    [8] Nikolina-Alexandra Xaxiri, Eleni Nikouli, Panagiotis Berillis, Konstantinos Ar. Kormas . Bacterial biofilm development during experimental degradation of Melicertus kerathurus exoskeleton in seawater. AIMS Microbiology, 2018, 4(3): 397-412. doi: 10.3934/microbiol.2018.3.397
    [9] Angel Valverde, María González-Tirante, Marisol Medina-Sierra, Raúl Rivas, Ignacio Santa-Regina, José M. Igual . Culturable bacterial diversity from the chestnut (Castanea sativa Mill.) phyllosphere and antagonism against the fungi causing the chestnut blight and ink diseases. AIMS Microbiology, 2017, 3(2): 293-314. doi: 10.3934/microbiol.2017.2.293
    [10] Mojtaba Moosavian, Mina Moradzadeh, Ataollah Ghadiri, Morteza Saki . Isolation and Identification of Legionella spp. in environmental water sources based on macrophage infectivity potentiator (mip) gene sequencing in southwest Iran. AIMS Microbiology, 2019, 5(3): 223-231. doi: 10.3934/microbiol.2019.3.223
  • We show that the quotient of Levy processes of jump-diffusion type has a fat-tailed distribution. An application is to price theory in economics, with the result that fat tails arise endogenously from modeling of price change based on an excess demand analysis resulting in a quotient of arbitrarily correlated demand and supply whether or not jump discontinuities are present. The assumption is that supply and demand are described by drift terms, Brownian (i.e., Gaussian) and compound Poisson jump processes. If P1dP/dt (the relative price change in an interval dt) is given by a suitable function of relative excess demand, (DS)/S (where D and S are demand and supply), then the distribution has tail behavior F(x)xζ for a power ζ that depends on the function G in P1dP/dt=G(D/S). For G(x)|x|1/q one has ζ=q. The empirical data for assets typically yields a value, ζ˜=3, or ζ[3,5] for many financial markets.

    Many theoretical explanations have been offered for the disparity between the tail behavior of the standard asset price equation and empirical data. This issue never arises if one models price dynamics using basic economics methodology, i.e., generalized Walrasian adjustment, rather than the usual starting point for classical finance which assumes a normal distribution of price changes. The function G is deterministic, and can be calibrated with a smaller data set. The results establish a simple link between the decay exponent of the density function and the price adjustment function, a feature that can improve methodology for risk assessment.

    The mathematical results can be applied to other problems involving the relative difference or quotient of Levy processes of jump-diffusion type.



    COVID-19 is the novel coronavirus that has spread among humans, mainly in China, since December 2019 [1]. It was first discovered in Wuhan, the capital city of the Hubei province [2]. The virus subsequently spread to other provinces in China. Cases of infection have also appeared in other countries [3,4]. COVID-19 is a respiratory virus that is transferred via contact with an infected person through droplets when a person coughs or sneezes, or through saliva droplets. The main clinical manifestations of the infection are fever, fatigue, respiratory symptoms (mainly dry cough), and emergence of dyspnea. Most patients have mild to moderate symptoms with good prognosis, but some patients can be in critical condition or die [1]. As of February 15, 2020, it has been reported that 68,500 people have confirmed with the COVID-19 infection and 1596 people have died in the mainland of China [5].

    To avoid the risk of a large-scale movement of population that can accelerate the spread of the disease, the Chinese government has implemented various measures. The government keeps track of people who had direct contact with the COVID-19 patients, known as contact tracing. From the last day of contact with a confirmed patient, contacts are monitored for signs of illness within 14 days. If contacts develop fever or other symptoms, they are isolated, tested, and hospitalized immediately to prevent further spread of the virus to others. Since Jan 23, 2020, metropolitan-wide quarantine and traffic restriction have been taken in Wuhan and several nearby cities. As the outbreak grows rapidly, China is imposing a more massive quarantine. All schools have postponed the start of the spring semester. Most companies have also postponed the starting date of work after the holiday of Spring Festival. A number of countries or regions have issued restrictions on the entry of Chinese citizens. A large number of domestic and international flights have been canceled. There is no doubt that the outbreak of this virus infection has severely affected people's life, economy and health. How long this situation will last and when the disease will be controlled is a great concern to everyone.

    In the last few decades, there were several major outbreaks of infectious diseases, such as atypical pneumonia (SARS) in 2003, H1N1 influenza in 2009, and H7N9 influenza in 2013. It is imperative to improve early predictive and warning capabilities for the disease endemic or pandemic. Mathematical models, combined with the prevalence data, have been used to study the dynamics, analyze the causes and key factors of the outbreaks, forecast the trend of disease spread, and provide optimal control strategies and measures. For example, Zhou and Ma [6] formulated a discrete mathematical model to investigate the transmission of SARS. Their simulation results agree with the data and indicate that early quarantine and a high quarantine rate are crucial to the control of SARS. Chowell et al. [7] and Lekone et al. [8] developed ordinary differential equations and stochastic SEIR models to study the dynamics of infectious disease and the effect of control interventions, respectively. Their models used the outbreak of Ebola in the Democratic Republic of Congo in 1995 as a case study. Very recently, some researchers have studied the spread of COVID-19 [9,10,11]. Most of them estimated the basic reproductive number R0, a key parameter to evaluate the potential of viral transmission [12,13,14,15,16]. Tang et al. [17,18] constructed a deterministic compartmental model and also estimated the basic reproductive number. They provided the predicted results under different degrees of control measures.

    In this study, we will use a discrete-time stochastic compartment model to study the dynamics of COVID-19 epidemics. The model includes the transmission of COVID-19 and the government's measures to control the disease spread. The model captures the epidemiological status of individuals in the clinical progression of the disease and the changes in each time interval. We make full use of the existing reported data and perform parameter estimation based on the data. Furthermore, we forecast the spread of the disease in the next period of time based on the parameter estimates and numerical simulation. We also use simulations to evaluate the risk of returning to work at different timing and provide suggestions on when people can start their routine work and life.

    Based on the development and epidemiological characteristics of COVID-19 infection, the SEIR model is appropriate to study the dynamic of this disease. The population is partitioned into subpopulations as susceptible (S(t)), exposed (E(t)), infected (I(t)), hospitalized (H(t)), and recovered (R(t)). Let N denote the total population size. The government takes measures to track and quarantine people who have close contact with confirmed cases. Thus, a fraction of the susceptible population is quarantined and identified as Sq(t) and a fraction of the exposed population is isolated and identified as Eq(t). We consider the discrete time point series T=0,1,2,,Tn as the time progression of the disease. Here the time step is chosen to be one day, i.e., h=1. At this timescale, the number of each compartment is dependent on the number in the previous day and the inflows and removals from other compartments during the day. Let Bij(t) be the number of individual transportation between compartments. We provide their detailed descriptions as follows:

    B11(t) is the number of susceptible individuals who become newly infected;

    B12(t) is the number of quarantined susceptible individuals who have contact with infected individuals but are not infected;

    B21(t) is the number of new cases with symptom onset;

    B31(t) is the number of new confirmed and admitted patients;

    B32(t) is the number of new death from infected individuals;

    B33(t) is the number of newly recovered from infected individuals;

    B41(t) is the number of people released from quarantine;

    B51(t) is the number of people admitted to hospital (also isolated);

    B61(t) is the number of newly recovered from hospitalized cases;

    B62(t) is the number of new death from hospitalized cases.

    The transition of an individual from one state to the next state can be considered as a stochastic process. The time length that an individual has been in a certain compartment obeys exponential distribution. If we assume that the parameter of the exponential distribution is λ(t), then the probability that individuals leave the current state in the time interval h is 1exp(λ(t)h). Further, the numbers of inflows and removals from other compartments during one day can be generated by a binomial distribution. The number of experiments in the binomial distribution is the number of individuals in the current compartment. The transmission of the disease is presumed to occur in the context of close contact between susceptible and infected persons. Assuming the transmission probability is β and the contact rate is c(t), then βc(t)I(t)N is the exponential rate that leads to the probability of individuals leaving the susceptible compartment. In view of contact tracing, we denote the quarantined proportion of individuals exposed to the virus is q. Based on the above assumptions and the stochastic SEIR model in [8], the discrete-time stochastic compartment model for COVID-19 infection is constructed as

    S(t+h)=S(t)B11(t)B12(t)+B41(t),E(t+h)=E(t)+(1q)B11(t)B21(t),I(t+h)=I(t)+B21(t)B31(t)B32(t)B33(t),Sq(t+h)=Sq(t)+B12(t)B41(t),Eq(t+h)=Eq(t)+qB11(t)B51(t),H(t+h)=H(t)+B31(t)+B51(t)B61(t)B62(t),R(t+h)=R(t)+B33(t)+B61(t), (2.1)

    where

    B11(t)Poi(S(t)P11(t)),  B12(t)Poi(S(t)P12(t)),B21(t)Bin(E(t),P21),  B31(t)Bin(I(t),P31),B32(t)Bin(I(t),P32),  B33(t)Bin(I(t),P33),B41(t)Bin(Sq(t),P41),  B51(t)Bin(Eq(t),P51),B61(t)Bin(H(t),P61),  B62(t)Bin(H(t),P32).

    The above random variables involve binomial distributions Bin(n,p) with the following probabilities:

    P11(t))=1exp[βc(t)I(t)Nh],  P12(t)=1exp[c(t)q(1β)I(t)Nh],P21=1exp(σh),  P31=1exp(δIh),P32=1exp(αh),  P33=1exp(γIh),P41=1exp(λh),  P51=1exp(δqh),P61=1exp(γHh).

    Note that the number of S(t) is approximately the total population N in China, which is a large number. The limit distribution of the binomial distribution is the poisson distribution. This is used instead for B11(t) and B12(t). The variation of these seven compartments and their relationship are illustrated in the diagram 1. The other parameters are summarized in Table 1.

    Table 1.  Estimates of the parameters and initial values of the model.
    Parameters Definition Baseline values Range Mean value (estimated)Std Source
    c0 Initial contact rate 31 (0, 50] 34.037 0.389 MCMC
    cu Finial contact rate 1 (0, 50] 0.933 0.0037 MCMC
    under control measures
    k Exponential decreasing rate 0.1 [0, 1] 0.144 0.0035 MCMC
    of the contact rate
    β Probability of transmission 0.095 [0, 1] 0.111 0.0015 MCMC
    q Quarantine rate of exposed individuals 0.4 [0, 1] 0.415 0.016 MCMC
    σ Transition rate of exposed individuals 1/7 [0, 1] 1/7 - [17]
    to the infected class
    λ Rate at which the quarantined uninfected 1/14 1/14 1/14 - [17]
    contacts were released
    δI Transition rate of symptomatic infected 0.1 [0.1] 0.304 0.001 MCMC
    individuals to the quarantined infected class
    δq Transition rate of quarantined exposed 0.42 [0, 1] 0.413 0.0116 MCMC
    individuals to the quarantined infected class
    γI Recovery rate of symptomatic infected 0.008 [0, 1] 0.0085 0.0001 MCMC
    individuals
    γH Recovery rate of quarantined infected 0.017 [0, 1] 0.018 0.0003 MCMC
    individuals
    α Disease-induced death rate 0.0027 [0, 1] 0.0027 0.0001 MCMC
    N The total population number - - 1.4×1010 - Data
    E(0) Initial exposed population 105 [1,200] 52 0.64 MCMC
    I(0) Initial symptomatic infected population 54 [1,100] 31 0.346 MCMC
    Sq(0) Initial quarantined susceptible population - - 734 - Data
    Eq(0) Initial quarantined exposed population 2 [1,100] 5 0.24 MCMC
    H(0) Initial quarantined infected population - - 41 - Data
    R(0) Initial recovered population - - 6 - Data

     | Show Table
    DownLoad: CSV

    Reducing exposure is one of the effective measures to control the spread of disease. The public raises awareness of precaution and takes fewer trips or wears a mask in the presence of a disease outbreak. The government has taken more control measures, e.g., blocking the entrances and exits of the city to stop the movement of population. In this model, the contact rate c(t) (i.e. the average number of contacts of a person per unit time) is assumed to be a piecewise function. It is a constant before the time point t at which control measures are taken. It is assumed to decrease gradually from c0 to cu. We assume that c(t) takes the following form:

    c(t)={c0,  tt,(c0cu)ek(tt)+cu,  t>t. (2.2)

    The city of Wuhan was blockaded on January 23, 2020. The actual data we used for data fitting start from January 11, 2020. Thus, we let t=12.

    We obtained the data from the National Health Commission of the People's Republic of China [5] and the Health Commission of the Hubei Province [19]. The data that will be used for parameter estimation include the cases in the mainland of China, such as the newly reported confirmed cases, the newly recovered cases, the new death cases, and the number of people released from medical observation. The total population N in this study is the population of China, approximately 1.4 billion. The cumulative number of reported confirmed cases can be used to compare with the model simulations. The total number of confirmed cases is the cumulative number of cases minus the cumulative number of recovered and death cases. As a complete reporting coverage of the cases has been available since January 11, 2020, the data used in our analysis are from January 11 to February 13, 2020. Since February 13, the Hubei province has carried out screening of the previous suspected cases and revised the diagnosis results, adding "clinical diagnosis" to the diagnosis classification. Although the number of clinical cases has been public as confirmed cases in Hubei, the clinical diagnosis is not included in the number of confirmed cases in this study.

    In general, Bayesian estimation or the maximum likelihood estimation method can be used to estimate unknown parameters in this type of stochastic models [8,20]. Because Bij(t), where (i,j) {(1, 1), (1, 2), (2, 1), (3, 1), (3, 2), (3, 3), (4, 1), (5, 1), (6, 1), (6, 2)}, are conditionally independent, the likelihood function can be the accumulation of probability densities of all variables, given by

    L(B11(t),B12(t),B21(t),B31(t),B32(t),B33(t),B41(t),B51(t),B61(t),B62(t)|Θ)=Tnt=0gi,j(Bij(t)|) (4.1)

    where gi,j is the binomial densities of Bij(t). The best scenario is that the value of transportation between compartments can be obtained. However, the exposed and some infected people are impossible to be identified. In the reported cases about COVID-19 infection, the data we obtained is the time series for B41,B61,B62, and the number of new cases, which is equal to the sum of B51 and B31. The remaining time series including B11,B12,B21,B32,B33 are hard to be determined. It is difficult to accomplish parameter estimates using Bayesian estimation or the maximum likelihood estimation due to such a data structure. Therefore, the Metropolis-Hastings (MH) algorithm will be used to estimate parameters in our model by carrying out the Markov Chain Monte Carlo (MCMC) procedure [21,22].

    In the model, we assume σ=1/7 and λ=1/14 because the quarantined individuals were isolated for 14 days and the incubation period of COVID-19 is about 7 days [17,23]. The initial conditions are the data on January 11, that is, H0=41, R0=6, and the sum of Sq and Eq is 739. We estimated the parameters from the mean values of 10,000 samples after a burn-in period of 5000 iterations. The estimated results of all parameters relevant to COVID-19 infection are displayed in Table 1. The model parameters fitted to the data of COVID-19 infection, including the newly reported confirmed cases and the total confirmed cases, are shown in Figure 2. In addition, stochastic simulations in Figure 3 show that the model can provide a good fit to the data of newly recovered cases and new death cases of the COVID-19 infection.

    Figure 1.  Flow diagram of the COVID-19 infection.
    Figure 2.  The data of newly reported confirmed cases and total number of confirmed cases of COVID-19 disease from January 11, 2020 to February 13, 2020. Stochastic fit was performed 100 times. The newly reported confirmed case data are shown in (a) and the total number of confirmed case data are shown in (b). The data and the fitted are represented by the deep blue and light blue curves, respectively.
    Figure 3.  The data of newly recovered cases and new death cases of COVID-19 disease from January 11 to February 13, 2020. Stochastic fit was performed 100 times. The newly recovered cases are shown in (a), death cases are shown in (b). The data and the fitted are represented by the deep blue and light blue curves, respectively.

    Comparing with the basic reproductive number, the effective reproductive number can be used to measure the number of secondary cases generated by one primary case in a population in which there is partial immunity or some intervention measures have been implemented [24,25,26,27,28]. It changes during the progress of the disease outbreak. As the population size is much larger than the resulting size of the outbreak, i.e., S(t)/N1, the effective or control reproductive number of our model is given by the following formula

    Rc(t)=βc(t)(1q)δI+α+γI.

    In contrast with the basic reproductive number that usually involves a constant contact rate, we use a time-varying contact rate c(t) in the effective reproductive number [17]. Using the parameter values estimated from the data fitting, the effective reproductive number Rc(t) can be computed numerically. The result shown in Figure 4 indicates that Rc(t) is large at the beginning of the disease outbreak. However, as control measures are implemented, the effective reproductive number decreases eventually to less than 1.

    Figure 4.  The effective reproductive number Rc(t) from January 11 to February 13, 2020. The gray line is the range due to disturbance in the stochastic model.

    An assumption of our model is that the contact rate is exponentially decreasing to cu due to various prevention measures. The estimates show that the maximum value of the contact rate is 34.03 at the beginning of disease spread. As control measures are implemented, the contact rate declines to 0.93. This indicates that susceptible people are well protected by these measures and there is a small chance of transmission by infected people. Using this assumption and estimated parameters, we can make predictions about the future trend of the disease spread. In Figure 5(a) and (b), the simulations show the dynamics of the number of newly confirmed cases and the total confirmed cases in the 350 days after January 11, 2020. The number of newly cases has begun to decrease slightly compared with before. Figure 5(a) also indicates that the number of newly confirmed cases will continue to decline in the future. The number of newly confirmed cases will decrease to 0 in about 102–119 days after January 11, i.e., April 22, 2020 to May 9, 2020. Figure 5(b) shows that the number of total confirmed cases is still increasing. It will reach the peak in about 46–48 days after January 11, that is, February 26, 2020 to February 28, 2020.

    Figure 5.  Predicted results of newly confirmed cases shown in (a) and total number of confirmed cases shown in (b) under control measures in the 350 days since January 11, 2020. The red dotted line in (a) is the time when the newly confirmed cases reach zero. The red dotted line in (b) is the time when the total confirmed cases reach the peak level.

    Up to now, the Chinese government has adopted very strict measures to control the spread of the disease. Most people are quarantined at home. Most companies have stopped working and all schools have postponed the opening date. How long will such measures last? If people start to work and students go to school as usual, the population flow will increase the risk of contact. Will this lead to a second outbreak of the disease? When does people's life can go back to normal? These are issues of great concern to everybody in the country. We consider the following four scenarios:

    (i). Suppose the time to return to work is March 1 but the protective measures are inadequate. The parameters in the model are assumed to be c(t)=3 when T>50;

    (ii). Suppose the time to return to work is March 1 and the protective measures are good. The parameters are chosen to be c(t)=1.5 when T>50;

    (iii). Suppose the time to return to work is March 20 but the protective measures are inadequate. The parameters are c(t)=3 when T>70;

    (iv). Suppose the time to return to work is March 20 and the protective measures are good. The parameters are c(t)=1.5 when T>70;

    Using the above assumptions and fixing other parameters according to data fitting, we numerically investigate the dynamics of the number of newly confirmed cases and total confirmed cases, which are shown in Figure 6. The simulations suggest that when people return to work early without sufficient protection it is likely to observe the increase of newly confirmed cases. The total number of confirmed cases will also increase by more than 2000. If people return to work on March 20, it will result in a small increase in the number of newly confirmed cases. However, it's going to decrease quickly. Thus, postponing the return to work would be of great help to control the disease transmission.

    Figure 6.  Predicted results of newly confirmed cases shown in (a) and total number of confirmed cases shown in (b) with four assumptions (see text) in the 350 days since January 11, 2020. The results are the mean of 50 stochastic simulations.

    In view of the randomness in the transmission of the COVID-19 infection, we construct a discrete-time stochastic compartment system to study the dynamic behavior of the disease outbreak, in which the population in each compartment is assumed to obey a binomial distribution. In the model, we further assume that the contact rate between susceptible and infected individuals decreases exponentially since the government has implemented strict control measures. This discrete system can make a good use of the newly reported data, including the number of newly reported confirmed cases, newly recovered cases, new death cases, and those quarantined that are released per day, to calibrate the model parameters. Because the diagnostic criterion for confirmed cases was revised on February 13, we use the data reported from January 11 to February 13, 2020 in this study. Although there are no data for the populations in the other groups, we estimated the parameters in the model using the MCMC procedure based on the exiting reported data.

    The maximum value of contact rate we estimated is 34.02 and then dropped to 0.93. It indicates that the control measures implemented are very effective. Figure 3 shows that the effective reproductive number of the infection declines from 6.96 to 0.47 over the time period of simulation. Some other groups have also estimated R0 of the COVID-19 infection [29,30]. See a summary of the estimates in the reference [31]. Some estimates of R0 are comparable to ours (before various control measures are taken) and some are smaller. This may be due to different estimation methods under different model assumptions. In any case, the results indicate that this novel coronavirus is highly contagious in the early stage (e.g., higher than the SARS coronavirus outbreak in 2003 [31]). However, the control measures implemented so far are shown to be efficient in lowering the effective reproductive number to below 1.

    The stochastic fits show that the reported data are less than those simulated in the initial phase, but the subsequent simulation agrees well with the reported data. Most of the data are in the region of the stochastic simulation. By the time of the submission of this manuscript, the number of newly reported confirmed cases has begun to decrease. Assuming that the current control measures adopted by the government and individuals are maintained, the predicted results shown in Figure 4 indicate that the peak of total confirmed cases will reach around late February of 2020, followed by a decrease. The newly confirmed cases will decline to zero in late April or early May of 2020. To study the timing of returning to work on the disease dynamics, we consider four different scenarios on the time to return to work and the strength of protective measures. We assume the contact rate is 3 if the resumption of work causes wide migration of people and the protective measures are not sufficient. In this case, the simulation shows that the infection will have a second outbreak if people return to work on March 1. The situation will be better if people can return to work after March 20, 2020.

    In conclusion, our simulation shows that the contact rate is a key factor on the control of the COVID-19 outbreak. When there is a sign of epidemic, we should raise awareness of self-protection and take all possible protective measures, such as wearing a mask and staying indoor to reduce the risk of getting infected. Although the number of new cases of infection is decreasing, there is still a possibility of future outbreaks if there are no adequate protective measures or people return to work early. The public should not relax their vigilance against the transmission of this highly contagious disease.

    This research was partially supported by the National Natural Science Foundation of China (grant numbers: 61772017 (ST), 11631012 (ST)) and by the Fundamental Research Funds for the Central Universities (grant numbers: 2018CBLZ001 (SH), GK201901008 (ST)). L. Rong is supported by the National Science Foundation (grant number: DMS-1758290).

    No conflict of interest.



    [1] M. Aoki, H. Yoshikawa, Reconstructing Macroeconomics: A Perspective from Statistical Physics and Combinatorial Stochastic Processes, Cambridge University Press, 2006.
    [2] L. Bachelier, Théorie de la spéculation, Ann. Sci. École Norm. Sup., 17 (1900), 21–86.
    [3] W. A. Brock, C. H. Hommes, Heterogeneous beliefs and routes to chaos in a simple asset pricing model, Journal of Economic Dynamics and Control, 22 (1998), 1235–1274. doi: 10.1016/S0165-1889(98)00011-6
    [4] J. P. Bouchaud, M. Potters, Theory of financial risks. From Statistical Physics to Risk Management, Cambridge University Press, 2000.
    [5] F. Black, M. Scholes, The pricing of options and corporate liabilities, J. Political Economy, 81 (1973), 637–654. doi: 10.1086/260062
    [6] G. Caginalp, D. Balevonich, Asset flow and momentum: Deterministic and stochastic equations, Phil. Trans. Royal Soc. A., 357 (1999), 2119–2113. doi: 10.1098/rsta.1999.0421
    [7] C. Caginalp, G. Caginalp, The quotient of normal random variables and application to asset price fat tails, Physica A, 499 (2018), 457–471. doi: 10.1016/j.physa.2018.02.077
    [8] C. Caginalp, G. Caginalp, Price equations with symmetric supply/demand; implications for fat tails, Economics Letters, 176 (2019), 79–82. doi: 10.1016/j.econlet.2018.12.037
    [9] G. Caginalp, G. B. Ermentrout, A kinetic thermodynamics approach to the psychology of fluctuations in financial markets, Appl. Math. Lett., 3 (1990), 17–19.
    [10] C. Carroll, The method of endogenous gridpoints for solving dynamic stochastic optimization problems, Economics Letters, 91 (2006), 312–320. doi: 10.1016/j.econlet.2005.09.013
    [11] N. Champagnat, M. Deaconu, A. Lejay, N. Navet, S. Boukherouaa, An empirical analysis of heavy-tails behavior of financial data: The case for power laws, HAL archives-ouvertes, 2013.
    [12] R. Cont, P. Tankov, Financial modeling with jump processes, Chapman & Hall/CRC, London, 2004.
    [13] J. Daníelsson, B. N. Jorgensen, G. Samorodnitsky, M. Sarma, C. G. de Vries, Fat tails, VaR and subadditivity, J. Econometrics, 172 (2013), 283–291. doi: 10.1016/j.jeconom.2012.08.011
    [14] E. Díaz-Francés, F. J. Rubio, On the existence of a normal approximation to the distribution of the ratio of two independent normal random variables, Stat. Papers, 54 (2013), 309–323. doi: 10.1007/s00362-012-0429-2
    [15] M. DeSantis, D. Swigon, Slow-fast analysis of a multi-group asset flow model with implications for the dynamics of wealth, PloS one, 13 (2018), e0207764. doi: 10.1371/journal.pone.0207764
    [16] C. Eom, T. Kaizoji, S. H. Kang, L. Pichl, Bitcoin and investor sentiment: statistical characteristics and predictability, Physica A, 514 (2019), 511–521. doi: 10.1016/j.physa.2018.09.063
    [17] C. Eom, T. Kaizoji, E. Scalas, Fat tails in financial return distributions revisited: Evidence from the Korean stock market, Physica A, 526 (2019), 121055. doi: 10.1016/j.physa.2019.121055
    [18] A. Erdelyi, Asymptotic Expansions, Dover Publications, 1956.
    [19] L. C. Evans, Partial Differential Equations, 2nd Ed., Am. Math. Soc., 2010.
    [20] G. W. Evans, S. Honkapohja, Stochastic gradient learning in the cobweb model, Economics Letters, 61 (1998), 333–337. doi: 10.1016/S0165-1765(98)00176-1
    [21] E. Fama, The Behavior of Stock-Market Prices, J. Business, 38 (1965), 34–105. doi: 10.1086/294743
    [22] J. P. Fouque, G. Papanicolaou, R. Sircar, K. Sølna, Multiscale stochastic volatility for equity, interest rate, and credit derivatives, Cambridge University Press, 2011.
    [23] X. Gabaix, Power laws in economics: An introduction. Journal of Economic Perspectives, 30 (2016), 185–206.
    [24] X. Gabaix, P. Gopikrishnan, V. Plerou, H. E. Stanley, Institutional investors and stock market volatility, The Quarterly Journal of Economics, 121 (2006), 461–504. doi: 10.1162/qjec.2006.121.2.461
    [25] S. Gjerstad, J. Dickhaut, Price formation in double auctions, Games and economic behavior, 22 (1998), 1–29. doi: 10.1006/game.1997.0576
    [26] S. Gjerstad, The competitive market paradox, Journal of Economic Dynamics and Control, 31 (2007), 1753–1780. doi: 10.1016/j.jedc.2006.07.001
    [27] S. Gjerstad, Price dynamics in an exchange economy, Economic Theory, 52 (2013), 461–500. doi: 10.1007/s00199-011-0651-5
    [28] R. Geary, The frequency distribution of the quotient of two normal variates, J. Royal Stat. Soc., 93 (1930), 442–446. doi: 10.2307/2342070
    [29] P. Gopikrishnan, M. Meyer, L. Amaral, H. E. Stanley, Inverse cubic law for the distribution of stock price variations, Eur. Phys. J. B, 3 (1998), 139–140.
    [30] J. Henderson, R. Quandt, Microeconomic theory - A mathematical approach, McGraw-Hill, New York, 1980.
    [31] J. Hirshleifer, A. Glazer, D. Hirshleifer, Price theory and applications: decisions, markets, and information, Cambridge Univ. Press, 2005.
    [32] J. C. Hull, Options, Futures and Other Derivatives, 9th Ed., Pearson, 2015.
    [33] T. D. Jeitschko, Equilibrium price paths in sequential auctions with stochastic supply, Economics Letters, 64 (1999), 67–72. doi: 10.1016/S0165-1765(99)00066-X
    [34] D. W. Jansen, C. G. De Vries, On the frequency of large stock returns: Putting booms and busts into perspective, The review of economics and statistics, (1991), 18–24.
    [35] K. Kawamura, The structure of bivariate Poisson distribution, Kodai Math Sem. Rep., 25 (1973), 246–256. doi: 10.2996/kmj/1138846776
    [36] M. Kemp, Extreme events - robust portfolio construction in the presence of fat tails, Wiley Finance, Hoboken, NJ, 2011.
    [37] P. J. Klenow, B. A. Malin, Microeconomic evidence on price-setting, Handbook of monetary economics, 3 (2010), 231–284. doi: 10.1016/B978-0-444-53238-1.00006-5
    [38] M. Kirchler, J. Huber, Fat tails and volatility clustering in experimental asset markets, Journal of Economic Dynamics and Control, 31 (2007), 1844–1874. doi: 10.1016/j.jedc.2007.01.009
    [39] T. Lux, The stable Paretian hypothesis and the frequency of large returns: an examination of major German stocks, Applied financial economics, 6 (1996), 463–475. doi: 10.1080/096031096333917
    [40] B. LeBaron, R. Samanta, Extreme value theory and fat tails in equity markets, 2005. Available at SSRN 873656.
    [41] B. Mandelbrot, Sur certain prix speculatifs: faits empirique et modele base sur les processes stables additifs de Paul Levy, Comptes Rendus de l'Academie de Sciences, 254 (1962), 3968–3970.
    [42] A. M. Mood, F. A. Graybill, D. C. Boes, Introduction to the theory of statistics, 3rd Ed., McGraw-Hill, 1973.
    [43] B. Mandelbrot, R. Hudson, The misbehavior of markets: A fractal view of financial turbulence, Basic Books, New York, 2007.
    [44] G. Marsaglia, Ratios of normal variables, J. Stat. Softw., 16 (2006), 1–10.
    [45] H. Merdan, M. Alisen, A mathematical model for asset pricing, App. Math. Comp., 218 (2011), 1449–1456. doi: 10.1016/j.amc.2011.06.028
    [46] P. Milgrom, Discovering Prices: Auction Design in Markets with Complex Constraints, Columbia University Press, 2017.
    [47] R. N. Mantegna, H. E. Stanley, Introduction to econophysics: correlations and complexity in finance, Cambridge University press, 1999.
    [48] T. Pham-Gia, N. Turkkan, E. Marchand, Density of the ratio of two normal random variables and applications, Communications in Statistics-Theory and Methods, 35 (2006), 1569–1591. doi: 10.1080/03610920600683689
    [49] C. Plott, K. Pogorelskiy, Call market experiments: efficiency and price discovery through multiple calls and emergent Newton adjustments, American Economic Journal: Microeconomics, 9 (2017), 1–41.
    [50] D. Porter, S. Rassenti, A. Roopnarine, V. Smith, Combinatorial auction design, Proceedings of the National Academy of Sciences, 100 (2003), 11153–11157. doi: 10.1073/pnas.1633736100
    [51] S. T. Rachev, Handbook of Heavy Tailed Distributions in Finance, Handbooks in Finance, Book 1, Elsevier, 2003.
    [52] M. Raberto, S. Cincotti, S. M. Focardi, M. Marchesi, Agent-based simulation of a financial market, Physica A, 299 (2001), 319–327. doi: 10.1016/S0378-4371(01)00312-0
    [53] S. T. Rachev, C. Menn, F. J. Fabozzi, Fat-tailed and skewed asset return distributions: implications for risk management, portfolio selection, and option pricing, John Wiley & Sons, 2005.
    [54] Z. Rachev, B. Rocheva-Iotovo, S. Stoyanov, Capturing fat tails, Risk, 23 (2010), 72.
    [55] A. W. Rathgeber, J. Stadler, S. Stöckl, Financial modelling applying multivariate Lévy processes: New insights into estimation and simulation, Physica A, 532 (2019), 121386. doi: 10.1016/j.physa.2019.121386
    [56] Y. Si, S. Nadarajah, X. Song, On the distribution of quotient of random variables conditioned to the positive quadrant, Communications in Statistics-Theory and Methods, 49 (2020), 2514–2528. doi: 10.1080/03610926.2019.1576893
    [57] N. Taleb, G. Daniel, The problem is beyond psychology: The real world is more random than regression analyses, International Journal of Forecasting, Forthcoming, 2011.
    [58] Y. Tong, The multivariate normal distribution, Springer-Verlag, New York, 1990.
    [59] D. Watson, M. Getz, Price theory and its uses, 5th Ed., University Press of America, 1981.
    [60] E. R. Weintraub, Microfoundations, Cambridge University Press, 1979.
  • This article has been cited by:

    1. Millicent R. Smith, Idupulapati M. Rao, 2021, 9780128191941, 384, 10.1016/B978-0-12-819194-1.00012-8
    2. A. H. Gunnabo, J. van Heerwaarden, R. Geurts, E. Wolde-meskel, T. Degefu, K. E. Giller, Phylogeography and Symbiotic Effectiveness of Rhizobia Nodulating Chickpea (Cicer arietinum L.) in Ethiopia, 2021, 81, 0095-3628, 703, 10.1007/s00248-020-01620-8
    3. Mark A. Uebersax, Karen A. Cichy, Francisco E. Gomez, Timothy G. Porch, Jim Heitholt, Juan M. Osorno, Kelvin Kamfwa, Sieglinde S. Snapp, Scott Bales, Dry beans ( Phaseolus vulgaris L.) as a vital component of sustainable agriculture and food security—A review , 2022, 2639-6181, 10.1002/leg3.155
    4. Ravinder K. Goyal, Autar K. Mattoo, Maria Augusta Schmidt, Rhizobial–Host Interactions and Symbiotic Nitrogen Fixation in Legume Crops Toward Agriculture Sustainability, 2021, 12, 1664-302X, 10.3389/fmicb.2021.669404
    5. Anna Szczerba, Agnieszka Płażek, Przemysław Kopeć, Magdalena Wójcik-Jagła, Franciszek Dubert, Effect of different Bradyrhizobium japonicum inoculants on physiological and agronomic traits of soybean (Glycine max (L.) Merr.) associated with different expression of nodulation genes, 2024, 24, 1471-2229, 10.1186/s12870-024-05911-x
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3214) PDF downloads(100) Cited by(0)

Figures and Tables

Figures(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog