
The global impact of the COVID-19 pandemic is widely recognized as a significant concern, with human flow playing a crucial role in its propagation. Consequently, recent research has focused on identifying and analyzing factors that can effectively regulate human flow. However, among the multiple factors that are expected to have an effect, few studies have investigated those that are particularly associated with human flow during the COVID-19 pandemic. In addition, few studies have investigated how regional characteristics and the number of vaccinations for these factors affect human flow. Furthermore, increasing the number of verified cases in countries and regions with insufficient reports is important to generalize conclusions. Therefore, in this study, a group-level analysis was conducted for Narashino City, Chiba Prefecture, Japan, using a human flow prediction model based on machine learning. High-importance groups were subdivided by regional characteristics and the number of vaccinations, and visual and correlation analyses were conducted at the factor level. The findings indicated that tree-based models, especially LightGBM, performed better in terms of prediction. In addition, the cumulative number of vaccinated individuals and the number of newly infected individuals are likely explanatory factors for changes in human flow. The analyses suggested a tendency to move with respect to the number of newly infected individuals in Japan or Tokyo, rather than the number of new infections in the area where they lived when vaccination had not started. With the implementation of vaccination, attention to the number of newly infected individuals in their residential areas may increase. However, after the spread of vaccination, the perception of infection risk may decrease. These findings can contribute to the proposal of new measures for efficiently controlling human flows and determining when to mitigate or reinforce specific measures.
Citation: Masaya Mori, Yuto Omae, Yohei Kakimoto, Makoto Sasaki, Jun Toyotani. Analyzing factors of daily travel distances in Japan during the COVID-19 pandemic[J]. Mathematical Biosciences and Engineering, 2024, 21(8): 6936-6974. doi: 10.3934/mbe.2024305
[1] | Dao Viet Thang, Nguyen Manh Hung, Nguyen Cao Khang, Le Thi Mai Oanh . Structural and multiferroic properties of (Sm, Mn) co-doped BiFeO3 materials. AIMS Materials Science, 2020, 7(2): 160-169. doi: 10.3934/matersci.2020.2.160 |
[2] | Xuan Luc Le, Nguyen Dang Phu, Nguyen Xuan Duong . Enhancement of ferroelectricity in perovskite BaTiO3 epitaxial thin films by sulfurization. AIMS Materials Science, 2024, 11(4): 802-814. doi: 10.3934/matersci.2024039 |
[3] | Z. Aboub, B. Daoudi, A. Boukraa . Theoretical study of Ni doping SrTiO3 using a density functional theory. AIMS Materials Science, 2020, 7(6): 902-910. doi: 10.3934/matersci.2020.6.902 |
[4] | Yu Shiratsuchi, Yoshinori Kotani, Saori Yoshida, Yasunori Yoshikawa, Kentaro Toyoki, Atsushi Kobane, Ryoichi Nakatani, Tetsuya Nakamura . Direct observations of ferromagnetic and antiferromagnetic domains in Pt/Co/Cr2O3/Pt perpendicular exchange biased film. AIMS Materials Science, 2015, 2(4): 484-496. doi: 10.3934/matersci.2015.4.484 |
[5] | Wolfgang Zúñiga-Mera, Sonia Gaona Jurado, Alejandra Isabel Guerrero Duymovic, Claudia Fernanda Villaquirán Raigoza, José Eduardo García . Effect of the incorporation of BiFeO3 on the structural, electrical and magnetic properties of the lead-free Bi0.5Na0.5TiO3. AIMS Materials Science, 2021, 8(5): 792-808. doi: 10.3934/matersci.2021048 |
[6] | Murat Aycibin, Naciye ECE . First-principles calculation of the electronic and optical properties of BiRhO3 compound. AIMS Materials Science, 2017, 4(4): 894-904. doi: 10.3934/matersci.2017.4.894 |
[7] | Evangelos Karagiannis, Dimitra Papadaki, Margarita N. Assimakopoulos . Circular self-cleaning building materials and fabrics using dual doped TiO2 nanomaterials. AIMS Materials Science, 2022, 9(4): 534-553. doi: 10.3934/matersci.2022032 |
[8] | Asima Adak (Maity), Soumya Mukherjee, Mahua Ghosh Chaudhuri, Siddhartha Mukherjee . Phase evolution, microscopic analysis, optical and dielectric property evaluation of Co-doped BaSnO3 by mechanical mix assisted solid state sintering method. AIMS Materials Science, 2016, 3(3): 1281-1293. doi: 10.3934/matersci.2016.3.1281 |
[9] | Yana Fajar Prakasa, Sumari Sumari, Aman Santoso, Muhammad Roy Asrori, Ririn Cahyanti . The performance of radar absorption of MnxFe3–xO4/rGO nanocomposites prepared from iron sand beach and coconut shell waste. AIMS Materials Science, 2023, 10(2): 227-248. doi: 10.3934/matersci.2023013 |
[10] | Leydi J. Cardenas F., Josep Ma. Chimenos, Luis C. Moreno A., Elaine C. Paris, Miryam R. Joya . Enhancing Co3O4 nanoparticles: Investigating the impact of nickel doping and high-temperature annealing on NiCo2O4/CoO heterostructures. AIMS Materials Science, 2023, 10(6): 1090-1104. doi: 10.3934/matersci.2023058 |
The global impact of the COVID-19 pandemic is widely recognized as a significant concern, with human flow playing a crucial role in its propagation. Consequently, recent research has focused on identifying and analyzing factors that can effectively regulate human flow. However, among the multiple factors that are expected to have an effect, few studies have investigated those that are particularly associated with human flow during the COVID-19 pandemic. In addition, few studies have investigated how regional characteristics and the number of vaccinations for these factors affect human flow. Furthermore, increasing the number of verified cases in countries and regions with insufficient reports is important to generalize conclusions. Therefore, in this study, a group-level analysis was conducted for Narashino City, Chiba Prefecture, Japan, using a human flow prediction model based on machine learning. High-importance groups were subdivided by regional characteristics and the number of vaccinations, and visual and correlation analyses were conducted at the factor level. The findings indicated that tree-based models, especially LightGBM, performed better in terms of prediction. In addition, the cumulative number of vaccinated individuals and the number of newly infected individuals are likely explanatory factors for changes in human flow. The analyses suggested a tendency to move with respect to the number of newly infected individuals in Japan or Tokyo, rather than the number of new infections in the area where they lived when vaccination had not started. With the implementation of vaccination, attention to the number of newly infected individuals in their residential areas may increase. However, after the spread of vaccination, the perception of infection risk may decrease. These findings can contribute to the proposal of new measures for efficiently controlling human flows and determining when to mitigate or reinforce specific measures.
Growth models have traditionally been used for representing the time course of microbial concentration [1,2,3]. Also, there have been applications for representing the time course of other variables, namely: Colony diameter [4,5,6,7]; tumor size [8]; weight of plant seedlings [9]; dry weight of wheat grains [10]; cumulative epicotyl emergence (germination) of oaks [11]; body weight of Boer goats [12], weight of pigs [13], size of Desinognathus [14], pressure variation [15], among others.
The quality of growth models fitting may be assessed in terms of its capability for representing the slope at the inflection point (k) and the duration of the lag phase (tlag), rather than the squared error (SSE) [16]. The lag phase duration indicates the time that the exponential phase begins. The inflection point characterizes the shape of the curve, and it allows to determine the maximum slope of the growth curve and the lag phase duration via the classical approach. In turn, the maximum growth rate is a measure of the growth speed during the exponential phase. Particularly, in the study of growth of undesired microorganisms on products, the lag phase duration indicates the initial time of the exponential phase of contaminants, and the maximum growth rate indicates its speed of growth [16,17].
One way to have growth models capable of generating the lag phase is through a differential equation, including a multiplicative term that suppresses the vector field for low values of the population size, i.e. during the lapse after the initial time. In the model proposed by Baranyi, Roberts and McClure [1], a time dependent adaptation function called adjustment function (αt) is incorporated into a logistic-type model, yielding a non-autonomous model and achieving satisfactory representation of the lag phase. However, autonomous differential equations also have a significant capability for representing the lag phase and the inflection point, some examples are: The Verhulst's logistic equation [18], the generalized Gompertz function [19], the Putter-Bertalanffy model [20], the generalized logistic model [19], the autonomous neoclassical model of Vadasz, Vadasz [21], and the model of Fujikawa, Kai, Morozumi [16]. Also, they allow performing equilibrium stability analysis [20]. Additional discussion is presented in section 2.
It is worth noticing that Tsoularis and Wallace [19] and Ohnishi, Yamakawa and Akamine [20] determined the analytical expression for the population at the inflection point, and assessed the effect of exponential parameter values on the existence or absence of the inflection point and lag phase, based on the curve of the growth rate versus population size. Tsoularis and Wallace [19] noticed that there is no logarithmic inflection point for certain parameter values. Also, they determined the limit values of the population at the inflection point for parameter values approaching zero or infinite. Ohnishi, Yamakawa and Akamine [20] observed that negative parameter values also lead to inflection point and lag phase, thus representing real life behavior. In contrast to these studies, in this work we consider the generalized Fujikawa's model which included the function fc=(1−Nmin/N)c that characterizes the original Fujikawa's model, the power law function Nα and a saturation term. We assess the effect of exponential parameters on the growth curve, considering separately the population N and the logarithmic population y=ln(N/Nto), with negative and positive parameter values. Indeed: ⅰ) we determine the ranges of parameter values that lead to existence or absence of inflection point, and those for the existence or absence of lag phase; ⅱ) we assess the effect of different combinations of positive and negative parameter ranges on the existence or absence of inflection point and lag phase; ⅲ) we determine the parameter ranges that lead to monotonic decreasing behavior of the slope of the growth rate curve; ⅳ) we propose an approach for the analytical determination of the lag phase duration. Numerical simulations illustrate the approach. In summary, the main contribution of this work is the determination of the capability of a generalized Fujikawa's model to capture all the features of real life monotonic growth behavior, including lag phase and inflection point.
The organization of the study is as follows. Section 2 presents some Preliminaries on the study of growth models. Section 3 presents the proposed model and the properties of the adaptation function. Section 4 presents the features of the regular growth model: Convergence and boundedness properties; determination of the regular inflection point and conditions for its existence; effect of model parameters on the existence of the inflection point and the lag phase; and properties of the regular inflection point, including population size of the inflection point under limiting values of the model parameters, and parameter values that lead to inflection point located at the mean value of the curve. Section 5 presents the features of the logarithmic growth model: Convergence and boundedness properties; determination of the logarithmic inflection point and conditions for its existence; effect of model parameters on the existence of the inflection point and the lag phase; and properties of the logarithmic inflection point, including population size of the logarithmic inflection point under limiting values of the model parameters, and parameter values that lead to inflection point located at the mean value of the curve. In section 6, an approach is proposed for determining the lag phase duration, for each the regular and the logarithmic growth curves. In section 7 the model is applied to experimental data, what includes the arrangement of the model and its properties in terms of the logarithmic measure, and a numerical fitting example. Finally, in section 8 the conclusions are drawn.
In microbial culture studies, a typical monotonic growth curve comprises the following phases: Lag phase, exponential growth phase, and stationary phase, although the lag phase may not appear. Thus, monotonic growth does not comprise decay phase or oscillations [17,21,22]. The lag phase involves the physiological adaptation of microbial cells to the new environment, resulting in low or negligible growth rate. The exponential growth phase involves exponential reproduction of microorganisms. The stationary phase comprises constant population of microorganisms [16,22,23,24]. The maximum growth rate occurs in the exponential growth phase, at the so-called inflection point [17,24,25].
In growth model studies for microbial count data (cfu/mL), the logarithmic measure defined as y=ln(N/Nto) or y=ln(N), is used as dependent variable for curve fitting and determination of the inflection point and lag phase duration [2,26,27].
The inflection point is characterized by the fact that the derivative of the dependent variable is maximum (the growth rate is maximum) [19,21]. The existence of the inflection point requires a convex shape in the growth rate versus population curve [17,25]. In turn, the existence of the lag phase requires the existence of the inflection point, and it is possible that the inflection point exists while the lag phase does not. In autonomous growth models the lag phase is generated if the initial population size (Nto) is sufficiently close to the unstable equilibrium point, thus yielding small value of the initial growth rate (dN/dt|t=to). In this case, the escaping of the population from the neighborhood of the unstable equilibrium point takes more time [17,21].
The lag phase duration can be quantified via the classical approach, or based on the third derivative (maximum acceleration of the growth rate) or the maximum curvature [18,24]. The classical approach is based on the intersection of the tangent line at the inflection point with the horizontal line that meets the initial point [17,24].
The model of Vadasz and Vadasz [17,21] comprises a quadratic term, which defines a stable and an unstable equilibrium point. Certain values of the model parameters, corresponding to a specific region, render the model capable of representing the inflection point and the lag phase. For this parameter region, the value of the stable equilibrium point is higher than the unstable one, both being positive. The stable equilibrium point is the carrying capacity. The state variable behaves monotonically, escaping from the unstable equilibrium point towards the stable one. In addition, the model showed a high capacity to represent the lag phase and the inflection point for experimental data.
The generalized logistic model proposed by Tsoularis and Wallace [19] is:
dNdt=rNα[1−(NNmax)β]γ | (2.1) |
Where r, α, β, γ, Nmax are positive constants, N is the population size, and Nmax is the limit value, usually called carrying capacity. This model combines the function Nα, α>0, with the logistic type model.
dNdt=r[1−(NNmax)β]γ | (2.2) |
Therein, model simulations with α=1.5 and α=3 indicated their capability for representing the lag phase. In addition, several growth models that are capable of representing the lag phase can be expressed in this form, namely: The Korf’s form, the hyperbolic form for regenerative growth, the Smith’s equation and the Richard model [19,25]. Moreover, the Smith’s equation has achieved satisfactory representation of the lag phase for experimental data [9].
The Putter-Bertalanffy model analyzed by [20], can be expressed as
dNdt=rN1−ar[1−(NNmax)γ] | (2.3) |
where r, Nmax and γ are positive constants, and N is the population size. It is a Savageau's model, which is in turn a particular case of the generalized logistic function analyzed by [19]. The possible ranges of (1−ar), either positive, negative or zero, lead to different possibilities of the existence of the inflection point and the lag phase.
The model of Fujikawa, Kai and Morozumi is [16,28]:
dNdt=rN(1−NNmax)(1−NminN)c | (2.4) |
where N is the population size, r, c, Nmin, Nmax are positive constants, Nmax>Nmin and Nto∈(Nmin,Nmax). This model was constructed by adding the function fc=(1−Nmin/N)c to the logistic model
dNdt=rN(1−NNmax) |
The function fc strongly suppresses the growth rate during the lag phase, whereas the term (1−N/Nmax) strongly suppresses the growth rate during the stationary phase. The fc function is not mechanistically inspired. Fujikawa, Kai and Morozumi fitted the new model and the traditional model of Baranyi and Roberts [29] for data of E. coli culture. The tlag of the new model was similar to that of the Baranyi’s model, although slightly less ascertained. In contrast, the slope (k) was more ascertained than that of Baranyi’s model. In the fitting, it was found that a small (Nto−Nmin) value is necessary for the estimated c being less dependent on initial population. To this end, the value Nmin=(1−1×10−6)Nto was used, yielding appropriate representation of the lag phase
We propose the following modification of the Fujikawa’s model:
dNdt=rNα(1−NNmax)γ(1−NminN)c | (3.1) |
where r, α, γ, c, Nmin, Nmax are constants, and r>0, Nmax>Nmin>0; N is the population size, and Nmax is the carrying capacity. In turn, the model (3.1) combines the term fc=(1−Nmin/N)c with the hyperlogistic function:
dNdt=rNα(1−NNmax)γ | (3.2) |
The difference of model (3.1) with the Fujikawa’s model (2.4) is the exponents α and γ. Indeed, Fujikawa’s model can be regarded as a simplified variant corresponding to α=1, γ=1. The model coefficients are assumed constant in order to facilitate the study of model properties. Despite this assumption, the model parameters can be defined as function of growth variables in further studies.
Remark 3.1. The term fc=(1−Nmin/N)c with c>0 slows the growth rate of the population size N during the lag phase, provided small positive (Nto–Nmin). In turn, this is related to the fact that dN/dt|N=Nto is proportional to (Nto–Nmin)c. In contrast, a c<0 value and a positive small (Nto–Nmin) value imply: i) dN/dt|N=Nto values are high, so that the initial time period of the N versus t curve features large slope; ii) (1/N)(dN/dt)|N=Nto values are large, so that the initial time period of the y versus t curve features large slope.
Remark 3.2. In the generalized logistic form of Tsoularis and Wallace [19], Eq (2.1), and in the hyperlogistic form (3.2), the suppression of the growth rate during the lag phase is due to the term Nα, with α>0, whereas in the proposed model (3.1), it is due to the term fc with c>0.
Remark 3.3. The c=0 value implies absence of the fc term in model (3.1), so that it reduces to the hyperlogistic model (3.2), which was studied by Tsoularis & Wallace [19].
Remark 3.4. The shape of the growth curve and the existence of the lag phase are related to the equilibrium points of the system. The term Nα, α>0 generates the equilibrium point Neq=0, whereas the term (1−N/Nmax)γ, γ>0 yields Neq=Nmax; and the term (1−Nmin/N)c, c>0 yields Neq=Nmin. For γ>0, c>0 and Nto>Nmin, the N versus t curve escapes from the equilibrium Neq=Nmin and approaches equilibrium Neq=Nmax. The lag phase occurring at N≈Nmin is related to a trajectory that begins close to Nmin, so that it takes more time near Nmin. Therefore, the existence of the lag phase occurring at N≈Nmin requires: i) the existence of the equilibrium point Neq=Nmin, what in turn implies c>0; ii) a Nto value such that Nto−Nmin is small positive.
At what follows, we state the definitions of the growth rate curve related to the mean between lower and higher populations. The mean between Nmin and Nmax is:
Nmn=Nmin+Nmax2 | (3.3) |
The logarithmic population y, and its minimum and maximum values are defined as:
y=ln(NNto) | (3.4) |
ymin=ln(NminNto) | (3.5) |
ymax=ln(NmaxNto) | (3.6) |
As a consequence of the above definitions, y|t=to=0. The mean between ymin and ymax is:
ym=ymin+ymax2 | (3.7) |
Substituting the definition of ymin and ymax, that is, Eqs (3.5) and (3.6), into ym, gives:
ym=ln(√NminNmaxNto) | (3.8) |
To obtain the value of N corresponding to y=ym, we use the definition y=ln(N/Nto), what gives:
Nmy=√NminNmax | (3.9) |
The adaptation function fc=(1−Nmin/N)c, c>0, has the following properties:
i)fc|N=Nmin=0;fc|N=Nmax=(1−NminNmax)c<1;fc∈(0,1)forN∈(Nmin,Nmax] | (3.10) |
ii)dfcdN=cNminN−2(1−NminN)c−1,hencedfc/dN>0forN∈(Nmin,Nmax],iii)dfcdc=ln(1−Nmin/N)eln(1−Nmin/N)c,so thatdfc/dc<0forN∈(Nmin,Nmax] | (3.11) |
iv)ddc(fc|N=Nmax)<0,so thatascincreases,fc|N=Nmaxdecreases | (3.12) |
Remark 3.5. Property ii implies that fc is monotonically increasing with respect to N, for N∈(Nmin,Nmax]
Remark 3.6. Properties ii, iii and iv imply that fc is lower for high c values and for small positive values of (N−Nmin). As a consequence, the dN/dt value is suppressed for small positive values of (Nto−Nmin).
Property iv follows from:
ddc(fc|N=Nmax)=ln(1−NminNmax)eln(1−Nmin/Nmax)c | (3.13) |
Moreover, fc exhibits an inflection point with the following features:
i)Ifc>1,thenfchas an inflection point atN⋆=(c+1)Nmin/2,so thatN⋆>Nmin | (3.14) |
ii)fc|N=N⋆=(c−1c+1)c,c>1 | (3.15) |
iii)dfcdN|N=N⋆=4c(c+1)(c−1)Nmin(c−1c+1)c,c>1 | (3.16) |
Remark 3.7 The inflection point of fc leads to a steep change of the fc versus N curve; and consequently a steep change of the dN/dt versus t curve for c>1.
Figure 1 shows the fc versus N curve for Nmin=(1−1×10−6)Nto, Nto=1 and several c values.
Proposition 4.1. Boundedness and convergence. Consider the growth model (3.1) with γ>0, real values of α and c, and Nmax>Nmin>0. If Nto∈(Nmin,Nmax), the population size N exhibits the following properties: i) N increases monotonically from N=Nto to Nmax, so that limt→∞N=Nmax; ii) N∈[Nto,Nmax)∀t≥to.
Proof. The above proposition follows from the values of the growth rate dN/dt: ⅰ) dN/dt=0 for N=Nmax; ⅱ) dN/dt>0 for N∈(Nmin,Nmax).
In this study, the regular inflection point (t#, N#) is defined as the point of the N versus t curve where dN/dt is maximum, being dN/dt|N=N# the maximum growth rate.
Proposition 4.2. Consider the model (3.1) with γ>0; real values of α and c; Nmax>Nmin>0; and Nto∈(Nmin,Nmax). With these conditions the regular inflection point population, N#, is provided by:
αN#−γNmax−N#+cNminN#(N#−Nmin)=0 | (4.1) |
Proof. As the regular inflection point corresponds to a maximum of the dN/dt versus t curve, then
d2N/dt2|N=N#=0,dN/dt|N=N#>0 | (4.2) |
The condition Nto∈(Nmin,Nmax) implies that N∈(Nmin,Nmax) ∀t≥to, according to Proposition 4.1. The condition dN/dt|N=N#>0 implies that N#≠Nmax. The above two results imply N#∈(Nmin,Nmax). This result and Eq (4.2) give
ddN(dNdt)|N=N#=0,N#∈(Nmin,Nmax) | (4.3) |
finally, applying condition (4.3) for the growth model (3.1), yields Eq (4.1).
Remark 4.1. The condition (4.3) implies that N# corresponds to a maximum of the dN/dt versus N curve.
Proposition 4.3. Consider growth model (3.1) with real values of α, c and γ, Nmax>Nmin>0 and Nto∈(Nmin,Nmax). If N#<Nto, then N>N#∀t≥to; that is, the N versus t curve exhibits no inflection point. If N#=Nto, then N=N# for t=to; that is, the regular inflection point occurs at t=to, N=Nto.
Proof. Since the N versus t curve is monotonic and N≥Nto ∀t≥to, then N>N#∀t≥to if N#<Nto.
Remark 4.2. Consider growth model (3.1), with γ>0 and real values of α and c, and Nmax>Nmin>0. The existence of the inflection point in the N versus t curve requires: i) the existence of a maximum in the dN/dt versus N curve, which is given by Eq (4.1), subject to N#∈(Nmin,Nmax), and ii) that Nto<N# where Nto>Nmin>0. The existence of the lag phase requires the existence of the inflection point and a small value of dN/dt|t=to. In turn, this implies a positive c value and a small positive (Nto−Nmin) value.
Proposition 4.4. Consider growth model (3.1) with real values of α, c, γ, and Nmax>Nmin>0 and Nto∈(Nmin,Nmax):
i) if
γ>0,andc∈(0,∞) | (4.4) |
holds, then the dN/dt versus N curve exhibits a maximum at N=N#, N#∈(Nmin,Nmax).
ii) the inflection point population, N#, satisfies:
αN−γNmax−N+cNminN(N−Nmin)>0forNmin<N<N#αN−γNmax−N+cNminN(N−Nmin)<0forN#<N≤Nmax | (4.5) |
αN−γNmax−N+cNminN(N−Nmin)=0forN=N# | (4.6) |
iii) if either
i)γ>0,α∈(0,∞)andc∈(−∞, −α(Nmax/Nmin−1)]orii)γ>0,α∈(−∞,0],andc∈(−∞,0] | (4.7) |
holds, the dN/dt versus t curve is monotonic decreasing for N∈(Nmin,Nmax).
Proof. Model (3.1) with Nmax>Nmin>0 subject to parameter regime (4.4) gives:
dNdt|N=Nmin=0;dNdt|N=Nmax=0;dNdt|N∈(Nmin,Nmax)>0 | (4.8) |
Therefore, the dN/dt versus N curve exhibits a maximum at N=N#, with N#∈(Nmin,Nmax).
In addition, we notice from Eq (4.1) that the inflection point population N# is provided by D2|N=N#=0, where
D2=αN+cNminN(N−Nmin)−γNmax−N | (4.9) |
The D2 term exhibits vertical asymptotes at N=Nmin and at N=Nmax. In addition, if parameter regime (4.4) holds, then
limN→N−maxD2=−∞ | (4.10) |
limN→N+minD2=+∞ | (4.11) |
so that there is at least one N value that belongs to range (Nmin,Nmax) and leads to D2=0. Then, there exists a maximum in the dN/dt versus N curve, at N=N#, N# provided by Eq (4.6), and N#∈(Nmin,Nmax).
For the second part, we recall that for model (3.1) with Nto∈(Nmin,Nmax), the inflection point of the N versus t curve satisfies Eq (4.2). As the inflection point corresponds to a maximum of the dN/dt versus t curve, then
ddt(dNdt)>0fort<t#ddt(dNdt)<0fort>t#ddt(dNdt)=0fort=t# | (4.12) |
according to [30]. Equivalently,
(dNdt)⋅ddN(dNdt)>0fort<t#(dNdt)⋅ddN(dNdt)<0fort>t#(dNdt)⋅ddN(dNdt)=0fort=t# | (4.13) |
As dN/dt>0 for N∈(Nmin,Nmax), real values of α, c, γ and Nmax>Nmin>0, then the conditions t<t#, t=t# and t>t# are equivalent to N<N#, N=N# and N>N#, respectively, for N∈(Nmin,Nmax). Therefore, expressions (4.13) lead to
ddN(dNdt)>0forNmin<N<N#ddN(dNdt)<0forN#<N<NmaxddN(dNdt)=0forN=N# | (4.14) |
For model (3.1), and N∈(Nmin,Nmax), we have
ddN(dNdt)=rNα(1−NNmax)γ(1−NminN)c(αN−γNmax−N+cNminN(N−Nmin)) | (4.15) |
This and expressions (4.14) lead to
αN−γNmax−N+cNminN(N−Nmin)>0forNmin<N<N#αN−γNmax−N+cNminN(N−Nmin)<0forN#<N<NmaxαN−γNmax−N+cNminN(N−Nmin)=0forN=N# | (4.16) |
This and condition N∈(Nmin,Nmax) result in second part of Proposition.
For the third part, we notice that if γ>0, and (4.7) holds, then
αN−γNmax−N+cNminN(N−Nmin)<0forN∈(Nmin,Nmax) | (4.17) |
This and Eq (4.15) lead to
ddN(dNdt)<0forN∈(Nmin,Nmax) | (4.18) |
As dN/dt>0 for N∈(Nmin,Nmax), then
(dNdt).ddN(dNdt)<0forN∈(Nmin,Nmax) | (4.19) |
thus d2N/dt2<0 for N∈(Nmin,Nmax). Equivalently, the dN/dt versus t curve is decreasing.
Remark 4.3. In first part of Proposition 4.4 not only the values belonging to the parameter regime (4.4) lead to existence of a maximum in the dN/dt versus N curve. Indeed, certain negative c values with γ>0 also lead to a maximum. The third part of Proposition 4.4 corresponds to the parameter conditions for the absence of a maximum in the dN/dt versus N curve, hence absence of the regular inflection point.
At what follows, the effect of model parameters on the existence of the inflection point and lag phase is assessed, on the basis of the growth rate versus population size. The effect of model parameters on the dN/dt vs N curve and the N vs t growth curve for the model (3.1) is illustrated in Figures 2–4. The following parameter values are used: Nmax=100, γ=1, Nmin=(1−1×10−6)Nto, Nto=1. As the advantage of the proposed model is its capacity to represent the lag phase, it is assumed that γ is positive, and its effect is not assessed. There are four different cases of the α, c values, yielding different possibilities of the existence of the lag phase and the inflection point:
● α>0, c>0: Pⅰ) dN/dt=0 for N=Nmin; Pⅱ) dN/dt=0 for N=Nmax; Pⅲ) dN/dt>0 for N∈(Nmin,Nmax). The first part of Proposition 4.4 applies, so that the inflection point occurs provided Nto<N# and the lag phase existence depends on the (Nto−Nmin) value. Simulations are shown in Figure 2.
● α>0 and c<0: Pⅰ) dN/dt=∞ for N=Nmin, that is, there is a vertical asymptote at N=Nmin; Pⅱ) dN/dt=0 for N=Nmax; Pⅲ) dN/dt>0 for N∈(Nmin,Nmax). Property Pi implies that dN/dt|N=Nto is overly high for low values of (Nto−Nmin), so that there is no lag phase. The occurrence of the inflection point is determined by Eq (4.6) subject to N#∈(Nmin,Nmax), provided Nto<N#. The third part of Proposition 4.4 indicates that for c≤−α(Nmax/Nmin−1) values, there is no inflection point. Moreover, simulations show that for some parameter values there is a local minimum and a local maximum in the range N∈(Nmin,Nmax). As a consequence, there may be inflection point but there is no lag phase. Simulation shown in Figure 3 confirm the above facts, and also shows that the inflection point is owed to the term Nα.
● α<0 and c>0: The properties Pi to Piv of case α>0, c>0 are fulfilled. The first part of Proposition 4.4 applies, so that the inflection point occurs provided Nto<N#, and the lag phase existence depends on the (Nto−Nmin) value. Simulation shown in Figure 4 confirms the above facts, and also show that: ⅰ) the effect of negative values of α on the lag phase duration is low; ⅱ) the maximum growth rate is overly low, so that the slope of the population at the exponential growth phase is overly low. This confirms the importance of the Nα term for the slope of the exponential growth phase.
● α<0 and c<0: Pⅰ) dN/dt=∞ for N=Nmin, that is, there is a vertical asymptote at N=Nmin; Pⅱ) dN/dt=0 for N=Nmax; Pⅲ) dN/dt>0 for N∈(Nmin,Nmax). Proposition 4.4 applies, so that the dN/dt versus N curve is monotonically decreasing for N∈(Nmin,Nmax). Thus, no inflection point occurs, and therefore no lag phase occurs. Simulations in Figure 2 confirm the above properties.
Intermediate cases are:
● α>0 and c=0: The term fc is absent, so that the model becomes dN/dt=rNα(1−N/Nmax)γ, and: Pⅰ) dN/dt>0 for N=Nmin; Pⅱ) dN/dt=0 for N=Nmax; Pⅲ) dN/dt>0 for N∈(Nmin,Nmax); Piv) dN/dt|N=Nto is positive and it is proportional to Nαto. Thus, the dN/dt versus N curve exhibits a maximum with N#∈(Nmin,Nmax), so that there is an inflection point provided Nto<N#, and it is owed to the term Nα with α>0. For general (not small) values of Nto, dN/dt|t=to is not small, so that there is no lag phase. Simulations are shown in Figure 3.
● α=0 and c<0: The term Nα is absent, and the properties are similar to those of case α<0 and c<0: Pⅰ) dN/dt=∞ for N=Nmin, that is, there is a vertical asymptote at N=Nmin; Pⅱ) dN/dt=0 for N=Nmax; Pⅲ) dN/dt>0 for N∈(Nmin,Nmax). According to Proposition 4.4, the dN/dt versus N curve exhibits monotonic decreasing behavior, so that no inflection point occurs, and therefore lag phase neither occurs, what is confirmed by simulation in Figure 3.
● α=0 and c>0: The term Nα is absent, and the properties Pi to Piv of case α<0 and c>0 are fulfilled. The first part of Proposition 4.4 applies, so that an inflection point occurs provided Nto<N# and the lag phase existence depends on the (Nto−Nmin) value. Although simulations are not shown, they indicate that the maximum occurs at low values of the population size (N# is low), and the maximum growth rate is overly low, so that the slope of the growth curve at the exponential growth phase is overly low.
Remark 4.4. The above cases of α and c values indicate that the c parameter has a significant influence on the lag phase duration but the influence of the α value is low. Concerning the slope at the exponential growth phase, the α value has a significant influence, whereas the influence of the c value is low.
Remark 4.5. The parameter region α<0 leads to low slope of the population at the exponential growth phase, whereas c<0 implies that no lag phase occurs. Thus, the condition c>0 is required for representing the lag phase, whereas α>0 is required for higher slopes of the population at the exponential growth phase. This is in accordance with the fact that the formulation of the adaptation function fc by Fujikawa, Kai and Morozumi [16] aimed at suppressing the growth rate during the lag phase, but without affecting the rest of the population growth behavior.
Proposition 4.5. Consider the model (3.1) with γ>0. The population at the inflection point, N#, is defined by:
ˉax2+ˉbx+ˉc=0,x=N#Nmin | (4.20) |
where
ˉa=α+γ;ˉb=(−1)(α(1+NmaxNmin)+γ−c);ˉc=(α−c)NmaxNmin | (4.21) |
And the discriminant is
ˉb2–4ˉaˉc=(α(1−NmaxNmin)+γ−c)2+4γcNmaxNmin | (4.22) |
so that no imaginary solutions are possible for γc≥0. If ˉa>0, c>0, and ˉb2–4ˉaˉc is non-negative, only the positive square root can lead to N#>Nmin.
Proof. The population at the regular inflection point is obtained by applying Eq (4.3) to the growth model (3.1), leading to
α(Nmax−N)(N–Nmin)−γN(N–Nmin)+cNmin(Nmax−N))|N=N#=0 | (4.23) |
Rearranging yields the quadratic Eq (4.20). In order to determine whether the negative square root term can lead to N#>Nmin for ˉa>0, we use Proposition presented in Appendix B. Using the definition of ˉa, ˉb, ˉc, we have: ˉc+ˉa+ˉb=(−1)c(Nmax/Nmin−1). Since (Nmax/Nmin)>1, then ˉc+ˉa+ˉb<0 for c>0, so that the condition ˉc+ˉa+ˉb>0 is not fulfilled for c>0. Consequently, only the positive square root term leads to N#/Nmin>1. End of the proof.
Proposition 4.6. Limit values of N#. Consider model (3.1), with α>0, γ>0, c≥0. The limit values of the regular inflection point, N#, are:
i)limγ→0N#=Nmaxii)limγ→∞N#=Nminiii)limc→0N#=Nmin2(α+γ)[α(1+NmaxNmin)+γ+|(−1)α(NmaxNmin−1)+γ|]iv)limc→∞N#=Nmaxv)limα→0N#=Nmin2γ[γ−c+√(c−γ)2+4cγNmaxNmin]vi)limα→∞N#=Nmax | (4.24) |
Proof. The above limit values are obtained by substituting γ→0, γ→∞, c→0, c→∞, α→0, α→∞ in Eq (4.23), and solving for N#.
Remark 4.6. The above proposition indicates that the limit values of N# for α→0 depends on the values of c and γ; whereas the limit value of N# for c→0 depends on the values of α, γ.
Proposition 4.7. Conditions for N#=Nto. Consider the model (3.1). The population at the regular inflection point, N#, satisfies: N#=Nto for
cγ+α(NmaxNmin−NtoNmin)=(NtoNmin−1)(NtoNmin−αγ+αNmaxNmin) | (4.25) |
Proof. Using N#=Nto in Eq (4.23) and arranging yields Eq (4.25).
Remark 4.7. The above proposition implies that in the N versus t curve, with Nto∈(Nmin,Nmax), the maximum (dN/dt) occurs at t=to, N=Nto, if the α, γ, c values are such that Eq (4.25) holds true.
Proposition 4.8. Conditions for N#=Nmn. Consider model (3.1). The population at the regular inflection point, N#, satisfies N#=Nmn for
cγ+α=12(1+γ−αγ+αNmaxNmin) | (4.26) |
Proof. Using N#=Nmn in Eq (4.23) and arranging yields Eq (4.26).
Remark 4.8. The above proposition implies that in the N versus t curve, with Nto∈(Nmin,Nmn), the maximum dN/dt occurs at N=Nmn where Nmn=(Nmin+Nmax)/2, if the α, γ, c values are such that Eq (4.26) holds true.
Remark 4.9. The determination of the lag time of the regular growth curve via the classical approach requires the determination of the inflection point, as shown in Appendix A. To this end, the population at the inflection point (N#) can be computed using expressions (4.20), (4.21), and the inflection time (t#) must be computed numerically.
Proposition 5.1. Boundedness and convergence. Consider growth model (3.1) with γ>0 and real values of α, c. If Nto∈(Nmin,Nmax), equivalently yto∈(ymin,ymax), then ymin<0<ymax, yto∈(ymin,ymax) and the logarithmic population y exhibits the following features: i) y increases monotonically from y=yto to y=ymax, so that limyt→∞=ymax; ii) y∈[yto,ymax)∀t≥to.
Proof. For y=ln(N/Nto) and the growth model (3.1), with γ>0 and real values of α, c, the specific growth rate dy/dt=(1/N)(dN/dt) fulfills the following properties: ⅰ) (dy/dt)=0 for N=Nmax; ⅱ) (dy/dt)>0 for N∈(Nmin,Nmax).
In this study, the logarithmic inflection point (t∗, N∗, y∗) is defined as the point of the y versus t curve where dy/dt is maximum, being [dy/dt]y∗ the maximum specific growth rate; and N∗ the value of N corresponding to y∗, ie N∗=Ntoey∗.
Proposition 5.2. Consider the growth model (3.1) with γ>0; real values of α, c; Nmax>Nmin>0; and Nto∈(Nmin,Nmax), equivalently yto∈(ymin,ymax): The logarithmic inflection point population, N∗, is provided by
α−1N∗−γNmax−N∗+cNminN∗(N∗−Nmin)=0 | (5.1) |
Proof. As the logarithmic inflection point corresponds to a maximum of the dy/dt=(1/N)(dN/dt) versus t curve, then
d2y/dt2|N=N∗=0,dy/dt|N=N∗>0. | (5.2) |
The condition Nto∈(Nmin,Nmax) implies that N∈(Nmin,Nmax)∀t≥t0, equivalently y∈(ymin,ymax) ∀t≥t0, according to Proposition 5.1. The condition dy/dt|N=N∗>0 implies that N∗≠Nmax (equivalently y∗≠ymax). The above two results imply that N∗∈(Nmin,Nmax), equivalently y∗∈(ymin,ymax). This result, Eq (5.2) and using y=ln(N/Nto), yield
ddN(1NdNdt)|N=N∗=0;N∗∈(Nmin,Nmax) | (5.3) |
Applying this condition to the growth model (3.1), yields Eq (5.1).
Remark 5.1. The condition (5.3) implies that N∗ provided by Eq (5.1) corresponds to a maximum of the (1/N)(dN/dt) versus N curve [21].
Proposition 5.3. Consider the logarithmic measure y=ln(N/Nto), being N provided by growth model (3.1) with real values of α, c, γ, and Nmax>Nmin>0; and Nto∈(Nmin,Nmax). If N∗<Nto, (equivalently y∗<yto) then N>N∗∀t≥to; that is, the y versus t curve exhibits no inflection point. If N∗=Nto, (equivalently y∗=yto), then N=N∗ for t=to; that is, the logarithmic inflection point occurs at t=to, N=Nto.
Proof. Recall that the N versus t curve is monotonic and N≥Nto∀t≥to. Then, N>N∗∀t≥to for N∗<Nto. Consequently, y>y∗∀t≥to for y∗<yto.
Remark 5.2. Consider the logarithmic measure y=ln(N/Nto) provided by growth model (3.1), with γ>0, real values of α, c, and Nmax>Nmin>0, and Nto∈(Nmin,Nmax). The existence of the inflection point of the y versus t curve requires: i) the existence of a maximum in the (1/N)(dN/dt) versus N curve which is given by Eq (5.1) subject to N∗∈(Nmin,Nmax), and ii) that Nto<N∗ where Nto>Nmin>0. The existence of the lag phase requires the existence of the inflection point and a small value of (1/N)(dN/dt)|t=to. In turn, this implies a positive c value and a small positive (Nto−Nmin) value.
Proposition 5.4. Consider the logarithmic measure y=ln(N/Nto), being N provided by growth model (3.1) with real values of α, c, γ, and Nmax>Nmin>0; and Nto∈(Nmin,Nmax):
i) If
γ>0,andc∈(0,∞) | (5.4) |
holds, then the (1/N)(dN/dt) versus N curve exhibits a maximum at N=N∗, N∗∈(Nmin,Nmax).
ii) The logarithmic inflection point population N∗, satisfies:
α−1N−γNmax−N+cNminN(N−Nmin)>0forNmin<N<N∗α−1N−γNmax−N+cNminN(N−Nmin)<0forN∗<N≤Nmax | (5.5) |
α−1N−γNmax−N+cNminN(N−Nmin)=0forN=N∗ | (5.6) |
iii) If either
γ>0,α>1,andc∈(−∞, −(α−1)(Nmax/Nmin−1)]orγ>0,α∈(−∞,1]andc∈(−∞,0] | (5.7) |
then (1/N)(dN/dt) versus t curve is monotonic decreasing for N∈(Nmin,Nmax).
Proof. Model (3.1) with Nmax>Nmin>0 subject to parameter regime (5.4) gives:
1NdNdt|N=Nmin=0;1NdNdt|N=Nmax=0;1NdNdt|N∈(Nmin,Nmax)>0 | (5.8) |
Therefore, the (1/N)(dN/dt) versus N curve exhibits a maximum at N=N∗, with N∗∈(Nmin,Nmax). In addition, we notice from Eq (5.1) that the inflection point population N∗ is provided by D2|N=N∗=0, where
D2=(α−1)N+cNminN(N−Nmin)−γNmax−N | (5.9) |
The D2 term exhibits vertical asymptotes at N=Nmin and at N=Nmax. In addition, if parameter regime (5.4) holds, then
limN→N−maxD2=−∞ | (5.10) |
limN→N+minD2=+∞ | (5.11) |
so that there is at least one N value that belongs to range (Nmin,Nmax) and leads to D2=0. Then, there exists a maximum in the (1/N)(dN/dt) versus N curve, at N=N∗, N∗ provided by Eq (5.1), and N∗∈(Nmin,Nmax).
For the second part, we recall that for the model (3.1) with Nto∈(Nmin,Nmax), the inflection point of the y versus t curve satisfies Eq (5.2). As the inflection point corresponds to a maximum of the dy/dt versus t curve, then
ddtdydt>0fort<t∗ddtdydt<0fort>t∗ddtdydt=0fort=t∗ | (5.12) |
according to [30]. As y=ln(N/Nto), then
d2ydt2=(dNdt).ddN(1NdNdt) | (5.13) |
Combining this with Eq (5.12) yields:
dNdt⋅ddN(1NdNdt)>0fort<t∗dNdt⋅ddN(1NdNdt)<0fort>t∗dNdt⋅ddN(1NdNdt)=0fort=t∗ | (5.14) |
As (dN/dt)>0 for N∈(Nmin,Nmax), real values of α, c, γ and Nmax>Nmin>0, then the conditions t<t∗, t=t∗ and t>t∗ are equivalent to N<N∗, N=N∗ and N>N∗, respectively, for N∈(Nmin,Nmax). Therefore, expression (5.14) leads to
ddN(1NdNdt)>0forNmin<N<N#ddN(1NdNdt)<0forN#<N<NmaxddN(1NdNdt)=0forN=N# | (5.15) |
For model (3.1), and N∈(Nmin,Nmax), we have
ddN(1NdNdt)=rNα−1(1−NNmax)γ(1−NminN)c[α−1N−γNmax−N+cNminN(N−Nmin)] | (5.16) |
This and expressions (5.15) lead to
α−1N−γNmax−N+cNminN(N−Nmin)>0forNmin<N<N∗α−1N−γNmax−N+cNminN(N−Nmin)<0forN∗<N<Nmaxα−1N−γNmax−N+cNminN(N−Nmin)=0forN=N∗ | (5.17) |
This and condition N∈(Nmin,Nmax) results in second part of Proposition.
For the third part, we notice that if (5.7) holds, then
α−1N−γNmax−N+cNminN(N−Nmin)<0forN∈(Nmin,Nmax) | (5.18) |
This and Eq (5.16) lead to
ddN(1NdNdt)<0forN∈(Nmin,Nmax) | (5.19) |
As (dN/dt)>0 for N∈(Nmin,Nmax), then
(dNdt).ddN(1NdNdt)<0forN∈(Nmin,Nmax) | (5.20) |
thus d2y/dt2<0 for N∈(Nmin,Nmax). Equivalently, the dy/dt versus t curve is decreasing.
Remark 5.3. In first part of Proposition 5.4, not only the parameter values belonging to parameter regime (5.4) lead to existence of a maximum in the (1/N)(dN/dt) versus N curve. Indeed, certain negative c values with γ>0 also lead to a maximum. The third part of Proposition 5.4 corresponds to absence of a maximum in the (1/N)(dN/dt) versus N curve, therefore absence of the logarithmic inflection point.
At what follows, the effect of model parameters on the existence of the inflection point and lag phase is assessed, on the basis of the specific growth rate versus population size curve.
The effect of model parameters on the (1/N)(dN/dt) vs N curve and the y vs t curve is illustrated in Figures 5–7. The following parameter values are used: Nmax=100, γ=1, Nmin=(1−1×10−6)Nto, Nto=1. As the advantage of the proposed model is its capacity to represent the lag phase, it is assumed that γ is positive and its effect is not assessed. There are four different cases of the α, c values, yielding different possibilities of the existence of the lag phase and the inflection point in the y versus t curve:
● α>1, c>0: Pⅰ) (1/N)(dN/dt)=0 for N=Nmin; Pⅱ) (1/N)(dN/dt)=0 for N=Nmax; Pⅲ) (1/N)(dN/dt)>0 for N∈(Nmin,Nmax). The first part of Proposition 5.4 applies, so that the (1/N)(dN/dt) versus N curve exhibits a maximum at N=N∗, N∗∈(Nmin,Nmax). Thus, the inflection point occurs provided Nto<N∗, and the lag phase existence depends on the (Nto−Nmin) value. Simulation is shown in Figure 5.
● α<1 and c>0: The properties Pi to Piv of case α>1, c>0 are fulfilled. The first part of Proposition 5.4 applies, so that the inflection point occurs provided Nto<N∗, and the lag phase existence depends on the (Nto−Nmin) value. Simulation is shown in Figure 6, confirming the above facts, and also shows that: ⅰ) the effect of α<1 values on the lag phase duration is low; ⅱ) the maximum growth rate is overly low, leading to low values of the slope of the population at the exponential growth phase. This confirms the importance of the Nα term for the slope of the exponential growth phase.
● α>1 and c<0: Pⅰ) (1/N)(dN/dt)=∞ for N=Nmin, that is, there is a vertical asymptote at N=Nmin; Pⅱ) (1/N)(dN/dt)=0 for N=Nmax; Pⅲ) (1/N)(dN/dt)>0 for N∈(Nmin,Nmax). Property Pi implies that (1/N)(dN/dt)|N=Nto is overly high for low values of (Nto−Nmin), so that there is no lag phase. The occurrence of inflection point is determined by Eq (5.6) subject to N∗∈(Nmin,Nmax), and provided Nto<N∗. The third part of Proposition 5.4 indicates that for c≤(−1)(α−1)(Nmax/Nmin−1) values there is no inflection point. Moreover, simulations show that for some parameter values there is a local minimum and a local maximum in the range N∈(Nmin,Nmax). As a consequence, there may be inflection point but there is no lag phase. Simulation in Figure 7 confirms the above facts.
● α<1 and c<0: Pⅰ) (1/N)(dN/dt)=∞ for N=Nmin, that is, there is a vertical asymptote at N=Nmin; Pⅱ) (1/N)(dN/dt)=0 for N=Nmax; Pⅲ) (1/N)(dN/dt)>0 for N∈(Nmin,Nmax). Proposition 5.4 applies, so that the (1/N)(dN/dt) versus N curve is monotonically decreasing for N∈(Nmin,Nmax). Thus, no inflection point occurs, and therefore lag phase neither occurs. Simulation in Figures 5 and 7 confirms the above facts.
Intermediate cases are:
● α>1 and c=0: The term fc is absent, so that the model becomes (1/N)(dN/dt)=rNα−1(1−N/Nmax)γ, and: Pⅰ) (1/N)(dN/dt)>0 for N=Nmin; Pⅱ) (1/N)(dN/dt)=0 for N=Nmax; Pⅲ) (1/N)(dN/dt)>0 for N∈(Nmin,Nmax); Piv) (1/N)(dN/dt)|N=Nto is positive and is proportional to Nα−1to. Thus, the (1/N)(dN/dt) versus N curve exhibits a maximum with N∗∈(Nmin,Nmax), so that the inflection point occurs provided Nto<N∗ and it is owed to the term Nα−1 with α>1. For general (not small) values of Nto, (1/N)(dN/dt)|t=to is not small, so that there is no lag phase. Simulations are shown in Figure 7.
● α=1 and c<0: The term Nα−1 is absent, and the properties are similar to those of case α<1 and c<0: Pⅰ) (1/N)(dN/dt)=∞ for N=Nmin, that is, there is a vertical asymptote at N=Nmin; Pⅱ) (1/N)(dN/dt)=0 for N=Nmax; Pⅲ) (1/N)(dN/dt)>0 for N∈(Nmin,Nmax). According to the third part of Proposition 5.4, the (1/N)(dN/dt) versus N curve exhibits monotonic decreasing behavior, so that no inflection point occurs, and therefore lag phase neither occurs.
● α=1 and c>0: The term Nα−1 is absent, and the properties Pi to Piv of case α<1 and c>0, are fulfilled. The first part of Proposition 5.4 applies, so that the inflection point occurs provided Nto<N∗; and the lag phase existence depends on the (Nto−Nmin) value. Simulation in Figure 5 confirms the above properties and: ⅰ) the maximum occurs at low values of the population size (N∗ is low); ⅱ) the maximum growth rate is overly low, so that the slope of the growth curve at the exponential growth phase is overly low.
Remark 5.4. The above cases of α and c values indicate that the c parameter has a significant influence on the lag phase duration but the influence of the α value is low. Concerning the slope at the exponential growth phase, the α value has a significant influence, whereas the influence of the c value is low.
Remark 5.5. The parameter region α<1 leads to low slope of the population at the exponential growth phase, whereas c<0 implies that no lag phase occurs. Thus, the condition α>1 is required for higher slopes of the population at the exponential growth phase, whereas c>0 is required for representing the lag phase.
Proposition 5.5. Consider y=ln(N/Nto), and the model (3.1) with γ>0. The population at the logarithmic inflection point, N∗, is defined by
ˉax2+ˉbx+ˉc=0,x=N∗Nmin | (5.21) |
where
ˉa=α−1+γ,ˉb=−((α−1)(1+NmaxNmin)+γ−c),ˉc=(α−1−c)NmaxNmin | (5.22) |
and the discriminant is
ˉb2–4ˉaˉc=((α−1)(1−NmaxNmin)+γ−c)2+4γcNmaxNmin | (5.23) |
so that no imaginary solutions are possible for γc≥0. If ˉa>0, c>0 and ˉb2–4ˉaˉc is non-negative, only the positive square root leads to N∗>Nmin.
Proof. The population at the logarithmic inflection point, N∗, is obtained by applying Eq (5.3) to the growth model (3.1), leading to
((α−1)(Nmax−N)(N−Nmin)−γN(N–Nmin)+cNmin(Nmax−N))|N=N∗=0 | (5.24) |
Rearranging yields the quadratic Eq (5.21). In order to determine whether the negative square root term can lead to N∗>Nmin for ˉa>0, c>0 we use Proposition presented in Appendix B. Using the definition of ˉa, ˉb, ˉc, we have ˉc+ˉa+ˉb=(−1)c(Nmax/Nmin−1). Since (Nmax/Nmin)>1, then ˉc+ˉa+ˉb<0 for c>0. Hence, the condition ˉc+ˉa+ˉb>0 is not fulfilled for c>0. Consequently, only the positive square root term leads to N∗/Nmin>1.
Proposition 5.6. Limit values of N∗. Consider the model (3.1) with α−1>0, γ>0, c≥0. The limit values of the population at the logarithmic inflection point, N∗, are:
i)limγ→0N∗=Nmaxii)limγ→∞N∗=Nminiii)limc→0N∗=Nmin2(α−1+γ)[(α−1)(1+NmaxNmin)+γ+|(−1)(α−1)(NmaxNmin−1)+γ|]iv)limc→∞N∗=Nmaxv)lim(α−1)→0N∗=Nmin2γ(γ−c+√(c−γ)2+4cγNmaxNmin)vi)lim(α−1)→∞N∗=Nmax | (5.25) |
Proof. The above limit values of N∗ are obtained by substituting γ→0, γ→∞, c→0, c→∞, (α−1)→0, (α−1)→∞ in Eq (5.24) and solving for N∗.
Remark 5.6. The limit values of y∗ can be obtained straightforwardly by using y∗=ln(N∗/Nto) and the definitions (3.5), (3.6).
Proposition 5.7. Conditions for N∗=Nto. Consider the model (3.1). The population at the logarithmic inflection point, N∗, satisfies: N∗=Nto for
cγ+α−1(NmaxNmin−NtoNmin)=(NtoNmin–1)(NtoNmin−α−1γ+α−1NmaxNmin) | (5.26) |
Proof. Using N∗=Nto in Eq (5.24) and arranging yields Eq (5.26).
Remark 5.7. The above proposition implies that in the y versus t curve, with Nto∈(Nmin,Nmax), the maximum dy/dt occurs at t=to, N=Nto, if the α, γ, c values are such that Eq (5.26) holds true.
Proposition 5.8. Conditions for N∗=Nmy (equivalently, y∗=ym). Consider model (3.1). The population at the logarithmic inflection point, N∗, satisfies N∗=Nmy for
cγ+α−1=1−α−1γ+α−1√NmaxNmin | (5.27) |
Proof. Using N∗=Nmy, Nmy=(NminNmax)1/2 in Eq (5.24) and arranging yields Eq (5.27).
Remark 5.8. The above proposition implies that in the y versus t curve, with Nto∈(Nmin,Nmy), the maximum dy/dt occurs at N=Nmy; where Nmy=(NminNmax)1/2, if the α, γ, c values are such that Eq (5.27) holds true.
Remark 5.9. The determination of the lag time of the logarithmic growth curve via the classical approach requires the determination of the inflection point, as shown in Appendix A. To this end, the population at the inflection point (N∗) can be computed using expressions (5.21), (5.22), and the inflection time (t∗) must be computed numerically.
In this section, we propose an approach for the analytical determination of the lag time, for each the regular and the logarithmic measures. Also, a numerical example is provided, using the proposed approach and the classical approach.
The classical measure of the lag time requires the determination of the inflection point, as shown in Appendix A. In turn, the regular and the logarithmic population at the inflection point (N# and N∗) can be computed using the expressions (4.20), (4.21), (5.21), (5.22). The inflection times (t#,t∗) must be computed numerically.
The proposed approach is based on an application of the classical approach to a simplified growth model whose lag time value is quite similar to that of the original model. The simplified model is derived from proposed model (3.1) so as to fulfill the following properties: Rⅰ) its lag time value is quite similar to that of the original N versus t curve provided by model (3.1); Rⅱ) it allows to obtain an analytical expression of the population N; Rⅲ) it allows to apply the classical tlag approach. Thus, given a model simplified on the basis of the above facts, the tlag value of the simplified model can be considered as a tlag measure of the original model.
Disregarding the saturation term (1−N/NNmax) simplifies the determination of the lag phase duration, as was used by Baranyi, Roberts and McClure [1]. For the proposed model (3.1), with given values of α, c, γ, remarks 11, 20 indicate that: i) the effect of parameter α on the lag phase and the lag time is low; ii) during the lag phase and the early exponential growth phase, N≈Nto, so that the effect of simplification (1−N/Nmax)γ=(1−Nto/Nmax)γ on the lag time is low. Thus, a simplified form that accomplishes requirements Ri to Riii is:
dNs2dt=rNk1s2(1−NtoNmax)γ(1−NminNs2)c | (6.1) |
withk1∈{c−k2,c,2,3} | (6.2) |
where k2 is a positive integer, and the k1 option that is the nearest to the α value must be chosen. The Ns2 model is obtained from the proposed model (3.1) with the simplification α=k1 and (1−N/Nmax)γ=(1−Nto/Nmax)γ. The Ns2 curve has no inflection point, and it does not converge towards an equilibrium point, but its lag phase and the early part of the exponential growth phase are similar to those of the model (3.1). In the particular case k1=c, the analytical solution of Eq (6.1) is:
11−c[(Ns2−Nmin)1−c−(Nto−Nmin)1−c]=r(1−NtoNmax)γ(t−to),c≠1 | (6.3) |
The application of the classical tlag approach is different for the regular and the logarithmic growth curves.
To apply the classical tlag approach on the Ns2 versus t curve, we notice that this curve has no inflection point. Thus, we choose other point, located in the exponential growth phase, whose N value, denoted as Nps2 equals the inflection point of the simplified N versus t curve provided by
dNs1dt=rNk1s1(1−Ns1Nmax)γ(1−NminNs1)c | (6.4) |
so that Nps2=N#s1. This curve involves the simplification α=k1, k1 satisfying (6.2). In accordance with the classical tlag approach, the intersection of the line that is tangent to the Ns2 versus t curve at Ns2=N#s1, with the line Ns2=Nto yields the tlag value. The N#s1 value can be obtained via Eqs (4.20) and (4.21), with α=k1. In the particular case k1=c, the N#s1 value is obtained by applying Eq (4.3) to model (6.4) with k1=c:
N#s1=cNmax+γNminc+γ | (6.5) |
The time for Ns2=N#s1 (denoted ts) is obtained from the analytical solution of Eq (6.1), and using Ns2=N#s1. In the particular case k1=c, ts is obtained using Eq (6.3):
ts=t0+1r(1−NtoNmax)γ(1−c)[(N#s1−Nmin)1−c−(Nto−Nmin)1−c] | (6.6) |
The tlag can be computed from
tlag=ts−(N#s1−Nto)dNs2dt|Ns2=N#s1 | (6.7) |
where (dNs2/dt) at Ns2=N#s1 can be computed using Eq (6.1).
In summary, the application of the proposed approach for the regular growth curve involves the following steps:
● Recall the parameters of model (3.1) or obtain them by fitting to experimental data
● Select the k1 option (6.2) that is the nearest to the α value
● Determine the analytical expression of the Ns2 curve by using Eq (6.1), and in the case k1=c use Eq (6.3)
● Determine N#s1 value using Eqs (4.20) and (4.21), with α=k1; in the particular case k1=c, use Eq (6.5)
● Compute the time for Ns2=N#s1 (that is, ts), using the analytical expression for Ns2.
● Compute the tlag value using Eq (6.7).
As a numerical example we consider model (3.1), with r=1, Nmin=(1−1×10−6)Nto, Nto=1, γ=1; α=0.7 and c values c=0.74, c=1.1, c=1.2. We apply the aforementioned procedure, choosing k1=c, and the results are shown in Figure 8. Upper left: The time course of N for c=0.74, c=1.1, c=1.2 (solid, dash-dotted and dotted lines), with the inflection points (black squares) and the lag time points (gray squares) determined via the classical approach. Each of the remaining figures comprises: The N time course (thin solid line) with the inflection point (black square) and the lag time point (gray square) determined via the classical approach; the Ns1 time course (dash-dot line); and the Ns2 time course (thick solid line) with the point Ns2=N#s1 (black point), and the lag time point (gray point) determined via the proposed approach. The c values used are c=0.74 (upper right), c=1.1 (lower left), and c=1.2 (lower right).
To apply the classical tlag approach on the ys2=ln(Ns2/Nto) versus t curve, we notice that this curve has no inflection point. Thus, we choose other point, located in the exponential growth phase, whose ys2 value, denoted as yps2, results from the inflection point of the simplified y versus t curve provided by:
ys1=ln(Ns1/Nto) | (6.8) |
dNs1dt=rNk1s1(1−Ns1Nmax)γ(1−NminNs1)c | (6.9) |
so that yps2=y∗s1. This curve involves the simplification α=k1, being k1 chosen from options in (6.2). In accordance with the classical tlag approach, the intersection of the line that is tangent to the ys2 versus t curve at ys2=y∗s1, with the line ys2=ys2|t=to=0 yields the tlag value.
The y∗s1 value and the corresponding N∗s1 value (which satisfies y∗s1=ln(N∗s1/Nto)) are obtained from Eqs (5.21) and (5.22) with α=k1. The time for ys2=y∗s1, equivalently Ns2=N∗s1, denoted as ts, is obtained from the analytical solution of Eq (6.1) using Ns2=N∗s1. In the particular case k1=c, the time ts is obtained using Eq (6.3) with Ns2=N∗s1:
ts=t0+1r(1−NtoNmax)γ(−c+1)[(N∗s1−Nmin)−c+1−(Nto−Nmin)−c+1] | (6.10) |
The tlag value can be computed from
tlag=ts−ln(N∗s1/Nto)dys2dt|Ns2=N∗s1 | (6.11) |
where dys2/dt|Ns2=N∗s1 can be computed using dys2/dt=(1/Ns2)dNs2/dt, being dNs2/dt provided by Eq (6.1) with Ns2=N∗s1.
In summary, the application of the proposed approach for he logarithmic growth curve involves the following steps:
● Recall the parameters of model model (3.1) or obtain them by fitting to experimental data
● Select the k1 option (6.2) that is the nearest to the α value
● Determine the analytical expression of the Ns2 curve by using Eq (6.1), and in the case k1=c use Eq (6.3).
● Determine the y∗s1=ln(N∗s1/Nto) value, using Eqs (5.21) and (5.22) with α=k1.
● Compute the time for ys2=y∗s1 (that is, ts), using the analytical expression for Ns2 with Ns2=N∗s1.
● Compute the tlag value, using Eq (6.11)
As a numerical example, we consider model (3.1), with r=1, Nmin=(1−1×10−6)Nto, Nto=1, γ=1; α=0.7, and c values c=0.74, c=1.1, c=1.2. We apply the aforementioned procedure, choosing k1=c, as shown in Figure 9. Upper left: Time course of y for c values c=0.74, c=1.1, c=1.2 (solid, dash-dotted and dotted lines), with the inflection points (black squares) and lag time points (gray squares) determined via the classical approach. Each of the remaining figures comprise: The y time course (thin solid line) with the inflection point (black square) and the lag time point (gray square) determined via the classical approach; the ys1=ln(Ns1/Nto) time course (dash-dot line); and the ys2=ln(Ns2/Nto) time course (thick solid line) with the point ys2=y∗s1 (black point) and the lag time point (gray point) determined via the proposed approach. The c values used are c=0.74 (upper right), c=1.1 (lower left), and c=1.2 (lower right).
In this section, the proposed model (3.1) is applied to experimental microbial growth curves.
The differential Eq (3.1), the expression for the population at the inflection point N∗ and the expression for the lag time (tlag) are rewritten in terms of the logarithmic measure y=ln(N/Na) as follows, being Na a user defined positive constant. Using model (3.1), the above definition of y instead of definition (3.4), and the definition
ymin=ln(Nmin/Na),ymax=ln(Nmax/Na), | (7.1) |
one obtains:
dydt=ˉre(α−1)y(1−e(y−ymax))γ(1−e−(y−ymin))cˉr=rNα−1a, | (7.2) |
and its parameters are ˉr, α, γ, c, ymin, ymax. Using definition (7.1), the term Nmax/Nmin appearing in definitions (5.22) can be expressed as
NmaxNmin=eymax−ymin | (7.3) |
Definition y=ln(N/Na) implies y∗=ln(N∗/Na). Using definition (7.1) and Eq (5.21), one further obtains y∗=ln(eyminx), being x obtained from Eq (5.21). The t∗ value can be computed numerically by using Eq (7.2), so as to obtain y=y∗. The lag time can be computed via the classical approach with
tlag=t∗+(yto−y∗)dydt|y=y∗ | (7.4) |
being (dy/dt)|y∗ obtained using Eq (7.2).
Remark 7.1. In summary, the model properties in terms of the logarithmic measure includes: i) the logarithmic population at the inflection point (y∗), which can be computed using expression y∗=ln(eyminx), being x obtained from Eqs (5.21) and (5.22) with (Nmax/Nmin)=eymax−ymin; ii) the inflection time (t∗), which can be computed numerically using Eq (7.2); and iii) the lag time tlag which can be computed by means of Eq (7.4), being (dy/dt)|y∗ obtained using Eq (7.2). In this way, the computation of y∗, t∗, tlag neither requires N∗ nor Na.
The model of Baranyi and Roberts [31] is:
y=Ln(N), | (7.5) |
y(t)=yo+μmaxAt−1mLn(1+emμmaxAt−1em(ymax−yo)), | (7.6) |
At=t+1vLn(e−vt+qo1+qo), | (7.7) |
its parameters are μmax, m, ymax, qo, v; and the lag time is tlag=(1/v)Ln(1+1/qo).
The Gompertz model is [32]:
y=ya+(ymax−ya)exp(−exp(1+e1μmaxymax−ya(λ−t))) | (7.8) |
its parameters are μmax, ymax, λ, ya; and the lag time is tlag=λ.
In order to assess the capability of the proposed model to simulate growth curves, it and also the model of Baranyi and Roberts and the Gompertz model were fitted to growth records of Salmonellae. The datasets were taken from [1], corresponding to constant temperature and pH. The proposed model, the inflection point and the lag time expressed in terms of the logarithmic measure y are stated in subsection 7.1, whereas the model of Baranyi and Roberts and the Gompertz model are stated in subsection 7.2. The value Na=1 is used for the proposed model, in order to compare the simulation of y=ln(N) for the three models.
Model simulation and parameter estimation were performed using Matlab software (the Math Works Inc., Natick, Mass.). The differential Eq (7.2) was numerically integrated using the ode45 routine. Parameter estimation was carried out by minimizing the sum of the squared deviations between the values calculated by the model and the experimental data. The R2 value was used to evaluate the goodness of fit (see [33]), whereas the Akaike's Information Criterion (AIC) allows to compare the model acceptability on the basis of both fitting capability and complexity, being the complexity related to number of parameters [34,35]:
AIC=NLn(RSSN)+2(Nprm+1)+2(Nprm+1)(Nprm+2)N−Nprm−2 | (7.9) |
where N is the number of data points, Nprm is the number of model parameters, and RSS is the sum of the squared deviations between the values calculated by the model and the experimental data. The model with the lowest AIC value is the most acceptable. The fitting results are:
Generalized Fujikawa′s model:ˉr=0.92614,γ=0.4745,α=0.9647,c=0.8210,ymin=7.7897,ymax=20.15,y∗=10.977,t∗=8.715h,tlag=3.492h,R2=0.9919,AIC=−11.35,Model of Baranyi and Roberts:μmax=0.5518,m=1.1051,ymax=20.31,qo=0.0004574,v=2.7169tlag=2.8305h,R2=0.9922,AIC=−16.4288,Gompertz model:μmax=0.69007,ymax=21.013,λ=3.7876h,Ya=7.5373R2=0.9883,AIC=−11.4761, |
There is a good agreement between model simulation and experimental measurements (Figure 10), with a relatively high and similar R2 value for all the models. The estimated α, γ, c, are positive, and α is close to 1.0. The lag time tlag values were not quite similar. In terms of the AIC value, the Baranyi model is the most acceptable. This is partially related to the number of parameters, indeed Nprm=6, Nprm=5, Nprm=4 for the proposed model, the Baranyi's model and the Gompertz model, respectively. This worsens the AIC value of the proposed model in comparison with that of Baranyi's model. The fitting of the three models exhibit disagreement for the lag phase points (first seven points). In view of the above limitations, the fitting of the models to the lag phase points was improved, although at the cost of lowering the overall R2 value; α=1 was used in the proposed model, what is a common value, see [16,28]; and the R2 for the lag phase (R2lag) was calculated. The fitting results are:
GeneralizedFujikawa′smodel(α=1):ˉr=0.61107,γ=1.1423,c=0.88021,ymin=7.7785,ymax=20.2143,y∗=13.866,t∗=14.94h,tlag=4.9822h,R2=0.98845,AIC=−4.9322,R2lag=0.9663ModelofBaranyiandRoberts:μmax=0.6114,m=0.92099,ymax=20.1474,qo=0.047274,v=0.69237tlag=4.4745h,R2=0.9886,AIC=−8.1242,R2lag=0.96475Gompertzmodel:μmax=0.69168,ymax=20.9808,λ=4.00h,Ya=7.3672R2=0.98695,AIC=−9.0126,R2lag=0.84902 |
The three models achieved satisfactory agreement between simulation and experimental data, with a high R2 value (Figure 11). In terms of the AIC value, the Gompertz model is the most acceptable, but it exhibits significantly flawed fitting of several data points. The Baranyi model exhibits the second best AIC value. Both the Baranyi's model and the proposed model achieved proper fitting of the lag phase points, and the highest R2lag was achieved by the proposed model.
The obtained results confirm the capability of the proposed model to simulate monotonic growth curves, mainly those involving lag phase.
The studied generalized form of the Fujikawa’s growth model is autonomous, exhibits monotonic behavior of the population size and is capable of generating the common features of monotonic growth of microorganisms along time, including growth with decreasing slope without inflection point and growth with inflection point and lag phase. Indeed, either a large or a small lag time can be combined with either a large or a small slope of the exponential growth phase. Thus, the model is highly flexible and is suitable for different applications. Although Fujikawa, Kai and Morozumi [16,28] only employed positive values of α, c, negative values lead to features of the y versus t curve that are possible in real situations.
A deep insight on the growth curves is obtained by analyzing the dN/dt versus N curve at N=Nmin, N=Nmax and for N∈(Nmin,Nmax), in addition to assessing the presence of a maximum (which corresponds to an inflection point) or a minimum. A minimum and a maximum occur for α>0 and c<0, what is not common in growth curves.
The study of the effect of α and c values on the existence of the regular inflection point and lag phase indicates that: ⅰ) if both are positive, the inflection point occurs provided adequate Nto; ⅱ) if either α or c is negative, the occurrence of the inflection point is possible; ⅲ) if both parameters are negative, there is no inflection point nor lag phase. The existence of lag phase depends on the value of the growth rate at initial time, what is proportional to (Nto−Nmin).
The α and c parameters also affect the lag time magnitude and the slope of the exponential growth phase: ⅰ) the α value has a significant influence on the slope at the exponential growth phase, whereas the influence of the c value is low; ⅱ) the c parameter has a significant influence on the lag phase duration but the influence of the α value is low. For adequate representation of a growth curve that exhibits lag phase and steep exponential growth phase, it is required that both the α and the c values be positive.
The limit values of the inflection point population N# confirms it great flexibility, due to the parameters α, c, γ.
The numerical fitting example indicated the capability of the proposed model for fitting to experimental data of microbial growth. Indeed, the R2 value is high and similar to that of the Baranyi model.
A Rincón was supported by Universidad Católica de Manizales.
The authors declare that there is no conflict of interests regarding the publications of this paper.
For the regular measure (N), the tangent line at the inflection point (t#,N#) is
N=dNdt|t=t#t+N#−dNdt|t=t#t# | (A.1) |
The lag phase duration (tlag) is obtained from the intersection of this line with the horizontal line N=Nto, what gives
tlag=t#–(N#−Nto)dN/dt|t=t# | (A.2) |
For the logarithmic measure y=ln(N/Nto), the tangent line at the inflection point (t∗,y∗) is
y=dydt|t∗t+y∗−dydt|t=t∗t∗ | (A.3) |
where
dydt|t=t∗=1NdNdt|N=N∗ | (A.4) |
The lag phase duration (tlag) is obtained from the intersection of this line with the horizontal line y=yto=0, what gives
tlag=t∗–y∗dy/dt|t=t∗=t∗–ln(N∗/Nto)dy/dt|N=N∗ | (A.5) |
Proposition A.1. Consider the quadratic algebraic equation
ˉax2+ˉbx+ˉc=0 | (A.6) |
being ˉb2–4ˉaˉc non-negative. If ˉa>0, the solution with negative square root requires ˉc+ˉa+ˉb>0 for obtaining x>1.
Proof Consider the quadratic Eq (A.6), with ˉa>0, and the solution comprising the negative square root term. Accomplishment of condition x>1 implies
−√ˉb2−4ˉaˉc>2ˉa+ˉb√ˉb2−4ˉaˉc<(−1)(2ˉa+ˉb) | (A.7) |
This expression implies that 2ˉa+ˉb<0. Taking squares to Eq (A.7) and arranging yields −ˉc<ˉa+ˉb. Therefore, the solution with negative square root term, requires the condition ˉc+ˉa+ˉb>0 for achieving x>1.
[1] |
V. Saladino, D. Algeri, V. Auriemma, The psychological and social impact of Covid-19: new perspectives of well-being, Front. Psychol., (2020), 2550. https://doi.org/10.3389/fpsyg.2020.577684 doi: 10.3389/fpsyg.2020.577684
![]() |
[2] |
S. Shanbehzadeh, M. Tavahomi, N. Zanjari, I. Ebrahimi-Takamjani, S. Amiri-Arimi, Physical and mental health complications post-COVID-19: Scoping review, J. Psychosom. Res., 147 (2021), 110525. https://doi.org/10.1016/j.jpsychores.2021.110525 doi: 10.1016/j.jpsychores.2021.110525
![]() |
[3] |
I. Ali, O. M. Alharbi, COVID-19: Disease, management, treatment, and social impact, Sci. Total Environ., 728 (2020), 138861. https://doi.org/10.1016/j.scitotenv.2020.138861 doi: 10.1016/j.scitotenv.2020.138861
![]() |
[4] |
S. Pokhrel, R. Chhetri, A literature review on impact of COVID-19 pandemic on teaching and learning, Higher Educ. Future, 8 (2021), 133–141. https://doi.org/10.1177/2347631120983481 doi: 10.1177/2347631120983481
![]() |
[5] |
A. D. Kaye, C. N. Okeagu, A. D. Pham, R. A. Silva, J. J. Hurley, B. L. Arron, et al., Economic impact of COVID-19 pandemic on healthcare facilities and systems: International perspectives, Best Pract. Res. Clin. Anaesthesiol., 35 (2021), 293–306. https://doi.org/10.1016/j.bpa.2020.11.009 doi: 10.1016/j.bpa.2020.11.009
![]() |
[6] |
S. Naseer, S. Khalid, S. Parveen, K. Abbass, H. Song, M. V. Achim, COVID-19 outbreak: Impact on global economy, Front. Public Health, 10 (2023), 1009393. https://doi.org/10.3389/fpubh.2022.1009393 doi: 10.3389/fpubh.2022.1009393
![]() |
[7] |
M. Coccia, Factors determining the diffusion of COVID-19 and suggested strategy to prevent future accelerated viral infectivity similar to COVID, Sci. Total Environ., 729 (2020), 138474. https://doi.org/10.1016/j.scitotenv.2020.138474 doi: 10.1016/j.scitotenv.2020.138474
![]() |
[8] |
M. Coccia, How do low wind speeds and high levels of air pollution support the spread of COVID-19?, Atmos. Pollut. Res., 12 (2021), 437–445. https://doi.org/10.1016/j.apr.2020.10.002 doi: 10.1016/j.apr.2020.10.002
![]() |
[9] |
M. Coccia, The effects of atmospheric stability with low wind speed and of air pollution on the accelerated transmission dynamics of COVID-19, Int. J. Environ. Stud., 78 (2021), 1–27. https://doi.org/10.1080/00207233.2020.1802937 doi: 10.1080/00207233.2020.1802937
![]() |
[10] |
E. Bontempi, M. Coccia, International trade as critical parameter of COVID-19 spread that outclasses demographic, economic, environmental, and pollution factors, Environ. Res., 201 (2021), 111514. https://doi.org/10.1016/j.envres.2021.111514 doi: 10.1016/j.envres.2021.111514
![]() |
[11] |
E. Bontempi, M. Coccia, S. Vergalli, A. Zanoletti, Can commercial trade represent the main indicator of the COVID-19 diffusion due to human-to-human interactions? A comparative analysis between Italy, France, and Spain, Environ. Res., 201 (2021), 111529. https://doi.org/10.1016/j.envres.2021.111529 doi: 10.1016/j.envres.2021.111529
![]() |
[12] |
Y. Diao, S. Kodera, D. Anzai, J. Gomez-Tames, E. A. Rashed, A. Hirata, Influence of population density, temperature, and absolute humidity on spread and decay durations of COVID-19: A comparative study of scenarios in China, England, Germany, and Japan, One Health, 12 (2021), 100203. https://doi.org/10.1016/j.onehlt.2020.100203 doi: 10.1016/j.onehlt.2020.100203
![]() |
[13] |
C. Magazzino, M. Mele, M. Coccia, A machine learning algorithm to analyse the effects of vaccination on COVID-19 mortality, Epidemiol. Infect., 150 (2022), e168. https://doi.org/10.1017/S0950268822001418 doi: 10.1017/S0950268822001418
![]() |
[14] |
A. Nunez-Delgado, E. Bontempi, M. Coccia, M. Kumar, J. L. Domingo, SARS-CoV-2 and other pathogenic microorganisms in the environment, Environ. Res., 201 (2021), 111606. https://doi.org/10.1016/j.envres.2021.111606 doi: 10.1016/j.envres.2021.111606
![]() |
[15] |
M. Coccia, COVID-19 pandemic over 2020 (with lockdowns) and 2021 (with vaccinations): similar effects for seasonality and environmental factors, Environ. Res., 208 (2022), 112711. https://doi.org/10.1016/j.envres.2022.112711 doi: 10.1016/j.envres.2022.112711
![]() |
[16] |
H. Kato, A. Takizawa, Human mobility and infection from Covid-19 in the Osaka metropolitan area, npj Urban Sustainability, 2 (2022), 20. https://doi.org/10.1038/s42949-022-00066-w doi: 10.1038/s42949-022-00066-w
![]() |
[17] |
Y. Nohara, T. Manabe, Impact of human mobility and networking on spread of COVID-19 at the time of the 1st and 2nd epidemic waves in Japan: An effective distance approach, PLoS One, 17 (2022), e0272996. https://doi.org/10.1371/journal.pone.0272996 doi: 10.1371/journal.pone.0272996
![]() |
[18] |
A. Lison, J. Persson, N. Banholzer, S. Feuerriegel, Estimating the effect of mobility on SARS-CoV-2 transmission during the first and second wave of the COVID-19 epidemic, Switzerland, March to December 2020, Eurosurveillance, 27 (2022), 2100374. https://doi.org/10.2807/1560-7917.ES.2022.27.10.2100374 doi: 10.2807/1560-7917.ES.2022.27.10.2100374
![]() |
[19] |
K. Hibiya, A. Shinzato, H. Iwata, T. Kinjo, M. Tateyama, K. Yamamoto, et al., Effect of voluntary human mobility restrictions on vector-borne diseases during the COVID-19 pandemic in Japan: A descriptive epidemiological study using a national database (2016 to 2021), PLoS One, 18 (2023), e0285107. https://doi.org/10.1371/journal.pone.0285107 doi: 10.1371/journal.pone.0285107
![]() |
[20] |
M. Zhang, S. Wang, T. Hu, X. Fu, X. Wang, Y. Hu, et al., Human mobility and COVID-19 transmission: a systematic review and future directions, Ann. GIS, 28 (2022), 501–514. https://doi.org/10.1080/19475683.2022.2041725 doi: 10.1080/19475683.2022.2041725
![]() |
[21] |
N. Askitas, K. Tatsiramos, B. Verheyden, Estimating worldwide effects of non-pharmaceutical interventions on COVID-19 incidence and population mobility patterns using a multiple-event study, Sci. Rep., 11 (2021), 1972. https://doi.org/10.1038/s41598-021-81442-x doi: 10.1038/s41598-021-81442-x
![]() |
[22] |
R. S. John, J. C. Miller, R. L. Muylaert, D. T. Hayman, High connectivity and human movement limits the impact of travel time on infectious disease transmission, J. R. Soc. Interface, 21 (2024), 20230425. https://doi.org/10.1098/rsif.2023.0425 doi: 10.1098/rsif.2023.0425
![]() |
[23] |
J. Li, C. Zhuang, W. Zou, A tale of lockdown policies on the transmission of COVID-19 within and between Chinese cities: A study based on heterogeneous treatment effect, Econ. Hum. Biol., 53 (2024), 101365. https://doi.org/10.1016/j.ehb.2024.101365 doi: 10.1016/j.ehb.2024.101365
![]() |
[24] |
M. G. Thompson, J. L. Burgess, A. L. Naleway, H. Tyner, S. K. Yoon, J. Meece, et al., Prevention and attenuation of Covid-19 with the BNT162b2 and mRNA-1273 vaccines, N. Engl. J. Med., 385 (2021), 320–329. https://doi.org/10.1056/NEJMoa2107058 doi: 10.1056/NEJMoa2107058
![]() |
[25] |
V. Hall, S. Foulkes, F. Insalata, P. Kirwan, A. Saei, A. Atti, et al., Protection against SARS-CoV-2 after Covid-19 vaccination and previous infection, N. Engl. J. Med., 386 (2022), 1207–1220. https://doi.org/10.1056/NEJMoa2118691 doi: 10.1056/NEJMoa2118691
![]() |
[26] |
S. J. Thomas, E. D. Moreira, N. Kitchin, J. Absalon, A. Gurtman, S. Lockhart, et al., Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine through 6 months, N. Engl. J. Med., 385 (2021), 1761–1773. https://doi.org/10.1056/NEJMoa2110345 doi: 10.1056/NEJMoa2110345
![]() |
[27] |
J. Sadoff, G. Gray, A. Vandebosch, V. C..c3..a1rdenas, G. Shukarev, B. Grinsztejn, et al., Safety and efficacy of the BNT162b2 mRNA Covid-19 vaccine through 6 months, N. Engl. J. Med., 384 (2021), 2187–2201. https://doi.org/10.1056/NEJMoa2101544 doi: 10.1056/NEJMoa2101544
![]() |
[28] |
Y. Z. Huang, C. C. Kuan, Vaccination to reduce severe COVID-19 and mortality in COVID-19 patients: a systematic review and meta-analysis, Eur. Rev. Med. Pharmacol. Sci., 26 (2022). https://doi.org/10.26355/eurrev_202203_28248 doi: 10.26355/eurrev_202203_28248
![]() |
[29] |
B. Trogen, A. Caplan, Risk compensation and COVID-19 vaccines, Ann. Intern. Med., 174 (2021), 858–859. https://doi.org/10.7326/M20-8251 doi: 10.7326/M20-8251
![]() |
[30] |
K. P. Iyengar, P. Ish, R. Botchu, V. K. Jain, R. Vaishya, Influence of the Peltzman effect on the recurrent COVID-19 waves in Europe, Postgrad. Med. J., 98 (2022), e110–e111. https://doi.org/10.1136/postgradmedj-2021-140234 doi: 10.1136/postgradmedj-2021-140234
![]() |
[31] |
J. Guo, C. Deng, F. Gu, Vaccinations, mobility and COVID-19 transmission, Int. J. Environ. Res. Public Health, 19 (2021), 97. https://doi.org/10.3390/ijerph19010097 doi: 10.3390/ijerph19010097
![]() |
[32] |
L. L. Liang, H. M. Le, C. Y. Wu, C. Y. Sher, A. McGuire, Human mobility increased with vaccine coverage and attenuated the protection of COVID-19 vaccination: A longitudinal study of 107 countries, J. Glob. Health, 13 (2023). https://doi.org/10.7189/jogh.13.06009 doi: 10.7189/jogh.13.06009
![]() |
[33] |
M. Coccia, Improving preparedness for next pandemics: Max level of COVID-19 vaccinations without social impositions to design effective health policy and avoid flawed democracies, Environ. Res., 213 (2022), 113566. https://doi.org/10.1016/j.envres.2022.113566 doi: 10.1016/j.envres.2022.113566
![]() |
[34] |
H. Barbosa, S. Hazarie, B. Dickinson, A. Bassolas, A. Frank, H. Kautz, et al., Uncovering the socioeconomic facets of human mobility, Sci. Rep., 11 (2021), 8616. https://doi.org/10.1038/s41598-021-87407-4 doi: 10.1038/s41598-021-87407-4
![]() |
[35] | C. Kang, S. Gao, X. Lin, Y. Xiao, Y. Yuan, Y. Liu, et al., Analyzing and geo-visualizing individual human mobility patterns using mobile call records, in Proceedings of the 18th International Conference on Geoinformatics, (2010), 1–7. https://doi.org/10.1109/GEOINFORMATICS.2010.5567857 |
[36] |
A. Paez, F. A. Lopez, T. Menezes, R. Cavalcanti, M. G. d. R. Pitta, A spatio‐temporal analysis of the environmental correlates of COVID‐19 incidence in Spain, Geogr. Anal., 53 (2021), 397–421. https://doi.org/10.1111/gean.12241 doi: 10.1111/gean.12241
![]() |
[37] |
T. Chowdhury, H. Chowdhury, E. Bontempi, M. Coccia, H. Masrur, S. M. Sait, et al., Are mega-events super spreaders of infectious diseases similar to COVID-19? A look into Tokyo 2020 Olympics and Paralympics to improve preparedness of next international events, Environ. Sci. Pollut. Res., 30 (2023), 10099–10109. https://doi.org/10.1007/s11356-022-22660-2 doi: 10.1007/s11356-022-22660-2
![]() |
[38] |
M. Murakami, K. Fujii, W. Naito, M. Kamo, M. Kitajima, T. Yasutaka, et al., COVID-19 infection risk assessment and management at the Tokyo 2020 Olympic and Paralympic Games: A scoping review, J. Infect. Public Health, 17 (2024), 18–26. https://doi.org/10.1016/j.jiph.2023.03.025 doi: 10.1016/j.jiph.2023.03.025
![]() |
[39] |
S. Hu, C. Xiong, M. Yang, H. Younes, W. Luo, L. Zhang, A big-data driven approach to analyzing and modeling human mobility trend under non-pharmaceutical interventions during COVID-19 pandemic, Transp. Res. Part C: Emerg. Technol., 124 (2021), 102955. https://doi.org/10.1016/j.trc.2020.102955 doi: 10.1016/j.trc.2020.102955
![]() |
[40] |
D. Nakamoto, S. Nojiri, C. Taguchi, Y. Kawakami, S. Miyazawa, M. Kuroki, et al., The impact of declaring the state of emergency on human mobility during COVID-19 pandemic in Japan, Clin. Epidemiol. Glob. Health, 17 (2022), 101149. https://doi.org/10.1016/j.cegh.2022.101149 doi: 10.1016/j.cegh.2022.101149
![]() |
[41] |
M. Chakraborty, M. S. Mahmud, T. J. Gates, S. Sinha, Analysis and prediction of human mobility in the United States during the early stages of the COVID-19 pandemic using regularized linear models, Transp. Res. Rec., 2677 (2023), 380–395. https://doi.org/10.1177/03611981211067794 doi: 10.1177/03611981211067794
![]() |
[42] |
M. U. G. Kraemer, C. H. Yang, B. Gutierrez, C. H. Wu, B. Klein, D. M. Pigott, et al., The effect of human mobility and control measures on the COVID-19 epidemic in China, Science, 368 (2020), 493–497. https://doi.org/10.1126/science.abb4218 doi: 10.1126/science.abb4218
![]() |
[43] |
H. Fang, L. Wang, Y. Yang, Human mobility restrictions and the spread of the novel coronavirus (2019-nCoV) in China, J. Public Econ., 191 (2020), 104272. https://doi.org/10.1016/j.jpubeco.2020.104272 doi: 10.1016/j.jpubeco.2020.104272
![]() |
[44] |
J. Y. L. Chan, S. M. H. Leow, K. T. Bea, W. K. Cheng, S. W. Phoong, Z. W. Hong, et al., Mitigating the multicollinearity problem and its machine learning approach: a review, Mathematics, 10 (2022), 1283. https://doi.org/10.3390/math10081283 doi: 10.3390/math10081283
![]() |
[45] |
C. Strobl, A. L. Boulesteix, T. Kneib, T. Augustin, A. Zeileis, Conditional variable importance for random forests, BMC Bioinf., 9 (2008), 1–11. https://doi.org/10.1186/1471-2105-9-307 doi: 10.1186/1471-2105-9-307
![]() |
[46] | Prime Minister's Office of Japan, Prime minister's office of Japan. Available from: https://japan.kantei.go.jp. |
[47] | Ministry of Health, Labour and Welfare of Japan, Welcome to ministry of health, labour and welfare. Available from: https://www.mhlw.go.jp/english. |
[48] | Tokyo Metropolitan Government, Home-tokyo metropolitan government. Available from: https://www.metro.tokyo.lg.jp/english/index.html. |
[49] | Japan Meteorological Agency, Japan meteorological agency. Available from: https://www.jma.go.jp/jma/indexe.html. |
[50] | Portal Site of Official Statistics of Japan, Portal site of official statistics of Japan. Available from: https://www.e-stat.go.jp/en. |
[51] |
T. Watanabe, T. Yabu, Japan's voluntary lockdown, PLoS One, 16 (2021), e0252468. https://doi.org/10.1371/journal.pone.0252468 doi: 10.1371/journal.pone.0252468
![]() |
[52] |
Q. He, Z. Zhang, Y. Xie, The impact of COVID-19 on Americans' attitudes toward China: Does local incidence rate matter?, Soc. Psychol. Q., 85 (2022), 84–107. https://doi.org/10.1177/01902725211072773 doi: 10.1177/01902725211072773
![]() |
[53] | X. Chen, X. Di, How the covid-19 pandemic influences human mobility? similarity analysis leveraging social media data, in Proceedings of the IEEE 25th International Conference on Intelligent Transportation Systems (ITSC), IEEE, (2022), 2955–2960. https://doi.org/10.1109/ITSC55140.2022.9922060 |
[54] |
J. Yuan, M. Li, G. Lv, Z. K. Lu, Monitoring transmissibility and mortality of COVID-19 in Europe, Int. J. Infect. Dis., 95 (2020), 311–315. https://doi.org/10.1016/j.ijid.2020.03.050 doi: 10.1016/j.ijid.2020.03.050
![]() |
[55] |
H. Ito, N. Kawazoe, Examining transportation mode changes during COVID-19 in Toyama, Japan, Reg. Stud. Reg. Sci., 10 (2023), 253–272. https://doi.org/10.1080/21681376.2023.2180425 doi: 10.1080/21681376.2023.2180425
![]() |
[56] |
Y. Hara, H. Yamaguchi, Japanese travel behavior trends and change under COVID-19 state-of-emergency declaration: Nationwide observation by mobile phone location data, Transp. Res. Interdiscip. Perspect., 9 (2021), 100288. https://doi.org/10.1016/j.trip.2020.100288 doi: 10.1016/j.trip.2020.100288
![]() |
[57] |
F. Benita, Human mobility behavior in COVID-19: A systematic literature review and bibliometric analysis, Sustainable Cities Soc., 70 (2021), 102916. https://doi.org/10.1016/j.scs.2021.102916 doi: 10.1016/j.scs.2021.102916
![]() |
[58] |
E. Thuillier, L. Moalic, S. Lamrous, A. Caminada, Clustering weekly patterns of human mobility through mobile phone data, IEEE Trans. Mobile Comput., 17 (2017), 817–830. https://doi.org/10.1109/TMC.2017.2742953 doi: 10.1109/TMC.2017.2742953
![]() |
[59] | T. Mahajan, G. Singh, G. Bruns, An experimental assessment of treatments for cyclical data, in Proceedings of the 2021 Computer Science Conference for CSU Undergraduates, Virtual, 6 (2021), 22. |
[60] |
T. Horanont, S. Phithakkitnukoon, T. W. Leong, Y. Sekimoto, R. Shibasaki, Weather effects on the patterns of people's everyday activities: a study using GPS traces of mobile phone users, PLoS One, 8 (2013), e81153. https://doi.org/10.1371/journal.pone.0081153 doi: 10.1371/journal.pone.0081153
![]() |
[61] |
M. Yuan, Y. Lin, Model selection and estimation in regression with grouped variables, J. R. Stat. Soc. Ser. B: Stat. Methodol., 68 (2006), 49–67. https://doi.org/10.1111/j.1467-9868.2005.00532.x doi: 10.1111/j.1467-9868.2005.00532.x
![]() |
[62] | J. Friedman, T. Hastie, R. Tibshirani, A note on the group lasso and a sparse group lasso, preprint, arXiv: 1001.0736. https://doi.org/10.48550/arXiv.1001.0736 |
[63] |
N. Simon, J. Friedman, T. Hastie, R. Tibshirani, A sparse-group lasso, J. Comput. Graph. Stat., 22 (2013), 231–245. https://doi.org/10.1080/10618600.2012.681250 doi: 10.1080/10618600.2012.681250
![]() |
[64] |
L. Breiman, Random forests, Mach. Learn., 45 (2001), 5–32. https://doi.org/10.1023/A:1010933404324 doi: 10.1023/A:1010933404324
![]() |
[65] |
F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, et al., Scikit-learn: Machine learning in Python, J. Mach. Learn. Res., 12 (2011), 2825–2830. https://doi.org/10.48550/arXiv.1201.0490 doi: 10.48550/arXiv.1201.0490
![]() |
[66] | G. Ke, Q. Meng, T. Finley, T. Wang, W. Chen, W. Ma, et al., Lightgbm: A highly efficient gradient boosting decision tree, Adv. Neural Inf. Process. Syst., 30 (2017). |
[67] | T. Akiba, S. Sano, T. Yanase, T. Ohta, M. Koyama, Optuna: A next-generation hyperparameter optimization framework, in Proceedings of the 25th ACM SIGKDD International Conference on Knowledge Discovery & Data Mining, (2019), 2623–2631. https://doi.org/10.1145/3292500.3330701 |
[68] |
B. Gregorutti, B. Michel, P. Saint-Pierre, Grouped variable importance with random forests and application to multiple functional data analysis, Comput. Stat. Data Anal., 90 (2015), 15–35. https://doi.org/10.1016/j.csda.2015.04.002 doi: 10.1016/j.csda.2015.04.002
![]() |
[69] | T. Hastie, R. Tibshirani, J. H. Friedman, The Elements of Statistical Learning: Data Mining, Inference, and Prediction, Springer, 2 (2009). https://doi.org/10.1007/978-0-387-21606-5 |
[70] |
S. Luo, T. Chen, Two derivative algorithms of gradient boosting decision tree for silicon content in blast furnace system prediction, IEEE Access, 8 (2020), 196112–196122. https://doi.org/10.1109/ACCESS.2020.3034566 doi: 10.1109/ACCESS.2020.3034566
![]() |
[71] |
X. Zeng, Length of stay prediction model of indoor patients based on light gradient boosting machine, Comput. Intell. Neurosci., 2022. https://doi.org/10.1155/2022/9517029 doi: 10.1155/2022/9517029
![]() |
[72] |
Y. Hu, Z. Sun, Y. Han, W. Li, L. Pei, Evaluate pavement skid resistance performance based on Bayesian-LightGBM using 3D surface macrotexture data, Materials, 15 (2022), 5275. https://doi.org/10.3390/ma15155275 doi: 10.3390/ma15155275
![]() |
[73] |
S. Okamoto, State of emergency and human mobility during the COVID-19 pandemic in Japan, J. Transp. Health, 26 (2022), 101405. https://doi.org/10.1016/j.jth.2022.101405 doi: 10.1016/j.jth.2022.101405
![]() |
[74] | S. Fukui, Long-term changes in human mobility responses to COVID-19-related information in Japan, medRxiv, (2022), 2022–08. https://doi.org/10.1101/2022.08.15.22278703 |
[75] |
K. Henk, F. Rosing, F. Wolff, S. B. Frenzel, R. van Dick, V. A. Erkens, et al., An examination and extension of the Peltzman effect during the Covid-19 pandemic, Curr. Res. Ecol. Soc. Psychol., 4 (2023), 100091. https://doi.org/10.1016/j.cresp.2023.100091 doi: 10.1016/j.cresp.2023.100091
![]() |
[76] |
S. Peltzman, The effects of automobile safety regulation, J. Polit. Econ., 83 (1975), 677–725. https://doi.org/10.1086/260352 doi: 10.1086/260352
![]() |
[77] |
R. O. Nanda, A. A. Nursetyo, A. L. Ramadona, M. A. Imron, A. Fuad, A. Setyawan, et al., Community mobility and COVID-19 dynamics in Jakarta, Indonesia, Int. J. Environ. Res. Public Health, 19 (2022), 6671. https://doi.org/10.3390/ijerph19116671 doi: 10.3390/ijerph19116671
![]() |
[78] |
M. Coccia, COVID-19 vaccination is not a sufficient public policy to face crisis management of next pandemic threats, Public Organ. Rev., 23 (2023), 1353–1367. https://doi.org/10.1007/s11115-022-00661-6 doi: 10.1007/s11115-022-00661-6
![]() |
[79] |
M. Shibamoto, S. Hayaki, Y. Ogisu, COVID-19 infection spread and human mobility, J. Jap. Int. Econ., 64 (2022), 101195. https://doi.org/10.1016/j.jjie.2022.101195 doi: 10.1016/j.jjie.2022.101195
![]() |
[80] |
H. Kwon, C. Koylu, Revealing associations between spatial time series trends of COVID-19 incidence and human mobility: an analysis of bidirectionality and spatiotemporal heterogeneity, Int. J. Health Geogr., 22 (2023), 33. https://doi.org/10.1186/s12942-023-00357-0 doi: 10.1186/s12942-023-00357-0
![]() |
[81] |
A. Serisier, S. Beale, Y. Boukari, S. Hoskins, V. Nguyen, T. Byrne, et al., A case-crossover study of the effect of vaccination on SARS-CoV-2 transmission relevant behaviours during a period of national lockdown in England and Wales, Vaccine, 41 (2023), 511–518. https://doi.org/10.1016/j.vaccine.2022.11.073 doi: 10.1016/j.vaccine.2022.11.073
![]() |
[82] |
K. L. Andrejko, J. R. Head, J. A. Lewnard, J. V. Remais, Longitudinal social contacts among school-aged children during the COVID-19 pandemic: the Bay Area Contacts among Kids (BACK) study, BMC Infect. Dis., 22 (2022), 242. https://doi.org/10.1186/s12879-022-07218-4 doi: 10.1186/s12879-022-07218-4
![]() |
[83] | Ministry of Health, Labour and Welfare of Japan, COVID-19 Vaccine Q & A. Available from: https://www.mhlw.go.jp/stf/seisakunitsuite/bunya/vaccine_qa_archive.html#effect_4. |
[84] | Nippon Television Network Corporation, Corporate data. Available from: https://www.ntv.co.jp/english/an/cd.html. |
[85] | TV Asahi Corporation, Corporate profile. Available from: https://company.tv-asahi.co.jp/e/profile/index.html. |
[86] | TBS HOLDINGS, INC., Corporate profile. Available from: https://www.tbsholdings.co.jp/en/about/corporate/companyprofile.html. |
[87] | TV TOKYO Corporation, Corporate profile. Available from: https://www.tv-tokyo.co.jp/kaisha/company/profile.html. |
[88] | Fuji Television Network, Inc., Corporate profile. Available from: https://www.fujitv.com/about/corporate_profile. |
[89] |
M. Coccia, Effects of strict containment policies on COVID-19 pandemic crisis: lessons to cope with next pandemic impacts, Environ. Sci. Pollut. Res., 30 (2023), 2020–2028. https://doi.org/10.1007/s11356-022-22024-w doi: 10.1007/s11356-022-22024-w
![]() |
[90] |
S. Lee, E. Ko, K. Jang, S. Kim, Understanding individual-level travel behavior changes due to COVID-19: Trip frequency, trip regularity, and trip distance, Cities, 135 (2023), 104223. https://doi.org/10.1016/j.cities.2023.104223 doi: 10.1016/j.cities.2023.104223
![]() |
[91] |
I. Benati, M. Coccia, Global analysis of timely COVID-19 vaccinations: improving governance to reinforce response policies for pandemic crises, Int. J. Health Gov., 27 (2022), 240–253. https://doi.org/10.1108/IJHG-07-2021-0072 doi: 10.1108/IJHG-07-2021-0072
![]() |
[92] |
M. Coccia, Optimal levels of vaccination to reduce COVID-19 infected individuals and deaths: A global analysis, Environ. Res., 204 (2022), 112314. https://doi.org/10.1016/j.envres.2021.112314 doi: 10.1016/j.envres.2021.112314
![]() |
[93] |
M. Coccia, I. Benati, Negative effects of high public debt on health systems facing pandemic crisis: Lessons from COVID-19 in Europe to prepare for future emergencies, AIMS Public Health, 11 (2024), 477–498. https://doi.org/10.3934/publichealth.2024024 doi: 10.3934/publichealth.2024024
![]() |
[94] |
M. Coccia, Preparedness of countries to face COVID-19 pandemic crisis: strategic positioning and factors supporting effective strategies of prevention of pandemic threats, Environ. Res., 203 (2022), 111678. https://doi.org/10.1016/j.envres.2021.111678 doi: 10.1016/j.envres.2021.111678
![]() |
[95] |
M. Coccia, Pandemic prevention: lessons from COVID-19, Encyclopedia, 1 (2021), 36. https://doi.org/10.3390/encyclopedia1020036 doi: 10.3390/encyclopedia1020036
![]() |
[96] |
M. Coccia, Sources, diffusion and prediction in COVID-19 pandemic: lessons learned to face next health emergency, AIMS Public Health, 10 (2023), 145. https://doi.org/10.3934/publichealth.2023012 doi: 10.3934/publichealth.2023012
![]() |
[97] |
I. Benati, M. Coccia, Effective contact tracing system minimizes COVID-19 related infections and deaths: policy lessons to reduce the impact of future pandemic diseases, J. Public Admin. Govern., 12 (2022). https://doi.org/10.5296/jpag.v12i2.19834 doi: 10.5296/jpag.v12i2.19834
![]() |
[98] |
S. Khalilpourazari, H. H. Doulabi, A. Ö. Çiftçioğlu, G. W. Weber, Gradient-based grey wolf optimizer with Gaussian walk: Application in modelling and prediction of the COVID-19 pandemic, Expert Syst. Appl., 177 (2021), 114920. https://doi.org/10.1016/j.eswa.2021.114920 doi: 10.1016/j.eswa.2021.114920
![]() |
[99] | E. Savku, N. Azevedo, G. Weber, Optimal control of stochastic hybrid models in the framework of regime switches, in Modeling, Dynamics, Optimization and Bioeconomics Ⅱ: DGS Ⅲ, Porto, Portugal, February 2014, and Bioeconomy Ⅶ, Berkeley, USA, March 2014-Selected Contributions 3, Springer, (2014), 371–387. https://doi.org/10.1007/978-3-319-55236-1_18 |
[100] |
A. Özmen, E. Kropat, G. W. Weber, Robust optimization in spline regression models for multi-model regulatory networks under polyhedral uncertainty, Optimization, 66 (2017), 2135–2155. https://doi.org/10.1080/02331934.2016.1209672 doi: 10.1080/02331934.2016.1209672
![]() |
[101] |
E. Kropat, A. Özmen, G. W. Weber, S. Meyer-Nieberg, O. Defterli, Fuzzy prediction strategies for gene-environment networks–Fuzzy regression analysis for two-modal regulatory systems, RAIRO-Oper. Res., 50 (2016), 413–435. https://doi.org/10.1051/ro/2015044 doi: 10.1051/ro/2015044
![]() |
[102] |
E. Kropat, G. W. Weber, Fuzzy target-environment networks and fuzzy-regression approaches, Numer. Algebra Control Optim., 8 (2018), 135–155. https://doi.org/10.3934/naco.2018008 doi: 10.3934/naco.2018008
![]() |
[103] |
E. Özceylan, T. Paksoy, Fuzzy multi-objective linear programming approach for optimising a closed-loop supply chain network, Int. J. Prod. Res., 51 (2013), 2443–2461. https://doi.org/10.1080/00207543.2012.740579 doi: 10.1080/00207543.2012.740579
![]() |
[104] |
P. Haimerl, T. Hartl, Modeling COVID-19 infection rates by regime-switching unobserved components models, Econometrics, 11 (2023), 10. https://doi.org/10.3390/econometrics11020010 doi: 10.3390/econometrics11020010
![]() |
[105] |
R. Ananda, L. Harsyiah, M. Alfian, Classification of perceptions of the Covid-19 vaccine using multivariate adaptive regression spline, J. Varian, 6 (2023), 137–148. https://doi.org/10.30812/varian.v6i2.2639 doi: 10.30812/varian.v6i2.2639
![]() |
[106] |
F. Chen, H. Nakanishi, Y. Sekizawa, S. Ochi, M. So, Investigating the causal effects of COVID-19 vaccination on the adoption of protective behaviors in Japan: Insights from a fuzzy regression discontinuity design, PLoS One, 19 (2024), e0305043. https://doi.org/10.1371/journal.pone.0305043 doi: 10.1371/journal.pone.0305043
![]() |
[107] |
A. Szmelter-Jarosz, J. Ghahremani-Nahr, H. Nozari, A neutrosophic fuzzy optimisation model for optimal sustainable closed-loop supply chain network during COVID-19, J. Risk Financ. Manage., 14 (2021), 519. https://doi.org/10.3390/jrfm14110519 doi: 10.3390/jrfm14110519
![]() |
1. | Soumya G. Nair, Jyotirmayee Satapathy, N. Pavan Kumar, Influence of synthesis, dopants, and structure on electrical properties of bismuth ferrite (BiFeO3), 2020, 126, 0947-8396, 10.1007/s00339-020-04027-x | |
2. | Dao Viet Thang, Nguyen Manh Hung, Nguyen Cao Khang, Le Thi Mai Oanh, Structural and multiferroic properties of (Sm, Mn) co-doped BiFeO3 materials, 2020, 7, 2372-0484, 160, 10.3934/matersci.2020.2.160 | |
3. | A.T. Apostolov, I.N. Apostolova, J.M. Wesselinowa, Magnetic field effect on the dielectric properties of rare earth doped multiferroic BiFeO3, 2020, 513, 03048853, 167101, 10.1016/j.jmmm.2020.167101 | |
4. | S. R. Dhanya, Jyotirmayee Satapathy, Pavan Kumar, A review on the structural and magnetic properties of differently doped bismuth-ferrite multiferroics, 2021, 27, 0889-311X, 178, 10.1080/0889311X.2021.2020262 | |
5. | Mekonnen Tefera Kebede, Sheela Devi, Babita Tripathi, Sunil Chauhan, Venus Dillu, Structural transition and enhanced magnetic, optical and photocatalytic properties of novel Ce–Ni co-doped BiFeO3 nanoparticles, 2022, 152, 13698001, 107086, 10.1016/j.mssp.2022.107086 | |
6. | Dao Viet Thang, Nguyen Manh Hung, Du Thi Xuan Thao, Le Thi Mai Oanh, Do Danh Bich, Nguyen Cao Khang, Van Quang Nguyen, Nguyen Van Minh, Structural, Electrical, and Magnetic Properties of Bi0.90Nd0.10Fe0.98TM0.02O3 (TM=Mn, Co, Ni, and Cu) Materials, 2019, 10, 1949-307X, 1, 10.1109/LMAG.2019.2896847 | |
7. | Ibtissem Ragbaoui, Samia Aydi, Souad Chkoundali, Marwa Enneffati, Abdelhedi Aydi, Effect of substitution on the structural, electrical properties, and dielectric relaxor behavior in lead-free BiFeO3-based ceramics, 2024, 14, 2046-2069, 1330, 10.1039/D3RA06962A |