Physica A: Statistical Mechanics and its Applications
Fractional calculus approach to the statistical characterization of random variables and vectors
Introduction
Fractional moments are very useful in dealing with random variables with power-law distributions, , , where is the distribution function. Indeed, in such cases, moments exist only if and integer order moments greater than diverge. Distributions of this type are encountered in a wide variety of contexts; see the extensive literature in Refs. [1], [2] where power-law statistics appear in the framework of anomalous diffusion in many fields of applied science. To name but a few examples, we mention the travel length distribution in human motion patterns and animal search processes [3], [4], [5], [6], fluctuations in plasma devices [7], [8], scaling laws in polymer physics and its applications to gene regulation models [9], [10], as well as the distribution of time scales in processes such as the motion of charge carriers in amorphous semiconductors [11], tracer dispersion in groundwater [12], or sticking times in turbulent flows [13].
The lack of moments is of course a great limitation for the characterization of such distributions, because many methods of statistical analysis fail, like for example, the representation of the characteristic function or the log-characteristic function by moments or by cumulants, respectively. Then, it is necessary to resort to other distribution properties and a good option is given by the so-called fractional moments which have been extensively studied in the seventies and have recently attracted new interest both in theoretical [14], [15], [16] and in experimental settings [17]. In particular, in Refs. [15], [16] the moment problem originally stated in terms of integer moments, has been extended to fractional moments where it has been shown that the knowledge of some fractional moments improves significantly the convergence speed of the maximum entropy method for non-negative densities–which is defined in –if a proper optimization procedure is used in selecting the order . This can be explained by the non-local nature of fractional moments that will be better highlighted in the following.
Fractional moments of a non-negative random variable are expressible by the Mellin transform of the density and this fact has been widely used in literature principally in the field of the algebra of random variables. That is, the Mellin transform is the principal mathematical tool to handle problems involving products and quotients of independent random variables. As the approach presented in this paper is based on the Mellin transform, a brief outline of previous results is given in the following.
The first author who developed a systematic method to express the density of the product of independent random variables by the Mellin transform was Epstein. As consequence, in his pioneering work [18] he recognizes that any density function can be obtained as Mellin convolution of two independent densities. It must be remarked that such a result was already obtained by Ref. [19] outside the framework of integral transforms. Many other studies on the use of Mellin transform extending Epstein’s work have been reported providing explicit analytic forms for products of independent random variables with assigned densities. In Ref. [20] examples on product and quotient of independent random variables with Rayleigh distribution and moments of Rice distribution are obtained by Mellin transform tables. An extensive use of such concepts and the state-of-art on the algebra of random variables may be found in Ref. [21].
Another research direction on Mellin applications in probability is represented by the use of special functions like Fox’ H-functions and the Meijer’s G-functions, due principally to Mathai and co-workers. Such functions are indeed representable as Mellin–Barnes integrals of the product of gamma functions and are therefore suited to represent statistics of products and quotients of independent random variables whose fractional moments are expressible as gamma or gamma related functions. Applications of special functions to statistics and probability theory can be found in Refs. [22], [23], [24], [25], [26], [27] and references therein.
Whether fractional moments exists and how they are correlated to local properties of the characteristic function has been investigated in Refs. [28], [29], [30], [31]. In particular, in Ref. [28] a relation is derived between the absolute moment of real order and integrals involving the Marchaud fractional derivative of the characteristic function. Other relevant studies on the existence of fractional (absolute) moments can also be found in Refs. [32], [33], [34].
All the cited works share the common issue of using the Mellin transform because it naturally coincides with moments of the type if the random variable has a density defined only in the positive domain. In this paper, following the previous results in Ref. [35], it is shown that applying the Mellin transform to the characteristic function, and not to the density, a sound representation of the statistics of general multivariate random variables is possible. In this way, exploiting the Hermitian nature of the characteristic function, a unique representation valid for every random variable and vectors, defined in the whole real range is presented. Further, by means of simple properties of fractional operators, a generalization of the well-known link between moments and the derivative of the characteristic function evaluated in zero is provided. Such a link is useful because it allows one to find existence criteria for fractional moments that can be derived in a simpler way with respect to Ref. [28] and clarifies the non-local nature of fractional moments. Moreover, this approach that starts from the characteristic function plainly applies to random variables defined in the whole domain without the need of partitioning. Such theoretical aspects go along with a better numerical treatment of the inverse Mellin transform involved in this approach as it will be stressed in the following.
In the application section, it will be shown that the proposed approach is extremely useful in the statistical characterization of raw data and in the representation of both random variables and vectors with small computational effort. Thus, the good numerical convergence makes the proposed approach a good and reliable tool also for data analysis. We want to stress that the representation proposed in this paper can be applied to any Fourier pair functions used in statistical analysis and, very recently, applications to the representation of stochastic processes have been proposed in Refs. [36], [37].
Section snippets
Probabilistic characterization of random variables
A brief outline of the basic tools of probability will serve to frame the problem in hand and to introduce symbols. Let be a probability space, where is a collection of events defining a -algebra of subsets on the sample space , and a probability measure that assigns a number to each event on and a random variable, where is a Borel set. is the cumulative distribution and is the probability density function (PDF). The Fourier transform of
Extension to multivariate random vectors
In this section, the concepts briefly outlined in the previous section will be derived in case of multivariate random vectors. Multivariate random vectors are fully characterized in probabilistic setting by the joint probability density function, or by its spectral counterpart, namely the joint characteristic function; alternatively also by joint moments or by joint cumulants of every order. Let be a probability space, , that is denotes a -dimensional random
Applications
In this section it will be shown how the previous results can be used to characterize probability distributions. The integral forms in the monovariate case, Eqs. , (14), and in the multivariate case, Eqs. (29), (31), will be discretized in a suitable form. A remarkable effect of having applied the Mellin transform to the characteristic function rather than to the density is the presence of the gamma function in the integrand. Indeed, by its fast decay with increasing imaginary part of the
Conclusions
In this paper, we provide new insights into the problem of the statistical representation of random variables by means of the Mellin transform and fractional calculus. We define a class of complex moments (comprehending fractional moments) that can represent every kind of joint PDF and joint characteristic function, including joint characteristic functions that are non-differentiable in zero, and densities with heavy tails, like the bivariate Cauchy distribution or other power-law
References (40)
- et al.
The random walk’s guide to anomalous diffusion: A fractional dynamics approach
Phys. Rep.
(2000) - et al.
On the use of fractional calculus for the probabilistic characterization of random variables
Probabilistic Engineering Mechanics
(2009) - et al.
Some multidimensional fractional integral operators involving a general class of polynomials
Journal of mathematical analysis and applications
(1995) - et al.
The restaurant at the end of the random walk: Recent developments in the description of anomalous transport by fractional dynamics
Journal of Physics A: Mathematical and General
(2004) - et al.
The scaling laws of human travel
Nature
(2006) - et al.
Understanding individual human mobility patterns
Nature
(2008) - et al.
Lévy strategies in intermittent search processes are advantageous
Proceedings of the National Academy of Sciences of the United States of America
(2008) - et al.
Lévy walk patterns in the foraging movements of spider monkeys (ateles geoffroyi)
Behavioral Ecology and Sociobiology
(2004) - et al.
Fractional kinetics for relaxation and superdiffusion in a magnetic field
Physics of Plasmas
(2002) Chaotic transport in zonal flows in analogous geophysical and plasma systems
Physics of Plasmas
(2000)
Statistical mechanics of polymer networks of any topology
Journal of Statistical Physics
Optimal target search on a fast-folding polymer chain with volume exchange
Physical Review Letters
Anomalous transit-time dispersion in amorphous solids
Physical Review B
The dynamical foundation of fractal stream chemistry: The origin of extremely long retention times
Geophys. Res. Lett.
Observation of anomalous diffusion and lévy flights in a two-dimensional rotating flow
Physical Review Letters
Maxentropic solution of fractional moment problems
Applied Mathematics and Computation
Hausdorff moment problem via fractional moments
Applied Stochastic Models in Business and Industry
Stieltjes moment problem via fractional moments
Applied Mathematics and Computation
Characterization of a benzoic acid modified glassy carbon electrode expressed quantitatively by new statistical parameters
Physica E: Low-dimensional and Nanostructures
Some applications of the mellin transform in statistics
The Annals of Mathematical Statistics
Cited by (41)
Harmonic transform-based density estimation method for uncertainty propagation and reliability analysis involving multi-modal distributions
2023, Mechanical Systems and Signal ProcessingGeneralized distribution reconstruction based on the inversion of characteristic function curve for structural reliability analysis
2023, Reliability Engineering and System SafetyCitation Excerpt :The basic ideas in the proposed method are illustrated in Fig. 1. Actually, the complex fractional moments, which are the extensions to integer moments and fractional moments and never be divergent, are applied to approximate the CF curve [35–40], and then the inverse Fourier transform is employed to obtain the probability density function (PDF) by inversion of CF. The rest of the paper is arranged as follows. In Section 2, the problem formulation is presented.
An improved complex fractional moment-based approach for the probabilistic characterization of random variables
2018, Probabilistic Engineering MechanicsCross-correlation and cross-power spectral density representation by complex spectral moments
2017, International Journal of Non-Linear MechanicsCitation Excerpt :The capabilities of the complex moments have been shown in several articles. Precisely, they are applied for the characterization of random variables [9–11], for the digital simulation of the random processes [12,13], for the approximated solution of the Fokker-Planck equation [14,15], Kolmogorov-Feller equation [16], for the evaluation of the stationary correlation function of fractional-order oscillator [17], and for a new way to perform the wavelet analysis [18]. The present paper aims to extends the application of the method described in [3] in order to provide a complete characterization of stochastic response processes when the problem is multivariate stochastic vector process.
Efficient reliability assessment of structural dynamic systems with unequal weighted quasi-Monte Carlo Simulation
2016, Computers and Structures