arxiv:math/ v3 [math.pr] 26 Sep 2005

Size: px
Start display at page:

Download "arxiv:math/ v3 [math.pr] 26 Sep 2005"

Transcription

1 ariv:math/ v3 [math.pr] 6 Sep 005 q WIENER AND RELATED PROCESSES. A BRYC PROCESSES CONTINUOUS TIME GENERALIZATION PAWEL J. SZABLOWSKI Abstract. We define two Markov processes. The finite dimensional distributions of the first one say = t t 0 depend on one parameter q 1, 1 > and of the second one say Y =Y t t R on two parameters q, α 1, 1 > 0,. The first one resembles Wiener process in the sense that for q = 1 it is Wiener process but also that for q < 1 and n 1 t n/ H n t/ t q, where H n n 0 are so called q Hermite polynomials, are martingales. It does not have however independent increments. The second one resemble Orstein- Ulehnbeck processes. For q = 1 it is a classical OU process. For q < 1 it is stationary with correlation function equal to exp α t s. On the way we deny Wesolowski s martingale characterization of Wiener process 1. Introduction As announced in abstract, we are going to define two continuous time families of Markov processes. One of them will resemble Wiener process and the other Orstein- Uhlenbeck OU process. They will be indexed apart from time parameter by additional parameter q 1, 1 > and in case q = 1 those processes will be classical Wiener and OU process. Of course for OU process there will additional parameter α > 0 responsible for covariance function of the process. For q 1, 1 both processes will assume values in a compact space < / 1 q, / 1 q >. One dimesional probabilities and transitional probabilities of these will be given explicitly. Moreover two families of polynomials orthogonal with respect to these measures will also be presented. Some properties of conditional expectations given the past of these processes including martingale properties will also be exposed. Thus quite detailed knowledge concerning these processes will be presented. The problem is that these processes will be introduced via discrete time one dimensional random fields Bryc processes that have been studied in detail recently see e.g. [1], [], [4], [5], [3]. That is why the first section will be dedicated to definition of Bryc processes and recollection of their basic properties. The interesting thing about Bryc and related processes is that their one dimensional densities belong to so called family of q Gaussian distributions introduced Bożejko et. al. in [10] using purely non-commutative probability arguments including special interpretation of parameter q. Hence these distributions and consequently Bryc processes can be viewed as the mysterious connection between commutative and non-commutative probability theories. Since 1997 there appeared Date: July, Mathematics Subject Classification. Primary 60J5, 60G44; Secondary 05A30, 60K40. Key words and phrases. q Wiener process, q Orstein-Uhlenbeck process, polynomial martingales property, q-gaussian distributions, quadratic harnesses. 1

2 PAWEL J. SZABLOWSKI couple of papers mostly in journals dedicated to Physics on the properties of q Gaussian distributions. See e.g. [7], [14], [1]. There is also a different path of research followed by W lodek Bryc and Jacek Weso lowski see e.g. [6], [13]. Their starting point is process with continuous time that satisfies several exactly 5 conditions on covariance function and on first and second conditional moments. Under these assumptions they prove that these processes are Markov and several properties of families of polynomials that orthogonalize transitional and one dimensional probabilities. They give several examples illustrating developed theory. One of the processes considered by them is so called q Brownian process. As far as one dimensional probabilities and transitional probabilities are concerned this process is identical with introduced in this paper q Wiener process. The point is that we derive this process under less assumptions on the first and second conditional moments we need only. The construction is also different. Our starting point is discrete time Bryc process. q Wiener process is obtained as continuous time transformation of continuous time generalization of the discrete time process. Besides we list many properties of q Wiener process that justify its name are sort of q analogies of well known martingale properties of Wiener process. Some of these properties can be derived from the definition of q Brownian process presented in [6]. They are not however stated explicitly. Recently W lodek Bryc and Jacek Weso lowski together with Wojtek Matysiak in [8] follow another more general path and considered so called continuous time quadratic harnesses. Vaguely speaking quadratic harnesses Z =Z t t>0 are such L processes with zero mean and covariance function EZ s Z t = min s, t that have respectively linear and quadratic first and second conditional moments of Z t with respect to field Z :, s > < u, for s < t < u. In their paper they show that each quadratic harness is characterized by 5 parameters apart from time parameter. They show also that under mild technical assumptions for each harness one can define a system of orthogonal polynomials {P n x, t} n 1 such that n : P n Z t, t, Z : t t>0 is a martingale. One of the processes that are defined in this paper, namely q Wiener process, although emerging from different context, turns out to be quadratic harnesses with 4 out of 5 parameters assuming value equal to 0, the family of orthogonal polynomials of its one dimensional distribution satisfies 3 term recursive fomulae given in paper of Bryc Matysiak and Weso lowski. They provide also an example of quadratic harness that is not an independent increment process. So far all the concrete examples given by these there authors concerned independent increment processes. Our approach is totally commutative, classical as the one presented in recent papers of Bryc and Weso lowski. q by no means is an operator. It is number parameter. Yet we are touching q series and special functions. The paper is organized as follows. As we mentioned in the first section we recall definition and summarize basic properties of Bryc processes. The next section is still dedicated mainly to certain properties of discrete time Bryc processes. In the third section first introduce continuous version of Bryc processes then we prove their existence, later parts of this section are devoted to presentation of the two processes that were mentioned earlier. We show their connection with emerging quadratic harness theory. We show here that Weso lowski s martingale characterization contained in [9] of Wiener process is not true. The last section contains proofs of the results from previous sections.

3 GENERALIZATIONS OF BRYC S PROCESSES 3 Although it is still not known if the Bryc process exist with parameter q > 1, we show in this paper that q Wiener and q Ortein-Uhlenbeck processes do not exist for q > 1.. Bryc processes By Bryc processes we mean square integrable random field = { n } n Z indexed by the integers, with non-degenerate covariance matrices and constant first two moments, that satisfy the following two sets of conditions :.1 a, b R; n Z : E n F n = a n 1 + n+1 + b, a.s. and A, B, C R; n Z : E n F n. =.3 A n 1 + n+1 + Bn 1 n+1 + D n 1 + n+1 + C, a.s., where F m := k : k m. Non-singularity of covariance matrices implies that all random variables n are non-degenerate and there is no loss of generality in assuming that E k = 0 and Ek = 1 for all k Z, which implies b = 0. It has been shown in [4] that.1 implies L -stationarity of. Since the case ρ := corr 0, 1 = 0 contains sequences of independent random variables which satisfy.1 and. but can have arbitrary distributions, we shall exclude it from the considerations. Observe that non-singularity of the covariance matrices implies ρ < 1. By Theorem 3.1 from [1] see also Theorems 1 and in [4], corr 0, k = ρ k and one-sided regressions are linear.4 E k..., 1, 0 = ρ k 0. Let us define algebras F m := k : k m, F m := k : k m, F m := k : k m, F m, j := k : k m k j. It turns out that parameters a, ρ, A, B, C are related to one another in the following way:.5 a = ρ 1 + ρ, C = 1 A ρ B, and.6 A In [5] it was shown that ρ + 1ρ + B = 1. Theorem 1. There do not exist standardized random fields = k k Z with nonsingular covariance matrices, satisfying.1 and. and ρ 0, unless D = 0 and either B = 0 or.6 holds with B B 1 B B 3, where { ρ 1 B 1 = ρ 1 + ρ { ρ 1 + ρ B 3 = }, B = 0, ρ 1 + ρ 1 ρ ρ /m ρ /m ρ 4 : m N Furthermore, if D = B = 0 or D = 0 and.6 holds with B B then the relevant fields exist and have the following distributions: ], }.

4 4 PAWEL J. SZABLOWSKI 1 If B = 0 and A = 1/ then k = RY k, R and Y k are independent, R 0, ER = 1 and all random variables Y k are symmetric with values ±1; If B = 0 and A 1/ then all random variables k are symmetric with values ±1; 3 If.6 holds then for all B 0, ρ / 1 + ρ all one-dimensional distributions of are equal to the uniquely determined symmetric and absolutely continuous distribution depending on B supported on a finite interval; 4 If.6 holds and B = ρ / 1 + ρ then k N 0, 1 for all k Z. Theorem. If.6 holds with D = 0 and B B 1 then there do not exist the relevant random fields satisfying.1-. with ρ 0 and such that sequence k k is uniformly integrable and Pr k 1 = 0 < 1. There also i.e. in [5] the following problem has been put up. Problem 1. Do there exist random fields satisfying.1-. with ρ 0, D = 0, and such that.6 holds with B B 3? Theorem 3 below, strengthens Theorem.3 from [] by stating that is always a Markov chain and relaxing assumptions for the uniquely determined case. Theorem 3. If is a standardized random field with non-singular covariance matrices, satisfying.1 and.1-., ρ 0, and such that D = B = 0 or D = 0 and.6 holds with B B, then is a stationary Markov chain. has uniquely determined distributions unless B = D = 0 and A = 1/. In [1] it has been shown regular B-processes are stationary random chains and that one can redefine parameters by introducing new parameter q < 1, 1 > and define parameters A and B as functions of q and ρ. It turns out that q = 1 is equivalent to the B = 0, case q 1, 1 > is equivalent to B B, the case B B 1 leads to undefined q and finally case B B 3 leads to q { 1/ρ n : n N }. It is known that the case q < 1 implies to nonprobabilistic solutions see e.g.[1]. Following.[1] we can express parameters A, B, C using parameters ρ and q : ρ 1 qρ.7a A = ρ + 11 qρ 4,.7b B = ρ 1 ρ 1 + q ρ + 11 qρ 4. Since we have.5 we can also express C as a function of q and ρ : 1 ρ C = 1 qρ 4. Conversely parameter q can be expressed by parameter B by the relationship: ρ 4 + B ρ + ρ q =. 1 + ρ B 4 ρ + 1 ρ 1 We can rephrase and summarize the results of [1] and [4] in the following way. Each Bryc process is characterized by two parameters ρ and q. For q outside the

5 GENERALIZATIONS OF BRYC S PROCESSES 5 set < 1, 1 > { 1/ρ /n : n N } Bryc processes do not exist. For 0 < ρ < 1 and q < 1, 1 > Bryc processes exist and all their finite dimensional distribution are uniquely determined and known. For 0 < ρ < 1 and q { 1/ρ /n : n N } it is not known if Bryc processes exist. Bryc sequences with 0 < ρ < 1 and q 1, 1 > we will denote regular B- processes. To complete recollection of basic properties of regular B-processes let us introduce so called q Hermite polynomials {H n x q} n 1 defined by the following recurrence:.9a.9b H 1 x q = 0, H 0 x q = 1, n 0 : xh n x q = H n+1 x q + [n] q H n 1 x q, where [0] q = 0, [1] q = 1, n [n] q = 1 + q q n 1. W. Bryc in [1] has shown that: n Z, k 1 : EH k n q F n 1 = ρ k H k n 1 q, a.s. and also that {H n x q} n 1 are orthogonal polynomials of distribution of n. The case q = 1 is equivalent to B = 0 and leads to distribution described above. We will concentrate thus on the case q 1, 1 >. It turns out then that this distribution has density for q 1, 1 and x 1 q, 1 q 1 q f H x q = π 4 1 qx 1 + q k 1 qx q k 1 q k+1, k=0 and.10 f H x 1 = 1 exp x, x R. π W. Bryc in [1] also found density of the conditional distribution n n 1 = y and later W. Bryc, W. Matysiak and P. J. Szab lowski in [3] found orthogonal polynomials of this conditional distribution. Namely it turns out that for q 1, 1 and x 1 q, 1 q we have 1 q f x y, q, ρ = π 4 1 qx.11 k=0 k=0 1 ρ q k 1 q k q k 1 qx q k 1 ρ q k 1 qρq k 1 + ρ q k xy + 1 qρ x + y q k, and that polynomials {p n x y, q, ρ} n 1 defined by.1a p 1 x y, q, ρ = 0, p 0 x y, q, ρ = 1,.1b n 0 : p n+1 x y, q = x ρyq n p n x y, q, ρ 1 ρ q n 1 [n] q p n 1 x y q, ρ, are orthogonal with respect to measure defined by the density.11 For q = 1 we have 1.13 f x y, 1, ρ = exp x ρy 1 ρ 1 ρ.

6 6 PAWEL J. SZABLOWSKI In this paper we are going to consider continuous time generalization of. Namely, more precisely, we are going to study the existence and basic properties of process Y = {Y t } t R, such that for every positive δ, random sequences are regular B-processes. δ n = Y nδ, n Z, 3. Auxiliary properties of Bryc processes The following two properties follow assumptions.1 and. and assumptions concerning parameters ρ and A, B, and C. Lemma 1. If is a Bryc process then: 3.1a 3.1b 3.1c E 4 n = + q, E n m = 1 + ρ n m 1 + q, E n n j n+k = ρ k+j E 4 n = ρ j+k + q. Proposition 1. If is a Bryc process then E n F n k, n+j is a linear function of n k and n+j. More precisely we have 3. E n F n k, n+j = ρk 1 ρ j 1 ρ j+k n k + ρj 1 ρ k 1 ρ n+j. j+k Proposition. If is Bryc process, then for n, j, k N, E n F n j, n+k is a linear function of n j, n+k and n j n+k. In particular 3.3 E n F n j, n+k = A 1 jk n j + A jk n+j + B jk n j n+j + C jk where: 3.4a 3.4b 3.4c 3.4d A 1 jk = ρj 1 ρ k 1 qρ k 1 qρ j+k 1 ρ j+k, A jk = ρk 1 ρ j 1 qρ j 1 qρ j+k 1 ρ j+k B jk = q + 1ρj+k 1 ρ j 1 ρ k 1 qρ j+k 1 ρ j+k, C jk = 1 ρ j 1 ρ k 1 qρ j+k. As a simple corollary of Proposition and its proof we have the following Lemma Lemma. Let { i } i Z Bryc process with parameters ρ and q. Let us define Z k = kj+m for fixed j N, m = 0,...,j 1 and k N. {Z k } k Z is also Bryc process with parameters ρ j and q. 4. Continuous time generalizations of regular B processes 4.1. Existence. From consideration { } of the previous section it follows that for each δ > 0 regular B-process n δ is L stationary with covariance function K n, m = ρ δ n m = E δ 0 δ n m n Z. But we have also for any integer k: ρ δ n = E δ 0 δ n = E δ/k 0 δ/k nk = ρ δ/k nk.

7 GENERALIZATIONS OF BRYC S PROCESSES 7 From this identity by standard methods we deduce that ρ δ = ρ δ, for certain positive parameter ρ < 1. Another words we deduce that if Y exists, then it is L stationary with covariance function K t, s = K s t = ρ s t. for some ρ 0, 1. Let us introduce new parameter We have then α = log 1 ρ > 0. K s, t = exp α s t. Existence of Y will be shown for two cases separately. Since for q = 1 we have normality of one dimensional and conditional distributions. Thus we deduce that process Y for q = 1 exists. It is in fact well known Orstein-Uhlenbeck process. Now consider fixed q 1, 1. Let S =< / 1 q, / 1 q >. First we will deduce the existence of the process Ỹ = Y t t Q. This follows from Kolmogorov s extension theorem. Since having natural ordering of Q we need only consistency of the family of finite dimensional distributions of Ỹ. This can be however easily shown by the following argument. Let us take finite set of numbers r 1 < r < < r n from Q. Let M denote the smallest common denominator of these numbers. Let us consider regular B process M with q and ˆρ = ρ 1/M or for parameter α : ˆα = α/m. Note that then n M = Y n/m for n N. Then joint distribution of Y r1,...,y rn is in fact a joint distribution of R M 1,...,R M n where numbers R 1,..., R n are defined by the relationships r i = R i /M. Since process M exists we have consistency since if { 1,..., k } {r 1,...,r n } then distribution of Y 1,..., Y k being equal to the distribution of T M 1,...,T M k with Tj defined by j = T j /M for j = 1,...,k is a projection of the distribution of Y r1,..., Y rn. Hence we deduce that process Ỹ with values in compact space S, exist. Now we use separation theorem and view Ỹ as separable modification of process Y itself. Hence Y exists. 4.. q OU processes. Process Y with parameters q, α will be called continuous time q, α OU-process. OU standing for Orstein-Uhlenbeck. Let us summarize what properties of Y can be deduced from the properties of discrete time regular B-processes presented in [1], [], [4], [5] and [3]. Theorem 4. Let Y be a continuous time q, α OU-process. Then: t R : Y t f H x q, s > t : Y s Y t = y f x y, q, e α s t, Y is a stationary Markov process with f x y, q, e α s t as its transfer function.

8 8 PAWEL J. SZABLOWSKI 3 Y is time symmetric. Moreover we have for any s R and δ, γ > 0 : EY s F s δ, s+γ = e αδ 1 e αγ Y 1 e αδ+γ s δ + e αγ 1 e αδ Y 1 e αδ+γ s+γ, and for s R and δ, γ > 0 : E Y s F s δ, s+γ = A1 δ, γy s δ + A δ, γy s+γ + B δ, γy s δ Y s+γ + C δ, γ, where: A 1 δ, γ = e αδ 1 e αγ 1 qe αγ 1 qe αδ+γ 1 e αδ+γ, A δ, γ = e αγ 1 e αδ 1 qe αδ 1 qe αδ+γ 1 e αδ+γ B δ, γ = q + 1e αδ+γ 1 e αδ 1 e αγ 1 qe αδ+γ 1 e αδ+γ, C δ, γ = 1 e αδ 1 e αγ 1 qe αδ+γ. 4 the following families of polynomials respectively given by.9 and.1 {H n x q} n 0 and { p x y, q, e α s t }, are orthogonal polynomials of n and 4.. That is in particular 5 n 1, t R : EH n Y t q = 0 n 1, s > t : E p n Y s Y t, q, e α s t F t = 0 a.s., n 1, s > t : EH n Y s q = e nα s t H n Y t q a.s., 6 Consequently n 1 : H n Y t q t R is a stationary random process with covariance function and K n s, t = e nα s t nα S n ω α = A n ω + n α, as its spectral density, where A n = EH n Y t q. Remark 1. Thus in particular {H n Y t q} t R more and more resembles white noise the greater n we choose. Remark. Note that following Theorem 1 and considerations of the section 4.1, for the existence of q OU process one needed only the following two conditions apart from linear independence for any m Z, n N and d > 0 of random variables from the set { Y md, Y m+1d,...,y m+nd } : n Z, d > 0 : E Y nd F n 1d, n+1d = e αd 1 + e αd Yn 1d + Y n+1d

9 GENERALIZATIONS OF BRYC S PROCESSES n Z, d > 0 : E Y nd F n 1d, n+1d =  Yn 1d + Y n+1d + ˆBY n 1d Y n+1d +Ĉ, where  = A 1 d, d = A d, d, ˆB = B d, d, Ĉ = C d, d. That is we need only symmetric and discrete for all increments d > 0 versions of condition defining E Ys F s δ, s+γ. Remark 3. Note that continuous time q OU process does not exist for q > 1. It is so because for discrete time Bryc process with parameters q, ρ presuming such process existed - compare Theorem 1 and Problem 1 the following relationship between parameters q and ρ must be satisfied for some integer n : ρ q n = 1 or equivalently ratio log q log ρ must be equal to some integer. However if q OU process existed parameter ρ would depend on time parameter t in the following way ρ = exp αt for some fixed positive α and consequently log q αt would have to be integer for all real t which is impossible q Wiener process. Let Y be given q, α OU process. Let us define: = 0; > 0 : = Y log /α. Process = 0 will be called q Wiener process. Let us also introduce the following filtration: F = t : t = Y t : t log /α = F log /α Following this definition and Theorem 4 of the previous section we have: Theorem 5. Let be a q Wiener process, then it has the following properties: 1, 0 : cov, = min,, > 0 : 1 f H x q 3 = y 1 x + y f y, q, 4 For all n 1 and 0 < we have E n/ H n, q F = n/ 4.6 H n q, a.s E n/ H n q = n/ 4.7 H n q. F Hence n 1 the pair Z n, F 0 is a martingale and the pair 0 U n, where Zn = n/ H n / q,, F, where Un = n/ H n 0 is the reverse martingale. In particular is a martingale and, F is a reversed martingale. >0 q

10 10 PAWEL J. SZABLOWSKI 4.8a 4.8b 5 We have also: 1 q + γ δ δ E F δ, +γ = γ δ + γ δ + δ δ + γ +γ a.s. E F δγ δ, +γ = δ + γ 1 q + γ + qδ + 1 q + qδ +γ γ + q + 1 δ +γ + δ + γ Corollary 1. Let be a q Wiener process and Y related to it q, α OU process. We have i For q 1, 1 both and Y processes do not allow continuous path modifications. ii F =, a.s. > > 0 : E 4 F = > > 0 : E Hence does not have independent increments. iii δγ var F δ, +γ = 1 q + γ + qδ 1 q + + q + 1 q, a.s 1 1 q +γ δ δ + γ + γ δ δ +γ. δ + γ Thus following [?] we deduce from iii that is a quadratic harness with parameters introduced in[?] θ = η = = sigma = 0 thus q commutation equation is very simple in this case, namely [x.y] = I. Remark 4. From assertion 4. of the Theorem 5 it follows that, F, 0, F are martingales and 0 /, F >0 / 1/, F are >0 reversed martingales. Thus if main result of Weso lowski s paper [9] was true we would deduce that is the Wiener process. We will show that Weso lowski made a mistake in his calculations. Following these martingale properties we have:for all 0 < < : E F = a.s., E F =. In order to simplify calculations we will assume that the following equalities are satisfied with probability 1. Hence E F = + or var F =. Further from reversed martingale property we have : E F =, hence E F = and E 1 F = 1, hence E F = +. Now mimicking way of reasoning of Weso lowski we have: E = E E F = E 4 + E = E E F = E 4 + E. Thus: E 4 + = E4 +..

11 GENERALIZATIONS OF BRYC S PROCESSES 11 Dividing both sides by we get: E = E Thus we see that E 4 = c for some constant c. We know that c = + q.see e.g. the proof of Corollary 1 and the fact that H 4 = q + q q + q and the fact that EH 4 = 0 and E =. Hence following Weso lowski we get E 3 = E4 = + q, E = + q +, E 3 = E 4 = + q. So far the Weso lowskis and our calculations yield the same partial results. Now: E 4 = E 4 4E3 + 6E 4E 3 + E4 = + q 4 + q q q + + q = + 3q q = + q + 1 q. In Weso lowski s paper it is wrongly written that E 4 = c + c 3 +, where c = E 4 /. Following Theorem 5 we get: Lemma 3. Let = 0 be a q Wiener process. For s R and 0 < <, we have: E 1 1 q sq k + 1 qs q k 1 F 4.11 = k=0 1 1 q sq k + 1 qs q k 1, a.s. k=0 Remark 5. Bryc and Weso lowski in [6] also introduced q Wiener processes. However the process introduced by them was a particular element of general class of q Meixner processes. To define q Meixner processes they needed 5 assumptions to be satisfied. Hence according to Bryc and Weso lowski q Wiener process is a L process that satisfies 5 the following conditions. Namely: Let t t 0 be L -process with zero expectation and for all 0 s < t < u 1 E t s = min t, s E t F s, u = u t u s s + t s u s t. 3 E t F s, u = A1 s, t, us+a s, t, uu+b s, t, u s u +D s, t, u+ E s, t, u s + F s, t, u s for some nonrandom functions A 1, A, B, C, D, E, F of s, t, u. 4 E t F s = s + t s 5 var t F u = t u u t. The procedure presented in this paper and in particular remark only conditions are needed to define q Wiener process in the class of L processes. Namely, conditions. and 3. for special s = t δ, u = t+γ and assuming A 1 = A and E = F = 0. To see this apply change of time scale defined by 4.5 to conditions 4.3 and 4.4 to obtain conditions 4.8a and 4.8b with γ = expnαdexpαd 1 and δ = exp n 1αdexpαd 1. All other properties mentioned in theorem 5 and assumed in the paper [6] that is conditions 1., 4., 5. follow these two plus

12 1 PAWEL J. SZABLOWSKI linear independence discussed in remark. In particular properties 4. and 5. follow assertion 4.of the Theorem 5 considered for n = 1 and. Bryc & Weso lowski indicate that in general that is if for these particular δ and γ we would have A 1 A and E + F 0, then conditions could not be deduced from. and 3. One of the examples is discussed by them so called q-poisson process. Above mentioned δ, γ, and A 1 = A when applied to time transform defined by 4.5 imply symmetry with respect to time direction of conditions imposed on first two conditional moments. This symmetry assumption thus turns out extremely important, and crucial. Remark 6. Following Remark 3 we deduce that q Wiener processes do not exist for q > Proofs of the results Proof of Lemma Remembering that H 4 x q = x 4 3+q+q x +1+q+q thus consequently x 4 = H 4 x q q + q H x q + + q and we see that En 4 = +q since EH 4 m = EH m = 0. Now using the fact that EH n F m = ρ H m for m < n we get: E n m = E m H n q + 1 = 1 + E meh n q F m = 1 + ρ n m EmH m = 1 + ρ n m E m 4 m To get 3.1c. we have: = 1 + ρ n m + q 1. E n n j n+k = E n n j E n+k F n = ρ k E 3 n n j = ρ k+j E 3 ne n j F n = ρ k+j E 4 n = ρ k+j + q. Proof of Proposition 1. For fixed natural numbers k, j let us denote Z i = E n+i F n k, n+j. Let Z =[Z i ] j 1 i= k+1. Now notice that for each coordinate of vector Z we have: Z i = E i F n k, n+j = EE i F i F n k, n+j = ρ 1 + ρ E i 1 F n k, n+j + E i+1 F n k, n+j ρ 1+ρ Z i 1 + Z i+1 if i k + 1 j 1 ρ = 1+ρ n k + ρ 1+ρ Z k+ if i = k 1. ρ 1+ρ n+j + ρ 1+ρ Z j if i = j 1 Hence we have a vector linear equation Z = AZ+ ρ 1 + ρ n k + n+j, where A = ρ 0 1+ρ... 0 ρ 1+ρ ρ 1+ρ.

13 GENERALIZATIONS OF BRYC S PROCESSES 13 Now notice that matrix I A where I denotes unity matrix is nonsingular. This is so because sum of the absolute values of elements in each row of the matrix is less than 1 which means that eigenvalues of the matrix I A are inside circle at center in 1 and radius less than 1, thus nonzero. Consequently each component of the vector Z is a linear function of n k and n+j. Having linearity of E i F n k, n+j with respect to n k and n+j we get 3.. Proof of Proposition. Let us denote m i,j = E n+i n+j F n l, n+k, for i, j = l, l + 1,...,k 1, k. Notice that using 3. we get m i,j = m j,i and that m l,j = n l E n+j F n l, n+k = ρl+j 1 ρ k j 1 ρ k+l n l + ρk j 1 ρ l+j 1 ρ k+l n l n+k and m k,j = n+k E n+j F n l, n+k = ρl+j 1 ρ k j 1 ρ k+l n l n+k + ρk j 1 ρ l+j 1 ρ k+l n+k m l, l = n 1, m k,k = n+k. Besides we have for i, j = j + 1,...,k 1 and i j m i,j = E n+i n+j F n l, n+k = E n+i E n+j F n+j F n l, n+k = ρ 1 + ρ m i,j 1 + m i,j+1 and if i = j m i,i = E n+i F n l, n+k = E E n+i F n+i Fn l, n+k = Am i 1,i 1 + m i+1,i+1 + Bm i 1,i+1 + C Notice also that we have in fact l + k unknowns and l + k linear equations. Moreover random variables m ij are well defined since conditional expectation is uniquely defined up to set of probability 1. Thus we get the main assertion of the proposition. Now we know that for some A j, B j, C j we have 5.1 E n F n j, n+k = A 1 jk n j + A jk n+k + B jk n j n+k + C jk. First thing to notice is that 1 = A 1 jk + A jk + ρj+k B jk + C jk. Secondly let us multiply 5.1 by n j, n+k and n j n+k and let us take expectation of both sides of obtained in that way equalities. In doing so we apply assertions of Lemma 1. In this way we will get three equations: 1 + ρ j 1 + q = A 1 jk 1 + ρ k 1 + q = A 1 jk + q + A jk + B jk ρ j+k + q + 1 A ρ j+k 1 + q 1 + ρ j+k 1 + q jk A jk ρj+k B jk + A jk + q + B jk ρ j+k + q + 1 A 1 jk A jk ρj+k B jk ρ j+k + q = A 1 jk ρj+k + q + A jk ρj+k + q + B jk 1 + ρ j+k 1 + q + ρ j+k 1 A 1 jk A jk ρj+k B jk. Solution of this system of equations is, as it can be easily, checked 3.4. Secondly consider discrete time random field Z = {Z k } k Z such that Z k = kj. Obviously

14 14 PAWEL J. SZABLOWSKI we have EZ k F k = E kj F k 1j k+1j = ρj 1 + ρ j k 1j + k+1j = ρ j 1 + ρ j Z k 1 + Z k+1, and E Z k F k = E kj F k 1j k+1j = A j k 1j + k+1j + B j k 1j k+1j + C j = A j Z k 1 + Z k+1 + Bj Z k 1 Z k+1 + C j. where A j = A 1 jj = A jj, B j = B jj. Thus Z is a Bryc random field with different parameters. Notice that one dimensional distributions of processes and Z are the same. Hence parameters q for both processes Z and are the same. On the other hand parameter ρ Z of the process Z is related to parameter ρ of the process by the following relationship ρ Z = EZ 0 Z 1 = E 0 j = ρ j. Thus applying formulae.7 we get 3.4. Proof of Theorem 4. The proof of this theorem is in fact contained in above mentioned papers. For the proof of 5. notice that from 4. it follows that n, m 1, s, t R : EH m Y t q H n Y s q = δ n m e nα s t EHn Y t q. Further we use spectral decomposition theorem. Proof of Theorem Suppose <, cov, = EY log /α Y log /α = e α log /α log /α = e log// = =.. & 3. are just assertion 1 of theorem Notice that / = Y log /α.we have using assertion 5 of Theorem 4: E n/ H n / F = n/ E H n Ylog /α Flog /α = n/ e αnlog /α log /α H n Ylog /α = n/ H n /. a.s. Similarly E n/ H n q = n/ E H n Ylog /α q F log /α = F n/ e αlog /α log /α H n Ylog /α q = n/ n Hn Ylog /α q = n/ H n q. 5. to get first part we proceed as follows: E Ylog /α F log δ/α, log+γ/α = δ/ 1 / + γ / + γ1 δ/ Y log δ/α + Y log+γ/α 1 δ/ + γ 1 δ/ + γ = γ δ + γ δ + δ δ + γ +γ

15 GENERALIZATIONS OF BRYC S PROCESSES 15. To get second part we have: E F δ,+γ = E Y δ 1 +γ 1 q +γ = Y 1 δ +γ 1 q δ +γ +γ 1 δ 1 q δ + Y 1 δ 1 q δ q δ +γ 1 δ +γ log /α F log δ/α, log+γ/α +γ 1 +γ +γ 1 δ 1 q δ +γ Keeping in mind that = Y log /α we get: log δ/α log+γ/α Y log δ/α Y log+γ/α + E F δ,+γ γ 1 q + γ = δ + γ 1 q + γ + qδ δ 1 δ 1 +γ. 1 q δ +γ δ 1 q + qδ + δ + γ 1 q + γ + qδ +γ q + 1δγ + δ + γ 1 q + γ + qδ δ +γ δγ + 1 q + γ + qδ. Proof of Corollary 1. Recall that from well known Levy s characterization of Wiener process it follows that a process = t t 0 with continuous a.s. trajectories is a Wiener process iff condition t, F t, t t, F t are martingales. Recall that H x q = x 1. From assertion 3. of the above Proposition it follows that for any q Wiener process with q < 1, both, F and, F are martingales. Hence iff they had continuos trajectories they would be ordinary Wiener processes i.e. q = 1. A contradiction. Since we obtain q, α OU process from q Wiener process by continuous argument transformation, we deduce that any q, α OU process with q < 1 cannot have continuous a.s. trajectories. To prove ii we calculate: E F = E F + = + =. E 4 F = E 4 F 4 E 3 F + 6 E F 4 3 E F + 4 = E 4 F 4 E 3 F = E 4 F 4 E 3 F

16 16 PAWEL J. SZABLOWSKI Now recalling that H 3 x q = x 3 +qx and H 4 x q = x 4 3+q+q x + 1+q+q we see that 3/ H 3 q = E 3/ H 3 q = E 3 F + qe F and H 4 q = E H t 4 q F F = E 4 F 3 + q + q E F q + q. Thus E 3 F = + q q = + q + 3 and E 4 F = 3 + q + q + 1+q+q q+q q + q = q + q q + q. So E 4 F = q + q q + q 1 + q + q q = 1 q + + q 1 + q + q = 1 q + + q + 1 q. If had independent increments then E 4 F would not depend on. Note that for q = 1 E 4 F = 3 = 3E F. iii We have to subtract 4.8b square of 4.8b. var F δ,+γ δγ = 1 q + γ + qδ + + δ γ 1 q + γ δ + γ 1 q + γ + qδ γ δ + γ + +γ δ 1 q + qδ δ δ + γ 1 q + γ + qδ + δ + γ q + 1δγ δ +γ δ + γ 1 q + γ + qδ γδ δ + γ. After some simplifications we get = var F δ,+γ δγ = 1 q + γ + qδ 1 qδγ + γ 1 qδγ δ + γ + δ 1 q + γ + qδ δ + γ + δ 1 q + γ + qδ +γ 1 qγδ δ + γ γ + δ 1 q + γ + qδ δ +γ δγ = 1 q + γ + qδγ + δ γ + δ + 1 q + γ δ + δ +γ δ + γ δ +γ δγ 1 q + γ + qδγ + δ γ + δ + 1 q +γ δ δ +γ + γ δ

17 GENERALIZATIONS OF BRYC S PROCESSES 17 Proof of Lemma 3. We have P t x = P Y log /α x = P Y log /α x/. Knowing that f H is a density of Y t we get immediately assertion 1. To get assertion we have P x = y = P x + y = y = P Y log /α x + y Y log /α = y = P Y log /α x+y Y log /α = y. Now we recall that f x y, q, e α s t is the density of Y t given Y s = s. For q < 1 to get assertion 3 first we multiply both sides of 4.6 by s n / [n] q and sum over n from 0 to. Changing order of summation and conditional expectation and applying formula 1.4 of [3] giving generating function of q Hermite polynomials, leads to For q = 1 we multiply both sides of 4.6 by s n /n! and after summing up with respect to n we apply well known formula for generating functions of Hermite polynomials. References [1] Bryc, W. 001 Stationary random fields with linear regressions. Annals of Probability 9, No. 1, [] Bryc, W. 001 Stationary Markov chains with linear regressions. Stochastic Processes and Applications 93, [3] Bryc, W., Matysiak, W., Szab lowski, P.J. 003 Probabilistic aspects of Al-Salam Chihara polynomials. Proceedings of the American Mathematical Society 133, [4] Matysiak, W., Szab lowski, P.J. 00 A few remarks on Bryc s paper on random fields with linear regressions. Annals of Probability 30, No. 3, [5] Matysiak, W., Szab lowski, P. J. 005, Bryc s Random Fields: The Existence and Distributions Analysis, Ariv:math.PR/math/ [6] Bryc, W., Weso lowski, J. 005, Conditional Moments of $q$-meixner Processes, Probab. Theory Rel. Fields 131, [7] Bryc, W., Bożejko, M. 005, On a Class of Free Levy Laws Related to a Regression Problem, ariv:math.oa/ [8] Bryc, W., Matysiak, W., Weso lowski, J. 005, Quadratic Harnesses, $q$-commutations, and Orthogonal Martingale Polynomials, ariv:math.pr/ [9] Weso lowski, Jacek 1990, A Martingale Characterization of the Wiener Process, Statistics & Probability Letters 10, [10] Bożejko, M. Kümmerer, B., Speicher, R. 1997, $q$-gaussian Processes: No-Commutative and Classical Aspects, Comm. Math. Phys. 1851, [11] Manuy, R., Yor, M. 005, Levy Bridges and Monsieur Jourdain, Stochastic Processes and Their Applications 115, [1] Anshelevich, M. 004, $q$-levy Processes, ariv:math:oa/ [13] Bryc, W., Weso lowski, Jacek, Bi-Poisson Process, arxiv.org:math/pr/ [14] Bryc, W. 001, Classical versions of Gaussian processes. Conditional moments and Bell s inequality, Comm. Math. Phys. 19, Acknowledgement 1. References Acknowledgement. The Author would like to thank W lodek Bryc and Jacek Weso lowski for their critical remarks concerning first version of the paper. Department of Mathematics and Information Sciences, Warsaw University of Technology, pl. Politechniki 1, Warszawa, Poland address: pszablowski@elka.pw.edu.pl

Quadratic harnesses, q-commutations, and orthogonal martingale polynomials. Quadratic harnesses. Previous research. Notation and Motivation

Quadratic harnesses, q-commutations, and orthogonal martingale polynomials. Quadratic harnesses. Previous research. Notation and Motivation Quadratic harnesses, q-commutations, and orthogonal martingale polynomials Quadratic harnesses Slide 1 W lodek Bryc University of Cincinnati August 2005 Based on joint research with W Matysiak and J Wesolowski

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS

PROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please

More information

X n = c n + c n,k Y k, (4.2)

X n = c n + c n,k Y k, (4.2) 4. Linear filtering. Wiener filter Assume (X n, Y n is a pair of random sequences in which the first component X n is referred a signal and the second an observation. The paths of the signal are unobservable

More information

Stochastic process for macro

Stochastic process for macro Stochastic process for macro Tianxiao Zheng SAIF 1. Stochastic process The state of a system {X t } evolves probabilistically in time. The joint probability distribution is given by Pr(X t1, t 1 ; X t2,

More information

Continued fractions for complex numbers and values of binary quadratic forms

Continued fractions for complex numbers and values of binary quadratic forms arxiv:110.3754v1 [math.nt] 18 Feb 011 Continued fractions for complex numbers and values of binary quadratic forms S.G. Dani and Arnaldo Nogueira February 1, 011 Abstract We describe various properties

More information

A matrix over a field F is a rectangular array of elements from F. The symbol

A matrix over a field F is a rectangular array of elements from F. The symbol Chapter MATRICES Matrix arithmetic A matrix over a field F is a rectangular array of elements from F The symbol M m n (F ) denotes the collection of all m n matrices over F Matrices will usually be denoted

More information

4 Sums of Independent Random Variables

4 Sums of Independent Random Variables 4 Sums of Independent Random Variables Standing Assumptions: Assume throughout this section that (,F,P) is a fixed probability space and that X 1, X 2, X 3,... are independent real-valued random variables

More information

An almost sure invariance principle for additive functionals of Markov chains

An almost sure invariance principle for additive functionals of Markov chains Statistics and Probability Letters 78 2008 854 860 www.elsevier.com/locate/stapro An almost sure invariance principle for additive functionals of Markov chains F. Rassoul-Agha a, T. Seppäläinen b, a Department

More information

ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION

ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION Annales Univ. Sci. Budapest., Sect. Comp. 33 (2010) 273-284 ON SUM OF SQUARES DECOMPOSITION FOR A BIQUADRATIC MATRIX FUNCTION L. László (Budapest, Hungary) Dedicated to Professor Ferenc Schipp on his 70th

More information

ULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR ORTHOGONAL POLYNOMIALS

ULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR ORTHOGONAL POLYNOMIALS ULTRASPHERICAL TYPE GENERATING FUNCTIONS FOR ORTHOGONAL POLYNOMIALS arxiv:083666v [mathpr] 8 Jan 009 Abstract We characterize, up to a conjecture, probability distributions of finite all order moments

More information

1 Last time: least-squares problems

1 Last time: least-squares problems MATH Linear algebra (Fall 07) Lecture Last time: least-squares problems Definition. If A is an m n matrix and b R m, then a least-squares solution to the linear system Ax = b is a vector x R n such that

More information

Monte-Carlo MMD-MA, Université Paris-Dauphine. Xiaolu Tan

Monte-Carlo MMD-MA, Université Paris-Dauphine. Xiaolu Tan Monte-Carlo MMD-MA, Université Paris-Dauphine Xiaolu Tan tan@ceremade.dauphine.fr Septembre 2015 Contents 1 Introduction 1 1.1 The principle.................................. 1 1.2 The error analysis

More information

Gaussian Processes. 1. Basic Notions

Gaussian Processes. 1. Basic Notions Gaussian Processes 1. Basic Notions Let T be a set, and X : {X } T a stochastic process, defined on a suitable probability space (Ω P), that is indexed by T. Definition 1.1. We say that X is a Gaussian

More information

Linearization coefficients for orthogonal polynomials. Michael Anshelevich

Linearization coefficients for orthogonal polynomials. Michael Anshelevich Linearization coefficients for orthogonal polynomials Michael Anshelevich February 26, 2003 P n = monic polynomials of degree n = 0, 1,.... {P n } = basis for the polynomials in 1 variable. Linearization

More information

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539

Brownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539 Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory

More information

MATH 205C: STATIONARY PHASE LEMMA

MATH 205C: STATIONARY PHASE LEMMA MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)

More information

arxiv: v2 [math.pr] 4 Sep 2017

arxiv: v2 [math.pr] 4 Sep 2017 arxiv:1708.08576v2 [math.pr] 4 Sep 2017 On the Speed of an Excited Asymmetric Random Walk Mike Cinkoske, Joe Jackson, Claire Plunkett September 5, 2017 Abstract An excited random walk is a non-markovian

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY LINEAR ALGEBRA ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,

More information

P (A G) dp G P (A G)

P (A G) dp G P (A G) First homework assignment. Due at 12:15 on 22 September 2016. Homework 1. We roll two dices. X is the result of one of them and Z the sum of the results. Find E [X Z. Homework 2. Let X be a r.v.. Assume

More information

Boolean Inner-Product Spaces and Boolean Matrices

Boolean Inner-Product Spaces and Boolean Matrices Boolean Inner-Product Spaces and Boolean Matrices Stan Gudder Department of Mathematics, University of Denver, Denver CO 80208 Frédéric Latrémolière Department of Mathematics, University of Denver, Denver

More information

Weak convergence and Brownian Motion. (telegram style notes) P.J.C. Spreij

Weak convergence and Brownian Motion. (telegram style notes) P.J.C. Spreij Weak convergence and Brownian Motion (telegram style notes) P.J.C. Spreij this version: December 8, 2006 1 The space C[0, ) In this section we summarize some facts concerning the space C[0, ) of real

More information

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3 Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................

More information

Lecture 12: Detailed balance and Eigenfunction methods

Lecture 12: Detailed balance and Eigenfunction methods Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),

More information

LARGE DEVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILED DEPENDENT RANDOM VECTORS*

LARGE DEVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILED DEPENDENT RANDOM VECTORS* LARGE EVIATION PROBABILITIES FOR SUMS OF HEAVY-TAILE EPENENT RANOM VECTORS* Adam Jakubowski Alexander V. Nagaev Alexander Zaigraev Nicholas Copernicus University Faculty of Mathematics and Computer Science

More information

A linear algebra proof of the fundamental theorem of algebra

A linear algebra proof of the fundamental theorem of algebra A linear algebra proof of the fundamental theorem of algebra Andrés E. Caicedo May 18, 2010 Abstract We present a recent proof due to Harm Derksen, that any linear operator in a complex finite dimensional

More information

A linear algebra proof of the fundamental theorem of algebra

A linear algebra proof of the fundamental theorem of algebra A linear algebra proof of the fundamental theorem of algebra Andrés E. Caicedo May 18, 2010 Abstract We present a recent proof due to Harm Derksen, that any linear operator in a complex finite dimensional

More information

Proofs for Large Sample Properties of Generalized Method of Moments Estimators

Proofs for Large Sample Properties of Generalized Method of Moments Estimators Proofs for Large Sample Properties of Generalized Method of Moments Estimators Lars Peter Hansen University of Chicago March 8, 2012 1 Introduction Econometrica did not publish many of the proofs in my

More information

Notes on Quadratic Extension Fields

Notes on Quadratic Extension Fields Notes on Quadratic Extension Fields 1 Standing notation Q denotes the field of rational numbers. R denotes the field of real numbers. F always denotes a subfield of R. The symbol k is always a positive

More information

a 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b 2.

a 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b 2. Chapter 1 LINEAR EQUATIONS 11 Introduction to linear equations A linear equation in n unknowns x 1, x,, x n is an equation of the form a 1 x 1 + a x + + a n x n = b, where a 1, a,, a n, b are given real

More information

arxiv: v4 [math.pr] 18 May 2017

arxiv: v4 [math.pr] 18 May 2017 ON STATIONARY MARKOV PROCESSES WITH POLYNOMIAL CONDITIONAL MOMENTS arxiv:1312.4887v4 [math.pr] 18 May 2017 PAWE L J. SZAB LOWSKI Abstract. We study a class of stationary Markov processes with marginal

More information

Eötvös Loránd University, Budapest. 13 May 2005

Eötvös Loránd University, Budapest. 13 May 2005 A NEW CLASS OF SCALE FREE RANDOM GRAPHS Zsolt Katona and Tamás F Móri Eötvös Loránd University, Budapest 13 May 005 Consider the following modification of the Barabási Albert random graph At every step

More information

1 Linear Difference Equations

1 Linear Difference Equations ARMA Handout Jialin Yu 1 Linear Difference Equations First order systems Let {ε t } t=1 denote an input sequence and {y t} t=1 sequence generated by denote an output y t = φy t 1 + ε t t = 1, 2,... with

More information

Lecture 7. Econ August 18

Lecture 7. Econ August 18 Lecture 7 Econ 2001 2015 August 18 Lecture 7 Outline First, the theorem of the maximum, an amazing result about continuity in optimization problems. Then, we start linear algebra, mostly looking at familiar

More information

2. Matrix Algebra and Random Vectors

2. Matrix Algebra and Random Vectors 2. Matrix Algebra and Random Vectors 2.1 Introduction Multivariate data can be conveniently display as array of numbers. In general, a rectangular array of numbers with, for instance, n rows and p columns

More information

. Find E(V ) and var(v ).

. Find E(V ) and var(v ). Math 6382/6383: Probability Models and Mathematical Statistics Sample Preliminary Exam Questions 1. A person tosses a fair coin until she obtains 2 heads in a row. She then tosses a fair die the same number

More information

SOMEWHAT STOCHASTIC MATRICES

SOMEWHAT STOCHASTIC MATRICES SOMEWHAT STOCHASTIC MATRICES BRANKO ĆURGUS AND ROBERT I. JEWETT Abstract. The standard theorem for stochastic matrices with positive entries is generalized to matrices with no sign restriction on the entries.

More information

Supermodular ordering of Poisson arrays

Supermodular ordering of Poisson arrays Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore

More information

Lecture Notes 1: Vector spaces

Lecture Notes 1: Vector spaces Optimization-based data analysis Fall 2017 Lecture Notes 1: Vector spaces In this chapter we review certain basic concepts of linear algebra, highlighting their application to signal processing. 1 Vector

More information

I forgot to mention last time: in the Ito formula for two standard processes, putting

I forgot to mention last time: in the Ito formula for two standard processes, putting I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy

More information

Kernel families of probability measures. Saskatoon, October 21, 2011

Kernel families of probability measures. Saskatoon, October 21, 2011 Kernel families of probability measures Saskatoon, October 21, 2011 Abstract The talk will compare two families of probability measures: exponential, and Cauchy-Stjelties families. The exponential families

More information

Generalized Gaussian Bridges of Prediction-Invertible Processes

Generalized Gaussian Bridges of Prediction-Invertible Processes Generalized Gaussian Bridges of Prediction-Invertible Processes Tommi Sottinen 1 and Adil Yazigi University of Vaasa, Finland Modern Stochastics: Theory and Applications III September 1, 212, Kyiv, Ukraine

More information

CAAM 454/554: Stationary Iterative Methods

CAAM 454/554: Stationary Iterative Methods CAAM 454/554: Stationary Iterative Methods Yin Zhang (draft) CAAM, Rice University, Houston, TX 77005 2007, Revised 2010 Abstract Stationary iterative methods for solving systems of linear equations are

More information

Maximum Likelihood Drift Estimation for Gaussian Process with Stationary Increments

Maximum Likelihood Drift Estimation for Gaussian Process with Stationary Increments Austrian Journal of Statistics April 27, Volume 46, 67 78. AJS http://www.ajs.or.at/ doi:.773/ajs.v46i3-4.672 Maximum Likelihood Drift Estimation for Gaussian Process with Stationary Increments Yuliya

More information

Notes on Blackwell s Comparison of Experiments Tilman Börgers, June 29, 2009

Notes on Blackwell s Comparison of Experiments Tilman Börgers, June 29, 2009 Notes on Blackwell s Comparison of Experiments Tilman Börgers, June 29, 2009 These notes are based on Chapter 12 of David Blackwell and M. A.Girshick, Theory of Games and Statistical Decisions, John Wiley

More information

Other properties of M M 1

Other properties of M M 1 Other properties of M M 1 Přemysl Bejda premyslbejda@gmail.com 2012 Contents 1 Reflected Lévy Process 2 Time dependent properties of M M 1 3 Waiting times and queue disciplines in M M 1 Contents 1 Reflected

More information

A Concise Course on Stochastic Partial Differential Equations

A Concise Course on Stochastic Partial Differential Equations A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY LINEAR ALGEBRA ELEMENTARY LINEAR ALGEBRA K. R. MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Second Online Version, December 1998 Comments to the author at krm@maths.uq.edu.au Contents 1 LINEAR EQUATIONS

More information

A CLT FOR MULTI-DIMENSIONAL MARTINGALE DIFFERENCES IN A LEXICOGRAPHIC ORDER GUY COHEN. Dedicated to the memory of Mikhail Gordin

A CLT FOR MULTI-DIMENSIONAL MARTINGALE DIFFERENCES IN A LEXICOGRAPHIC ORDER GUY COHEN. Dedicated to the memory of Mikhail Gordin A CLT FOR MULTI-DIMENSIONAL MARTINGALE DIFFERENCES IN A LEXICOGRAPHIC ORDER GUY COHEN Dedicated to the memory of Mikhail Gordin Abstract. We prove a central limit theorem for a square-integrable ergodic

More information

For a stochastic process {Y t : t = 0, ±1, ±2, ±3, }, the mean function is defined by (2.2.1) ± 2..., γ t,

For a stochastic process {Y t : t = 0, ±1, ±2, ±3, }, the mean function is defined by (2.2.1) ± 2..., γ t, CHAPTER 2 FUNDAMENTAL CONCEPTS This chapter describes the fundamental concepts in the theory of time series models. In particular, we introduce the concepts of stochastic processes, mean and covariance

More information

RATE-OPTIMAL GRAPHON ESTIMATION. By Chao Gao, Yu Lu and Harrison H. Zhou Yale University

RATE-OPTIMAL GRAPHON ESTIMATION. By Chao Gao, Yu Lu and Harrison H. Zhou Yale University Submitted to the Annals of Statistics arxiv: arxiv:0000.0000 RATE-OPTIMAL GRAPHON ESTIMATION By Chao Gao, Yu Lu and Harrison H. Zhou Yale University Network analysis is becoming one of the most active

More information

MATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5

MATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5 MATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5.. The Arzela-Ascoli Theorem.. The Riemann mapping theorem Let X be a metric space, and let F be a family of continuous complex-valued functions on X. We have

More information

Stein s Method and Characteristic Functions

Stein s Method and Characteristic Functions Stein s Method and Characteristic Functions Alexander Tikhomirov Komi Science Center of Ural Division of RAS, Syktyvkar, Russia; Singapore, NUS, 18-29 May 2015 Workshop New Directions in Stein s method

More information

A PRIMER ON SESQUILINEAR FORMS

A PRIMER ON SESQUILINEAR FORMS A PRIMER ON SESQUILINEAR FORMS BRIAN OSSERMAN This is an alternative presentation of most of the material from 8., 8.2, 8.3, 8.4, 8.5 and 8.8 of Artin s book. Any terminology (such as sesquilinear form

More information

The properties of L p -GMM estimators

The properties of L p -GMM estimators The properties of L p -GMM estimators Robert de Jong and Chirok Han Michigan State University February 2000 Abstract This paper considers Generalized Method of Moment-type estimators for which a criterion

More information

Harmonic Analysis. 1. Hermite Polynomials in Dimension One. Recall that if L 2 ([0 2π] ), then we can write as

Harmonic Analysis. 1. Hermite Polynomials in Dimension One. Recall that if L 2 ([0 2π] ), then we can write as Harmonic Analysis Recall that if L 2 ([0 2π] ), then we can write as () Z e ˆ (3.) F:L where the convergence takes place in L 2 ([0 2π] ) and ˆ is the th Fourier coefficient of ; that is, ˆ : (2π) [02π]

More information

Kaczmarz algorithm in Hilbert space

Kaczmarz algorithm in Hilbert space STUDIA MATHEMATICA 169 (2) (2005) Kaczmarz algorithm in Hilbert space by Rainis Haller (Tartu) and Ryszard Szwarc (Wrocław) Abstract The aim of the Kaczmarz algorithm is to reconstruct an element in a

More information

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2016 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing

More information

Stochastic Realization of Binary Exchangeable Processes

Stochastic Realization of Binary Exchangeable Processes Stochastic Realization of Binary Exchangeable Processes Lorenzo Finesso and Cecilia Prosdocimi Abstract A discrete time stochastic process is called exchangeable if its n-dimensional distributions are,

More information

On Least Squares Linear Regression Without Second Moment

On Least Squares Linear Regression Without Second Moment On Least Squares Linear Regression Without Second Moment BY RAJESHWARI MAJUMDAR University of Connecticut If \ and ] are real valued random variables such that the first moments of \, ], and \] exist and

More information

We simply compute: for v = x i e i, bilinearity of B implies that Q B (v) = B(v, v) is given by xi x j B(e i, e j ) =

We simply compute: for v = x i e i, bilinearity of B implies that Q B (v) = B(v, v) is given by xi x j B(e i, e j ) = Math 395. Quadratic spaces over R 1. Algebraic preliminaries Let V be a vector space over a field F. Recall that a quadratic form on V is a map Q : V F such that Q(cv) = c 2 Q(v) for all v V and c F, and

More information

SAMPLE CHAPTERS UNESCO EOLSS STATIONARY PROCESSES. E ξ t (1) K.Grill Institut für Statistik und Wahrscheinlichkeitstheorie, TU Wien, Austria

SAMPLE CHAPTERS UNESCO EOLSS STATIONARY PROCESSES. E ξ t (1) K.Grill Institut für Statistik und Wahrscheinlichkeitstheorie, TU Wien, Austria STATIONARY PROCESSES K.Grill Institut für Statistik und Wahrscheinlichkeitstheorie, TU Wien, Austria Keywords: Stationary process, weak sense stationary process, strict sense stationary process, correlation

More information

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing

More information

1. Introduction to commutative rings and fields

1. Introduction to commutative rings and fields 1. Introduction to commutative rings and fields Very informally speaking, a commutative ring is a set in which we can add, subtract and multiply elements so that the usual laws hold. A field is a commutative

More information

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim

More information

UC Berkeley Department of Electrical Engineering and Computer Sciences. EECS 126: Probability and Random Processes

UC Berkeley Department of Electrical Engineering and Computer Sciences. EECS 126: Probability and Random Processes UC Berkeley Department of Electrical Engineering and Computer Sciences EECS 6: Probability and Random Processes Problem Set 3 Spring 9 Self-Graded Scores Due: February 8, 9 Submit your self-graded scores

More information

ON KRONECKER PRODUCTS OF CHARACTERS OF THE SYMMETRIC GROUPS WITH FEW COMPONENTS

ON KRONECKER PRODUCTS OF CHARACTERS OF THE SYMMETRIC GROUPS WITH FEW COMPONENTS ON KRONECKER PRODUCTS OF CHARACTERS OF THE SYMMETRIC GROUPS WITH FEW COMPONENTS C. BESSENRODT AND S. VAN WILLIGENBURG Abstract. Confirming a conjecture made by Bessenrodt and Kleshchev in 1999, we classify

More information

What is A + B? What is A B? What is AB? What is BA? What is A 2? and B = QUESTION 2. What is the reduced row echelon matrix of A =

What is A + B? What is A B? What is AB? What is BA? What is A 2? and B = QUESTION 2. What is the reduced row echelon matrix of A = STUDENT S COMPANIONS IN BASIC MATH: THE ELEVENTH Matrix Reloaded by Block Buster Presumably you know the first part of matrix story, including its basic operations (addition and multiplication) and row

More information

7 Curvature of a connection

7 Curvature of a connection [under construction] 7 Curvature of a connection 7.1 Theorema Egregium Consider the derivation equations for a hypersurface in R n+1. We are mostly interested in the case n = 2, but shall start from the

More information

IEOR 4701: Stochastic Models in Financial Engineering. Summer 2007, Professor Whitt. SOLUTIONS to Homework Assignment 9: Brownian motion

IEOR 4701: Stochastic Models in Financial Engineering. Summer 2007, Professor Whitt. SOLUTIONS to Homework Assignment 9: Brownian motion IEOR 471: Stochastic Models in Financial Engineering Summer 27, Professor Whitt SOLUTIONS to Homework Assignment 9: Brownian motion In Ross, read Sections 1.1-1.3 and 1.6. (The total required reading there

More information

Large Sample Properties of Estimators in the Classical Linear Regression Model

Large Sample Properties of Estimators in the Classical Linear Regression Model Large Sample Properties of Estimators in the Classical Linear Regression Model 7 October 004 A. Statement of the classical linear regression model The classical linear regression model can be written in

More information

Notes on Linear Algebra and Matrix Theory

Notes on Linear Algebra and Matrix Theory Massimo Franceschet featuring Enrico Bozzo Scalar product The scalar product (a.k.a. dot product or inner product) of two real vectors x = (x 1,..., x n ) and y = (y 1,..., y n ) is not a vector but a

More information

Krzysztof Burdzy University of Washington. = X(Y (t)), t 0}

Krzysztof Burdzy University of Washington. = X(Y (t)), t 0} VARIATION OF ITERATED BROWNIAN MOTION Krzysztof Burdzy University of Washington 1. Introduction and main results. Suppose that X 1, X 2 and Y are independent standard Brownian motions starting from 0 and

More information

Chebyshev coordinates and Salem numbers

Chebyshev coordinates and Salem numbers Chebyshev coordinates and Salem numbers S.Capparelli and A. Del Fra arxiv:181.11869v1 [math.co] 31 Dec 018 January 1, 019 Abstract By expressing polynomials in the basis of Chebyshev polynomials, certain

More information

Foundations of Matrix Analysis

Foundations of Matrix Analysis 1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the

More information

Detailed Proof of The PerronFrobenius Theorem

Detailed Proof of The PerronFrobenius Theorem Detailed Proof of The PerronFrobenius Theorem Arseny M Shur Ural Federal University October 30, 2016 1 Introduction This famous theorem has numerous applications, but to apply it you should understand

More information

Independence of some multiple Poisson stochastic integrals with variable-sign kernels

Independence of some multiple Poisson stochastic integrals with variable-sign kernels Independence of some multiple Poisson stochastic integrals with variable-sign kernels Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological

More information

A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University

A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University Lecture 19 Modeling Topics plan: Modeling (linear/non- linear least squares) Bayesian inference Bayesian approaches to spectral esbmabon;

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY LINEAR ALGEBRA ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Second Online Version, December 998 Comments to the author at krm@mathsuqeduau All contents copyright c 99 Keith

More information

CHUN-HUA GUO. Key words. matrix equations, minimal nonnegative solution, Markov chains, cyclic reduction, iterative methods, convergence rate

CHUN-HUA GUO. Key words. matrix equations, minimal nonnegative solution, Markov chains, cyclic reduction, iterative methods, convergence rate CONVERGENCE ANALYSIS OF THE LATOUCHE-RAMASWAMI ALGORITHM FOR NULL RECURRENT QUASI-BIRTH-DEATH PROCESSES CHUN-HUA GUO Abstract The minimal nonnegative solution G of the matrix equation G = A 0 + A 1 G +

More information

This operation is - associative A + (B + C) = (A + B) + C; - commutative A + B = B + A; - has a neutral element O + A = A, here O is the null matrix

This operation is - associative A + (B + C) = (A + B) + C; - commutative A + B = B + A; - has a neutral element O + A = A, here O is the null matrix 1 Matrix Algebra Reading [SB] 81-85, pp 153-180 11 Matrix Operations 1 Addition a 11 a 12 a 1n a 21 a 22 a 2n a m1 a m2 a mn + b 11 b 12 b 1n b 21 b 22 b 2n b m1 b m2 b mn a 11 + b 11 a 12 + b 12 a 1n

More information

SMSTC (2007/08) Probability.

SMSTC (2007/08) Probability. SMSTC (27/8) Probability www.smstc.ac.uk Contents 12 Markov chains in continuous time 12 1 12.1 Markov property and the Kolmogorov equations.................... 12 2 12.1.1 Finite state space.................................

More information

Elementary linear algebra

Elementary linear algebra Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The

More information

Analysis-3 lecture schemes

Analysis-3 lecture schemes Analysis-3 lecture schemes (with Homeworks) 1 Csörgő István November, 2015 1 A jegyzet az ELTE Informatikai Kar 2015. évi Jegyzetpályázatának támogatásával készült Contents 1. Lesson 1 4 1.1. The Space

More information

Mathematical Methods wk 2: Linear Operators

Mathematical Methods wk 2: Linear Operators John Magorrian, magog@thphysoxacuk These are work-in-progress notes for the second-year course on mathematical methods The most up-to-date version is available from http://www-thphysphysicsoxacuk/people/johnmagorrian/mm

More information

A NOTE ON FREE QUANTUM GROUPS

A NOTE ON FREE QUANTUM GROUPS A NOTE ON FREE QUANTUM GROUPS TEODOR BANICA Abstract. We study the free complexification operation for compact quantum groups, G G c. We prove that, with suitable definitions, this induces a one-to-one

More information

Characterizations of free Meixner distributions

Characterizations of free Meixner distributions Characterizations of free Meixner distributions Texas A&M University August 18, 2009 Definition via Jacobi parameters. β, γ, b, c R, 1 + γ, 1 + c 0. Tridiagonal matrix {(β, b, b,...), (1 + γ, 1 + c, 1

More information

Some functional (Hölderian) limit theorems and their applications (II)

Some functional (Hölderian) limit theorems and their applications (II) Some functional (Hölderian) limit theorems and their applications (II) Alfredas Račkauskas Vilnius University Outils Statistiques et Probabilistes pour la Finance Université de Rouen June 1 5, Rouen (Rouen

More information

STAT 7032 Probability Spring Wlodek Bryc

STAT 7032 Probability Spring Wlodek Bryc STAT 7032 Probability Spring 2018 Wlodek Bryc Created: Friday, Jan 2, 2014 Revised for Spring 2018 Printed: January 9, 2018 File: Grad-Prob-2018.TEX Department of Mathematical Sciences, University of Cincinnati,

More information

arxiv: v1 [math.pr] 1 Jan 2016

arxiv: v1 [math.pr] 1 Jan 2016 A CHARACTERIZATION OF THE NORMAL DISTRIBUTION BY THE INDEPENDENCE OF A PAIR OF RANDOM VECTORS arxiv:1601.00078v1 [math.pr] 1 Jan 2016 WIKTOR EJSMONT Abstract. Kagan and Shalaevski [11] have shown that

More information

SEPARABLE TERM STRUCTURES AND THE MAXIMAL DEGREE PROBLEM. 1. Introduction This paper discusses arbitrage-free separable term structure (STS) models

SEPARABLE TERM STRUCTURES AND THE MAXIMAL DEGREE PROBLEM. 1. Introduction This paper discusses arbitrage-free separable term structure (STS) models SEPARABLE TERM STRUCTURES AND THE MAXIMAL DEGREE PROBLEM DAMIR FILIPOVIĆ Abstract. This paper discusses separable term structure diffusion models in an arbitrage-free environment. Using general consistency

More information

ACO Comprehensive Exam October 14 and 15, 2013

ACO Comprehensive Exam October 14 and 15, 2013 1. Computability, Complexity and Algorithms (a) Let G be the complete graph on n vertices, and let c : V (G) V (G) [0, ) be a symmetric cost function. Consider the following closest point heuristic for

More information

Lecture 13: Simple Linear Regression in Matrix Format. 1 Expectations and Variances with Vectors and Matrices

Lecture 13: Simple Linear Regression in Matrix Format. 1 Expectations and Variances with Vectors and Matrices Lecture 3: Simple Linear Regression in Matrix Format To move beyond simple regression we need to use matrix algebra We ll start by re-expressing simple linear regression in matrix form Linear algebra is

More information

SOLUTION OF GENERALIZED LINEAR VECTOR EQUATIONS IN IDEMPOTENT ALGEBRA

SOLUTION OF GENERALIZED LINEAR VECTOR EQUATIONS IN IDEMPOTENT ALGEBRA , pp. 23 36, 2006 Vestnik S.-Peterburgskogo Universiteta. Matematika UDC 519.63 SOLUTION OF GENERALIZED LINEAR VECTOR EQUATIONS IN IDEMPOTENT ALGEBRA N. K. Krivulin The problem on the solutions of homogeneous

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY LINEAR ALGEBRA ELEMENTARY LINEAR ALGEBRA K. R. MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND Corrected Version, 7th April 013 Comments to the author at keithmatt@gmail.com Chapter 1 LINEAR EQUATIONS 1.1

More information

Meixner matrix ensembles

Meixner matrix ensembles Meixner matrix ensembles W lodek Bryc 1 Cincinnati April 12, 2011 1 Based on joint work with Gerard Letac W lodek Bryc (Cincinnati) Meixner matrix ensembles April 12, 2011 1 / 29 Outline of talk Random

More information

Spectral inequalities and equalities involving products of matrices

Spectral inequalities and equalities involving products of matrices Spectral inequalities and equalities involving products of matrices Chi-Kwong Li 1 Department of Mathematics, College of William & Mary, Williamsburg, Virginia 23187 (ckli@math.wm.edu) Yiu-Tung Poon Department

More information

STAT 100C: Linear models

STAT 100C: Linear models STAT 100C: Linear models Arash A. Amini June 9, 2018 1 / 56 Table of Contents Multiple linear regression Linear model setup Estimation of β Geometric interpretation Estimation of σ 2 Hat matrix Gram matrix

More information

7 Poisson random measures

7 Poisson random measures Advanced Probability M03) 48 7 Poisson random measures 71 Construction and basic properties For λ 0, ) we say that a random variable X in Z + is Poisson of parameter λ and write X Poiλ) if PX n) e λ λ

More information

Unit Roots in White Noise?!

Unit Roots in White Noise?! Unit Roots in White Noise?! A.Onatski and H. Uhlig September 26, 2008 Abstract We show that the empirical distribution of the roots of the vector auto-regression of order n fitted to T observations of

More information