1 Introduction
Stochastic processes with random time-change has become a well-established and highly ramified branch of the modern theory of stochastic processes which attracts more and more attention due to various applications in financial, biological, ecological, physical, technical and other fields of research.
Rich class of models is provided by time-changed Poisson processes, of which the most intensively studied are two fractional extensions of the Poisson process, namely, the space-fractional and time-fractional Poisson processes obtained by choosing a stable subordinator or its inverse in the role of time correspondingly. We refer, for example, to papers [8, 19, 21, 1, 2, 16, 20], among many others (see also references therein). In particular, in paper [21] a general class of time-changed Poisson processes $N({H^{\psi }}(t))$, $t\gt 0$, was introduced and studied, where $N(t)$ is a Poisson process and ${H^{\psi }}(t)$ is an arbitrary subordinator with the Laplace exponent ψ, independent of $N(t)$, and their distributional properties, hitting times and governing equations were presented (see, also [8]). In paper [5] Poisson processes time-changed by general inverse subordinators were studied, the governing equations for their marginal distributions were presented and some other properties were described. The Poisson process itself, being in a sense a core object concerning applicability to count data and simple tractability, however, as a reverse side of its simplicity, is a rather restrictive model. Therefore, it has been quite natural to search for its extensions and generalizations with some new useful features and properties needed for applications. For example, time-changed processes $N({H^{\psi }}(t))$ allow for jumps of arbitrary size and other interesting properties ([21, 8]). Another extension is provided by generalized counting processes introduced in [7]: such processes perform k kinds of jumps of amplitude $1,2,\dots ,k$ with rates ${\lambda _{1}},\dots ,{\lambda _{k}}$, for a fixed $k\ge 1$. Note that these processes and their time-changed versions are of interest in various applications, in particular, in risk theory (see, e.g., [13]).
Recently, generalized counting processes, their time-changed versions and fractional extensions have been intensively investigated, including Poisson and Skellam processes of order k, Pólya–Aeppli process of order k, Bell–Touchard process, Poisson-logarithmic process, etc. and their fractional extensions. We refer to papers [7, 15, 6, 10, 12, 11, 13], to mention only few, see also references therein.
In the present paper we consider several particular models of generalized counting processes time-changed by a general inverse subordinator, characterize their distributions and present governing equations for them, which are obtained following the technique presented in [5]. The equations are given in terms of the generalized Caputo–Djrbashian derivatives, which are called also convolution-type derivatives with respect to Bernštein functions. These generalized derivatives were introduced in [14] and [22] and have been widely used to describe various advanced stochastic models.
The convolution-type derivatives allow to study the properties of subordinators and their inverses in the unifying manner ([22, 17]), in particular, the governing equations for their densities can be given in terms of convolution-type derivatives. These properties will be the key tools for our study: by considering several models of time-changed processes, we elucidate the technique leading from the equations for the densities of inverse subordinators and their Laplace transforms to the equations for probabilities of processes time-changed by inverse subordinators and equations for some related functions.
Stable subordinators and their inverses are the most widely used tools for time-change and lead to so-called fractional in time/space or space-time fractional versions of original processes and corresponding fractional governing equations for probabilities of processes. In this paper we consider a general framework of double time-change by means of an arbitrary subordinator and an inverse to another subordinator, with a natural connection to generalized fractional calculus. To the best of our knowledge, the generalized counting processes with such general double time-change have not been considered in the literature. We also present a particular model of time-changed generalized counting process obtained by using the compound Poisson-Gamma subordinator. The state probabilities for this model can be given by the closed form expressions involving some special functions. In this way we complement the existing studies of time-changed processes.
The paper is organized as follows. Section 2 collects the basic definitions and facts on the convolution-type derivatives needed in our study. We also consider, as a warmup example, the time-changed Poisson process $N({Y^{f}}(t))$, $t\gt 0$, where ${Y^{f}}$ is the inverse to the subordinator with Laplace exponent f, independent of $N(t)$. We partly extend the results from [5] for this process. This simplest case allows to demonstrate very transparently the technique underpinning the derivations of the corresponding results for more general models considered in the next sections. Namely, in Sections 3 and 4 we study the processes ${N^{\psi ,f}}(t)=N\left({H^{\psi }}({Y^{f}}(t))\right)$, $t\gt 0$, and ${M^{\psi ,f}}(t)=M\left({H^{\psi }}({Y^{f}}(t))\right)$, $t\gt 0$, where M is the generalized counting process. Time-change is done by the independent subordinator ${H^{\psi }}$ and inverse subordinator ${Y^{f}}$ related to Bernštein functions ψ and f correspondingly. We present the differential-difference equations for probabilities of the processes ${N^{\psi ,f}}$ and ${M^{\psi ,f}}$, which are given in terms of the convolution-type derivative with respect to function f and the difference operator related to function ψ. Expressions for the probabilities of these processes are also presented in terms of the Laplace transform of the process ${Y^{f}}$. Some other properties of these processes are described and a comparison with related results in the literature is given. As examples, we consider the processes ${N^{\psi ,f}}$ and ${M^{\psi ,f}}$, where ${H^{\psi }}$ is a compound Poisson-Gamma process.
2 Preliminaries
Generalized fractional derivatives
We review the main definitions and some facts on the generalized fractional derivatives, which will be important for our study. (For more details, see, e.g., [17, 22].)
Let $f(x)$ be a Bernštein function,
${\overline{\nu }_{f}}(ds)$ is a nonnegative measure on $(0,\infty )$ (the Lévy measure for $f(x)$) such that
The generalized Caputo–Djrbashian (C-D) derivative, or convolution-type derivative, with respect to the Bernštein function f is defined on the space of absolutely continuous functions as follows ([22], Definition 2.4):
where ${\nu _{f}}(s)=a+{\overline{\nu }_{f}}(s,\infty )$ is the tail of the Lévy measure ${\overline{\nu }_{f}}(s)$ of the function f.
(1)
\[ f(x)=a+bx+{\int _{0}^{\infty }}\left(1-{e^{-xs}}\right){\overline{\nu }_{f}}(ds),\hspace{1em}x\gt 0,\hspace{2.5pt}a,b\ge 0,\](2)
\[ {\mathcal{D}_{t}^{f}}u(t)=b\frac{d}{dt}u(t)+{\int _{0}^{t}}\frac{\partial }{\partial t}u(t-s){\nu _{f}}(s)ds,\]In the case where $f(x)={x^{\alpha }}$, $x\gt 0$, $\alpha \in (0,1)$, the derivative (2) reduces to the fractional C-D derivative
For the Laplace transform of the derivative (2) the following relation holds ([22], Lemma 2.5):
\[ \mathcal{L}\left[{\mathcal{D}_{t}^{f}}u\right](s)=f(s)\mathcal{L}\left[u\right](s)-\frac{f(s)}{s}u(0),\hspace{1em}s\gt {s_{0}},\]
for u such that $|u(t)|\le \mathsf{M}{e^{{s_{0}}t}}$, M and ${s_{0}}$ being some constants. Similarly to the C-D fractional derivative, the convolution-type derivative can be alternatively defined via its Laplace transform.The generalization of the classical Riemann–Liouville (R-L) fractional derivative is introduced in [22] by means of another convolution-type derivative with respect to f given by the formula
The derivatives ${\mathcal{D}_{t}^{f}}$ and ${\mathbb{D}_{t}^{f}}$ are related as follows (see, [22], Proposition 2.7):
(4)
\[ {\mathbb{D}_{t}^{f}}u(t)=b\frac{d}{dt}u(t)+\frac{d}{dt}{\int _{0}^{t}}u(t-s){\nu _{f}}(s)ds.\]Bernštein functions are associated in a natural way with subordinators.
Let ${H^{f}}(t)$, $t\ge 0$, be a subordinator, that is, nondecreasing Lévy process. Its Laplace transform is of the form:
where the function f, called the Laplace exponent, is a Bernštein function. Consider a subordinator ${H^{f}}$ with the Laplace exponent f given by (1), and let ${Y^{f}}$ be its inverse process defined as
It was shown in [22] that the distribution of the inverse process ${Y^{f}}$ has a density ${\ell _{f}}(t,x)=\mathrm{P}\{{Y^{f}}(t)\in dx\}/dx$ provided that the following condition holds.
Condition I.
${\overline{\nu }_{f}}(0,\infty )=\infty $ and the tail ${\nu _{f}}(s)=a+{\overline{\nu }_{f}}(s,\infty )$ is absolutely continuous.
The Laplace transform of the density with respect to t has the form ([22])
The density ${\ell _{f}}(t,u)$ of the inverse process ${Y^{f}}$ satisfies the equation ([22], Theorem 4.1)
subject to
(7)
\[\begin{array}{r@{\hskip10.0pt}c@{\hskip10.0pt}l}& & \displaystyle {\mathbb{D}_{t}^{f}}{\ell _{f}}(t,u)=-\frac{\partial }{\partial u}{\ell _{f}}(t,u),\hspace{0.1667em}t\gt 0,\\ {} & & \displaystyle 0\lt u\lt \infty ,\hspace{2.5pt}\mathrm{if}\hspace{2.5pt}b=0,\hspace{2.5pt}\mathrm{and}\hspace{2.5pt}0\lt u\lt t/b,\hspace{2.5pt}\mathrm{if}\hspace{2.5pt}b\gt 0,\end{array}\]The space Laplace transform of the density ${\ell _{f}}(t,x)$,
is an eigenfunction of the operator ${\mathcal{D}_{t}^{f}}$, that is, satisfies the equation
with $\tilde{{\ell _{f}}}(0,\lambda )=1$ (see [5, 14, 17]).
(9)
\[ \tilde{{\ell _{f}}}(t,\lambda )={\int _{0}^{\infty }}{e^{-\lambda x}}{\ell _{f}}(t,x)dx=\mathsf{E}{e^{-\lambda {Y^{f}}(t)}},\hspace{1em}t\ge 0,\hspace{2.5pt}\lambda \gt 0,\](10)
\[ {\mathcal{D}_{t}^{f}}\tilde{{\ell _{f}}}(t,\lambda )=-\lambda \tilde{{\ell _{f}}}(t,\lambda ),\hspace{1em}t\gt 0,\]In the case where $f(x)={x^{\alpha }}$, $x\gt 0$, $\alpha \in (0,1)$, we have $\tilde{{\ell _{f}}}(t,\lambda )={\mathcal{E}_{\alpha }}(-\lambda {t^{\alpha }})$, where ${\mathcal{E}_{\alpha }}(\cdot )$ is the Mittag-Leffler function
Equations (7)–(10) are important for the study of the processes time-changed by inverse subordinators as can be seen in what follows. We first consider the Poisson process with an inverse subordinator as a simplest example which demonstrates the technique to be applied further to more general models.
Poisson process time-changed by an inverse subordinator
Let $N(t)$ be the Poisson process with intensity λ, and ${Y^{f}}(t)$ be the inverse subordinator. Consider the process ${N^{f}}(t)=N({Y^{f}}(t))$, $t\gt 0$.
Theorem 1.
Let Condition I hold. Then the marginal distributions ${p_{k}^{f}}(t)=P\left\{N({Y^{f}}(t))=k\right\}$, $k=0,1,\dots \hspace{0.1667em}$, satisfy the differential equation
with ${p_{k}^{f}}(0)=1$ for $k=0$, ${p_{k}^{f}}(0)=0$ for $k\ge 1$, and
where ${\tilde{\ell }_{f}}(t,\lambda )={\textstyle\int _{0}^{\infty }}{e^{-\lambda u}}{\ell _{f}}(t,u)du$ is the Laplace transform of the density of the inverse subordinator ${Y^{f}}(t)$.
(12)
\[ {\mathcal{D}_{t}^{f}}{p_{k}^{f}}(t)=\lambda \left({p_{k}^{f}}(t)-{p_{k-1}^{f}}(t)\right),\hspace{1em}k\ge 0,\hspace{2.5pt}t\gt 0\](13)
\[ {p_{k}^{f}}(t)=\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}{\tilde{\ell }_{f}}(t,\lambda ),\]Proof.
Proof of equation (12). For the probabilities ${p_{x}^{f}}(t)$ we have
We take the generalized R-L convolution-type derivative ${\mathbb{D}_{t}^{f}}$ given by (4) and use the equations (7)–(8) for the density ${\ell ^{f}}(t,u)$ of the inverse subordinator. We obtain
Using the relation (5) between convolution derivatives of C-D and R-L types, we have
and note also that
From (17), taking into account (18)–(19), we finally obtain
Proof of equation (13):
(16)
\[\begin{aligned}{}{p_{x}^{f}}(t)& =\mathsf{P}\left\{N\left({Y^{f}}(t)\right)=x\right\}={\int _{0}^{\infty }}{p_{x}}(u){\ell _{f}}(t,u)du\\ {} & ={\int _{0}^{\infty }}\frac{{e^{-\lambda u}}{(\lambda u)^{x}}}{x!}{\ell _{f}}(t,u)du,\hspace{1em}x=0,1,2,\dots .\end{aligned}\](17)
\[\begin{aligned}{}{\mathbb{D}_{t}^{f}}{p_{x}^{f}}(t)& ={\int _{0}^{\infty }}{p_{x}}(u){\mathbb{D}_{t}^{f}}{\ell _{f}}(t,u)du=-{\int _{0}^{\infty }}{p_{x}}(u)\frac{\partial }{\partial u}{\ell _{f}}(t,u)du\\ {} & ={\int _{0}^{\infty }}{\ell _{f}}(t,u)\frac{\partial }{\partial u}{p_{x}}(u)du-{p_{x}}(u){\ell _{f}}(t,u){\big|_{u=0}^{\infty }}\\ {} & ={\int _{0}^{\infty }}{\ell _{f}}(t,u)(-\lambda [{p_{x}}(u)-{p_{x-1}}(u)])du+{p_{x}}(0){\ell _{f}}(t,0)\\ {} & =-\lambda \left[{p_{x}^{f}}(t)-{p_{x-1}^{f}}(t)\right]+{p_{x}}(0){\nu _{f}}(t).\end{aligned}\](18)
\[ {\mathcal{D}_{t}^{f}}{p_{x}^{f}}(t)={\mathbb{D}_{t}^{f}}{p_{x}^{f}}(t)-{\nu _{f}}(t){p_{x}^{f}}(0),\](19)
\[ {p_{x}^{f}}(0)={\int _{0}^{\infty }}{p_{x}}(u){\ell _{f}}(0,u)du={\int _{0}^{\infty }}{p_{x}}(u)\delta (u)du={p_{x}}(0)=1.\]
\[\begin{aligned}{}\mathsf{P}\left\{N({Y^{f}}(t))=k\right\}& ={\int _{0}^{\infty }}\mathsf{P}(N(s)=k)\mathsf{P}({Y^{f}}(t)\in ds)\\ {} & ={\int _{0}^{\infty }}\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}{e^{-\lambda s}}\mathsf{P}\left\{{Y^{f}}(t)\in ds\right\}\\ {} & =\mathsf{E}\left[\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}\exp \left\{-\lambda {Y^{f}}(t)\right\}\right]\\ {} & =\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}\mathsf{E}\exp \left(-\lambda {Y^{f}}(t)\right)=\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}{\tilde{\ell }_{f}}(t,\lambda ).\end{aligned}\]
We used above that the distribution of the Poisson process N with the rate λ can be written as
For the probability generating function we have
\[\begin{aligned}{}{G^{f}}(u,t)& =\mathsf{E}{u^{N({Y^{f}}(t))}}={\int _{0}^{\infty }}{e^{-s\lambda (1-u)}}{\ell _{f}}(t,s)ds\\ {} & =\mathsf{E}{e^{-\lambda (1-u){Y^{f}}(t)}}={\tilde{\ell }_{f}}(t,\lambda (1-u)),\end{aligned}\]
therefore, we obtain formula (14). Now we use the fact, that ${\tilde{\ell }_{f}}(t,\lambda )$ is an eigenfunction of the operator ${\mathcal{D}_{t}^{f}}$ (see (10)), from which we conclude that ${G^{f}}(u,t)$, being given by ${\tilde{\ell }_{f}}(t,\lambda (1-u))$, satisfies equation (15). □3 Models of time-changed Poisson processes
Consider the time-changed Poisson processes ${N^{\psi }}(t)=N\left({H^{\psi }}(t)\right)$, $t\gt 0$, where $N(t)$ is the Poisson process with intensity $\lambda \gt 0$ and ${H^{\psi }}(t)$ is the subordinator with Bernštein function $\psi (u)$ (with $a=b=0$, see formula (1)), independent of $N(t)$. This class of processes was introduced and studied in [21] and called by the authors Poisson processes with Bernštein intertimes. It was shown in [21] that the distributions of ${N^{\psi }}(t)$, $t\gt 0$, can be presented as
and satisfy the difference-differential equations
Note that since ψ is a Bernštein function, it is in ${C^{\infty }}$, and the integrals in (21) with respect to its Lévy measure ${\overline{\nu }_{\psi }}$ are well defined and represent (up to the sign) the derivatives of ψ, so that (21) can be also written as
with the usual initial conditions ${p_{0}^{\psi }}(0)=1$, ${p_{k}^{\psi }}(0)=0$ for $k\ge 1$. The last equation can be represented in the form
where B is the backshift operator, $B{p_{k}^{f}}(t)={p_{k-1}^{f}}(t)$, and it is supposed that ${p_{-1}}(t)=0$.
(20)
\[ {p_{k}^{\psi }}(t)=\mathsf{P}\left\{{N^{\psi }}(t)=k\right\}=\frac{{(-1)^{k}}}{k!}\frac{{d^{k}}}{d{t^{k}}}{e^{-t\psi (\lambda u)}}{|_{u=1}},\](21)
\[ \frac{d}{dt}{p_{k}^{\psi }}(t)=-\psi (\lambda ){p_{k}^{\psi }}(t)+{\sum \limits_{m=1}^{k}}\frac{{\lambda ^{m}}}{m!}{p_{k-m}^{\psi }}(t){\int _{0}^{\infty }}{e^{-s\lambda }}{s^{m}}{\overline{\nu }_{\psi }}(ds),\hspace{1em}k\ge 0,\hspace{2.5pt}t\gt 0.\](22)
\[ \frac{d}{dt}{p_{k}^{\psi }}(t)=-\psi (\lambda ){p_{k}^{\psi }}(t)-{\sum \limits_{m=1}^{k}}{\psi ^{(m)}}(\lambda )\frac{{(-\lambda )^{m}}}{m!}{p_{k-m}^{\psi }}(t),\hspace{1em}k\ge 0,\hspace{2.5pt}t\gt 0,\](23)
\[ \frac{d}{dt}{p_{k}^{\psi }}(t)=-\psi (\lambda (I-B)){p_{k}^{\psi }}(t),\hspace{1em}k\ge 0,\hspace{2.5pt}t\gt 0,\]The probability generating function of the process ${N^{\psi }}$ is of the form ([21])
We refer for more detail on these processes to [21].
Consider the process ${N^{\psi }}$ time-changed by an independent inverse subordinator ${Y^{f}}$. We can state the following result.
Theorem 2.
The process ${N^{\psi ,f}}(t)=N\left({H^{\psi }}({Y^{f}}(t))\right)$, $t\gt 0$, has probability distribution function
and the probabilities ${p_{k}^{\psi ,f}}$ satisfy the equation
with the usual initial conditions.
(25)
\[ {p_{k}^{\psi ,f}}(t)=\mathsf{P}\left\{{N^{\psi ,f}}(t)=k\right\}=\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}{\tilde{\ell }_{f}}(t,\psi (\lambda )),\]
The probability generating function of the process ${N^{\psi ,f}}$ has the form
and satisfies the equation
with ${G^{\psi ,f}}(u,0)=1$.
Proof.
To prove equation (25), we perform the calculations as follows:
Using the relation between generalized derivatives of C-D and R-L types, we can write
\[\begin{aligned}{}{p_{k}^{\psi ,f}}(t)& =\mathsf{P}\left\{N\left({H^{\psi }}({Y^{f}}(t))\right)=k\right\}\\ {} & ={\int _{0}^{\infty }}\mathsf{P}\left\{N(s)=k\right\}\mathsf{P}\left\{{H^{\psi }}({Y^{f}}(t))\in ds\right\}\\ {} & ={\int _{0}^{\infty }}\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}{e^{-\lambda s}}\mathsf{P}\left\{{H^{\psi }}({Y^{f}}(t))\in ds\right\}\\ {} & =\mathsf{E}\left[\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}\exp \left\{-\lambda {H^{\psi }}({Y^{f}}(t))\right\}\right]\\ {} & =\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}\mathsf{E}\exp \left\{-\lambda {H^{\psi }}({Y^{f}}(t))\right\}.\end{aligned}\]
Now it is left to note that
\[ \mathsf{E}\exp \left\{-\lambda {H^{\psi }}({Y^{f}}(t))\right\}=\mathsf{E}\exp \left(-\psi (\lambda ){Y^{f}}(t)\right)={\tilde{\ell }_{f}}(t,\psi (\lambda )).\]
Let us state equation (26). We have
\[ {p_{k}^{\psi ,f}}(t)=\mathsf{P}\left\{N\left({H^{\psi }}({Y^{f}}(t))\right)=k\right\}={\int _{0}^{\infty }}{p_{k}^{\psi }}(u){\ell _{f}}(t,u)du,\hspace{1em}k=0,1,2,\dots .\]
We repeat reasoning in the same lines as in the proof of Theorem 1. We take the generalized R-L derivative ${\mathbb{D}_{t}^{f}}$ and use equations (22)–(23):
(29)
\[\begin{aligned}{}{\mathbb{D}_{t}^{f}}{p_{k}^{\psi ,f}}(t)& ={\int _{0}^{\infty }}{p_{k}^{\psi }}(u){\mathbb{D}_{t}^{f}}{\ell _{f}}(t,u)du=-{\int _{0}^{\infty }}{p_{k}^{\psi }}(u)\frac{d}{du}{\ell _{f}}(t,u)du\\ {} & ={\int _{0}^{\infty }}{\ell _{f}}(t,u)\frac{d}{du}{p_{k}^{\psi }}(u)du-{p_{k}^{\psi }}(u){\ell _{f}}(t,u){\big|_{0}^{\infty }}\\ {} & ={\int _{0}^{\infty }}{\ell _{f}}(t,u)\left(-\psi (\lambda (I-B))\right){p_{k}^{\psi }}(u)du+{p_{k}^{\psi }}(0){\ell _{f}}(t,0)\\ {} & ={\int _{0}^{\infty }}{\ell _{f}}(t,u)\Big(-\Big(\psi (\lambda ){p_{k}^{\psi }}(u)\\ {} & \hspace{1em}+{\sum \limits_{m=1}^{k}}{\psi ^{(m)}}(\lambda )\frac{{(-\lambda )^{m}}}{m!}{p_{k-m}^{\psi }}(u)\Big)\Big)du+{p_{k}^{\psi }}(0){\nu _{f}}(t)\\ {} & =-\psi (\lambda ){p_{k}^{\psi ,f}}(t)-{\sum \limits_{m=1}^{k}}{\psi ^{(m)}}(\lambda )\frac{{(-\lambda )^{m}}}{m!}{p_{k-m}^{\psi ,f}}(t)+{p_{k}^{\psi }}(0){\nu _{f}}(t)\\ {} & =-\psi (\lambda (I-B)){p_{k}^{\psi ,f}}(t)+{p_{k}^{\psi }}(0){\nu _{f}}(t).\end{aligned}\]
\[ {\mathcal{D}_{t}^{f}}{p_{k}^{\psi ,f}}(t)={\mathbb{D}_{t}^{f}}{p_{k}^{\psi ,f}}(t)-{\nu _{f}}(t){p_{k}^{\psi ,f}}(0);\]
and we have (see (8))
\[ {p_{k}^{\psi ,f}}(0)={\int _{0}^{\infty }}{p_{k}^{\psi }}(u){\ell _{f}}(0,u)du={\int _{0}^{\infty }}{p_{k}^{\psi }}(u)\delta (u)du={p_{k}^{\psi }}(0)=1.\]
Therefore, we obtain
Using the expression for ${G^{\psi }}(u,t)$ given by (24), we calculate ${G^{\psi ,f}}$:
\[\begin{aligned}{}{G^{\psi ,f}}(u,t)& =\mathsf{E}{u^{N\left({H^{\psi }}({Y^{f}}(t))\right)}}={\int _{0}^{\infty }}{G^{\psi }}(u,s){\ell _{f}}(t,s)ds=\\ {} & ={\int _{0}^{\infty }}{e^{-s\psi (\lambda (1-u))}}{\ell _{f}}(t,s)ds=\mathsf{E}{e^{-\psi (\lambda (1-u)){Y^{f}}(t)}}=\\ {} & ={\tilde{\ell }_{f}}(t,\psi (\lambda (1-u)));\end{aligned}\]
and since ${\tilde{\ell }_{f}}(t,\psi (\lambda (1-u)))$ is an eigenfunction of ${\mathcal{D}_{t}^{f}}$ with the eigenvalue $\psi (\lambda (1-u))$ (see (10)), it follows that ${G^{\psi ,f}}(u,t)$ satisfies equation (28). □Remark 2.
We notice at once that equations (26) for the probabilities of the process ${N^{\psi ,f}}(t)=N({H^{\psi }}({Y^{f}}(t)))$ mimic the corresponding equations for the process ${N^{\psi }}(t)=N({H^{\psi }}(t))$, that is, the process before the time-change by an inverse subordinator, only the ordinary derivative in time is replaced with the generalized fractional derivative. This can be anticipated quite straightforwardly, in view of the technique used for the proof, and holds also for other models, like those below and in the next section. On the other side, the equation (26) tells us, that for the probabilities of the processes with double time-change like $N\left({H^{\psi }}({Y^{f}}(t))\right)$, the action (in time) of the operator ${\mathcal{D}_{t}^{f}}$ (which is related to ${Y^{f}}$) is equal to the action (in space) of the operator $-\psi \left(\lambda \left(I-B\right)\right)$ (which is related to ${H^{\psi }}$, and depends also on the outer process N). Some more insight on these operators can be found within the approach applied in the paper [1]. Let ${Y^{f}}={Y^{\beta }}$ be the inverse stable subordinator with a parameter $\beta \in (0,1]$, that is, $f(\lambda )={\lambda ^{\beta }}$, then equation (26) becomes
where ${\mathcal{D}_{t}^{\beta }}$ is the Caputo–Djrbashian fractional derivative (3). If we suppose furthermore that ${H^{\psi }}={H^{\alpha }}$ is the stable subordinator with a parameter $\alpha \in (0,1]$, that is, $\psi (\lambda )={\lambda ^{\alpha }}$, then we come to the equation
Equations (30) and (31) were stated in paper [1], where the particular representation of the operator $\psi \left(\lambda \left(I-B\right)\right)$ was used, and the equations for the probabilities of time-changed processes were stated using the interplay between the operator $\psi \left(\lambda \left(I-B\right)\right)$ and the fractional operator ${\mathcal{D}_{t}^{\beta }}$ in the Fourier domain. The approach in [1] can be extended for the case where we deal with the generalized fractional operator ${\mathcal{D}_{t}^{f}}$, as soon as we know that its eigenfunction is given by ${\tilde{\ell }_{f}}(t,\lambda )$ (see (10)). We refer for more detail on this approach to [1].
(30)
\[ {\mathcal{D}_{t}^{\beta }}{p_{k}^{\psi ,\beta }}(t)=-\psi \left(\lambda \left(I-B\right)\right){p_{k}^{\psi ,\beta }}(t),\](31)
\[ {\mathcal{D}_{t}^{\beta }}{p_{k}^{\alpha ,\beta }}(t)=-{\lambda ^{\alpha }}{\left(I-B\right)^{\alpha }}{p_{k}^{\alpha ,\beta }}(t).\]Remark 3.
One of the important properties of the Poisson process is that it is a renewal process having i.i.d. waiting times (or inter-arrival times) ${J_{n}}$ with exponential distribution. Fractional in time Poisson process $N({Y^{\beta }}(t))$, $t\gt 0$, where ${Y^{\beta }}$ is the inverse stable subordinator with a parameter $\beta \in (0,1]$, keeps the renewal property with distribution of waiting times represented by means of the Mittag-Leffler function,
where λ is a rate of the Poisson process N, the Mittag-Leffler function is given by (11). Moreover, this property is retained by the process $N({Y^{f}})$, where ${Y^{f}}$ is a general inverse subordinator: $N({Y^{f}})$ is a renewal process with distribution of waiting times represented as
(see, e.g., [18, 14]).
Theorem 1 (equation (12)) shows that the distributions of the process $N({Y^{f}})$ solve an analogue of the Kolmogorov forward equation for the Poisson process, with the ordinary derivative in time replaced with the convolution-type derivative (or the C-D fractional derivative for the case of process $N({Y^{\beta }})$).
Equation (26) for distributions of processes ${N^{\psi ,f}}$ appears in the form which is substantially different from the case of the processes N and $N({Y^{f}})$: this equation contain additionally the probabilities ${p_{k-j}^{\psi ,f}}(t)$ for $2\le j\le k$.
For the processes with double time-change ${N^{\psi ,f}}(t)=N({H^{\psi }}({Y^{f}}(t)))$, $t\gt 0$, the renewal property does not hold in general. In particular, in [8] it was shown that for the space-time fractional Poisson processes $N({H^{\alpha }}({Y^{\beta }}(t)))$ (with ${H^{\alpha }}$ and ${Y^{\beta }}$ being stable and inverse stable subordinators) the renewal property does not hold, opposite to the case of time-fractional Possion process $N({Y^{\beta }}(t))$.
The waiting times of k-th event of the processes $N({H^{\alpha }}({Y^{\beta }}(t)))$ were investigated in [1]. Namely, for
Following [1], for the general case of process ${N^{\psi ,f}}$, we can consider the waiting times
For the case of space-time fractional Poisson process $N({H^{\alpha }}({Y^{\beta }}(t)))$ we have in (33) ${\tilde{\ell }_{f}}(t,\psi (\lambda ))={\mathcal{E}_{\beta }}(-{t^{\beta }}{\lambda ^{\alpha }})$ and the derivation of formula (32) relies on particular properties of the Mittag-Leffler function and the simple form of the Bernštein function of a stable subordinator (see the proof of Theorem 2 in [1]). It would be interesting to obtain expressions for other particular cases of the process ${N^{\psi ,f}}$. We address this problem to further research.
\[ {T_{k}^{\alpha ,\beta }}=\inf \{s\ge 0\hspace{0.1667em}:\hspace{0.1667em}N({H^{\alpha }}({Y^{\beta }}(t)))\gt k\}\]
the density of distribution was derived in the form ([1], Theorem 2)
(32)
\[ \mathsf{P}\{{T_{k}^{\alpha ,\beta }}\in dt\}/dt=\frac{\beta k}{\alpha t}\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}{\mathcal{E}_{\beta }}(-{t^{\beta }}{\lambda ^{\alpha }})=\frac{\beta k}{\alpha t}{p_{k}^{\alpha \beta }}(t).\]
\[ {T_{k}^{\psi ,f}}=\inf \{s\ge 0\hspace{0.1667em}:\hspace{0.1667em}N({H^{\psi }}({Y^{f}}(t)))\gt k\}\]
and write their distribution in the form
\[\begin{array}{r}\displaystyle \mathsf{P}\{{T_{k}^{\psi ,f}}\le t\}=\mathsf{P}\{N({H^{\psi }}({Y^{f}}(t)))\ge k\}={\sum \limits_{m=k}^{\infty }}{p_{m}^{\psi ,f}}(t)\\ {} \displaystyle ={\sum \limits_{m=k}^{\infty }}\frac{{(-\lambda )^{m}}}{m!}{\partial _{\lambda }^{m}}{\tilde{\ell }_{f}}(t,\psi (\lambda )),\end{array}\]
where for the last equality formula (25) was used. We note that the part of calculations in the proof of Theorem 2 in [1] still holds for the general case of the process ${N^{\psi ,f}}$ and allows to simplify the sum in the formula above to the expression
(33)
\[ \mathsf{P}\{{T_{k}^{\psi ,f}}\le t\}=\frac{{(-\lambda )^{k}}}{(k-1)!}{\partial _{\lambda }^{k-1}}\Big(\frac{{\tilde{\ell }_{f}}(t,\psi (\lambda ))}{\lambda }\Big).\]Consider a generalization of counting processes with Bernštein intertimes, which was introduced in the paper [8]:
\[ {N^{{\psi _{1}},{\psi _{2}},\dots ,{\psi _{n}}}}(t)=N\left({\sum \limits_{j=1}^{n}}{H^{{\psi _{j}}}}(t)\right),\hspace{1em}t\ge 0,\]
where ${H^{{\psi _{i}}}},i=1,\dots n$, are n independent subordinators with Bernštein functions ${\psi _{i}}$, independent of the Poisson process N. In [8] the following result was stated.Proposition 1 ([8]).
The distribution of the subordinated process ${N^{{\psi _{1}},{\psi _{2}},\dots ,{\psi _{n}}}}(t)$, $t\ge 0$, is the solution to the Cauchy problem
with the usual initial conditions
Consider the process ${N^{{\psi _{1}},{\psi _{2}},\dots ,{\psi _{n}},f}}(t)={N^{{\psi _{1}},{\psi _{2}},\dots ,{\psi _{n}}}}({Y^{f}}(t))$, $t\ge 0$.
Theorem 3.
The process ${N^{{\psi _{1}},{\psi _{2}},\dots ,{\psi _{n}},f}}$ has probability distribution function
\[ {p_{k}^{{\psi _{1}},{\psi _{2}},\dots ,{\psi _{n}},f}}(t)=\frac{{\left(-\lambda {\partial _{\lambda }}\right)^{k}}}{k!}{\tilde{\ell }_{f}}(t,{\sum \limits_{j=1}^{n}}{\psi _{j}}(\lambda )),\]
and the probabilities ${p_{k}^{{\psi _{1}},{\psi _{2}},\dots ,{\psi _{n}},f}}(t)$ satisfy the equation
with the usual initial conditions.
Example 1.
Consider the time-changed process
where ${N_{1}}(t)$ is the Poisson process with intensity ${\lambda _{1}}$, and
is the compound Poisson-Gamma subordinator with parameters λ, α, β, that is, with the Laplace exponent
\[ {\psi _{GN}}(u)=\lambda {\beta ^{\alpha }}\left({\beta ^{-\alpha }}-{(\beta +u)^{-\alpha }}\right).\]
In the case when $\alpha =1$ we have the compound Poisson-exponential subordinator, which we will denote as ${E_{N}}(t)$, and the corresponding time-changed process as ${N^{E}}(t)$. The detailed study of the processes ${N^{GN}}(t)$ and ${N^{E}}(t)$ was presented in [3, 4].We now consider the processes
The next theorem is obtained as a corollary of Theorem 2.
Theorem 4.
The process ${N^{GN,f}}(t)$ has the probability distribution function
\[ {p_{k}^{GN,f}}(t)=\frac{{\left(-{\lambda _{1}}{\partial _{{\lambda _{1}}}}\right)^{k}}}{k!}{\tilde{\ell }_{f}}(t,{\psi _{GN}}({\lambda _{1}})),\]
and ${p_{k}^{GN,f}}(t)$ satisfy the equation
\[\begin{aligned}{}{\mathcal{D}_{t}^{f}}{p_{k}^{GN,f}}(t)& =\left(\frac{\lambda {\beta ^{\alpha }}}{{({\lambda _{1}}+\beta )^{\alpha }}}-\lambda \right){p_{k}^{GN,f}}(t)\\ {} & \hspace{1em}+\frac{\lambda {\beta ^{\alpha }}}{{({\lambda _{1}}+\beta )^{\alpha }}}{\sum \limits_{m=1}^{k}}\frac{{\lambda _{1}^{m}}}{{\left({\lambda _{1}}+\beta \right)^{m}}}\frac{\Gamma (m+\alpha )}{m!\Gamma (\alpha )}{p_{k-m}^{GN,f}}(t).\end{aligned}\]
The process ${N^{E,f}}(t)={N_{1}}({E_{N}}({Y^{f}}(t)))$ has the probability distribution function
\[ {p_{k}^{E,f}}(t)=\frac{{\left(-{\lambda _{1}}{\partial _{{\lambda _{1}}}}\right)^{k}}}{k!}{\tilde{\ell }_{f}}(t,{\psi _{E}}({\lambda _{1}})),\]
where ${\psi _{E}}(u)=\frac{\lambda u}{\beta +u}$, and ${p_{k}^{E,f}}(t)$ satisfy the equation
4 Models of time-changed generalized counting processes
Consider the generalized counting process (GCP) $M(t)$, $t\ge 0$, which was introduced in [7]. The probabilities ${\tilde{p}_{n}}(t)=\mathsf{P}\left\{M(t)=n\right\}$ depend on k parameters ${\lambda _{1}},\dots ,{\lambda _{k}}$ and are given by the formula
where $\Omega (k,n)=\left\{({x_{1}},\dots ,{x_{k}})\hspace{0.1667em}:\hspace{0.1667em}{\textstyle\sum _{j=1}^{k}}j{x_{j}}=n,{x_{j}}\in {N_{0}}\right\}$, $\Lambda ={\textstyle\sum _{j=1}^{k}}{\lambda _{j}}$.
(34)
\[ {\tilde{p}_{n}}(t)={\sum \limits_{\Omega (k,n)}^{}}{\prod \limits_{j=1}^{k}}\frac{{\left({\lambda _{j}}t\right)^{{x_{j}}}}}{{x_{j}}!}{e^{-\Lambda t}},\hspace{1em}n\ge 0,\]GCP performs k kinds of jumps of amplitude $1,2,\dots ,k$ with rates ${\lambda _{1}},\dots ,{\lambda _{k}}$. Note that GCP comprises as particular cases such important for applications models as the Poisson process of order k and Pólya–Aeppli process of order k (see, e.g., [11, 10]).
The probability generating function of GCP is given by ([12])
The probabilities (34) satisfy
where ${z_{k}}={\textstyle\sum _{j=1}^{k}}{x_{j}}$ (see [13]).
\[ \frac{d{\tilde{p}_{n}}(t)}{dt}=-\Lambda {\tilde{p}_{n}}(t)+{\sum \limits_{j=1}^{\min \{n,k\}}}{\lambda _{j}}{\tilde{p}_{n-k}}(t),\hspace{1em}n\ge 0,\]
with the usual initial condition. The probabilities ${\tilde{p}_{n}}(t)$ can be also written as
(36)
\[ {\tilde{p}_{n}}(t)={\sum \limits_{\Omega (k,n)}^{}}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}{\left(-{\partial _{\Lambda }}\right)^{{z_{k}}}}{e^{-\Lambda t}},\hspace{1em}n\ge 0,\]Recently various models of time-changed GCP and fractional GCP were studied, in particular, with time-change given by a Lévy subordinators, inverse subordinators, including the cases of some specific subordinators, and also the fractional extensions of processes of the form $M({H^{\psi }}(t))$ (see, for example, [7, 10, 12, 11, 13], and references therein).
Following the lines of the previos section, we consider the time-changed process ${\mathsf{M}^{\psi ,f}}(t)=M({H^{\psi }}({Y^{f}}(t))$, that is, with double time-change by an independent subordinator ${H^{\psi }}$ and an inverse subordinator ${Y^{f}}$, which are independent of M. To the best of our knowledge, such general case has not been presented in the literature. In the next theorem we characterize its probabilities ${\tilde{p}_{n}^{\psi ,f}}(t)=\mathsf{P}\left\{{\mathsf{M}^{\psi ,f}}(t)=n\right\}$ and probability generating function.
Theorem 5.
The process ${\mathsf{M}^{\psi ,f}}$ has the probability distribution function
and ${\tilde{p}_{n}^{\psi ,f}}(t)$ satisfy the equation
with the initial conditions
and satisfies the equation
with ${\tilde{G}^{\psi ,f}}(u,0)=1$.
(37)
\[ {\tilde{p}_{n}^{\psi ,f}}(t)={\sum \limits_{\Omega (k,n)}^{}}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}{\left(-{\partial _{\Lambda }}\right)^{{z_{k}}}}{\tilde{\ell }_{f}}(t,\psi (\Lambda )),\hspace{1em}n\ge 0,\](38)
\[\begin{array}{r}\displaystyle {\mathcal{D}_{t}^{f}}{\tilde{p}_{n}^{\psi ,f}}(t)=-\psi (\Lambda ){\tilde{p}_{n}^{\psi ,f}}(t)-{\sum \limits_{m=1}^{n}}{\sum \limits_{\Omega (k,m)}^{}}{\psi ^{({z_{k}})}}\left(\Lambda \right){\prod \limits_{j=1}^{k}}\frac{{(-{\lambda _{j}})^{{x_{j}}}}}{{x_{j}}!}{\tilde{p}_{n-m}^{\psi ,f}}(t),\\ {} \displaystyle n\ge 0,\hspace{0.1667em}t\gt 0,\end{array}\]
\[ {\tilde{p}_{n}^{\psi ,f}}(0)=\left\{\begin{array}{l}1,\hspace{1em}n=0,\hspace{1em}\\ {} 0,\hspace{1em}n\ge 1.\hspace{1em}\end{array}\right.\]
The probability generating function of the process ${\mathsf{M}^{\psi ,f}}$ is of the form
(39)
\[ {\tilde{G}^{\psi ,f}}(u,t)={\tilde{\ell }_{f}}\Big(t,\psi \Big({\sum \limits_{j=1}^{k}}{\lambda _{j}}(1-{u^{j}})\Big)\Big),\hspace{1em}|u|\lt 1,\](40)
\[ {\mathcal{D}_{t}^{f}}{\tilde{G}^{\psi ,f}}(u,t)=-\psi \Big({\sum \limits_{j=1}^{k}}{\lambda _{j}}(1-{u^{j}})\Big){\tilde{G}^{\psi ,f}}(u,t)\]Proof.
For calculating ${\tilde{p}_{n}^{\psi ,f}}(t)$ we use the formula (36):
\[\begin{aligned}{}{\tilde{p}_{n}^{\psi ,f}}(t)& ={\int _{0}^{\infty }}{\tilde{p}_{n}}(s)P\left\{{H^{\psi }}({Y^{f}}(t))\in ds\right\}=\\ {} & ={\int _{0}^{\infty }}\sum \limits_{\Omega (k,n)}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}{\left(-{\partial _{\Lambda }}\right)^{{z_{k}}}}{e^{-\Lambda t}}P\left\{{H^{\psi }}({Y^{f}}(t))\in ds\right\}=\\ {} & =\sum \limits_{\Omega (k,n)}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}\mathsf{E}{\left(-{\partial _{\Lambda }}\right)^{{z_{k}}}}\exp \left\{-\Lambda {H^{\psi }}({Y^{f}}(t))\right\}=\\ {} & =\sum \limits_{\Omega (k,n)}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}{\left(-{\partial _{\Lambda }}\right)^{{z_{k}}}}\mathsf{E}\exp \left\{-\Lambda {H^{\psi }}({Y^{f}}(t))\right\}.\end{aligned}\]
This implies formula (37), since
\[ \mathsf{E}\exp \left\{-\Lambda {H^{\psi }}({Y^{f}}(t))\right\}=\mathsf{E}\exp \left\{-\psi (\Lambda ){Y^{f}}(t)\right\}={\tilde{\ell }_{f}}(t,\psi (\Lambda )).\]
To derive equation (38), we follow the same lines as in the proof of Theorem 2, but now, instead of equation (22) (or (23)) for probabilities of the process ${N^{\psi }}=N({H^{\psi }})$, we use the following equation for the probabilities ${\tilde{p}_{n}^{\psi }}(t)$ of the process ${M^{\psi }}=M({H^{\psi }})$ (see [13]):
Thus, we obtain
\[\begin{aligned}{}{\mathbb{D}_{t}^{f}}{\tilde{p}_{n}^{\psi ,f}}(t)& =\psi (\Lambda ){\tilde{p}_{n}^{\psi ,f}}(t)-{\sum \limits_{m=1}^{n}}\sum \limits_{\Omega (k,m)}{\psi ^{{z_{k}}}}(\Lambda ){\prod \limits_{j=1}^{k}}\frac{{(-{\lambda _{j}})^{{x_{j}}}}}{{x_{j}}!}{\tilde{p}_{n-m}^{\psi ,f}}(t)\\ {} & \hspace{1em}+{\tilde{p}_{n}^{\psi ,f}}(0){\nu _{f}}(t),\end{aligned}\]
and then apply the same reasonings as those after formula (29) to come to (38).Using the expression for the probability generating function ${\tilde{G}^{\psi }}(u,t)$ of the time-changed process ${M^{\psi }}(t)=M({H^{\psi }}(t))$ (see [12]),
\[ {\tilde{G}^{\psi }}(u,t)=\exp \Big\{-t\psi \Big({\sum \limits_{j=1}^{k}}{\lambda _{j}}(1-{u^{j}})\Big)\Big\},\]
we calculate ${\tilde{G}^{\psi ,f}}$ as
\[\begin{aligned}{}{\tilde{G}^{\psi ,f}}(u,t)& ={\int _{0}^{\infty }}{\tilde{G}^{\psi }}(u,s){\ell _{f}}(t,s)ds=\\ {} & ={\int _{0}^{\infty }}\exp \Big\{-t\psi \Big({\sum \limits_{j=1}^{k}}{\lambda _{j}}(1-{u^{j}})\Big)\Big\}{\ell _{f}}(t,s)ds=\\ {} & =\mathsf{E}{e^{-\psi \left({\textstyle\textstyle\sum _{j=1}^{k}}{\lambda _{j}}(1-{u^{j}})\right){Y^{f}}(t)}}={\tilde{\ell }_{f}}\Big(t,\psi \Big({\sum \limits_{j=1}^{k}}{\lambda _{j}}(1-{u^{j}})\Big)\Big),\end{aligned}\]
that is, (39) holds. Equation (40) follows in view of (39) and (10). □Example 2.
As a continuation of Example 1, consider GCP with time-change given by a compound Poisson-Gamma subordinator: ${M^{GN}}(t)=M({G_{N}}(t))$.
Probabilities of this process can be calculated and expressed in terms of special functions.
We will use the fact that the transition probability measure of the process ${G_{N}}(t)$ can be written in the closed form ([3])
therefore, the probability law of ${G_{N}}(t)$ has atom ${e^{-\lambda t}}$ at zero, that is, has a discrete part $P\left\{{G_{N}}(t)=0\right\}={e^{-\lambda t}}$, and the density of the absolutely continuous part is expressed in terms of the Wright function, more precisely, by means of its particular case defined as
(41)
\[\begin{aligned}{}P\left\{{G_{N}}(t)\in ds\right\}& ={e^{-\lambda t}}{\delta _{\{0\}}}(ds)+{\sum \limits_{n=1}^{\infty }}{e^{-\lambda t}}\frac{{\left(\lambda t{\beta ^{\alpha }}\right)^{n}}}{n!\Gamma (\alpha n)}{s^{\alpha n-1}}{e^{-\beta s}}ds\\ {} & ={e^{-\lambda t}}{\delta _{\{0\}}}(ds)+{e^{-\lambda t-\beta s}}\frac{1}{s}\Phi \left(\alpha ,0,\lambda t{(\beta s)^{\alpha }}\right)ds,\end{aligned}\]Using formula (41), we obtain:
\[\begin{aligned}{}{\tilde{p}_{n}^{GN}}(t)& =\mathsf{P}\left\{M({G_{N}}(t))=n\right\}={\int _{0}^{\infty }}\mathsf{P}(M(s)=n)\mathsf{P}({G_{N}}(t)\in ds)\\ {} & ={\int _{0}^{\infty }}{\sum \limits_{\Omega (k,n)}^{}}{\prod \limits_{j=1}^{k}}\frac{{\left({\lambda _{j}}s\right)^{{x_{j}}}}}{{x_{j}}!}{e^{-\Lambda s}}\\ {} & \hspace{1em}\times \Big\{{e^{-\lambda t}}{\delta _{\{0\}}}(ds)+{e^{-\lambda t-\beta s}}\frac{1}{s}\Phi \left(\alpha ,0,\lambda t{(\beta s)^{\alpha }}\right)\Big\}ds\\ {} & ={e^{-\lambda t}}{\int _{0}^{\infty }}{\sum \limits_{\Omega (k,n)}^{}}{\prod \limits_{j=1}^{k}}\frac{{\left({\lambda _{j}}s\right)^{{x_{j}}}}}{{x_{j}}!}{e^{-\Lambda s}}{e^{-\beta s}}\frac{1}{s}\Phi \left(\alpha ,0,\lambda t{(\beta s)^{\alpha }}\right)ds\\ {} & ={e^{-\lambda t}}{\sum \limits_{\Omega (k,n)}^{}}{\int _{0}^{\infty }}{\prod \limits_{j=1}^{k}}\frac{{\left({\lambda _{j}}s\right)^{{x_{j}}}}}{{x_{j}}!}{e^{-(\Lambda +\beta )s}}\frac{1}{s}{\sum \limits_{l=1}^{\infty }}\frac{{(\lambda t{\beta ^{\alpha }})^{l}}}{l!\Gamma (\alpha l)}{s^{\alpha l-1}}ds\\ {} & ={e^{-\lambda t}}{\sum \limits_{\Omega (k,n)}^{}}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}{\sum \limits_{l=1}^{\infty }}\frac{{(\lambda t{\beta ^{\alpha }})^{l}}}{l!\Gamma (\alpha l)}{\int _{0}^{\infty }}{e^{-(\Lambda +\beta )s}}{s^{{\textstyle\textstyle\sum _{j=1}^{k}}{x_{j}}+\alpha l-1}}ds\\ {} & ={e^{-\lambda t}}{\sum \limits_{\Omega (k,n)}^{}}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}{\sum \limits_{l=1}^{\infty }}\frac{{(\lambda t{\beta ^{\alpha }})^{l}}}{l!\Gamma (\alpha l)}\frac{\Gamma \left({z_{k}}+\alpha l\right)}{{\left(\Lambda +\beta \right)^{{z_{k}}+\alpha l}}}\\ {} & ={e^{-\lambda t}}{\sum \limits_{\Omega (k,n)}^{}}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}\frac{1}{{\left(\Lambda +\beta \right)^{{z_{k}}}}}{\sum \limits_{l=1}^{\infty }}\frac{{(\lambda t{\beta ^{\alpha }})^{l}}}{l!{\left(\Lambda +\beta \right)^{\alpha l}}}\frac{\Gamma \left({z_{k}}+\alpha l\right)}{\Gamma (\alpha l)},\end{aligned}\]
which can be written in the form
\[ {\tilde{p}_{n}^{GN}}(t)={e^{-\lambda t}}{\sum \limits_{\Omega (k,n)}^{}}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}\frac{1}{{\left(\Lambda +\beta \right)^{{z_{k}}}}}{\hspace{0.1667em}_{1}}{\Psi _{1}}\left(({z_{k}},\alpha ),(0,\alpha ),\frac{\lambda t{\beta ^{\alpha }}}{{\left(\Lambda +\beta \right)^{\alpha }}}\right),\]
where ${_{p}}{\Psi _{q}}$ is the generalized Wright function
\[ {_{p}}{\Psi _{q}}(({a_{i}},{\alpha _{i}}),({b_{j}},{\beta _{j}}),z)={\sum \limits_{k=0}^{\infty }}\frac{{\textstyle\prod \limits_{i=1}^{p}}\Gamma ({a_{i}}+{\alpha _{i}}k)}{{\textstyle\prod \limits_{j=1}^{q}}\Gamma ({b_{j}}+{\beta _{j}}k)}\frac{{z^{k}}}{k!}\]
defined for $z\in C$, ${a_{i}},{b_{i}}\in C$, ${\alpha _{i}},{\beta _{i}}\in R$, ${\alpha _{i}},{\beta _{i}}\ne 0$ and $\textstyle\sum {\alpha _{i}}-\textstyle\sum {\beta _{i}}\gt -1$ (see, e.g., [9]).If $\alpha =1$, that is, ${G_{N}}(t)={E_{N}}(t)$, then probabilities take the form
\[ {\tilde{p}_{n}^{E}}(t)={e^{-\lambda t}}{\sum \limits_{\Omega (k,n)}^{}}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}\frac{1}{{\left(\Lambda +\beta \right)^{{z_{k}}}}}{\mathcal{E}_{1,2}^{{z_{k}}+1}}\left(\frac{\lambda t\beta }{\Lambda +\beta }\right),\]
where ${\mathcal{E}_{\rho ,\delta }^{\gamma }}$ is the three-parameter generalized Mittag-Leffler function, which is defined as
\[ {\mathcal{E}_{\rho ,\delta }^{\gamma }}(z)={\sum \limits_{k=0}^{\infty }}\frac{\Gamma (\gamma +k)}{\Gamma (\gamma )}\frac{{z^{k}}}{k!\Gamma (\rho k+\delta )},\hspace{1em}z\in C,\hspace{2.5pt}\rho ,\delta ,\gamma \in C,\]
with $\hspace{2.5pt}\text{Re}(\rho )\gt 0$, $\text{Re}(\delta )\gt 0$, $\text{Re}(\gamma )\gt 0$ (see, e.g., [9]). For $n=0$ we have
\[ {\tilde{p}_{0}^{GN}}(t)={e^{-\lambda t}}+{e^{-\lambda t}}{\sum \limits_{l=1}^{\infty }}\frac{{\left(\lambda t{\beta ^{\alpha }}\right)^{l}}}{{\left(\Lambda +\beta \right)^{\alpha l}}}\frac{\Gamma (\alpha l)}{l!\Gamma (\alpha l)}=\exp \Big\{-\lambda t\Big(1-\frac{{\beta ^{\alpha }}}{{\left(\Lambda +\beta \right)^{\alpha }}}\Big)\Big\},\]
and, if $\alpha =1$,
Consider now the process ${M^{GN,f}}(t)=M({G_{N}}({Y^{f}}(t))$, $t\ge 0$. The next theorem follows from Theorem 5.
Theorem 6.
The process ${M^{GN,f}}$ has the probability distribution function given by formula (37) with $\psi ={\psi _{GN}}$ and the probabilities ${\tilde{p}_{n}^{GN,f}}(t)$ satisfy the equation
\[\begin{aligned}{}{\mathcal{D}_{t}^{f}}{\tilde{p}_{n}^{GN,f}}(t)& =\left(\frac{\lambda {\beta ^{\alpha }}}{{\left(\Lambda +\beta \right)^{\alpha }}}-\lambda \right){\tilde{p}_{n}^{GN,f}}(t)+\\ {} & \hspace{1em}+\frac{{\left(\lambda \beta \right)^{\alpha }}}{{\left(\Lambda +\beta \right)^{\alpha }}\Gamma (\alpha )}{\sum \limits_{m=1}^{n}}\sum \limits_{\Omega (k,m)}\frac{\Gamma \left({z_{k}}+\alpha \right)}{{\left(\Lambda +\beta \right)^{{z_{k}}}}}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}{\tilde{p}_{n-m}^{GN,f}}(t).\end{aligned}\]
The process ${M^{E,f}}$ has the probability distribution function given by (37) with $\psi ={\psi _{EN}}$ and the probabilities ${\tilde{p}_{n}^{E,f}}(t)$ satisfy the equation
\[\begin{aligned}{}{\mathcal{D}_{t}^{f}}{\tilde{p}_{n}^{E,f}}(t)& =-\lambda \frac{\Lambda }{\Lambda +\beta }{\tilde{p}_{n}^{E,f}}(t)+\\ {} & \hspace{1em}+\frac{\lambda \beta }{\Lambda +\beta }{\sum \limits_{m=1}^{n}}\sum \limits_{\Omega (k,m)}\frac{\Gamma \left({z_{k}}+1\right)}{{\left(\Lambda +\beta \right)^{{z_{k}}}}}{\prod \limits_{j=1}^{k}}\frac{{\lambda _{j}^{{x_{j}}}}}{{x_{j}}!}{\tilde{p}_{n-m}^{E,f}}(t).\end{aligned}\]