A set-indexed fractional Brownian motion

Erick Herbin Dassault Aviation 78 quai Marcel Dassault 92552 Saint-Cloud Cedex France erick.herbin@dassault-aviation.fr Ely MerzbachCorresponding author. Tel.: +972-2-9941168 Dept. of Mathematics Bar Ilan University 52900 Ramat-Gan Israel merzbach@macs.biu.ac.il

June 2004

Abstract
We define and prove the existence of a fractional Brownian motion indexed by a collection of closed subsets of a measure space. This process is a generalization of the set-indexed Brownian motion, when the condition of independance is relaxed. Relations with the Lévy fractional Brownian motion and with the fractional Brownian sheet are studied. We prove stationarity of the increments and a property of self-similarity with respect to the action of solid motions. Regularity conditions are exhibited. Finally, behavior of the set-indexed fractional Brownian motion along increasing paths is analysed.
AMS classification: 62 G 05, 60 G 15, 60 G 17, 60 G 18.
Keywords: fractional Brownian motion, Gaussian processes, stationarity, self-similarity, set-indexed processes.

1 Introduction

Recently, different extensions of the fractional Brownian motion were introduced and widely used to describe complex or chaotic phenomena in several fields of sciences. Here we define a new very natural set-indexed generalization of the fractional Brownian motion. The set-indexed fractional Brownian motion studied here seems to be well-adapted to modelize problems in applied mathematics (see [11).
The fractional Brownian motion was defined by B. B. Mandelbrot and J. W. Van Ness, and extended essentially into two directions. One is generally called the multifractional Brownian motion, replacing the index parameter of self-similarity (called also the Hurst parameter) by a real measurable function (see [2, [10). The other one are multiparameter fractional Brownian motions in which the set of the indices is a subset of the Euclidean space R N   (see [8and [12, [14for trajectory properties).
At least two types of multiparameter fractional Brownian motions were introduced. One is called the Lévy fractional Brownian motion because it extends the Lévy Brownian motion and the other is called the fractional Brownian sheet because it can be seen as an extension of the Brownian sheet. We refer to [4for a survey on all these processes. Moreover, in [4, some multiparameter extensions of the multifractional Brownian motion are well studied.
The set-indexed fractional Brownian motion introduced here is a simple extension of the set-indexed Brownian motion (called also the white noise) and possesses the main properties required for a fractional Brownian motion.
Moreover, by choosing the class of sets parametrizing the process, we get great flexibility and possibilities to obtain particular types of fractional Brownian motion.
In this paper, we prove that our definition of set-indexed fractional Brownian motion satisfies both self-similarity and a condition of stationarity. Moreover such a process is a time-changed classical fractional Brownian motion on each increasing path (flow).
Conversely, we compute the covariance function of any self-similar and stationary set-indexed process. For Gaussian processes, we get an extension of our set-indexed fractional Brownian motion.
This paper is divided as follow :
In the next section, we present the general framework needed for set-indexed processes. We define the concept of set-indexed fractional Brownian motion (SifBm). We prove the existence of this process showing that its covariance function is positive definite. Moreover, we compare our definition to previous definitions given in some particular cases and our definition seems to be quite natural and satisfactory. The two fractal properties which are stationarity and self-similarity are studied in section 3. As it will be see, stationarity of increments can be defined in different non equivalent ways. In section 4, we discuss the possibility of finding a characterization of the set-indexed fractional Brownian motion by the two main properties: stationarity and self-similarity.
Section 5 deals with the problem of continuity. In fact, we show that the SifBm is continuous when the set-indexed Brownian motion is also continuous. Finally in the last section, we study SifBm on increasing paths.

2 Framework and definition

2.1 Indexing collection, set-indexed processes

Let T   be a locally compact complete separable metric and measure space with metric d   and measure m   . All processes will be indexed by a class A   of compact connected subsets of T   .
In what follows, for any class of sets D   , the class of finite unions of sets from D   will be denoted by D ( u )   . In the terminology of [7, we assume that A   is an indexing collection:
Definition 2.1 A nonempty class A   of compact, connected subsets of T   is called an indexing collection if it satisfies the following:
  • 1. , T A   , and A A   if A   or T   .
  • 2. A   is closed under arbitrary intersections and if A , B A   are nonempty, then A B   is nonempty. If ( A i )   is an increasing sequence in A   then i A i ¯ A   .
  • 3. The σ   -algebra generated by A   , σ ( A ) =   , the collection of all Borel sets of T   .
  • 4. Separability from above There exists an increasing sequence of finite subclasses A n = { A 1 n , . . . , A k n n }   of A   closed under intersections and satisfying , T A n   and a sequence of functions g n : A A n ( u )   such that
    • (a) g n   preserves arbitrary intersections and finite unions (i.e.
      g n ( A A A ) = A A g n ( A )   for any A A   , and if i = 1 k A i = j = 1 m A j   , then i = 1 k g n ( A i ) = j = 1 m g n ( A j )   ),
    • (b) for each A A   , A ( g n ( A ) )   ,
    • (c) g n ( A ) g m ( A )   if n m   ,
    • (d) for each A A   , A = n g n ( A )   ,
    • (e) if A , A A   then for every n   , g n ( A ) A A   , and if A A n   then g n ( A ) A A n   .
    • (f ) g n ( ) = n   .
  • 5. Every countable intersection of sets in A ( u )   may be expressed as the closure of a countable union of sets in A   .
(Note: `   ' indicates strict inclusion and ` ( ) ¯   ' and` ( )   ' denote respectively the closure and the interior of a set.)
We shall define the semi-algebra C   to be the class of all subsets of T   of the form C = A \ B , A A , B A ( u ) .   C   is closed under intersections and any set in C ( u )   may be expressed as a finite disjoint union of sets in C   . Note that if B = i = 1 k A i A ( u )   , without loss of generality we can require that for each i   , A i j i A j   . Such a representation of B A ( u )   will be called extremal. If C = A \ B , A A , B A ( u )   , then the representation of C   is called extremal if that of B   is. Unless otherwise stated, it will always be assumed that all representations of sets in A ( u )   and C   are extremal.
Numerous examples of topological spaces T   and indexing collections A   satisfying the preceding assumptions may be found in [7. In particular, our framework generalizes the usual multiparameter setting: if T = R + N   , then the class A = { [ 0 , t ] : t R + N }   satisfies all the assumptions. More generally, we can allow A   to consist of all the lower layers of R + N   : a set A   is a lower layer if [ 0 , t ] A , t A   .
Now, let ( Ω , , P )   be any complete probability space. A filtration (indexed by A   ) is a class { A : A A }   of complete sub- σ   -fields of   which satisfies the following conditions:
  • A , B A   , A B   , if A B   .
  • Monotone outer-continuity: A i = A i   for any decreasing sequence ( A i )   in A   .
Definition 2.2 A ( A   -indexed) stochastic process X = { X A : A A }   is a collection of random variables indexed by A   with X = 0   , and is said to be adapted if X A   is A   -measurable, for every A A   . X   is said to be integrable (square integrable) if E [ | X A | ] <   ( E [ ( X A ) 2 ] <   ) for every A A   .
Remark 2.3 Any multiparameter process X ~   can be considered as a set-indexed process X   , setting T = R + N   , A = { [ 0 , t ] ; t R + N }   and X ~ t = X [ 0 , t ]   .
Definition 2.4 A ( A   -indexed) stochastic process X   is additive if it has an (almost sure) additive extension to C ( u )   : i.e., X = 0   and if C , C 1 , C 2 C ( u )   with C = C 1 C 2   and C 1 C 2 =   , then almost surely X C 1 + X C 2 = X C .  
In addition to assuming that X = 0   , to avoid technicalities we will generally assume as well that X = 0   , where : = A A A   .
Definition 2.5 An additive process X   is increasing if X C 0 C C   and if for any decreasing sequence ( A n )   in A   , X ( n A n ) = lim n X ( A n )   .
It is observed in Corollary 1.4.11 of [7that an increasing process in fact defines a measure on   for each ω Ω   .
Definition 2.6 Let Λ   be a non-negative increasing function defined on A   with Λ = 0   . We say that an A   -indexed additive process X   is a Brownian motion with variance measure Λ   if X = 0   , and if for disjoint sets C 1 , . . . , C n C , X C 1 , . . . , X C n   are independent mean-zero Gaussian random variables with variances Λ C 1 , . . . Λ C n   , respectively.

2.2 Increments of a set-indexed process

The notion of increments for a set-indexed process X = { X U ; U A }   is not as simple as in the case of real indices, where it is only the difference between values of the process.
In the case of multiparameter processes, we use to define the increment between s , t R + N   by
Δ X s , t = r { 0 , 1 } N ( 1 ) N i r i X [ s i + r i ( t i s i ) ] i (1)
which is different from the simple difference X t X s   (see [4).
In the case of set-indexed processes, the increments are defined from the collection of subsets C   .
For all C = U \ 1 i n U i   , we define the increment of the process X   on C   by
Δ X C = X U i = 1 n X U U i + i < j X U ( U i U j ) + ( 1 ) n X U ( 1 i n U i ) (2)
According to remark  2.3 , this expression, applied to the multiparameter case, gives the definition ( 1 ) of the increments.
In the following, it would be important to consider the particular increments C 0 = { C = U \ V ; U , V A }   . Moreover the definition of the increment process Δ X   can be extended to C ( u )   , the finite unions of elements of C   . Particularly, for all U , V A   , Δ X U V   is well-defined.
Remark 2.7 The process Δ X   could be seen as an extension of the process X   for the set of indices C   . For all A A C   , we have Δ X A = X A   (because A = A \   ).
In the case of an additive process, the definition of the increment Δ X C   coincides with the additive extension of X C   to C C   . However, if X   is not additive, which is the case of a direct definition of the process for the set of indices C   , in general Δ X C X C   for C C   . For this reason, we use a different notation for the increments of X   .
Remark 2.8 If X = { X U ; U A }   is Gaussian, then Δ X = { Δ X C ; C C }   is clearly Gaussian.

2.3 Definition of sifBm

Recall that the fractional Brownian motion is defined to be a mean-zero Gaussian process such that s , t R + ; E [ X t X s ] 2 = | t s | 2 H   The natural set-indexed extension of this process is to substitute the term | t s | 2 H   with d ( U , V ) 2 H   , where d   is some distance between two subsets of T   .
In this paper, we consider the choice of d ( U , V ) = m ( U V )   , where   is the symmetric difference between two sets and m   is a measure on T   .
Lemma 2.9 Let m   be a finite measure on T   . For all α ( 0 , 2 ]   , the function ( U , V ) m ( U ) α + m ( V ) α m ( U V ) α   is positive definite.

Proof

For all measurable subset U   of T   such that m ( U ) < +   , we define the elementary function f = 1 U   . Finite linear combinations of elementary functions are called simple functions. It is well-known that simple functions are dense in L 2 ( m )   .
Moreover, for all U , V T   , we have
1 U V = 1 ( U \ V ) ( V \ U )
= 1 U ( 1 1 V ) + 1 V ( 1 1 U )
= 1 U + 1 V 21 U . 1 V
= ( 1 U 1 V ) 2
= | 1 U 1 V |
Then we only have to show that the function
L 2 ( m ) × L 2 ( m ) R
( f , g ) m ( f 2 ) α + m ( g 2 ) α m ( | f g | 2 ) α
is positive definite.
Let f 1 , f 2 , . . . , f n L 2 ( m )   and u 1 , u 2 , . . . , u n R   . We have to show that
i = 1 n j = 1 n { m ( f i 2 ) α + m ( f j 2 ) α m ( | f i f j | 2 ) α } u i u j 0 (3)
Setting u 0 = i = 1 n u i   and f 0 = 1   , we can write
i = 1 n j = 1 n { m ( f i 2 ) α + m ( f j 2 ) α m ( | f i f j | 2 ) α } u i u j = i = 0 n j = 0 n m ( | f i f j | 2 ) α u i u j (4)
But, for all λ > 0   , we have
i = 0 n j = 0 n e λ m ( | f i f j | 2 ) α u i u j = i = 0 n j = 0 n ( e λ m ( | f i f j | 2 ) α 1 ) u i u j
= λ i = 0 n j = 0 n m ( | f i f j | 2 ) α u i u j + o ( λ ) (5)
Then, ( 3 ) is equivalent to
i = 0 n j = 0 n e λ m ( | f i f j | 2 ) α u i u j 0 (6)
in the neighborhood of λ = 0   .
In L 2 ( m )   , let us define the bilinear form < f , g > = m ( f . g )   . If we identify the elements of L 2 ( m )   that are almost everywhere equal, then L 2 ( m )   is a complete separable metric space with this scalar product.
Let us show that there exists a random variable X   taking its values in L 2 ( m )   , such that
f L 2 ( m ) ; E [ e i < f , X > ] = e λ f 2 α (7)
Consider a α   -stable real random variable Y = ( d ) S α ( ( cos π α 2 ) 1 α , 1 , 0 )   where 0 < α 2   and α 1   (see [13, prop. 1.2.12, p 15).
As the function f e 1 2 f 2   is positive definite (see remark  2.10 ), by a theorem of Bochner-Minlos, there exists a random variable G   , such that f L 2 ( m ) ; E [ e i < f , G > ] = e 1 2 f 2 .   Moreover, we can suppose that G   is independant from Y   .
For all f L 2 ( m )   , we compute
E [ e i < f , Y 1 / 2 G > ] = E { E [ e i < f , Y 1 / 2 G > | Y ] }
= E { E [ e i Y 1 / 2 < f , G > | Y ] }
= E [ e Y 2 f 2 ]
As E [ e γ Y ] = e γ α   for all γ > 0   , we get E [ e i < f , Y 1 / 2 G > ] = e 2 α f 2 α   Then X = { 2 λ 1 / α Y 1 / 2 G if α 1 2 λ G if α = 1   satisfies ( 7 ), and f e λ f 2 α   is non-negative definite. That proves ( 6 ) and the result follows.  
Remark 2.10 For all family ( f 1 , . . . , f n )   of L 2 ( m )   provided with the scalar product previously defined, there exist p n   and a family ( x 1 , . . . , x n )   of R p   such that
i , j ; x i x j R p = f i f j L 2 ( m ) (8)
To show this result, let us consider an orthonormal basis ( e 1 , . . . , e p )   of V e c t ( f i ) 1 i n   , and the canonical basis ( ε 1 , . . . , ε p )   of R p   . For all i   , there exists a family ( λ 1 i , . . . , λ p i )   s.t. f i = k λ k i . e k   . Then the vectors ( x 1 , . . . , x n )   defined by x i = k λ k i . ε k   satisfy ( 8 ).   This remark allows to show directly ( 3 ) using lemma 2.10.8 in [13Since the existence of a mean-zero Gaussian process is equivalent to the positive definite property of its covariance function, we can define
Definition 2.11 A mean-zero Gaussian process B H = { B U H ; U A }   satisfying
E [ B U H B V H ] = 1 2 [ m ( U ) 2 H + m ( V ) 2 H m ( U V ) 2 H ] (9)
where H ( 0 , 1 )   , is called a set-indexed fractional Brownian motion (sifBm).
H   is the index of self-similarity of the process.
Remark 2.12 If H = 1 2   , the process B 1 2   is the well known set-indexed Brownian motion. Indeed, let us compute the covariance function of this process E [ B U 1 2 B V 1 2 ] = 1 2 [ m ( U ) + m ( V ) m ( U V ) ]   As
m ( U V ) = m ( U \ V ) + m ( V \ U )
= m ( U \ U V ) + m ( V \ U V )
= m ( U ) + m ( V ) 2 m ( U V )
we have E [ B U 1 2 B V 1 2 ] = m ( U V )   which is the covariance function of the set-indexed Brownian motion.
Remark 2.13 In the case of T = R +   and A = { [ 0 , t ] ; t R + }   , the process B H   is the classical fBm. Indeed, the covariance function is E [ B [ 0 , s ] H B [ 0 , t ] H ] = 1 2 [ s 2 H + t 2 H | t s | 2 H ]   which is the covariance function of the fractional Brownian motion.
Remark 2.14 In the case of T = R + N   and A   is the set of rectangles of the form [ 0 , t ]   , the process B H   can be seen as a multiparameter process. Then it is interesting to compare it with the other known multiparameter fractional Brownian motions (see [4).
  • the Lévy fractional Brownian motion is a mean-zero Gaussian process X = { X t ; t R + N }   such that s , t R + N ; E [ X s X t ] = 1 2 [ s 2 H + t 2 H t s 2 H ]   where H ( 0 , 1 )   .
  • the fractional Brownian sheet is a mean-zero Gaussian process X = { X t ; t R + N }   such that s , t R + N ; E [ X s X t ] = 1 2 N i = 1 [ s i 2 H i + t i 2 H i | t i s i | 2 H i ]   where H = ( H 1 , . . . , H N ) ( 0 , 1 ) N   .
As for all s , t R + N   , E [ B [ 0 , s ] H B [ 0 , t ] H ] = 1 2 [ N i = 1 s i 2 H + N i = 1 t i 2 H ( i s i + i t i 2 i s i t i ) 2 H ]   we can see that B H   is different from these two processes.
This fact will be also shown in the next sections in the study of properties of sifbm, and its restriction on flows. It is therefore natural to wonder if the Levy fBm and the fBs can have set-indexed extension. The answer seems to be negative.
Actually the definition of the sheet is strongly associated with the Euclidean structure of R N   . Therefore it is incompatible with a set-indexed viewpoint.
Moreover the Levy fBm can be seen as a simple one parameter process where the increment between two points only depends from distance between them.

3 Fractal properties

The fractional Brownian motion has two important properties which make it the most natural fractal process:
  • its increments are stationary h R + ; ( X t + h X h ) t R + = ( d ) ( X t X 0 ) t R +  
  • it is self-similar a R + ; ( X a t ) t R + = ( d ) ( a H X t ) t R +  
Moreover, the fBm is the only Gaussian process which has these two properties.
In this section, we show that in some sense these properties still hold for the set-indexed fractional Brownian motion. Moreover they characterize the covariance function of the process between two sets U   and V   such that U V   .

3.1 Stationarity of the increments

Stationarity of increments of a set-indexed process can be defined in various ways. The set-indexed Brownian motion satisfies all of them, but the different extensions of fractional Brownian motion do not.
In the case of T = R + N   and A   is the collection of rectangles, the classical definition of stationarity of increments can be studied.
A process X = { X [ 0 , t ] ; t R + N }   is said to have stationary increments against translations if for all τ R + N   , the two processes { Δ X [ τ , t + τ ] ; t R + N }   and { Δ X [ 0 , t ] ; t R + N }   have the same law.
Both Lévy fractional Brownian motion and fractional Brownian sheet satisfy this property of stationarity (see [4).
Remark 3.1 This definition is weaker than stationarity of increments against isometries of R + N   , i.e. for all g G ( R N )   , { Δ X g ( [ 0 , t ] ) ; t R + N } = ( d ) { Δ X [ 0 , t ] ; t R + N }   where G ( R N )   is the group of isometries of R N   preserving C   .
On the contrary to stationarity against translations, the context of set-indexed processes imposes additional assumptions to the strict context of multiparameter processes. Actually as the image of C C   by any isometry of R N   does not necessarily belong to C   , a stability assumption is needed for the definition to make sense.
However in the strict context of multiparameter processes, this assumption does not need to be considered. The Lévy fractional Brownian motion satisfy this property of increment stationarity in the strong sense (see [13).
However in general, there is no reason that the SifBm possesses the stationarity increments property against translations. In some particular cases, we can show directly using the next lemma that this property is not satisfied.
Lemma 3.2 Let B H = { B U H ; U A }   be a SifBm of index H ( 0 , 1 )   . For all C = U \ ( 1 i n U i ) C   , where i { 1 , . . . , n } ; U i U   , we have
E [ Δ B C H ] 2 = k = 1 n ( 1 ) k i 1 < < i k m ( U [ p { i 1 , . . . , i k } U p ] ) 2 H + 1 2 k , l ( 1 ) k + l i 1 < < i k j 1 < < j l m ( [ p { i 1 , . . . , i k } U p ] [ p { j 1 , . . . , j l } U p ] ) 2 H (10)

Proof

By definition, Δ B C H = B U H + k = 1 n ( 1 ) k i 1 < < i k B p { i 1 , . . . , i k } U p H   Then,
E [ Δ B C H ] 2 = E [ B U H ] 2 + 2 k ( 1 ) k i 1 < < i k E [ B U H . B p { i 1 , . . . , i k } U p H ] + k , l ( 1 ) k + l i 1 < < i k j 1 < < j l E [ B p { i 1 , . . . , i k } U p H . B p { j 1 , . . . , j l } U p H ]  
and, using the covariance function of X  
E [ Δ B C H ] 2 = m ( U ) 2 H + k ( 1 ) k i 1 < < i k { m ( U ) 2 H + m ( p { i 1 , . . . , i k } U p ) 2 H m ( U [ p { i 1 , . . . , i k } U p ] ) 2 H } + 1 2 k , l ( 1 ) k + l i 1 < < i k j 1 < < j l { m ( p { i 1 , . . . , i k } U p ) 2 H + m ( p { j 1 , . . . , j l } U p ) 2 H m ( [ p { i 1 , . . . , i k } U p ] [ p { j 1 , . . . , j l } U p ] ) 2 H } (11)
Let us consider the two following terms in expression ( 11 ):
  • term in m ( U ) 2 H   m ( U ) 2 H + k ( 1 ) k i 1 < < i k m ( U ) 2 H = k = 0 n C n k m ( U ) 2 H = 0  
  • term in m ( p { i 1 , . . . , i k } U p ) 2 H  
    k ( 1 ) k i 1 < < i k m ( p { i 1 , . . . , i k } U p ) 2 H + k , l ( 1 ) k + l i 1 < < i k j 1 < < j l m ( p { i 1 , . . . , i k } U p ) 2 H k ( 1 ) k i 1 < < i k m ( p { i 1 , . . . , i k } U p ) 2 H l = 1 n ( 1 ) l C n l
    therefore this term is equal to k ( 1 ) k i 1 < < i k m ( p { i 1 , . . . , i k } U p ) 2 H l = 0 n ( 1 ) l C n l = 0  
The two other terms of expression ( 11 ) give the result.   The main idea to define a set-indexed extension of the fractional Brownian motion, was to extend s , t R + ; E [ X t X s ] 2 = | t s | 2 H   in U , V A ; E [ X U X V ] 2 = m ( U V ) 2 H   However, it should be more interesting to get
C C ; E [ Δ X C ] 2 = m ( C ) 2 H (12)
According to lemma  3.2 , the set-indexed fractional Brownian motion B H   satisfies U , V A ; E [ Δ B U \ V H ] 2 = m ( U \ V ) 2 H   but the property ( 12 ) does not hold.
Moreover, we will see that there is no set-indexed process satisfying ( 12 ) for H 1 2   (theorem  4.3 ).
Proposition 3.3 Let X = { X t ; t R + N }   be a fractional Brownian sheet of constant parameter H   in every axis.
For all a , b R + N   such that a b   (i.e. i = 1 , . . . , N ; a i < b i   ), we have E [ Δ X [ a , b ] ] 2 = m ( [ a , b ] ) 2 H   and consequently, For all a , b , a , b R + N   , such that m ( [ a , b ] ) = m ( [ a , b ] )   , we have Δ X [ a , b ] = ( d ) Δ X [ a , b ]  

Proof

For all a , b R + N   with a b   , as X   has stationary increments against translations, we have
E [ Δ X [ a , b ] ] 2 = E [ Δ X [ 0 , b a ] ] 2
= E [ X b a ] 2
= N i = 1 | b i a i | 2 H
= m ( [ a , b ] ) 2 H
As Δ X   is a mean-zero Gaussian process, the result follows.  
Remark 3.4
  • If a fractional Brownian sheet of index H = ( H 1 , . . . , H N )   satisfies the property of proposition  3.3 , we have H 1 = = H N   .
  • The Levy fractional Brownian motion does not satisfy this property.
Definition 3.5 A set-indexed process X   is said to be C 0   -increment stationary if for all C , C C 0   such that m ( C ) = m ( C )   , we have Δ X C = ( d ) Δ X C   .
In the case of one-parameter fractional Brownian motion, the increment stationarity property gives an equality between laws of some increment processes.
However, in the set-indexed case, definition  3.5 does not tell anything about correlation between increments.
In section  4 , we see that a stronger property is not worth to be considered.
Proposition 3.6 The set-indexed fractional Brownian motion B H   is C 0   -increment stationary.

Proof

For all C C 0   , there exist U , V A   where V U   , such that C = U \ V   .
Then Δ B C H = B U H B U V H = B U H B V H   . We compute
E [ Δ B C H ] 2 = E [ B U H B V H ] 2
= m ( U V ) 2 H
= m ( C ) 2 H
Thus, as Δ B H   is a Gaussian process, for all C , C C 0   such that m ( C ) = m ( C )   , we have Δ B C H = ( d ) Δ B C H   .  
Remark 3.7 In the proof of proposition  3.6 , we saw that C C 0 ; E [ Δ B C H ] 2 = m ( C ) 2 H   However, in general E [ Δ B C H . Δ B C H ] 1 2 [ m ( C ) 2 H + m ( C ) 2 H m ( C C ) 2 H ]   In fact, it can be shown that for C = U \ V   where V U   , and C = U \ V   where V U   , and U , V , U , V A   E [ Δ B C H . Δ B C H ] = 1 2 [ m ( U V ) 2 H + m ( V U ) 2 H m ( U U ) 2 H m ( V V ) 2 H ]  

3.2 Self-similarity

To study a set-indexed version of the notion of self-similarity for a set-indexed process, we need some assumptions about the set A   .
We suppose that A   is provided with the operation of a group G   that can be extended satifying
U , V A , g G ; g . ( U V ) = g . U g . V (13)
g . ( U \ V ) = g . U \ g . V
and assume there exists a function μ : G R +  
U A , g G ; m ( g . U ) = μ ( g ) . m ( U ) (14)
Remark 3.8 We can see easily that μ   is an group-homomorphism.
Example 3.9 In the case of T = R N   and A = { [ 0 , t ] ; t R N }   , we can consider the multiplication by the elements of R +   g R + , t R N ; g . [ 0 , t ] = [ 0 , g . t ]   Moreover, g R + , t R N ; m ( g . [ 0 , t ] ) = g N m ( [ 0 , t ] )  
The following result will be useful in the next section.
Lemma 3.10 Under the assumptions about the group G   , the cardinal of G   is not finite.

Proof

If the function μ   is not constant, there exists g G   such that μ ( g ) > 1   (take g ~   s.t. μ ( g ) 1   and then g = g ~   or g = g ~ 1   ).
For all integer n   , we have μ ( g n ) = [ μ ( g ) ] n   . If G   is finite, the set { g n ; n N }   is finite, which is in conflict with lim n μ ( g n ) =   .  
Definition 3.11 A set-indexed process X = { X U ; U A }   is said to be self-similar of index H   , if for all g G   ,
{ X g . U ; U A } = ( d ) { μ ( g ) H . X U ; U A } (15)
Proposition 3.12 The set-indexed fractional Brownian motion B H   is self-similar of index H   , for the operation of the group G   .

Proof

Let g   be an element of the group G   . For all U , V A   , we have
E [ B g . U H B g . V H ] = 1 2 [ m ( g . U ) 2 H + m ( g . V ) 2 H m ( g . U g . V ) 2 H ]
As g . ( U V ) = g . U g . V   , we get
E [ B g . U H B g . V H ] = μ ( g ) 2 H 2 [ m ( U ) 2 H + m ( V ) 2 H m ( U V ) 2 H ]
= μ ( g ) 2 H E [ B U H B V H ]
Therefore, the two mean-zero Gaussian processes { B g . U H ; U A }   and { μ ( g ) H . B U H ; U A }   have the same law.  

4 Pseudo-characterisation of SifBm

Recall that the fractional Brownian motion is the only mean-zero Gaussian process which is self-similar and has stationary increments. In the same way, the only multiparameter mean-zero Gaussian process which is self-similar and whose increments are stationary in the strong sense (under isometries of R N   ), is the Levy fractional Brownian motion ([13).
In the case of the set-indexed processes, there is not such a characterisation.
However, the two properties of self-similarity and stationarity of increments characterise the covariance function of the process between all U   and V   such that U V   .
Proposition 4.1 Let X = { X U ; U A }   be a set-indexed process satisfying the two following properties
  • 1. self-similarity of index H  
  • 2. C 0   -increment stationarity
Then, the covariance function between two subsets U   and V   such that U V   is
E [ X U X V ] = K [ m ( U ) 2 H + m ( V ) 2 H m ( V \ U ) 2 H ] (16)

Proof

Let U 0   be a non m   -null fixed element of A   . For all U , V A   such that U V   , there exists g G   such that m ( g . U 0 ) = m ( V \ U )   , i.e. μ ( g ) . m ( U 0 ) = m ( V \ U )   .
Then using C 0   -increment stationarity property, we have
E [ X V X U ] 2 = E [ Δ X V \ U ] 2
= E [ Δ X g . U 0 ] 2
As g . U 0 A   , we have Δ X g . U 0 = X g . U 0   and by self-similarity
E [ X V X U ] 2 = E [ X g . U 0 ] 2
= [ μ ( g ) ] 2 H E [ X U 0 ] 2
= [ m ( V \ U ) ] 2 H E [ X U 0 ] 2 m ( U 0 ) 2 H (17)
The result follows from ( 17 ).  
Remark 4.2 The proposition  4.1 shows that our set-indexed extension of the fractional Brownian motion is very natural provided the two properties of self-similarity and stationarity of the increments.
However, if there exists a mean-zero Gaussian process with covariance function
E [ X U X V ] = 1 2 [ m ( U ) 2 H + m ( V ) 2 H m ( V \ U ) 2 H m ( U \ V ) 2 H ] (18)
it satisfies proposition  4.1 as well.
To determine completely the covariance function of a self-similar, C 0   -increments, set-indexed process, we need assumptions about E [ Δ X U \ V . Δ X V \ U ]   , where U , V A   .
For all U , V A   , we have Δ X U \ V = X U X U V   and Δ X V \ U = X V X U V   .
Then,
E [ Δ X U \ V . Δ X V \ U ] = E [ X U . X V ] E [ X U . X U V ] E [ X V . X U V ] + E [ X U V ] 2
= E [ X U . X V ] 1 2 [ m ( U ) 2 H + m ( V ) 2 H m ( V \ U ) 2 H m ( U \ V ) 2 H ]
Particularly, assuming the independance of Δ X U \ V   and Δ X V \ U   , is equivalent to ( 18 ), provided that such a process X   exists.
The property of C 0   -increment stationarity seems too weak to characterize completely the covariance of a self-similar process. It can be tempting to define a process which would have a stronger kind of increment stationarity. For instance, does it exist a self-similar process which satisfy E [ Δ X C ] 2 = m ( C ) 2 H   for all C C   ?
Theorem 4.3 The only Gaussian set-indexed process X = { X U ; U A }   such that
C C ; E [ Δ X C ] 2 = K . m ( C ) 2 H (19)
where K > 0   and H ( 0 , 1 )   , is the set-indexed Brownian motion.

Proof

Let X   be a set-indexed Gaussian process satisfying ( 19 ).
First, we can see that X   is C 0   -increment stationary. Moreover, for all U A   and g G   , we have E [ X g . U ] 2 = K . μ ( g ) 2 H . m ( U ) 2 H = μ ( g ) 2 H . E [ X U ] 2   . Then, as X   is Gaussian, we conclude that X   is self-similar. Proposition  4.1 implies that for all U , V A   , such that U V   ,
E [ X U . X V ] = K [ m ( U ) 2 H + m ( V ) 2 H m ( V \ U ) 2 H ] (20)
For all U 1   and U 2   in A   , let us consider U A   such that U 1 U   and U 2 U   . The subset of T   , C = U \ ( U 1 U 2 )   belongs to C   and Δ X C = X U X U 1 X U 2 + X U 1 U 2   . Then we have
E [ Δ X C ] 2 = E [ X U ] 2 + E [ X U 1 ] 2 + E [ X U 2 ] 2 + E [ X U 1 U 2 ] 2
2 E [ X U . X U 1 ] 2 E [ X U . X U 2 ] + 2 E [ X U . X U 1 U 2 ]
+ 2 E [ X U 1 . X U 2 ] 2 E [ X U 1 . X U 1 U 2 ] 2 E [ X U 2 . X U 1 U 2 ]
Hence
2 E [ X U 1 . X U 2 ] = E [ Δ X C ] 2 E [ X U ] 2 E [ X U 1 ] 2 E [ X U 2 ] 2 E [ X U 1 U 2 ] 2
+ 2 E [ X U . X U 1 ] + 2 E [ X U . X U 2 ] 2 E [ X U . X U 1 U 2 ]
+ 2 E [ X U 1 . X U 1 U 2 ] + 2 E [ X U 2 . X U 1 U 2 ]
Using ( 20 ), we get
2 E [ X U 1 . X U 2 ] = K { m ( U \ ( U 1 U 2 ) ) 2 H m ( U \ U 1 ) 2 H m ( U \ U 2 ) 2 H + m ( U \ ( U 1 U 2 ) ) 2 H }
+ K [ m ( U 1 ) 2 H + m ( U 2 ) 2 H m ( U 1 \ U 2 ) 2 H m ( U 2 \ U 1 ) 2 H ] (21)
Taking V A   such that U ⊊︀ V   , we get an expression of 2 . E [ X U 1 . X U 2 ]   different from ( 21 ) if H 1 2   .  
Corollary 4.4 There exists no set-indexed process which is H   -self-similar (for H 1 2   ) and whose increments satisfy one of the following
  • 1. C   -increment stationarity C , C C ; m ( C ) = m ( C ) Δ X C = ( d ) Δ X C  
  • 2. for all function f : C C   , such that C C ; m ( f ( C ) ) = m ( C )   { Δ X f ( C ) ; C C } = ( d ) { Δ X C ; C C }  

Proof

First, we can see easily that the second property implies C   -increment stationarity.
Then we only need to consider the first property.
Let U 0   be a fixed element of A   . For all C C   , there exists g G   such that m ( g . U 0 ) = m ( C )   , i.e. μ ( g ) . m ( U 0 ) = m ( C )   . Therefore, by C   -increment stationarity,
E [ Δ X C ] 2 = E [ X g . U 0 ] 2
= μ ( g ) 2 H E [ X U 0 ] 2
= m ( C ) 2 H E [ X U 0 ] 2 m ( U 0 ) 2 H
By theorem  4.3 , the result follows.  

5 Continuity of the SifBm

The results about the existence of a continuous version of set-indexed processes are not as simple as processes indexed by R +   . Even in the simple case of the set-indexed Brownian motion, if the collection A   is too rich, there does not exist any version which is continuous on the whole A   (see [1).
To study the continuity of a set-indexed process X   , we have to consider the behavior of | X U X V |   when U   and V   are close. In order to do this, we provide A   with some distance. In the classical case of Gaussian processes, we used to consider the canonical distance d 2 ( U , V ) = E [ X U X V ] 2   for U , V A   (see [1,[3). Let us mention two other distances that are also classical:
  • The measure m   on T   induces the pseudo-metric d m   on A   U , V A ; d m ( U , V ) = m ( U V )  
  • we recall the definition of the Hausdorff metric d H a u s   on K \   , the nonempty compact subsets of T   U , V K \ ; d H a u s ( U , V ) = inf { ε > 0 : U V ε and V U ε }  
The notion of continuity depends of the distance considered. However, if ( A , d m )   (resp. ( A , d H a u s )   ) is compact, then d   -continuity and d m   -continuity (resp. d H a u s   -continuity) are equivalent (see [1).
For any function x : A R   , define x A = sup A A | x ( A ) |   and let B ( A ) = { x : A R ; x A < } .   Let C ( A )   denote the class of functions in B ( A )   which are d H a u s   -continuous on A   .
Then, studying d   -continuity, we consider the balls ( U , ε ) = { V A ; d ( U , V ) < ε }   , and the metric entropy D ( ; A , d )   of ( A , d )   , which gives the smallest number D ( ε ; A , d )   of balls of radius ε > 0   required to cover A   .
If ( A , d )   is totally bounded, i.e. D ( ε ; A , d )   is finite for all ε > 0   , then under the assumption 0 1 ln D ( ε ; A , d ) . d ε <   Dudley's theorem states that the process X   has a continuous modification (see [1, [3, [9).
Theorem 5.1 Let B H = { B U H ; U A }   be a set-indexed fractional Brownian motion.
The two following statements are equivalent
  • 1. B H   is almost surely continuous on A   .
  • 2. the set-indexed Brownian motion W   is almost surely continuous on A   .

Proof

As
U , V A ; E [ B U H B V H ] 2 = ( E [ W U W V ] 2 ) 2 H (22)
the two canonical pseudo-metrics associated to the set-indexed fractional Brownian motion B H   and the set-indexed Brownian motion W   , are equivalent.
Then, Dudley's theorem gives the result.   A simple consequence of this result is that the sifBm has a continuous modification on rectangles of R + N   .

6 SifBm on increasing paths

The notion of flows is the key to reducing the proofs of many of theorems on characterization and weak convergence to a one-dimensional problem (see [7).

6.1 Generality on flows

In general, A ( u )   is not closed under countable intersections, so we will occasionally require a larger class A ~ ( u )   , which is the class of countable intersections of sets in A ( u )   : i.e. U A ~ ( u )   if there exists a sequence ( U n ) n N   in A ( u )   such that n U n = U   .
Definition 6.1 Let S = [ a , b ] R   . An increasing function f : S A ~ ( u )   is called a flow.
  • A flow f   is right-continuous if f ( s ) = v > s f ( v ) , s [ a , b ) ,   and f ( b ) = u < b f ( u ) ¯   .
  • A flow f   is continuous if it is right-continuous and f ( s ) = u < s f ( u ) ¯ s ( a , b ) .  
  • A flow f   is simple if there exists a finite sequence ( t 0 , . . . t n )   with a = t 0 . . . t n = b   and flows f i : [ t i 1 , t i ] A   , i = 1 , . . . , n   such that for s [ t i 1 , t i ]   , f ( s ) = f i ( s ) j = 1 i 1 f j ( t j )   .
Any process X   indexed by A   can be projected by a simple flow f   onto a process indexed by a subset of R   :
Definition 6.2 Let X   be an A   -indexed process and f   a simple flow on S = [ a , b ]   . Then the S   -indexed process X f   is defined as follows:
X f ( s ) : = X f ( s ) , s S .   X f   is called the projection of X   along f   .
In the case that X   can be extended to an additive process on A ~ ( u )   , X   can be projected by any arbitrary flow according to the preceding definition.
The following lemma shows the importance of the concept of flows for set-indexed processes.
Let S ( A )   denote the class of simple continuous flows defined on [ 0 , 1 ]   .
Lemma 6.3 ([6) The finite dimensional distributions of an (additive) A   -indexed process X   determine and are determined by the finite dimensional distributions of the class { X f : f S ( A ) }   .
In this section, we study the set-indexed fractional Brownian motion on flows.

6.2 Levy fractional Brownian motion and fractional Brownian sheet on flows

As a preliminary, let us study the classical cases of Levy fractional Brownian motion and fractional Brownian sheet.
Let us consider T = R + N   and A = { [ 0 , t ] ; t R + N }   . We can associate to any flow f   , an increasing function f ~ : [ 0 , 1 ] R + N   such that t [ 0 , 1 ] ; f ( t ) = [ 0 , f ~ ( t ) ]  
  • If X   is a Levy fBm, s , t [ 0 , 1 ] ; E [ X f ~ ( t ) X f ~ ( s ) ] 2 = f ~ ( t ) f ~ ( s ) 2 H   Then, if f ~ ( t ) = t . α   , where α R + N   , X f ~   is a classical fractional Brownian motion, otherwise it is not.
  • If X   is a fBs, s , t [ 0 , 1 ] ; E [ X f ~ ( s ) . X f ~ ( t ) ] = N i = 1 1 2 [ f ~ i ( s ) 2 H + f ~ i ( t ) 2 H | f ~ i ( t ) f ~ i ( s ) | 2 H ]   Then, if the function f ~   is a line parallel to one axis of R + N   , the process X f ~   is a fractional Brownian motion. However, if f ~ ( t ) = t . α   , where α R + N   , s , t [ 0 , 1 ] ; E [ X f ~ ( s ) . X f ~ ( t ) ] = [ s 2 H + t 2 H | t s | 2 H ] N N i = 1 α i 2 H 2   which is not a fBm.
In the two cases of the classical multiparameter extensions of the fractional Brownian motion, we saw that the projection of the process along a flow, is not in general a real-indexed fBm.

6.3 SifBm on flows is a standard fBm

Our definition for a set-indexed fractional Brownian motion is also justified by the following proposition.
Proposition 6.4 Let X   be a set-indexed fractional Brownian motion, and f   be a flow on [ 0 , 1 ]   . Then the process X f = { X f ( t ) ; t [ 0 , 1 ] }   is a time-changed fractional Brownian motion.

Proof

The process X f   is clearly a mean zero Gaussian process indexed by [ 0 , 1 ]   .
Moreover, its covariance function can be computed
E [ X f ( s ) X f ( t ) ] = 1 2 { m [ f ( s ) ] 2 H + m [ f ( t ) ] 2 H m [ f ( s ) f ( t ) ] 2 H }
For all s t   , we have f ( s ) f ( t )   and then
E [ X f ( s ) X f ( t ) ] = 1 2 { m [ f ( s ) ] 2 H + m [ f ( t ) ] 2 H m [ f ( t ) \ f ( s ) ] 2 H }
= 1 2 { m [ f ( s ) ] 2 H + m [ f ( t ) ] 2 H ( m [ f ( t ) ] m [ f ( s ) ] ) 2 H }
The function θ : [ 0 , 1 ] R +   such that for all t [ 0 , 1 ]   , θ ( t ) = m [ f ( t ) ]   is clearly increasing. Thus it defines a time change and we have
s , t [ 0 , 1 ] ; E [ X f ( s ) X f ( t ) ] = 1 2 { θ ( s ) 2 H + θ ( t ) 2 H | θ ( t ) θ ( s ) | 2 H } (23)
Then { X f θ 1 ( t ) ; t R + }   is a classical fractional Brownian motion.   Proposition  6.4 allows to identify the self-similarity index of the SifBm, as the Hölder exponent of the projection along any flow.
Let us recall the definition of the two classical Hölder exponents of a stochastic process X   at t 0 R +   :
  • the pointwise Hölder exponent α X ( t 0 ) = sup { α : limsup ρ 0 sup s , t B ( t 0 , ρ ) | X t X s | ρ α < }  
  • the local Hölder exponent α ~ X ( t 0 ) = sup { α : limsup ρ 0 sup s , t B ( t 0 , ρ ) | X t X s | | t s | α < }  
Corollary 6.5 Let X   be a set-indexed fractional Brownian motion with self-similarity index H   . The pointwise and local Hölder exponents of the projection X f   along any flows f   at t 0 [ 0 , 1 ]   , satisfy almost surely
α X f ( t 0 ) = { α θ ( t 0 ) . H if α θ ( t 0 ) < 1 H otherwise
α ~ X f ( t 0 ) = { α ~ θ ( t 0 ) . H if α ~ θ ( t 0 ) < 1 H otherwise
where θ   is the real function such that θ ( t ) = m [ f ( t ) ]   ( t [ 0 , 1 ]   ), and α θ ( t 0 )   (resp. α ~ θ ( t 0 )   ) is the pointwise (resp. local) Hölder exponent of θ   at t 0   .

Proof

Let f   be a flow and X f   the projection of X   along f   . By ( 23 ), we have
s , t [ 0 , 1 ] ; E [ X t f X s f ] 2 = | θ ( t ) θ ( s ) | 2 H (24)
  • If θ   is differentiable on [ 0 , 1 ]   , for all t 0 [ 0 , 1 ]   and ρ > 0   ,
    s , t B ( t 0 , ρ ) ; | θ ( t ) θ ( s ) | K . | t s |
    as ρ   tends to 0   . Then,
    s , t B ( t 0 , ρ ) ; E [ X t f X s f ] 2 K . | t s | 2 H (25)
    In [5, we see that equation ( 25 ) implies P { t 0 [ 0 , 1 ] ; α X f ( t 0 ) = α ~ X f ( t 0 ) = H } = 1  
  • If θ   is not differentiable in t 0 [ 0 , 1 ]   (i.e. if α ~ θ ( t 0 ) < 1   ), α < α ~ θ ( t 0 ) ; limsup ρ 0 sup s , t B ( t 0 , ρ ) E [ X t f X s f ] 2 | t s | 2 α . H = 0   and α > α ~ θ ( t 0 ) ; limsup ρ 0 sup s , t B ( t 0 , ρ ) E [ X t f X s f ] 2 | t s | 2 α . H = +   then, α ~ X f ( t 0 ) = α ~ θ ( t 0 ) . H   almost surely (see [5).
    In the same way, we get α X f ( t 0 ) = α θ ( t 0 ) . H   almost surely.
 
Remark 6.6 As a set-indexed process on a flow only depends on its covariance between subsets U   and V   such that U V   , the result stated in proposition  6.4 still holds for the set-indexed mean-zero Gaussian process defined by ( 18 ).
More generally, it holds for all process which is self-similar and C 0   -increment stationary. Therefore this result gives a supplementary legitimity that such a process could be called a set-indexed fractional Brownian motion.

7 Concluding remarks

The preceeding discussions permit the following general definition
Definition 7.1 A set-indexed Gaussian process which is self-similar of index H ( 0 , 1 )   and C 0   -increment stationary, is called general set-indexed fractional Brownian motion.
Corollary 7.2 Let X = { X U ; U A }   be a general set-indexed fractional Brownian motion. Then,
  • 1. For all H ( 0 , 1 )   , such a process exists (Definition  2.11 ),
  • 2. for any flow f   , X f   is a time-changed real-indexed fractional Brownian motion,
  • 3. the covariance function between two subsets U   and V   such that U V   is E [ X U X V ] = K [ m ( U ) 2 H + m ( V ) 2 H m ( V \ U ) 2 H ] .  
Acknowledgement The second author wishes to thank the first author for his kind invitation, and Marco Dozzi for preliminary discussions during a stay at Nancy. References

  1. R. J. Adler, An Introduction to Continuity, Extrema, and Related Topics for General Gaussian Processes, IMS Lect. Notes, Monograph Series, vol 12, Hayward, California, 1990.
  2. A. Benassi, S. Jaffard and D. Roux, Elliptic Gaussian random processes, Rev. Mat. Ibe., vol. 13, 19-89, 1998.
  3. R. M. Dudley, Sample functions of the Gaussian process, Annals of Prob., vol. 1, No. 1, 66-103, 1973.
  4. E. Herbin, From N   -parameter fractional Brownian motions to N   -parameter multifractional Brownian motions, to appear in Rocky Mountains J. Math., 2005.
  5. E. Herbin and J. Lévy-Véhel, Fine analysis of the regularity of Gaussian processes: Stochastic 2   -microlocal analysis, preprint, 2004.
  6. G. Ivanoff, Set-indexed processes: distributions and weak convergence, in: Topics in Spatial Stochastic Processes, Lecture Notes in Mathematics, 1802, 85-126, Springer, 2003.
  7. G. Ivanoff and E. Merzbach, Set-Indexed Martingales, Chapman & Hall/CRC, 2000.
  8. A. Kamont, On the fractional anisotropic Wiener field, Prob. and Math. Stat., vol. 16, fasc. 1, 85-98, 1996.
  9. D. Khoshnevisan, Multiparameter Processes: an introduction to random fields, Springer, 2002.
  10. R. Peltier and J. Lévy-Véhel, Multifractional Brownian motion : definition and preliminary results, Rapport de recherche INRIA 2645, 1995.
  11. B. Pesquet-Popescu, Modélisation bidimensionnelle de processus non stationnaires et application à l'étude du fond sous-marin, thèse de l'ENS Cachan, 1998.
  12. M. Talagrand, Hausdorff measure of trajectories of multiparameter fractional Brownian motion, Annals of Prob., vol. 23, No 2, 767-775, 1995.
  13. G. Samorodnitsky and M. S. Taqqu, Stable non-Gaussian random processes, Chapman & Hall/CRC, 1994.
  14. Y. Xiao, Hausdorff measure of the graph of fractional Brownian motion, Math. Proc. Camb. Phil. Soc., vol. 122, 565-576, 1997.