Gompertz function
A Gompertz curve or Gompertz function, named after Benjamin Gompertz, is a sigmoid function. It is a type of mathematical model for a time series, where growth is slowest at the start and end of a time period. The right-hand or future value asymptote of the function is approached much more gradually by the curve than the left-hand or lower valued asymptote, in contrast to the simple logistic function in which both asymptotes are approached by the curve symmetrically. It is a special case of the generalised logistic function.
|
|
|
Formula
where
- a is an asymptote, since
- b, c are positive numbers
- b sets the displacement along the x-axis (translates the graph to the left or right)
- c sets the growth rate (y scaling)
- e is Euler's Number (e = 2.71828...).
Derivation
The function curve can be derived from a Gompertz law of mortality, which states the rate of mortality (decay) falls exponentially with current size. Mathematically,
where
- is the rate of growth
- k is an arbitrary constant.
Example uses
Examples of uses for Gompertz curves include:
- Mobile phone uptake, where costs were initially high (so uptake was slow), followed by a period of rapid growth, followed by a slowing of uptake as saturation was reached.[1]
- Population in a confined space, as birth rates first increase and then slow as resource limits are reached.[2]
- Modeling of growth of tumors
- Modelling market impact in finance.[3]
Growth of tumors
In the 1960s A.K. Laird[4] for the first time successfully used the Gompertz curve to fit data of growth of tumors. In fact, tumors are cellular populations growing in a confined space where the availability of nutrients is limited. Denoting the tumor size as X(t) it is useful to write the Gompertz Curve as follows:
where:
- X(0) is the tumor size at the starting observation time;
- K is the carrying capacity, i.e. the maximum size that can be reached with the available nutrients. In fact it is:
independently on X(0)>0. Note that, in absence of therapies etc.. usually it is X(0)<K, whereas, in presence of therapies, it may be X(0)>K;
- α is a constant related to the proliferative ability of the cells.
- log() refers to the natural log.
It is easy to verify that the dynamics of X(t) is governed by the Gompertz differential equation:
i.e. is of the form:
where F(X) is the instantaneous proliferation rate of the cellular population, whose decreasing nature is due to the competition for the nutrients due to the increase of the cellular population, similarly to the logistic growth rate. However, there is a fundamental difference: in the logistic case the proliferation rate for small cellular population is finite:
whereas in the Gompertz case the proliferation rate is unbounded:
As noticed by Steel[5] and by Wheldon,[6] the proliferation rate of the cellular population is ultimately bounded by the cell division time. Thus, this might be an evidence that the Gompertz equation is not good to model the growth of small tumors. Moreover, more recently it has been noticed[7] that, including the interaction with immune system, Gompertz and other laws characterized by unbounded F(0) would preclude the possibility of immune surveillance.
Gompertz growth and logistic growth
The Gompertz differential equation
is the limiting case of the generalized logistic differential equation
(where is a positive real number) since
- .
In addition, there is an inflection point in the graph of the generalized logistic function when
and one in the graph of the Gompertz function when
- .
Gomp-ex law of growth
Based on the above considerations, Wheldon[6] proposed a mathematical model of tumor growth, called the Gomp-Ex model, that slightly modifies the Gompertz law. In the Gomp-Ex model it is assumed that initially there is no competition for resources, so that the cellular population expands following the exponential law. However, there is a critical size threshold such that for the growth follows the Gompertz Law:
so that:
Here there are some numerical estimates[6] for :
- for human tumors
- for murine (mouse) tumors
See also
References
- ↑ Islam, Towhidul; Fiebig, Denzil G.; Meade, Nigel (2002), "Modelling multinational telecommunications demand with limited data", International Journal of Forecasting, 18 (4): 605–624, doi:10.1016/S0169-2070(02)00073-0.
- ↑ Zwietering, M. H.; Jongenburger, I.; Rombout, F. M.; van 't Riet, K. (1990), "Modeling of the Bacterial Growth Curve", Applied and Environmental Microbiology, 56 (6): 1875–1881.
- ↑ Caravelli, F.; Sindoni, L.; Caccioli, F.; Ududec, C. (2015), Optimal leverage trajectories in presence of market impact, arXiv:1510.05123.
- ↑ Laird A. K. (1964). "Dynamics of tumor growth". Br J of Cancer. 18 (3): 490–502. doi:10.1038/bjc.1964.55.
- ↑ Steel, G.G. (1977). Growth Kinetics of Tumors. Oxford: Clarendon Press. ISBN 0-19-857388-X.
- 1 2 3 Wheldon, T.E. (1988). Mathematical Models in Cancer Research. Bristol: Adam Hilger. ISBN 0-85274-291-6.
- ↑ d'Onofrio A. (2005). "A general framework for modeling tumor-immune system competition and immunotherapy: Mathematical analysis and biomedical inferences". Physica D. 208 (3–4): 220–235. doi:10.1016/j.physd.2005.06.032.