Mellin ’ s Transform and Application to Some Time Series Models

This paper uses the Mellin transform to establish the means, variances, skewness, and kurtosis of fuzzy numbers and applied them to the random coefficient autoregressive (RCA) time series models. We also give a close form expression to the moment generating function related to fuzzy numbers. It is shown that the results of the proposed time series models are consistent with those of the conventional time series models and that the developed concepts are straightforward and easily implemented.


Introduction
The use of fuzzy set theory as a methodology for modeling and analyzing certain financial problems is of particular interest to a number of researchers due to fuzzy set theory's ability to quantitatively and qualitatively model those problems which involve vagueness and imprecision.Fuzzy time series models provide a new avenue to deal with subjectivity observed in most financial time series models.Most of the fuzzy financial models developed so far have, generally, been confined to modeling parameters through some form of defuzzification or linear type of fuzzy numbers such as Trapezoidal Fuzzy Number (Tr.F.N.) or Triangular Fuzzy Number (T.F.N.).
There are some connections between fuzzy numbers based on convolution principle and some integral transformations such as the Mellin transform.The Mellin transform and its inverse are related to the two-sided Laplace transform.The Mellin transform is used often in applied sciences such as physics, engineering, and computer science because of its scale invariance property.In fact, this scale invariance property is analogous to the Fourier transform's shift invariance property.In particular, this transform is used for the analysis of linear time-invariant systems.This transform is considered as a transformation from the time-domain in which inputs and outputs are functions of time.Thus, the Mellin transform has significant applications in probability theory, Markov chains, renewal theory, and time series.The objective of this paper is to use the Mellin transform to establish the means, variances, skewness, and kurtosis of fuzzy numbers and then to apply to the random coefficient autoregressive (RCA) time series models.We also provide a close form expression to the moment generating function related to fuzzy numbers.
We summarize the preliminaries and notations in Section 1.The remainder of the paper is organized as follows.Section 2 introduces the Mellin transform to obtain statistical moments of any order.The useful Mellin transforms for fuzzy numbers are derived in Section 3. We also define the moment generating function of fuzzy numbers and provide some illustrative examples.In Section 4, we discussed the fuzzy RCA model under different innovation terms.Finally, concluding remarks are given in Section 5.
1.1.Preliminaries and Notation.The purpose of this section is to recall some concepts which will be used throughout the rest of the paper.We quote definitions and properties of fuzzy sets theory with relevant operations.Most of these related definitions and properties may be found in the following two references [1,2].In the sequel we will denote the family of fuzzy numbers by F. For any  ∈ F, we use the notation () = [ 1 (),  2 ()] for -level sets of .Definition 3. The crisp set of elements that belong to the fuzzy set  at least to the degree  is called the -level set (or cut); that is, Alternatively, a fuzzy set  is convex if its every -level sets is a convex set.Definition 5.  fuzzy set , which is both convex and normal, is defined to be a fuzzy number on the universal set R. Definition 6.The standard Gaussian fuzzy number can be represented by  = ( : , , ) with membership function as where  is called the mean and  the standard deviation.Alternatively, defining the interval of confidence at level  as () = [ 1 (),  2 ()], and setting exp[−(1/2)(( − )/) 2 ] = , we get  1 () =  − √(−2 ln ) and  2 () =  + √(−2 ln ).Thus, we get the -cut representation of the standard gaussian fuzzy number as Definition 7. A Triangular Fuzzy Number (T.F.N.) can be represented completely by a triplet  = ( 1 ,  2 ,  3 ), where Alternatively, defining the interval of confidence at level  as () = [ 1 (),  2 ()] and setting ( −  1 )/( As in Definition 1, we obtain the parametric representation of the Tr.F.N.  = ( 1 ,  2 ,  3 ,  4 ), where  1 <  2 <  3 <  4 ∈ , as In Section 2, we discuss the use of the Mellin transform to find moments of fuzzy numbers, and nonzero weighted central possibilistic moments.

Probability Density Function from Membership Function
In order to add fuzziness to data, Kaufmann and Gupta [2] obtained a membership function from a probability density function by way of a simple linear transformation, which indicates a possible conversion of a membership function into a density function.The conversion of a membership function into a probability density function can be made by one of two linear transformations, which will be elaborated in the next Section.The conversion of the membership function of a fuzzy number (), into a probability density function can be achieved by using one of the two linear transformations presented in Sheen [3].
(a) Proportional probability density function (ppdf), (b) Uniform probability density function (updf): where ℎ  and ℎ  are the conversion constants which ensure that the area under the continuous probability function is equals to one.Sheen [3] shows, that for the proportional conversion, the domain of the triangular ppdf is the same as the converted T.F.N., while its height changes from 1 to ℎ  = 2/( 3 −  1 ).In the case of the updf function, the height of the triangular updf is given by ℎ  = √2/( 3 −  1 ), while its vertexes are expressed as: . For a Tr.F.N., the height of the trapezoidal ppdf is given by ℎ  = 2/[( 4 +  3 ) − ( 2 +  1 )].Meanwhile, the height and the vertexes, respectively, of the trapezoidal updf are expressed as: ).It is found that the conversion of a fuzzy number into its corresponding proportional distribution is computationally straightforward.

The Mellin Transform and Ranking of Fuzzy Numbers.
The Mellin transform   () of a function (), where  is positive, is defined in Sheen [3] as follows: The Mellin transform has a unique one-to-one correspondence; that is, () ↔   ().
The function () is called the inverse the Mellin transformation of   ().Under certain conditions, inverse Mellin transformation may be represented as an integral.The tables of both the Mellin transforms and inverse Mellin transforms are given by Erdelyi [4].The Mellin transforms are virtually two-sided Laplace transforms and can be expressed either as exponential Fourier transforms or as combinations of Laplace transforms.For this reason, the Mellin transforms are found in several works on Fourier and Laplace transforms.The two most important sources are Doetsch [5] and Titchmarsh [6].For more application on the Mellin transform, we can refer to González-Gaxiola and Santiago [7] and Frontczak and Schöbel [8] for details.
The moments of a distribution represent the expected values of the power of a random variable with a () distribution.In other words, the Mellin transform   () = [ −1 ] provides an alternative method to establish a series of moments of a distribution if () is viewed as a probability density function.Comparing the first two moments of a distribution with the Mellin transform allows the mean, variance, and other moments to be expressed as follows; that is, Thus,

Moments Properties of Some of Fuzzy Numbers
In this Section, on the line of Sheen [3], we discuss moment of a T.F.N fuzzy numbers and derive expression for the kurtosis and skewness.
Theorem 9.For a T.F.N.,  one has the following easy to prove results.In what follows, one will assume that ).Then, we have the following results: , ) . ( Proof.In order to find ℎ  , the conversion constants which ensure that the area under the continuous probability function is equal to one, we proceed as follows: set solving for ℎ  ; we have The Mellin transform   () of the T.F.N., where  is positive, is as follows: To find the variance, we proceed as follows: thus, To compute the Skewness and Kurtosis, we have Thus, and the kurtosis is given by To find the skewness, we will need Thus, ) . (23)

ISRN Applied Mathematics
Remark 10.It is important to note here that, for the case of a Tr.F.N., we have the following easy to prove results: ) . ( From the previous results, we can easily compute the skewness and kurtosis of .
Example 11.Consider the following symmetric quadratic fuzzy number given by Then, we have the following easy to prove results (see Figure 1): ) , )) (26)

Moment Generating Function.
If random variable  has a m.g.f., then the th order moment of its distribution is equal to the value of the th order derivative of its m.g.f. for  = 0.
In what follows, we will discuss the mgf of a T.F.N. and Tr.F.N and discuss some of theirs properties: Figure 1: Nonlinear quadratic membership function.
The th moments can easily be computed as The expected value is given by We can easily define the centered moment around the mean  is () .
(30) Thus, we have and the variance is )) ; (32) thus, we have Var(), once we substitute for value of : Along similar lines we can also find the skewness and kurtosis as well.
Remark 12.The centered moments of a Tr.F.N is given by the following formula: (34)

Autoregressive Time Series Models
Theorem 13.Consider the following AR(1) Models, where ).Under the stationary conditions, || < 1, we have the following moments in a fuzzy sense: Proof.Consider The variance is computed as follows: ISRN Applied Mathematics 9 thus, and the fourth moment is The skewness and kurtosis of the process follow from expressions for Var(  ), ( 3  ) and ( 4  ), respectively, and are omitted.
Remark 14.In Theorem 9, if we assume that the error term can be model by a Tr.F.N., given by   = [, 2, 3, 4], then we can easily establish the following results: where   (2) = (5/2),   (3) = (20/3) 2 ,   (4) = (75/4) 3 ,   (5) = (826/15) 4 , and   (6) = (335/2) 5 : , 4.1.Random Coefficient Volatility Models.Consider the class of random coefficient autoregressive (RCA) models defined by allowing random additive perturbations of the autoregressive (AR) coefficients of ordinary AR models.That is, we assume that the process   is given by where the parameters   , ( = 2, . . ., ), are assumed to be known,   and   () are zero mean square integrable independent processes, and the variances are denoted by  2  and  2  .  (), ( = 1, 2, . ..)   are independent of   and  − and may be thought of as incorporating structural changes.In order to motivate nonlinear forecasts for nonlinear models, we consider a class of estimating functions of the form   = ∑  =2  − ℎ  as in Thavaneswaran and Heyde [9], where That is, the estimating function method can be used to obtain a nonlinear forecast for a nonlinear models by considering a class of elementary martingale estimating functions generated by nonlinear functions of the observations.Using a similar argument, we could also obtain forecasts for various classes of GARCH models; see Thavaneswaran and Heyde [9] for details.The main message is that RCA models could be used to improve the forecasting performance of stochastic volatility models.Appadoo et al. [11] derive the kurtosis of the correlated RCA model as well as the normal GARCH model under the assumption that the errors are correlated.Appadoo et al. [12] consider some volatility models with quadratic GARCH innovations and derive the kurtosis of the process.

RCA Models.
Random coefficient autoregressive time series were introduced by Nicholls and Quinn [10] and some of their properties have been studied recently by Thavaneswaran et al. [13].RCA models exhibiting long memory properties have been considered in Leipus and Sugailis [14].A sequence of random variables {  } is called an RCA(1) time series if it satisfies the equations where  denotes the set of integers and and for an AR(1) process  () reduces to 3.
(c) The autocorrelation of  2  is given by   2  = ( 2 +  2  )  and for an AR (1) process it turns out to be    . (47)

Conclusion
The authors feel that the approach developed in this paper may be applicable in other areas as well.There are several applications of the Mellin transformation in statistical inference.For example, Samanta [15] used the Mellin transformation (through a modified Bessel function tabulated in Abramowitz and Stegun [16]) to obtain the minimum variance unbiased estimation (MVUE) of a probability function.
The method is based on finding solution of a certain integral function.The results obtained using fuzzy assumptions in this paper are more flexible, and the methodology proposed in this paper may also be applicable for other time series models.The evidence in favor of the fuzzy approach supported by the Mellin transpose highlights the advantage of the model developed in this paper.

Definition 1 .
Fuzzy set  in  ⊂ R, the set of real numbers, is a set of ordered pairs  = {(, () :  ∈ }, where () is the membership function or grade of membership, or degree of compatibility or degree of truth of  ∈  which maps  ∈  on the real interval [0, 1].Definition 2. If Sup () = 1,  ∈ R, then the fuzzy set  is called a normal fuzzy set in R.