arxiv: v1 [math.pr] 23 Jan 2018

Transcription

1 TRANSFER PRINCIPLE FOR nt ORDER FRACTIONAL BROWNIAN MOTION WIT APPLICATIONS TO PREDICTION AND EQUIVALENCE IN LAW TOMMI SOTTINEN arxiv: v1 [math.pr 3 Jan 18 Department of Mathematics and Statistics, University of Vaasa, P.O. Box 7, FIN-6511 Vaasa, FINLAND LAURI VIITASAARI Department of Mathematics and Statistics, University of elsinki, elsinki, P.O. Box 68, FIN-14 University of elsinki, FINLAND Abstract. The nth order fractional Brownian motion was introduced by Perrin et al. [13. It is the (upto a multiplicative constant) unique self-similar Gaussian process with urst index (n 1, n), having nth order stationary increments. We provide a transfer principle for the nth order fractional Brownian motion, i.e., we construct a Brownian motion from the nthe order fractional Brownian motion and then represent the nthe order fractional Brownian motion by using the Brownian motion in a nonanticipative way so that the filtrations of the nthe order fractional Brownian motion and the associated Brownian motion coincide. By using this transfer principle, we provide the prediction formula for the nthe order fractional Brownian motion and also a representation formula for all the Gaussian processes that are equivalent in law to the nth order fractional Brownian motion. 1. Introduction The fractional Brownian motion is probably the most well-known generalisation of the Brownian motion. It is a centred Gaussian process B, depending on the urst parameter (, 1), that is -self-similar and has stationary increments. In the centred Gaussian case these two properties characterize the process completely (upto a multiplicative constant). Nowadays the properties and stochastic analysis with respect to the fractional addresses: tommi.sottinen@iki.fi, lauri.viitasaari@iki.fi. Date: January 4, Mathematics Subject Classification. Primary: 6G; Secondary: 6G15, 6G5, 6G35, 699. Key words and phrases. fractional Brownian motion, stochastic analysis, transfer principle, prediction, equivalence in law. 1

2 SOTTINEN AND VIITASAARI Brownian motion are studied extensively and widely understood. We refer to Mishura [9 for details of fractional Brownian motion and its stochastic analysis. The nth order fractional Brownian motion was introduced by Perrin et al. [13. Their motivation was to extend the self-similarity index of the standard fractional Brownian motion beyond the limitation < 1. That is, the nth order fractional Brownian motion B with urst parameter (n 1, n) is -self-similar, and in the case n 1 one recovers the fractional Brownian motion. Moreover, the nth order fractional Brownian motion is n-stationary, meaning that its nth differences are stationary. As in the case of the fractional Brownian motion, these two properties: -self-similarity and n-stationarity characterize the process in the centred Gaussian world (upto a multiplicative constant). Perrin et al. [13 defined the nth order fractional Brownian motion by using the socalled Mandelbrot Van Ness representation [8 of the fractional Brownian motion. owever, while the kernel in the Mandelbrot Van Ness representation is rather simple, the representation requires the knowledge to the infinite past of the generating Brownian motion, which is very problematic in many applications. In this article, we propose to define the nth order fractional Brownian motion by using the Molchan Golosov representation [11 of the fractional Brownian motion. As such, we obtain a compact interval representation with the minor cost of having a slightly more complicated kernel. In addition, we show that the filtrations of the nth order fractional Brownian motion and the standard Brownian motion generating it coincide. This is very important for filtering and prediction. In addition, we provide a transfer principle for the nth order fractional Brownian motion which can be used to develop stochastic calculus in an elementary and simple way by using the stochastic calculus with respect to the standard Brownian motion. The rest of the paper is organised as follows. In Section we recall some elementary definitions and preliminaries. We also recall the transfer principle with respect to the fractional Brownian motion. In section 3 we present our main results, which we apply to the equivalence of laws problem and prediction problem in Section 4. We end the paper with a simulation study in Section 5 and some concluding remarks in Section 6.. Preliminaries We begin by recalling the concept of -self-similarity. Definition.1. Let X (X t ) t be a stochastic process and >. The stochastic process X is called -self-similar if for any a > we have X at law a X t, t, where law denotes the equality of the finite dimensional distributions. Definition.. Denote an increment of size l of a function g by l g(x) g(x + l) g(x).

3 nt ORDER FBM 3 For k 1, the kth order increment of size l is defined recursively by (k) l g(x) (k 1) where the th-order increment is () l g(x) g(x). l g(x + l) (k 1) l g(x). Example.1. The nd order increment of size l, () l g(x), is given by () l g(x) g(x + l) g(x + l) + g(x) which corresponds to the usual second order increment. Definition.3. A stochastic process X is called n-stationary if for every l > the process (k) l X is non-stationary for k, 1,..., n 1, but the process l X is stationary. Remark.1. Note that the definition of n-stationary is simply a continuous analogue of the classical notion of difference stationarity, used e.g. in series analysis and ARIMA-processes [5. Recall that the fractional Brownian motion B with urst index (, 1) is a centred Gaussian process having the covariance function r (t, s) 1 [ t + s t s, t, s R. (.1) It is easy to check that the process B is not stationary, but has stationary increments, i.e., the fractional Brownian motion is 1-stationary. It is also easy to check that the process B is -self-similar, and that the valid range for is (, 1. The case 1 is obviously a degenerate one: B (t) tb (1). The fractional Brownian motion B is connected to the standard Brownian motion W via integral representations. We begin by recalling the Mandelbrot Van Ness representation [8 of the fractional Brownian motion B on the real line with urst index (, 1): B (t) where W is a standard Brownian motion and the kernel is g (t, u) 1 Γ( + 1 ) g (t, u) dw (u), t R, (.) [ (t u) 1 ( u) 1 + ere x + : x : max(x, ) and Γ denotes the Gamma function: Γ(x) t x 1 e t dt. Remark.. We remark that (.) defines the fractional Brownian motion on the whole real line, but the filtrations of the fractional Brownian motion B and the generating Brownian motion W coincides only for t >. This is problematic for many applications..

4 4 SOTTINEN AND VIITASAARI For many practical applications it is useful to represent a process as an integral with respect to a Brownian motion on a compact interval. In the case of the fractional Brownian motion, one can use the so-called invertible Molchan-Golosov representation (see, for example, [1, 14). Proposition.1 ([1, 11). The Molchan Golosov representation of the fractional Brownian motion B with (, 1) is where k (t, s) d [ ( t s with normalising constant B (t) k (t, s) dw (s), t, (.3) ) 1 ( (t s) 1 1 ) t s 1 z 3 (z s) 1 dz, (.4) s d Γ( 3 ) Γ( + 1 )Γ( ). Moreover, the representation (.3) is invertible, i.e., the Brownian motion W is constructed from the fractional Brownian motion B by where k ( 1) (t, u) d with normalising constant W (t) k ( 1) (t, u) db (u), (.5) [ ( ) t 1 ( 1 (t s) 1 ) t s 1 z 3 (z s) 1 dz, s s d ere B denotes the Beta function: B ( 1, + 1 Γ( + 1 )Γ( ) ) ( + 1 B(a, b) Γ(a)Γ(b) Γ(a + b). ) Γ( 1 ). Remark.3. We emphasise the invertible role of the Molchan Golosov representation. We have that the filtrations F W (t) σ{w (s), s t} and F B (t) σ{b (s), s t} coincide. This is of paramount importance in, for example, prediction and filtering. Remark.4. The Wiener integral in (.5) can be understood as a Riemann Stieltjes integral, see [1. More generally it can be understood by using a so-called transfer principle, as will be explained the the next subsection below.

5 nt ORDER FBM 5.1. Transfer principle for the fractional Brownian motion. In this subsection we recall the Wiener integral and the transfer principle for the fractional Brownian motion. For this purposes we consider the fractional Brownian motion B on some compact interval [, T. Definition.4 (Isonormal process). The isonormal process associated with the fractional Brownian motion B (B (t), t [, T ) with the urst index (, 1) is the centred Gaussian family (B (h), h ), where the ilbert space ([, T ) is generated by the covariance r given by (.1) as follows: (i) indicators 1 t : 1 [,t), t T, belong to. (ii) is endowed with the inner product 1 t, 1 s : r (t, s), and the centred Gaussian family is then defined by the covariance [ E B (h)b ( h) h, h. Definition.4 states that B (h) is the image of h in the isometry that extends the relation B (1 t ) : B (t) linearly. This gives rise to the definition of Wiener integral with respect to the fractional Brownian motion. Definition.5 (Wiener integral). B (h) is the Wiener integral of the element h with respect to B, and it is denoted by T h(t) db (t). Remark.5. Due to the completion under the inner product, it may happen that the space is not a space of functions, but contains distributions. Indeed, it was shown by Pipiras and Taqqu [14 that the space is a proper function space only in the case 1, while it contains distributions for > 1. We also note that in the case of a standard Brownian motion, that is 1, we have 1 L ([, T ). The kernel representation of the fractional Brownian motion can be used to provide a transfer principle for Wiener integrals. In other words, the Wiener integrals with respect to the fractional Brownian motion can be transferred to a Wiener integrals with respect to the corresponding Brownian motion. To do this, we define a dual operator k on L ([, T ) associated with the kernel k by extending the relation linearly and closing it in. k 1 t k (t, ) (.6) Remark.6. For general (, 1) and step functions f we have (k f)(t) k (T, t)f(t) + T t [f(u) f(t) k (u, t) u du. (.7)

6 6 SOTTINEN AND VIITASAARI If, moreover, ( 1, 1), then (.7) can be simplified into (k f)(t) since for > 1 we have k (t, t). T t f(u) k (u, t) u Remark.7. The dual operator k depends on the interval [, T, even though the kernel k does not. This is the reason we have to consider the ilbert space ([, T ). Theorem.1 ([1). Let k be given by (.6). Then k and L ([, T ). Moreover, for any f we have T f(u) db (u) T du, (k f)(u) dw (u). 3. The nth order fractional Brownian motion provides an isometry between Perrin et al. [13 defined the nth order fractional Brownian motion by using the Mandelbrot Van Ness [8 representation of the fractional Brownian motion B with urst index (, 1): Definition 3.1 ([13). Let W (W (u), u R) be the two-sided standard Brownian motion and let Γ be the gamma function. The nth order fractional Brownian motion with urst index (n 1, n) is defined as where B (t) g (t, u) 1 Γ( + 1 ) ( 1 g (t, u) dw (u), t R, (3.1) [ (t u) 1 ( u) 1 + ) ( n 3 ) ( u) n+ 1 + t (n 1). (3.) (n 1)! Remark 3.1. In the case n 1, the classical fractional Brownian motion is recovered, as Definition 3.1 reduces to the Mandelbrot Van Ness representation of the fractional Brownian motion. Remark 3.. Some properties of the nth order fractional Brownian motion provided in [13 are (i) B (ii) B (iii) B is n-stationary, is -self-similar, has n 1 s continuously differentiable paths. In particular, d dt B (t) B(n 1) 1 (t). (3.3) Properties (i) and (ii) are unique to the nth order fractional Brownian motion in the class of centred Gaussian processes, i.e., they can be used as a qualitative definition. Property (iii) follows from the properties (i) and (ii).

7 nt ORDER FBM 7 The covariance function r (t, s) E[B (t)b (s) of the nth order fractional Brownian motion is r (t, s) ( 1) C t s (n 1) ( ( 1) j j j1 ) [ ( ) t j ( s ) j s + t s t, where and For the variance we have ( ) α j α(α 1) (α (j 1)), j! C 1 Γ( + 1) sin(π). [ ( ) V B (t) C 1 t. n 1 The following theorem provides an invertible Volterra representation on compact interval for the nth order fractional Brownian motion with respect to a standard Brownian motion generated from it. Theorem 3.1. Let n 1 be an integer and let (n 1, n). Let W (W (u), u ) be a one-sided standard Brownian motion. Define a sequence of Volterra kernels recursively as Then B k (1) (t, u) k (t, u), (3.4) (t, u) k (n 1) 1 (s, u) ds. (3.5) (t) u (t, u) dw (u), t, (3.6) defines an nth order fractional Brownian motion. Moreover, the Brownian motion W can be recovered from B by W (t) k ( 1) n+1 (t, u) d dn 1 du B n 1 (u). (3.7) In particular, the filtrations of W and B coincide: F W (t) F B (t) for all t. Proof. The proof is by induction. The case n 1 corresponds to the case of the fractional Brownian motion and follows from Proposition.1.

8 8 SOTTINEN AND VIITASAARI Assume now that the claim is valid for some k n 1. For k n we can, by square integrability of the kernels, apply stochastic Fubini theorem to have B (t) (t, u) dw (u) [ u [ s k (n 1) 1 (s, u) ds dw (u) k (n 1) 1 (s, u) dw (u) ds B (n 1) 1 (s) ds. Thus, (3.6) defines the nth order fractional Brownian motion by (3.3) together with the fact that B (). Furthermore, the claim (3.7) follows directly from (3.6) together with (3.3). Finally, the equivalence of filtrations follows from Proposition.1 together with the observation that as B () for all and n 1, the filtrations of B(n 1) 1 and B coincide by (3.3). Remark 3.3. As n increases, so does the smoothness of the paths of B. Since W is not smooth, the representation (3.6) implies that the kernels (t, u) have to become increasingly smooth in t as n increases. This is also obvious from (3.4) (3.5). Remark 3.4. Let us consider writing the inversion formula (3.7) as W (t) k ( n) (t, u) db (u). (3.8) Since B is smooth for n > 1, and W is not, the kernel k( n) (t, u) in (3.8) must be non-smooth in t. Actually it is a Schwarz kernel, i.e., a proper distribution. For example, if n, then we have k ( ) (t, u) k( 1) k( 1) 1 1 (t, u)δ(t u) (t, u), u where δ is the Dirac s delta function at point, and the partial derivative k( 1) 1 u (t, u) has to be understood in the sense of distributions. With the help of Theorem 3.1 we can provide the transfer principle for the nth order fractional Brownian motion. In what follows, the Wiener integral with respect to the nth order fractional Brownian motion is defined in the spirit of Definition.4 as follows. Define a ilbert space ([, T ) such that: (i) indicators 1 t : 1 [,t), t T, belong to (ii). is endowed with the inner product 1 t, 1 s is given by (3.4). covariance r : r (t, s), where the

9 Then B space nt ORDER FBM 9 (B (t), t [, T ) is the isonormal process associated with the ilbert ([, T ) and the Wiener integral T f(s) db (s) B (f) of f Gaussian random variable, and for f, g with respect to the nth order fractional Brownian motion B [ T E f(s) db (s) T we have g(s) db (s) f, g. is a centred The following provides a transfer principle for the nth order fractional Brownian motion, in the spirit of Theorem.1. Theorem 3. (Transfer principle). Let B be the nth order fractional Brownian motion on [, T with n 1 and (n 1, n). Define an operator : L ([, T ) by linearly extending 1 t (t, ). Then provides an isometry between we have T f(u) db (u) T and L ([, T ). Moreover, for any f ( f ) (u) dw (u). Furthermore, if n, then L ([, T ), and for any f L ([, T ) we have ( f ) (u) T f(t)k (n 1) 1 (t, u) dt. (3.9) Proof. The first part follows from the similar arguments as in the general case (see [18). For the reader s convenience we present the main arguments. Assume first that f is an elementary function of form f(t) a k 1 Ak k1 for some disjoint intervals A k (t k 1, t k. Then the claim follows by the very definition of the operator and Wiener integral with respect to B together with representation (3.6), and this shows that provides an isometry between and L ([, T ). ence can be viewed as a closure of elementary functions with respect to f f L ([,T ) which proves the claim. In order to complete the proof, we need to verify (3.9). For this, assume again that f is an elementary function. Then it is straightforward to check that (see also [18, Example 4.1) ( ) T f (u) f(t) k (t, u) dt, t

10 1 SOTTINEN AND VIITASAARI and thus, thanks to (3.4) (3.5), we have (3.9) for all elementary f. Let now f L ([, T ) and take a sequence f j of elementary functions such that f j f L ([,T ). Then (f j, j N) is a Cauchy sequence in L ([, T ). We have and thus f j ( f j ) (u) T f j L T ( T T T C T,,n f j L ([,T ). f j (t)k (n 1) 1 (t, u) dt, f j (t)k (n 1) 1 (t, u) dt ) du ( k (n 1) 1 (t, u) ) dt du fj L ([,T ) The same argument applied to f j f i shows that (f j, j N) is also a Cauchy sequence in. Thus L ([, T ) and (3.9) holds. Remark 3.5. We stress that for n 1 the equation (3.9) does not hold. On the other hand, in this case we also have L ([, T ) (1) for > 1, while for < 1 we have (1) L ([, T ). Finally, we also remark that, by the characterisation of the space (1) (1) in [14 one can characterise the space simply by using Fubini s theorem and (3.4) (3.5). Remark 3.6. The transfer principle provided in Theorem 3. extends in a straightforward manner to multiple Wiener integrals. Moreover, the transfer principle can be used as a simple approach to stochastic analysis and Malliavin calculus with respect to the nth order fractional Brownian motion. For the details on the topic, we refer to [ Applications 4.1. Equivalence in law. We first investigate the equivalence of law problem. For the treatment of the problem in the classical fractional Brownian motion or sheet case, see [16, 17. Gaussian processes are either singular or equivalent in law. By itsuda s representation theorem [6 a Gaussian process W ( W (t), t [, T ), is equivalent in law to a Brownian motion W (W (t), t [, T ) if and only if there there exists a Volterra kernel b L ([, T ) and a function a L ([, T ) such that W (t) W (t) s b(s, u)dw (u)ds + ere the Brownian motion W is constructed from W as W (t) W (t) s b (s, u)(d W (u) a(u)du)ds, a(s) ds. (4.1) a(s) ds.

11 nt ORDER FBM 11 where b is the unique resolvent Volterra kernel of b solving the equation s b (t, u)b(u, s)du b (t, s) + b(t, s) s b(t, u)b (u, s)du. The resolvent kernel can be constructed by using Neumann series, see [15 for details. The log-likelihood ratio of model W over W is l(t) log d P F t dp [ s b(s, u) dw (u) + a(s) dw (s) 1 [ s b(s, u) dw (u) + a(s) ds. (4.) B Consider then a Gaussian process ( B (t), t [, T ). This process is equivalent to the nth order fractional Brownian motion if and only if the process W (t) k ( 1) n+1 (t, u) d dn 1 du B n 1 is equivalent to a Brownian motion. From (4.1) and (3.6) it follows that B (t) (t, s)d W (s) B (t) s (t, s) b(s, u) dw (u) ds + Let us collect the discussion above as a theorem: B (u). (4.3) (t, s)a(s) ds. (4.4) Theorem 4.1. A Gaussian process is equivalent in law to an nth order fractional Brownian motion on [, T if and only if there exists a Volterra kernel b L ([, T ) and a function a L ([, T ) such that admits the representation (4.4), where W is B B constructed from by (4.3) and W is a Brownian motion connected to W via (4.1). The log-likelihood ratio of B over B is given by (4.). Proof. By (3.7), the process motion if and only if the process W (t) B is equivalent to the nth order fractional Brownian k ( 1) n+1 (t, u) d dn 1 du B n 1 is equivalent to a Brownian motion. From (4.1) and (3.6) it follows that B (t) (t, s)d W (s) B (t) s (t, s) b(s, u) dw (u) ds + (u). (4.5) (t, s)a(s) ds (4.6)

12 1 SOTTINEN AND VIITASAARI which concludes the proof. B Remark 4.1. For transformations B on [, T we have A(t) (t) B (t) + A(t) that are equivalent in law to (t, s)a(s) ds for some a L ([, T ). Consequently, A is + 1 d j A() for all j < n. Otherwise, in principle, the function A can be filtered out with dt j probability one given continuous data on any interval [, T with T >. In particular, s fractionally differentiable with it follows that the drift α in signal X(t) B (t) + αt for n, can be completely determined from continuous observations on any interval [, T. Indeed, α d dt X(). 4.. Prediction of the nthe order fractional Brownian motion. The equivalence of filtrations generated by the Brownian motion and the nth order fractional Brownian motion is of uttermost important in prediction. Indeed, this guarantees that the prediction of the nth order fractional Brownian motion can be done by using the same approach as taken in [19. Theorem 4.. The regular conditional distribution of the nth order fractional Brownian motion B (B (t), t [, T ) conditioned on the information F B (u) (v), v u}, is Gaussian process with random mean ˆB ( u) given by σ{b u [ ˆB (t u) B (u) + and a deterministic covariance ˆr (, u) given by ˆr (t, s u) r (t, s) u (t, v) k (u, v) dw (v) (4.7) (t, v)k (s, v) dv. (4.8) Proof. By the Gaussian correlation theorem (see Janson [7) the conditional law is Gaussian with mean [ t E B (t) FB (u) and a covariance function [( (t, s) ˆr (t, s u) E B B (t) ˆB (t u) ) ( B ) (s) ˆB (s u) FB (u). We start by proving equation (4.7). Since the nth order fractional Brownian motion admits the representation (3.6) and the filtrations of B W coincide, the prediction mean of B prediction mean under the observations F W (u): ˆB (t s) E [ B E given observations F B (t) FB [ B (t) FW (u) (u) and the Brownian motion (u) is the same as the

13 nt ORDER FBM 13 By using (3.6) and the independence of Brownian increments we obtain ˆB (t u) E [ [ u E [ u E u (t, v) dw (v) FW (u) (t, v) dw (v) + u (t, v) dw (v) FW (u) [ (t, v) dw (v) FW (u) (t, v) dw (v). This formula can also be written as ˆB proving (4.7). (t u) u u (t, v) dw (v) (u, v) dw (v) + u B (u) + u [ [ + E u (t, v) dw (v) (t, v) k (u, v) dw (v) (t, v) k (u, v) dw (v), It remains to prove (4.8). Proceeding similarly, the conditional covariance can be calculated as ˆr (t, s u) [( E [ E B u min(t,s) u ) ( (t) ˆB (t u) (t, v) dw (v) s u (t, v)k (t, v) dv. ) B (s) ˆB (s u) FB (s, v) dw (v) FW (u) (u) Since r (t, s) min(t,s) (t, u)k (s, u) du, formula (4.8) follows from this. This concludes the proof. Remark 4.. By using a the transfer principle, it is possible to write (4.7) as ˆB (t u) B where Ψ (t, u, ) is a Schwarz kernel. (u) + u Ψ (t, u, v) db (v),

14 14 SOTTINEN AND VIITASAARI Example 4.1. The prediction mean for a transformation f(b (t)) given the observation (u) can be calculated as F B E [ f(b (t)) F B 1 πˆr (t, t u) (u) f(x) exp 1 ( x ˆB n (t u) ˆr (t, t u) ) dx. (4.9) (If the transformation f in injective, then (4.9) is also the prediction under the observations σ{f(b (v)), v u}.) More generally, for f(b (t)) f(b (t 1),..., B (t k)) we have where [ E f(b (t)) F B (u) (π) k/ ˆr (t u) 1/ f(x) R { k exp 1 } (x ˆB n (t u)) ˆr (t u)( 1) (x ˆB n (t u)) dx, (4.1) ˆr (t u) [ˆr (t i, t j u) k i,j1, ˆB [ (t) ˆB (t i) k i1. Remark 4.3. The no-information and full-information asymptotics of the conditional covariance (4.8) can be computed similarly as in [19, Proposition 3. and Proposition 3.3. Actually, one can simply use the result for the standard fractional Brownian motion, [19, Proposition 3. and Proposition 3.3 and then apply (3.4) (3.5) together with the induction to obtain asymptotic expansions. Indeed, this follows directly from the Fubini s theorem. We leave the details to the reader. 5. Simulations In this section we illustrate the smoothening effect by simulating the paths of nth order fractional Brownian motion for different values of the urst index and for values of n 1,, 3, 4. In all of the simulations we have used interval [, 1 with grid points, and the path of the fractional Brownian motion is simulated by using the fast Fourier transform. In figures 1 5 we have generated one sample path of the fractional Brownian motion B with different values of, and the corresponding sample paths of B (), B(3), and B(4), by using (3.3). Overall, the smoothening effect is clearly visible from the pictures for all values {.1,.5,.5,.75,.9}. For all different values of, the paths of the processes B (3) and B (4) look rather smooth. This is clear also from the theoretical point of view, as B (3) is already twice continuously differentiable. The differences of the value can be seen obviously from the paths of B itself, but also from the paths of B (). Indeed, the

15 nt ORDER FBM 15 paths of B () corresponding to values.1 and.5 in figures 1 and are not clearly smooth. From theoretical point of view, the paths are barely continuously differentiable, as the derivative is very rough. In contrast, the paths of B () corresponding to values.75 and.9 are almost twice continuously differentiable, as the path of the B itself is almost differentiable in the sense that the ölder index is close to 1 corresponding to the smooth case. This can be seen also from figures 4 and 5. Indeed, comparing B (3).1 in Figure 1(C) to B().9 in Figure 5(B) the smoothness looks very similar. This is not surprising from theoretical point of view, as the regularity index of B (3).1 is.1 (meaning B (3).1 is twice continuously differentiable and the second derivative is.1-ölder) while the regularity index of B ().9 is 1.9. We also stress that in all of the Figures 1 5 it seems that the values seem to get smaller as n increases. This phenomena is again supported by the theory. Indeed, first of all we are integrating over [, t for t 1 which means that the extrema points gets smaller in absolute value. In addition, the small deviation probabilities gets higher (see [, 3, 4) as n increases. That is, for larger value of n the process is smoother which implies that the process stays in a small ɛ-ball around starting point B () with higher probability. 6. Concluding remarks Perrin et al. [13 defined the nth order fractional Brownian motion by using the Mandelbrot Van Ness representation (.), which lead to the formula (3.1). In this article we have defined the nth order fractional Brownian motion B of urst index (n 1, n) by using the Molchan Golosov representation of the fractional Brownian motion. In addition, we have provided the transfer principle for B that can be used to develop stochastic calculus with respect to B in a relatively simple manner. In addition, our compact interval representation is very useful, e.g. for simulations. Finally, for filtering problems it is of utmost important that the filtrations of B and the corresponding Brownian motion W coincide. We have shown here that by using the Molchan Golosov representation, this property is inherited directly from the same property of the underlying fractional Brownian motion. We have also used this fact to study prediction law of B and the equivalence of law problem that is important for maximum likelihood estimation. References [1 E. Alòs, O. Mazet, and D. Nualart, Stochastic calculus with respect to Gaussian processes, Ann. Probab., 9 (1), pp [ F. Aurzada, Path regularity of Gaussian processes via small deviations, Probab. Math. Stat., 31 (11), pp [3 F. Aurzada, F. Gao, T. Kühn, W. Li, and Q.-M. Shao, Small deviations for a family of smooth Gaussian processes, J. Theor. Prob., 6 (13), pp [4 F. Aurzada, I. Ibragimov, M. Lifshits, and J. van Zanten, Small deviations of smooth stationary Gaussian processes, Theory Probab. Appl., 53 (9), pp

16 3. order fbm 4. order fbm fbm. order fbm 16 SOTTINEN AND VIITASAARI (a) 1. order (b). order (c) 3. order (d) 4. order. Figure order fractional Brownian motion with.1. [5 P. Brockwell and R. Davis, Time Series: Theory and Methods, vol., Springer, [6 M. itsuda, Representation of Gaussian processes equivalent to Wiener process, Osaka J. Math., 5 (1968), pp [7 S. Janson, Gaussian ilbert spaces, vol. 19 of Cambridge Tracts in Mathematics, Cambridge University Press, Cambridge, [8 B. B. Mandelbrot and J. W. Van Ness, Fractional Brownian motions, fractional noises and applications, SIAM Rev., 1 (1968), pp [9 Y. S. Mishura, Stochastic calculus for fractional Brownian motion and related processes, vol. 199 of Lecture Notes in Mathematics, Springer-Verlag, Berlin, 8. [1 G. M. Molchan, istorical comments related to fractional Brownian motion, in Theory and applications of long-range dependence, Birkhäuser Boston, Boston, MA, 3, pp [11 G. M. Molchan and J. I. Golosov, Gaussian stationary processes with asymptotically a power spectrum, Dokl. Akad. Nauk SSSR, 184 (1969), pp [1 I. Norros, E. Valkeila, and J. Virtamo, An elementary approach to a Girsanov formula and other analytical results on fractional Brownian motions, Bernoulli, 5 (1999), pp

17 3. order fbm 4. order fbm fbm. order fbm nt ORDER FBM (a) 1. order (b). order (c) 3. order (d) 4. order. Figure order fractional Brownian motion with.5. [13 E. Perrin, R. arba, C. Berzin-Joseph, I. Iribarren, and A. Bonami, nth-order fractional Brownian motion and fractional Gaussian noises, IEEE Transactions on Signal Processing, 49 (1), pp [14 V. Pipiras and M. S. Taqqu, Are classes of deterministic integrands for fractional Brownian motion on an interval complete?, Bernoulli, 7 (1), pp [15 F. Smithies, Integral Equations, Cambridge Tracts in Mathematics, Cambridge University Press, [16 T. Sottinen, On Gaussian processes equivalent in law to fractional Brownian motion, J. Theoret. Probab., 17 (4), pp [17 T. Sottinen and C. A. Tudor, On the equivalence of multiparameter Gaussian processes, J. Theoret. Probab., 19 (6), pp [18 T. Sottinen and L. Viitasaari, Stochastic analysis of Gaussian processes via Fredholm representation, International Journal of Stochastic Analysis, (16), p. DOI:1.1155/16/ [19 T. Sottinen and L. Viitasaari, Prediction law of fractional Brownian motion, Stat. Probab. Lett., 19 (17), pp

18 3. order Bm 4. order Bm Bm. order Bm 18 SOTTINEN AND VIITASAARI (a) 1. order (b). order (c) 3. order (d) 4. order. Figure order standard Brownian motion.

19 3. order fbm 4. order fbm fbm. order fbm nt ORDER FBM (a) 1. order (b). order (c) 3. order (d) 4. order. Figure order fractional Brownian motion with.75.

20 3. order fbm 4. order fbm fbm. order fbm SOTTINEN AND VIITASAARI (a) 1. order (b). order (c) 3. order (d) 4. order. Figure order fractional Brownian motion with.9.