1. Introduction
Fractional calculus is presently applied to a lot of scientific fields. Despite the problem of defining fractional derivatives being quite old (see, for instance, [
1,
2]), it has mainly been developed in recent times (see [
3]). Due to its versatility in describing slower or also different time scales, fractional derivatives and fractional-order differential equations are very often used in applications, so that also different books have been written on the argument (see, for instance, [
4,
5,
6]). The main generalization of the classical Cauchy problems to the fractional order is achieved via the so-called Caputo-fractional derivative, introduced by Michele Caputo in [
7]. In such paper, the fractional derivative is used to study the Q-factor of some non-ferromagnetic solids, thus being introduced in an applicative context. From such moment, fractional calculus has been used to address a lot of different models: from epidemics [
8] to osmosis [
9], from neurophysiology [
10] to viscoelasticity [
11] and many others [
12].
Here we focus on fractional-order population growth models. A first model of population growth can be achieved by modifying the classical Malthus model by introducing a fractional-order derivative in place of the classical one (see [
12,
13]). As a second step, one could ask for a fractional-order generalization of a Gompertz model. Gompertz model are quite popular growth model. Such models take into account a time-varying birth rate, which describes the fact that a person’s resistance to death decreases with age. Such models have been used in particular to model cancer growth, starting from [
14] and then used to describe a single species growth (see for instance [
15]). For this and other reasons, Gompertz curves have been widely studied. For instance, knowing that some species of cancer evolved following a Gompertz law, optimal control of it has become necessary (see for instance [
16]). At the same time, stochastic models became necessary to describe eventual environmental (and thus unpredictable) effects (see [
17,
18,
19,
20] and many others).
Concerning fractional-order Gompertz models, the first one has been introduced in [
21], but it is not achieved by simply substituting the fractional derivative in place of the classical one. To understand how the fractional-order model is introduced, let us recall that the classical Gompertz curve
can be defined as the solution of the non-linear Cauchy problem
where
are dynamical parameters and
is the initial population density. It is also well known that the solution is given by
and then the model admits a carrying capacity
If we define
, the Cauchy problem (
1) can be rewritten as
In [
21], the fractional-order Gompertz model is achieved by substituting the Caputo-fractional derivative in place of the classical one only in the linear equation. In particular, the function
is defined as the solution of the fractional Cauchy problem
where
is the fractional Caputo derivative of order
, and then defining the fractional Gompertz curve as
. In such case, we have
where
is the Mittag-Leffler function (defined in Equation (
8)).
In [
22] another type of fractional-order Gompertz model has been introduced. To understand how such model is defined, let us recall that for the classical model we have
and then we can rewrite the first equation of Equations (
2) as
In [
22], they use the Caputo-fractional derivative with respect to another function, as defined in [
23], to define the improved fractional Gompertz curve
the solution of
given by
In this paper, we aim to define a class of stochastic Gompertz models that generalize the two proposed fractional Gompertz curves. To do this, we first need to investigate some results related to a class of stochastic linear fractional-integral equations, concerning in particular the existence of Gaussian solutions. Such equations generalize the Caputo-fractional stochastic differential equations studied for instance in [
24,
25].
In particular this approach leads to a method of construction for general fractional growth models with noise that preserves normal or log-normal one-dimensional distributions. The preservation of such laws permits recognition in some macroscopical observable functions (the mean in the normal case and the median in the log-normal case) of the original growth models. Thus, these stochastic models work as a noisy perturbation of the original deterministic ones. This procedure could not be achieved by using the classical tools of fractionalization via time-change (see for instance [
26,
27,
28,
29]) for different reasons. For instance, if we apply a time-change to the stochastic Gompertz model, since the stochastic differential equation that drives the model is non-linear, its mean does not solve Equation (
1) with the Caputo derivative in place of the classical one. However, such time-changed process can be still seen as the exponential of a time-changed Ornstein-Uhlenbeck process (in the sense of [
28]), but, being the latter not a Gaussian one, the time-changed Gompertz model is not log-normal and its median does not coincide with the function
given in Equation (
4), despite the mean of the time-changed Ornstein-Uhlenbeck process is still solution of Equation (
3). Our new approach overtakes such problems, giving then some log-normal or normal processes whose dynamics are given by perturbation of the deterministic ones.
The paper is structured as follows:
In
Section 2 we give some basic definitions and preliminaries on fractional calculus;
In
Section 3 we study a class of linear fractional-integral stochastic equations: we will need them to define the stochastic models for fractional Gompertz curves. In particular, we focus on existence and almost surely uniqueness of Gaussian solutions. Moreover, since they are obtained via a Picard approximation method, we also give an estimate of the speed of convergence of the method in terms of the distribution of the maximum of the chosen noise.
In
Section 4 we give some examples on possible choices of noise. In particular in
Section 4.4 we show that such fractional-integral equations are indeed a generalization of the fractional stochastic differential equations discussed in [
24,
25].
In
Section 5 we use the results from the previous sections to introduce stochastic models for fractional Gompertz curves. In particular in
Section 5.1 we give some generalities on the classical stochastic Gompertz model, while in
Section 5.2 and
Section 5.3 we give a stochastic version of the fractional Gompertz curve introduced in [
21] and of the improved fractional Gompertz curve introduced in [
22]. Finally, in
Section 5.4 we construct a new fractional Gompertz model obtained by merging the approach of the previous two models and we describe a stochastic counterpart for it.
2. Some Preliminaries on Fractional Calculus
Concerning the main properties of fractional integrals and derivatives, we refer to [
30]. Let us give the following definition of the fractional-integral.
Definition 1. Given the fractional-integral of order ν is defined asfor any suitable function . It is easy to see that for instance, for any
the fractional-integral
is defined. Moreover, for any
it holds
. It is also interesting to notice that the fractional-integral
is a convolution operator. Indeed if we define the kernel
, then, for any function
where ∗ is the convolution product and
f is extended to the whole real line by setting
for any
. If
, then the convolution kernel
is singular, but still in
for any
. Therefore, while for any
, one only needs
f to be in
, it is not enough if
.
Now we can define the Riemann-Liouville derivative.
Definition 2. Given , the Riemann-Liouville fractional derivative of order ν is defined asfor any suitable function . From the definition of
, one easily obtains that
for any
. Thus, by the semigroup property of the fractional-integral and the fact that
is the classical integral, we have
for any suitable function
. In particular we have that
is the left inverse of
and thus, vice versa,
is the right inverse of
for any
.
However, we also have another fractional derivative.
Definition 3. Given , the Caputo-fractional derivative of order ν is defined asfor any suitable function . The class of functions for which
is defined is smaller than the one for which
is: indeed, one has at least to ask that
f is absolutely continuous. Moreover, we have that
hence, working as before, we have
for any suitable function
. We can conclude that
is the left inverse of
, and then
is the right inverse of
. There is also a relation between Riemann-Liouville and Caputo derivative:
From now on we will denote
. This relation lets us also define the Caputo-fractional derivative for any Riemann-Liouville derivable function, hence for a much wider class of functions. Concerning Caputo derivatives, we can define fractional Cauchy problems by using them. Indeed, under suitable assumptions, the fractional Cauchy problem
is well-posed. In particular, the relaxation problem
admits as unique solution the function
where
is the Mittag-Leffler function, defined as
which is a generalization of the exponential function (observe that if
,
).
We need also to introduce fractional calculus with respect to other functions. Riemann-Liouville type fractional derivative of a function with respect of another function were introduced to deal with Leibniz rule and chain rule for fractional derivatives (see, for instance, [
31,
32]). For this part, we mainly refer to [
23]. Let us first give the definition of fractional-integral with respect to another function.
Definition 4. Given and an increasing function such that for any , the fractional-integral with respect to Ψ
is given byfor any suitable function . Observe that if
, we achieve the classical fractional-integral. Let us now define the Riemann-Liouville type fractional derivative.
Definition 5. Given and an increasing function such that for any , the Riemann-Liouville fractional derivative with respect to Ψ
is given byfor any suitable function . Observe that for
, we achieve the classical Riemann-Liouville fractional derivative. Moreover, we have in this case
Let us also give the definition of the Caputo type fractional derivative.
Definition 6. Given and an increasing function such that for any , the Caputo-fractional derivative with respect to Ψ
is given byfor any suitable function . In [
23] (Theorem 3) the following relation is shown
Using this relation, one can extend the definition of Caputo-fractional derivative of a function with respect to another function to the whole class of the Riemann-Liouville derivable (with respect to
) functions. Moreover, under suitable assumptions, the following fractional Cauchy problem is well-posed
In the spirit of [
31,
32], let us show a chain rule for Caputo-fractional derivatives of a function with respect to another function.
Proposition 1. Let g be a Caputo-derivable function and Ψ
be an increasing function in such that for any and . Define . Then Proof. First, let us observe that
Deriving both sides of this relation and dividing by
, by using Equation (
9), we have
Finally, by substituting
and
in place of
f and
g and using Equation (
10) we conclude the proof. □
This proposition leads us to easily give the solution for the relaxation equation
whenever
. Indeed, if we define
as the solution of the relaxation equation for the Caputo derivative, hence
, and
, we have, by the previous proposition
thus
is the solution of Equation (
11).
3. Stochastic Linear fractional-integral Equations with Constant Coefficients and Gaussian Solutions
From now on let us fix a complete filtered space .
In this section, we want to study existence and uniqueness of solutions of stochastic linear fractional-integral equations in the form
where
,
, and
is a given
-adapted Gaussian process. From now on, as shorthand notation, let us denote
where
for some
or
, and
Remark 1. Obviously, for any and , we have .
Moreover, let us denote for any , where with .
3.1. The fractional-integral of a Gaussian Process
First, one could ask if the fractional-integral of a Gaussian process is still a Gaussian process. Concerning this problem, we have the following Lemma.
Lemma 1. Let for some time interval J and define for . Then . Moreover, if for some , then .
Proof. Let us consider
and recall that
. Fix
and observe that
is well-defined and continuous (in
t). We need to show that it is a
-adapted Gaussian process. Let us define
which is well-defined as Riemann integral since, for fixed
,
is continuous in
. To show that
is
-adapted, let us observe that, by definition of Riemann integral,
for any
where
and
, with
. Hence we have that almost surely
Since
Z is
-adapted and with a.s. continuous paths, it is progressively measurable and then, for any
and
,
is
-measurable and thus, being
,
-measurable. Hence the variable
is
-measurable for any
and so it is its limit as
, concluding that for any
,
is
-measurable. Now let us consider, for
,
. Let us observe that for fixed
we have, for
,
which is a
function. Hence, we have, by Lebesgue dominated convergence theorem, that
thus also
(being a.s. limit of
-measurable r.v.) is
-measurable.
Now let us show that
is a Gaussian process. Let us fix
,
,
and let us consider the random variable
given by
As before, if we define, for fixed
and
,
and
for
, we have that, by definition of Riemann integral,
for any
. Hence we have, for any
,
Since is a Gaussian process the random variable is Gaussian for any . Hence is almost surely limit of Gaussian random variables, hence it must be Gaussian.
As before, if we consider
, we have that for
hence
is almost surely limit of Gaussian random variables and must be Gaussian itself. The arbitrariness of
and
gives us the fact that
is a Gaussian process.
Finally, suppose that
and let us consider
. Suppose, without loss of generality, that
and set
. Hence we have for
concluding the proof. □
Let us remark that fractionally integrated Gauss-Markov processes have been also studied in [
33].
3.2. Compatibility between Fractionally Integrated Gaussian Processes
Now we want to study the behavior of a fractionally integrated Gaussian process with respect to other Gaussian processes. To do this let us first give the following shorthand notation.
Definition 7. Let and be two -adapted Gaussian processes with a.s. continuous paths. We say that and are compatible if for any , any and any the random variable is still a Gaussian random variable. This obviously implies that . Let us denote, for any , It is obvious that if and are independent -adapted Gaussian processes with a.s. continuous paths, then and are compatible.
Now let us show the following Lemma.
Lemma 2. Let such that . Then, setting for , .
Proof. Let us consider
and recall that
. Fix
and observe that
by the previous Lemma and that for
,
is continuous. Thus, we have that
is continuous for any
. Moreover, since
and
G are
-adapted,
is
-adapted. Now we need to show the compatibility property.
Let us fix
,
and
and let us define the random variables
Let us work on the second one. Fix
. Thus, by recalling the definition of
and
for
and
given in the previous lemma, we have that
hence we have that almost surely
where on the RHS we have Gaussian random variables since
and
are compatible. Hence
is a Gaussian random variable. Moreover, if we define
for
, one has that
almost surely, thus
is a Gaussian random variable and then
and
G are compatible. □
Remark 2. Obviously, for any and , we have . Moreover, for any and , we also have .
3.3. Main Result
Now we are ready to show an existence and uniqueness result in
for the solution of Equation (
12) in the fashion of [
6] (Theorem
).
Theorem 1. For any , , and , Equation (12) admits a unique solution . Moreover, if , then . Proof. Let us consider
and recall that
. Fix
and define the operator
as
where
is the Banach space
of the continuous functions equipped with the Bielecki norm
for some
, which is equivalent to the classical
norm.
Let us show that
is well-posed, i.e.,
is continuous. To do this consider
and suppose, without loss of generality, that
. Then, we can set
. We have
hence, being
continuous, sending
, we have
and then
.
Now consider
, choose
and set
p such that
. We have
and then
Taking the maximum, we have
Thus, we can choose
big enough to have
With this choice of
, we have that
is a contraction and thus admits a unique fixed point (see [
34], Theorem
): let us denote it as
.
Moreover, let us consider the sequence (for fixed
)
This sequence is such that
in
by contraction theorem (see [
34]).
Now let us define a stochastic operator
. For
, let us define it as
for any stochastic process such that
is continuous for any
, while for
let us complete it as we wish, since
is a null set.
We can re-interpret our sequence as a sequence of stochastic processes given by
Now let us observe that
is a (degenerate)
-adapted Gaussian process with a.s. continuous paths. For
, we have that
which obviously belongs to
by Remark 2. Let us suppose that
. By using Remark 2 and Lemmas 1 and 2, we have that
. Hence we have that for any
,
.
Now, we have that
where the limit is in the a.s. sense, thus it is easy to see that
(a.s. continuity of the paths follows from the continuity of
for
, since
are fixed points of
). Finally, a.s. uniqueness follows from the fact that
are contractions for
, hence their fixed point is unique.
Now, if
G is a.s.
-Hölder continuous, let us define
and let us recall that
. Consider
and
and observe that, from (
14), we have
where
is such that
(that exists since we have chosen
). In particular, we have that for any
,
. Almost surely
-Hölder continuity of the paths of
thus follows from the fact that, for any
,
. □
Remark 3. Let us observe that the fractional-integral operator is a compact Hilbert-Schmidt operator in (see, for instance [35]) if . Indeed, the integral kernel (where is the indicator function of the interval ) is such that if and only if . In such case, one can use the structure of the equation to show that there exists a unique Gaussian solution. Setting for instance and , we have for fixed where I is the identity operator; hence we have In such a way, for , one has the characterization of the solution Y as and then Y is given by a linear operator applied to a Gaussian process, hence it is Gaussian.
3.4. Speed of Convergence
We could also investigate the speed of convergence of the sequence
defined in Equation (
15) to
. The following proposition is an easy consequence of the contraction theorem.
Proposition 2. Consider for some , and . Moreover, consider solution of Equation (12). Set p such that and and fix such that . Finally, define . Then, for the sequence defined in Equation (15) we have almost surely. As a consequence it holds Proof. Fix
(where the set
A is defined in Equation (
13)). By contraction theorem (see [
34]) we have, since
L is the Lipschitz constant of
,
Now let us recall that for any function
thus, we have
Now let us recall that
and
, thus we have
Since
, Equation (
16) holds. □
3.5. The Mean of
Let us introduce another class of Gaussian processes
We want to investigate the mean of
, solution of (
12), when
. We have the following result.
Proposition 3. Fix and let us suppose that is in , where is solution of (12) in J. Then, is solution of the fractional Cauchy problem Proof. First, let us observe that
hence
. Now, let us notice that
Hence we can use Fubini’s theorem to achieve
Rearranging the equation and applying
on both sides we have
Since on the RHS we have a
function, we can differentiate both terms and use (
6) and (
7) to achieve
□
Remark 4. It is not difficult to show that if , is Riemann-Liouville derivable and is in , then is solution of the Cauchy problem The proof of such result is analogous to the previous one.
6. Conclusions
In this paper, we have given some methods to construct stochastic fractional Gompertz models by using stochastic linear fractional equations with Gaussian driving processes. The choice of a Gaussian driving process is linked to the necessity (in the first class of models introduced in
Section 5.2) to preserve the lognormality of the Gompertz model. In
Section 5.3 and
Section 5.4 we obtain stochastic fractional Gompertz models with Gaussian one-dimensional law. These were actually only exemplifications. Indeed, one can use the construction method given in
Section 5.3 and
Section 5.4 to obtain Gaussian stochastic models for general growth models of the form
where
for some growth rate
.
One could also try to substitute the operator
in place of
in Equation (
12). In such a case one could show, by similar arguments, the existence and uniqueness of a Gaussian solution and then use the construction given in
Section 5.2 to obtain a log-normal stochastic model for (
28).
Concerning possible applications, it has been already observed in [
21,
22] that fractional Gompertz models are more appropriate than classical ones to describe some phenomena such as tumor growth (concerning the model in
Section 5.2), dark fermentation and other fermentation phenomena (concerning the model in
Section 5.3). In this paper we provided a method to introduce noise (due to eventual unpredictable variables in the environment) in such a way that a macroscopic observable function still preserves such laws. Concerning the choice of the noise, it depends on the autocorrelation one wants to introduce in the model. For instance, if one wants to introduce a long-range (or short-range) correlated noise, one could use a fractional Brownian motion as driving Gaussian process, while if a delta-correlated noise is needed one could use a classical Brownian motion as driving process.
Finally, we want to recall that our aim was to introduce some construction methods that could lead to log-normal or normal stochastic models for general fractional growth processes (as the ones in Equation (
28)) with a general Gaussian noise, in order to provide a wide range of models that could be possibly useful in future applications.