Parameter identification for gompertz and logistic dynamic equations

In this paper, we generalize and compare Gompertz and Logistic dynamic equations in order to describe the growth patterns of bacteria and tumor. First of all, we introduce two types of Gompertz equations, where the first type 4-paramater and 3-parameter Gompertz curves do not include the logarithm of the number of individuals, and then we derive 4-parameter and 3-parameter Logistic equations. We notice that Logistic curves are better in modeling bacteria whereas the growth pattern of tumor is described better by Gompertz curves. Increasing the number of parameters of Logistic curves give favorable results for bacteria while decreasing the number of parameters of Gompertz curves for tumor improves the curve fitting. Moreover, our results overshadow some of the existing results in the literature.


Introduction
Most of the growth curves are described by linear, power, parabolic, power-exponential, logistic, log-logistic, von Bertalanffy, Gompertz, and Richards curves; see [1], [2], [3], [4], [5], and [6] for the tumor, [7] for the human fetus, [8] for the human life. A recent research article [8] related with a human life modeled by Gompertz and Mirror Gompertz differential equations are literature, see [9] is given by where k is the proportionality constant and K is the carrying capacity. In this study, mathematical modeling is applied to the Pseudomonas putida and mammary tumor datas given in [10,11], respectively. Note that Pseudomonas putida is a bacterium found in most soil and water habitats, and is significant to the environment due to its complex metabolism and ability to control pollution, [12] and [13]. We model their growth patterns by continuous and discrete Gompertz and Logistic curves. To achieve our goal, we derive 4-parameter and 3-parameter Gompertz and Logistic dynamic equations. We first propose two types of Gompertz dynamic equations: The first type Gompertz dynamic equations are motivated by [14]. We contribute two first type continuous Gompertz curves to the literature. All of the discrete Gompertz curves in this type are new. 4-parameter second type Gompertz dynamic equations are motivated by [2] in which only 3-parameter discrete Gompertz curves are considered. 3-parameter second type continuous Gompertz are investigated earlier in [10]. Inspired by [15], we come up with 4-parameter Logistic dynamic equations while 3-parameter Logistic dynamic equations are constructed earlier in [16]. 4-parameter Logistic discrete curves are new. To establish both dynamic equations, we use the variation of constant formulas together with the circle dot multiplication and the circle minus substraction on time scales. We refer readers to [17] and [18] by Bohner and Peterson for the theory of time scales calculus.
The parameters of these models are estimated by NonlinearModelFit function of Wolfram Mathematica 11.0 applying Monte Carlo simulation and our comparison is based on outputs following from the p-values of parameters, adjusted R-squared, and RMSE (root mean square error), RRMSE (Relative Root Mean Square Error), MAPE (Mean Absolute Percent Error), MAE (mean absolute error), U1 (Theil inequality coefficient, Theil's U1), and U2 (Theil inequality coefficient, Theil's U2). We use the Mathematica 11 program for the goodness of fit test of the models. Having at least three small values of each determined statistical criterion, the p value less than 0.05 for each parameter, and adjusted R-squared value close to 1 show better performance in terms of goodness of fit.
Outline of this paper is as follows: In Section 2, we introduce the time scales calculus together with some preliminary results. Sections 3 and 4 are related with first and second type Gompertz dynamic equations. In each section we obtain 4-parameter and 3-parameter continuous and discrete Gompertz curves. In Section 6, Logistic dynamic equations are introduced and we explicitly calculate 4-parameter and 3-parameter continuous and discrete Logistic curves. In the last section, we discuss how Gompertz and Logistic curves fit the growth of Pseudomonas putida and mammary tumor and include our conclusion.

Preliminary results
A time scale, T, is an arbitrary nonempty closed subset of the real numbers R. The theory of time scales is to introduce a new calculus so that we can unify the continuous and discrete analysis. Here, we give basic definitions and some essential results without proofs. Nevertheless, we mainly refer readers two books [17] and [18] by Bohner and Peterson and the manuscript [16] by Akin-Bohner and Bohner.
The forward jump operator σ on T is defined as sðtÞ≔inffs > t : s 2 Tg 2 T; for all t 2 T. For this definition we also have sð;Þ ¼ sup T. The backward jump operator r on T is defined by rðtÞ≔ sup fs < t : s 2 Tg 2 T; for all t 2 T. Here, we have rð; The set T k is derived from T: If T has left-scattered maximum m, then T k ¼ T À fmg. Otherwise, T k ¼ T. The following notations are also useful: f σ (t) = f(σ(t)). Note that t 2 ½t 0 ; 1Þ T ¼ ½t 0 ; 1Þ \ T.
Assume f : T7 !R and let t 2 T k , then we define f Δ (t) to be the number (provided it exists) with the property that given any � > 0, there is a neighborhood U of t such that j ½f ðsðtÞÞ À f ðsÞ� À f D ðtÞ½sðtÞ À s� j� � j sðtÞ À s j; for all s 2 U. f Δ (t) is called the delta derivative of f(t) at t. Note that the delta-derivative turns out to be the usual derivative when T ¼ R while it is the forward difference operator when If f is differentiable and t is right-dense, then If f ; g : T7 !R are differentiable at t 2 T k , then the product fg : T7 !R is also differentiable at t with If f is continuous at each right-dense point t 2 T and whenever t 2 T is left-dense lim s!t À f ðsÞ exists as a finite number, then we say that f : T7 !R is rd-continuous. A function F : T k 7 !R is called an antiderivative of f : T7 !R provided F Δ (t) = f(t) holds for all t 2 T k . In this case, we define the integral of f by Z t a f ðsÞDs ¼ FðtÞ À FðaÞ for t 2 T: ð5Þ The set of all regressive and rdcontinuous functions is denoted by R. Note that �p = −p, p � q = p + q and α � p = αp for the continuous case. If p is regressive, then we define the exponential function by where x h ðzÞ ¼ 1 h Logð1 þ hzÞ; h > 0 is the cylinder transformation such that ξ 0 (z) = z. If p : T k 7 !R is rd-continuous and regressive, then the exponential function e p (t, t 0 ) is the unique solution of the IVP on T for each fixed t 0 2 T k : For data analysis we need to calculate exponential functions for a regressive constant β, see Table 2.4 in [17]. We use the following properties of the exponential function e p (t, s), t; s 2 T. 5. e p (t, s)e q (t, s) = e p�q (t, s); In addition, two of the useful formulas for a circle dot are where p is a regressive function and α is a constant, see [16].
The followings are the variation of constants formulas, see Theorems 2.74 and 2.77 in [17]. The equation is called regressive if x Δ = p(t)x is regressive (i.e., p is regressive) and f : T ! R is rdcontinuous. Theorem 0.2. Suppose (12) is regressive. Let t 0 2 T and x 0 2 R. The unique solution of the IVP

First type gompertz dynamic equations
In this section, we will introduce Gompertz dynamic curves motivated by the 4-parameter Gompertz curve given in [19] for the growth curve analyses of bacterial counts. Here, K can be found as the growth rate coefficient, t 0 is the initial time, A + B is the carrying capacity of the environment for the population. To explain the carrying capacity notion we can say that every environment has its own limits, therefore it is impossible for species to grow up infinitely. Thus, the number of the population should be finite. In order to obtain the Gompertz model in the continuous case, we differentiate Eq (13) and obtain In addition, note that we have on ½t 0 ; 1Þ T , where we use Theorem 0.1 and (10). Since e �K (t, t 0 ) = e −K(t−t 0 ) for t 2 R, and (14) holds, then we obtain Motivated by the calculation above, we have the following initial value problem modeling 4-parameter Gompertz curve on time scales. Theorem 0.4. The initial value problem has the solution of the form where K is the growth rate and t 0 is the initial time, ω 0 is the value of the function at the initial time and B is the coefficient that has an impact on carrying capacity. Proof. We notice that the positivity of K implies the positivity of e �K by Theorem 0.1. Since 1 þ mð�e �K Þ ¼ 1 1þme �K > 0, we have the positively regressivity of �e �K . Since 1 + μ(K��e �K ) = (1 + μ(�e �K )) K > 0 by (11), the dynamic equation in the IVP (16) is regressive. Therefore, we apply Theorem 0.3 and obtain the unique solution for t 2 ½t 0 ; where we use (14) and Theorem 0.1. Example 0.5. Let T ¼ R. Then the continuous Gompertz curve is obtained from (17) for t 2 ½t 0 ; 1Þ R by using Eqs (8) and (15). This is compatible with the continuous Gompertz growth curve (13) (9), (14) yields for t 2 ½t 0 ; 1Þ Z . Thus, the discrete Gompertz growth curve again follows from (17) for t 2 ½t 0 ; 1Þ Z . Motivated by the first variation of constant formula, Theorem 0.2, we derive another Gompertz curve on time scales.
Theorem 0.7. The initial value problem has the solution of the form for t 2 ½t 0 ; 1Þ T , where K is the decay rate coefficient and regeressive, t 0 is the initial time, ω 0 is the value of the function at the initial time and B is the coefficient that has an impact on carrying capacity. Proof. Since K is regressive, e �K is also regressive by (7). The dynamic equation in the IVP (21) is regressive. Then, in order to obtain the unique solution (22) we apply Theorem 0.3 for where we use (14) and Theorem 0.1. Example 0.8. Let T ¼ R. Then the alternative continuous Gompertz curve is obtained from (22) for t 2 R. It is worth to mention that Eq (23) is a new Gompertz curve in the continuous case.
Example 0.9. Let T ¼ Z. Then using (19), we have again follows from (22) for t 2 ½t 0 ; 1Þ Z . The Gompertz growth curve (23) is given as 4-parameter Gompertz growth curve in [19]. From this point of view, the 3-parameter Gompertz growth curve on time scales is obtained from (17) when B = 0, and so the 3-parameter continuous and discrete Gompertz curves are for t 2 ½t 0 ; 1Þ R and for t 2 ½t 0 ; respectively. From (22) when B = 0, the alternative 3-parameter continuous and discrete Gompertz curves and are gained to the literature, respectively.

The zwietering modification of gompertz growth curve
The Gompertz growth curve is reparameterized in order to model the bacteria growth population in food and is stated as in [14] for t 2 R, where K z is the absolute growth rate at time T, so called lag time, which is interpreted as the time between when a microbial population is transferred to a new habitat recovers and when a considerable cell division occurs.
In order to find a corresponding dynamic model, we rewrite (30) as so that we can use the property of circle dot (10) for the unification the continuous and discrete cases. Therefore, using (15) and (14) yield the dynamic Zwietering Modification Gompertz curve Therefore, we claim that (32) is the solution of the IVP Since (30) is the continuous modified Gompertz growth curve, the discrete modified Gompertz growth curve follows from (32).
Example 0.10. Let T ¼ Z. Then the discrete Zwietering modification of Gompertz growth curve is obtained from (19) for t 2 ½0; 1Þ Z .

Gompertz-laird growth curve
This model is mainly used for the modeling of tumor growth. The Laird re-parameterization prevails even today as the most frequently fitted Gompertz version in cancer research, and is now also commonly fitted to growth data in other fields such as those of domestic (e.g. poultry and livestock, marine (e.g. molluscs, fish, and dolphins) animals. The continuous Gompertz-Laird growth curve is given by for t 2 ½0; 1Þ R in [14], which is equivalent to for t 2 ½0; 1Þ R , where the parameter L describes the initial specific growth rate that is not a notion that measures the relative growth or absolute growth. More precisely, we can say that the absolute growth rate at t = 0 is ω 0 .L Thus, the term L can be described as division of the initial absolute growth rate with the initial value. Similarly, by using (15) we obtain the Gompertz-Laird growth curve on time scales as for t 2 ½0; 1Þ T and so (36) is the solution of the IVP Since (35) is the continuous Gompertz-Laird growth curve, the following example gives the discrete Gompertz-Laird growth curve.
Example 0.11. Let T ¼ Z. Then we obtain as the discrete Gompertz-Laird growth curve for t 2 ½t 0 ; 1Þ Z , where we use again (19). If L = Km in (36), the Zweifel and Lasker re-parametrization dynamic equation is obtained for studying fish growth. Moreover, the continuous and discrete curves are given in (35)

Second type gompertz dynamic equations
It is clear that (13) is not a Gompertz model when the dependent variable is log-transformed. In this subsection, we will derive Gompertz dynamic equations involving logarithmic functions. This idea of the derivation of Gompertz dynamic equations is inspired from the Gompertz difference equation where a is taken as the growth rate and b is taken as the exponential rate of growth deceleration, which was firstly given by Bassukas et. al. [3]. The equivalent form of (39) is given by and so the continuous version becomes Unifying (40) and (41), we end up by ½ln GðtÞ� D ¼ a þ ðb À 1Þ ln GðtÞ; t 2 T: By the second variation of constants formula, Theorem 0.3, we get the alternative dynamic equation ½ln GðtÞ� D ¼ a À ðb À 1Þ ln GðsðtÞÞ; t 2 T: ð43Þ Notice that Eq (42) turns out to be (39) when T ¼ Z while (43) turns out to be for a nonzero constant b and when T ¼ R; we obtain the following Gompertz differential equation from the Gompertz dynamic Eq (43) as which is equivalent to Gompertz differential Eq (1) when a = 0, b − 1 = β, and G = x in (45). In [8], Gompertz differential Eq (41) with a = 0 is called the Mirror Gompertz differential equation, and is equivalent to (2) when a = 0, b − 1 = β, and G = 1 − x in (41). From now on, we take a = α and b − 1 = β in (42) and (43) and assume where g 0 is a real number. The following theorems yield the second type Gompertz dynamic curves. Theorem 0.12. Suppose that β is regressive and α > 0. Then the solution of the IVP (42)-(46) is given by for t 2 ½t 0 ; 1Þ T . Proof. If ln G(t) is taken as u(t), then the IVP (42)-(46) becomes u Δ (t) = α + βu(t), u(t 0 ) = g 0 , t 2 ½t 0 ; 1Þ T . Then by Theorems 0.1 and 0.2, we obtain for t 2 ½t 0 ; 1Þ T .
and GðtÞ ¼ exp ðe �b ðt; t 0 Þ½g 0 þ a R t t 0 e �b ðt 0 ; tÞDt�Þ respectively for t 2 ½t 0 ; 1Þ Z and here we use (9). In (51) by taking t 0 = 0, β = b − 1 and α = a Equation 3.1 is obtained in [11]. Both continuous Gompertz growth curves (18) and (50) From these observations, we conclude the 3-parameter first type continuous Gompertz curve and the second type continuous Gompertz curve are identical. However, since such an intimate relation among the discrete curves cannot be observed, considering first and second type discrete Gompertz curves contributes to the literature.

Logistic dynamic equations
In this section, we derive 4-parameter and 3-parameter Logistic continuous and discrete curves from Logistic dynamic equations.

3-Parameter logistic dynamic curves
Since there are two versions of linear equations there are two Logistic dynamic equations and respectively, see [16]. By using the definition of circle minus (6), Logistic dynamic Eqs (53) and (54) turn out to be the typical Logistic differential Eq (3) under certain conditions on p and f when T ¼ R. In [16], it is shown that the solutions of (53) and (54) subject to are given by and see Theorem 4.2 in [16]. Here, we assume that p is regressive, f is a rd-continuous function. We now calculate continuous and discrete Logistic solutions in order to compare their data fitting. In population dynamics, one often assumes that there exists a constant N 6 ¼ 0 such that p (t) = Nf(t) for all t 2 T. Example 0.16. Let T ¼ R, t 0 = 1, f = α and p = β, where α and β are constants. Then (56) and (57) turn out to be and respectively. Example 0.17. Let T ¼ Z, t 0 = 1, f = α and p = β, where α and β are constants. (56) and (57) turn out to be and respectively.

4-Parameter logistic dynamic curves
The 4-parameter Logistic curve is introduced and discussed in [15], where k, f, and s are positive constants and b is a real where we use (9). Let b = 1 in (67). Note that p ¼ k f o and this means that Eq (66) turns out to be ð71Þ If T ¼ R; then we obtain (3) from (71). One of the logistic dynamic equations is (54). By taking (54) into account and using the definition of minus circle, we get This implies that By the simple useful formula, we have Solving the above equation for L Δ , we get If we take L = ω, p = k and f 1 p ¼ 1 f in (72), then we obtain Eq (71). At this point the following question arises: Is it possible to find an alternative 4-parameter logistic dynamic equation which turns out to be the equivalent form of the 3-parameter logistic dynamic equation? To find the answer of this question consider Eq (65) where e k is replaced arrive at the unique solution as follows: If we take b = 1 in (75), then we get q ¼ o s f : Then Eq (74) turns out to be Furthermore, if T ¼ R; (79) turns out to be the logistic differential Eq (3). Note that logistic dynamic Eq (79) is equal neither (53) nor (54). Therefore, (78) with b = 1 is different than (60) and (61). It means that obtaining 3-parameter logistic curves from 4-parameter logistic curves yields new 3-parameter logistic discrete curves. The following example consists of two new 3-parameter logistic discrete curves. Example 0.22. If we let T ¼ Z and b = 1 in (70) and (78), then we obtain and respectively.

Goodness-of-fit test for gompertz and logistic curves and conclusion
The main aim of this study for the statistical analysis of Gompertz and Logistic curves is to determine whether their equations are able to model Pseudomonas Putita and tumor data sets given in [10] and [11]. In order to achieve our goal, Here, e t shows the error component, y t the original time series values,ŷ t the estimated values of the time series, and T the number of observations of the series. The criteria to determine an equation showing better performance in terms of goodness of fit is to have statistically significant coefficients; in other words, meaningful p-values for each coefficient, the larger adjusted R-squared value and smaller errors (see S1 and S2 Figs). The coefficient estimates of these models are obtained by the Mathematica 11.0 and the Wolfram Language uses "Conjuga-teGradient", "Gradient", "LevenbergMarquardt", "Newton", "NMinimize", and "QuasiNewton" methods.
When we concentrate attention on the Logistic type growth curves, from S1 Fig, it is apparent that all of the Logistic type equations are successful in modeling the bacteria data set. In addition, growth curves (70), (78), (80), (81) are the new 4-parameter discrete Logistic curves that are obtained in this study. According the results in S1 Fig, among the Logistic type growth curves, 4-parameter Logistic growth curves are better in modeling when it is compared with 3parameter ones.
Moreover, we can infer that the 4-parameter continuous Logistic curve (63) and (70), (78) model better than Eq (18). Thus, Pseudomonas Putita data is modeled by Logistic type equations better than Gompertz type equations. In addition, Eq (63) highlighted with orange, (70) and (78) highlighted with green in S1 Fig [14]. On the other hand, 3-parameter discrete Gompertz growth curve (29) also models the tumor data set and this equation is developed in this study as well. In addition, 4-parameter continuous and discrete second type Gompertz curves (51) and (52) are also new. 2-parameter version 51a0 of (51) was studied in [8] as Mirror Gompertz curve. Nevertheless, its discrete version 52α0 is a new contribution to the literature. At this point, we declare that 3-parameter Gompertz curves are more successful in modeling tumor data than 4-parameter Gompertz curves. Although the errors of the Logistic type curves are smaller than the errors of Gompertz type curves, all of their parameters are not statistically significant. Therefore, the Gompertz-Liard curve (35) is the best curve in modeling when it is compared with all the other curves and so one can say that Gompertz type curves have better fitting than Logistic type curves for tumor data.
As a result, Logistic curves are better in modeling bacteria data whereas tumor data is modeled better by Gompertz curves. Increasing the number of parameters of Logistic curves give favorable results for bacteria data while decreasing the number of the parameters of Gompertz curves for tumor data turns out to be reasonable. S3 Fig gives us the curve fittings of 4-parameter discrete Logistic curve (70) and 3-parameter discrete Gompertz curve (29) for bacteria data set and also shows the importance of the number of parameters.