SUPPLEMENT TO ROP: MATRIX RECOVERY VIA RANK-ONE PROJECTIONS 1

Size: px
Start display at page:

Download "SUPPLEMENT TO ROP: MATRIX RECOVERY VIA RANK-ONE PROJECTIONS 1"

Transcription

1 SUPPLEMENT TO ROP: MATRIX RECOVERY VIA RANK-ONE PROJECTIONS BY T. TONY CAI and ANRU ZHANG Deartment of Statstcs The Wharton School Unversty of Pennsylvana We rove n ths sulement the techncal lemmas used n the roofs of the man results. Proof of Lemma.. Suose X : R! R n s gven by (.6), we consder rank- matrces: A = e ( ) e ( )T, A = () () k () k k ()T k Here e ( ),e ( ) are the - and - dmensonal vectors wth frst entry and others, resectvely. Then we have ka k F = ka k F =. EkX (A )k = E ( () ) ( () ) = n VarkX (A )k = Var ( () ) ( () ) = nvar By Chebyshev s nequalty, we have for all t>, =8n. (.3) P kx (A )k tn ale On the other hand, kx (A )k =k () k k () k + ( = k () k k () k ( ) 8 n(t ) ()T () k () k ) ( ( ) ()T () k () k ) The research was suorted n art by NSF FRG Grant DMS , NSF Grant DMS-898 and NIH Grant R CA

2 By Lemma n [8], we know P Hence, ( ) / ale ex( /4), P ( ) / ale ex( /4) (.4) P kx (A )k /4 ale ex( /4) + ex( /4). Combnng (.3) and (.4), we can see (.5) C /C kx (A )k /( nka k F ) kx (A )k /( nka k F ) = s kx (A )k kx (A )k r /4 tn holds wth robablty at least e /4 e /4 8 n(t ). Proof of Lemma 3.. Frst, the common used defnton for sub-gaussan dstrbuton of random varable X nclude the followng two. (.6) 9c, C >, such that P ( X t) ale C ex( ct ) (.7) 9c >, such that Ee tx ale ex(c t /) Suose P s fnte, then we have Ee tx = X k= t k (k)! EXk ale X k= t k (k)! k (k )!! = X k= (t ) k k k! namely X s sub-gaussan. Now suose X s sub-gaussan, then EX k =(k) = kc c k Z Z P ( X >t)t k (ct ) k dt ale kc e ct d(ct )= k!c c k whch mles that P ale max(c, )/c s fnte. Z t k e ct dt =ex(( t) /) max(c, ) k ale (k )!! c Proof of Lemma 7.. Wthout confuson, we smly use to reresent P n the roof. Note that we can multly A by a scale wthout loss of generalty. So we assume throughout the roof that kak F =. We ll rove ths lemma by stes. Frst, we show an nequalty on the even moments of T A ; next, we gve a bound on the moment generaton functon of T A. Fnally, we gve the desred tal bound.

3 . Ste : Even moments of T A. Assume that x =(x,,x ),y =(y,,y ) are two random..d. standard normal dstrbuted vectors. Based on the defnton of P n (3.), we know 3 (.8) E k = E j k ale k Ex k = k Eyj k ; E k = E j k = Ex k = Eyj k = Consder the exanson of E( T A ) k, where the non-zero terms can be wrtten as ky Y Y A l,j l E s E t j j. l= Here s + + s = t + t = k. By(.8), ths term can be bounded as j= ky l= Y A l,j l E s Y j= E t j j ale ky l= Y A l,j l 4k Ex s Y j= Ey t j j The rght hand sde s exact the term n the exanson of 4k E(x T A abs y) k, where A abs s the the element-wse absolute value of A. Therefore, we have (.9) E[ T A ] k ale 4k E[x T A abs y] k. Now we suose A abs has sngular value decomoston A abs = X a u v T = Udag(a)V T where U, V are orthogonal and a =(a,,a ) s the sngular value vector of A abs. A well-known fact s that P a = ka abs k F = kak F. Snce x, y are standard normal dstrbuted, we can see that x T A abs y and x T dag(a)y has the same dstrbuton. So Next, we note X E[x T A abs y] k = E[ a x y ] k z = X a y x q P j= a j y j,

4 4 q P then z s standard normal dstrbuted and ndeendent of j= a j y j snce a y A =. q P j= a j y j and z gven y,,y s always standard normal dstrbuted. For nteger k, v k E[x T A abs y] k u X X =E z t = E z k E( a jyj ) k j= =(k )!! = ale X k,k,,k, k + +k=k X k,k,,k, k + +k=k =((k )!!) =((k )!!) ( a j y j Together wth (.9), we have (.) E[ X k,k,,k, k + +k=k k! Q j= (k j)! k! Q j= k j! X k! Q j= k j! Y Y k,k,,k, k + +k=k X a ) k =((k T A ] k ale 4k ((k a k j= Y E(a y ) k Y (k )!! a k ((k + k + k ) )!! k! Q j= k j! Y (a ) k )!!) kak k F )!!) kak k F. Log-moment generaton functon of T A. By the bound of the even moments of T A, we can also gve the estmate of odd moments, for nteger k, ale E T A k+ ale q E[ T A ] k E[ T A ] k+ ale 4k+ (k )!! (k + )!! ale 4k+ (k + )! Also, E T A k ale 4k ((k )!!) ale 4k (k)!

5 So for all k ale t<,, E T A k ale k k!. Denote µ = E T A, then for X Ee t T A t k =+ k! E T A k ale +tµ+ For k= <t<, we have Ee t T A = + tµ + X t k 4k k= X k= 5 X t k k =+tµ+ t 4 t k= t k+ (k + )! E T A k+ ale+tµ + t 4 t 4 ale +tµ + t 4 Hence, we have for all / <t</, Ee t T A ale +tµ + t 4 t. Note that log( + x) ale +x for all (.) log E ex(t( T A µ)) = log E ex(t t <x<, wehave T A ) tµ ale t 4 t. for all / <t</. 3. The tal bound of kx (A)k /n. Fnally, we estmate the tal bound of kx (A)k/n. Note that kx (A)k /n (j)t A (j) A /n, j= based on (.), the logarthm of moment generatng functon of kx (A)k /n satsfes t log E ex (t(kx (A)k /n µ)) = n log E ex n ( T A µ) ale t 4 /n t /n By the roof of Lemma n [8], we know kx (A)k /n bound, P (kx (A)k /n µ (x/n + x/n)) ale ex( x) P (kx (A)k /n µ ale (x/n + x/n)) ale ex( x) µ has the tal Fnally we set = x/n, by Lemma 7., /(3 4 ) ale µ ale, we fnsh the roof of Lemma.

6 6 Proof of Lemma 7.. Snce P s symmetrc and of varance, we have E = E j =, E = E j =, E 4 = E j 4 ale 3 4 P for all, j. Thenby some exansons and calculatons, E T A = X,j A j j A = X,j E A j By the frst art n the roof of Lemma 7., wehave By Hölder s nequalty, E E T A ale T A 4 ale 9 8 PkAk 4 F q E T A = kak F j = X,j A j = kak F whch gves the rght of the orgnal nequalty. For the left, note that So E E Proof of Lemma 7.3. T A ale E T A T A /3 E T A 4 /3 s (E T A ) 3 kak F E T A 4 3 P 4. We frst rove the sub-gaussan art of the lemma. The moment generatng functon of z (t ) and z ( ale t</( )) satsfy Ee t z = Ee tz = ale + Z Z ex(t )dp ( X )=+ t ex(t ale +t ex(t /) Z ale +tex t / ale ex(t /) +t ale ex t /+ t ale + =+ Z Z /( ))d ex( ( ex(t )dp ( X )=+ 4t ex t /( ) t Z Z (/( ) t) d P ( X )d ex(t ) ) /( ))d P ( X )d ex(t )

7 Then the moment generatng functon of kzk /n and kzk /n satsfes Ee tkzk/n = Ee tkz n k/n ale ex t /(n)+ t Ee tkzk /n t/n n t = + /( ale ex ) t/n /( ) t/n Hence for C, P (kzk /n C) ale E ex(tkzk /n) ale ex t ex (tc) =ex@ n t + n( C) 7 /(n)+t C! n( C) A For C>, we can set t = n(c P (kzk /n C) ale ex,then ) n(c! ) Now we consder the tal bound of kzk /n. For C>4, P (kzk /n C) = E ex(tkzk /n) t C/n t(c 4 ) =ex ex(tc). t/n We set t = C 4 4C /n, P kzk /n C ale ex Fnally we consder kzk, n(c 4 ) n(c 4 ) 8 C( ale ex t/n) 8 C P (kzk ale C log n ) ale n ex( C log n /( )) = n (C / ) Next, we consder the Gaussan art of the lemma. The bound of kzk s already gven by Lemma 5. n [7]. For kzk, we can see E z =, Z E z = xe x /() dx = / Hence, E(kzk /n) = /, Var(kzk /n) = Var( z )/n = ( / ) /n. By Chebyshev nequalty, P (kzk /n ) ale P kzk /n / ( / ) ale Var(kzk /n) ( = + / / ) ( / )n

8 8 For the bound of kzk,wehave P (kzk log n ) ale ale n P ( z ale log n ) log n ex( 4 log n) = n log n Proof of Lemma 7.4. Agan wthout confuson, we smly use to reresent P n the roof. The roof also requres some knowledge of moment generaton functon and "-net method. We ll rove by stes. Moment Generaton Functon of a T X (z)b. Suose a R, b R are fxed unt vectors. In order to handle the oerator norm of X (z), we frst consder a T X (z)b. Note that (.) a T X (z)b = z a T () ()T b Denote X = a T (), Y = b T (),then{x }, {Y } are two ndeendent sets of..d. sub-gaussan samles. Moreover by (), () are..d. from symmetrc dstrbuton P, one can show E(X ) k = E( X j a j j ) k =, X E(X ) k = E( a j j ) k = ale j= X k + k =k k! k! k! X k + k =k Y k! k! k! a k X = k E( a x ) k ale k (k )!! Y E(x k ) k! a k E( k )!

9 d Here x N(, ). Smlarly, E(Y ) k =, E(Y ) k ale k (k )!!. Then for t < /, (.3) X t k E(X Y ) k X t k ( k (k )!!) E ex(tx Y )= ale k! (k)! k= k= X (t ) k (k )!! X / = k = (t ) k ( ) k k! k = k= t 4 Now for fxed z R n, the logarthm of the moment generatng functon of a T X (z)b satsfes log E ex(ta T X (z)b) = ale ale k= log E ex(tz X Y ) ale t z 4 ( t z 4 ) ale t kzk 4 ( t kzk 4 ) t kzk 4 ( t kzk ) for any t < /(kzk ). Here we used the fact that log( x) = X x ale X x = x x 9 log( t z 4 ) for ale x<. Tal Bound of a T X (z)b By the roof of Lemma n [8], we know for fxed z R n, a T X (z)b has tal bound: for x> and fxed a, b, wehave P a T X (z)b kzk x + kzk x ale ex( x); P a T X (z)b ale kzk x + kzk x ale ex( Set x = C( + ), we have (.4) P a T X (z)b kzk C( + )+kzk C( + ) aleex( C( + )). x).

10 For convenence, We denote (.5) T = kzk (C( + )) + kzk C( + ). "-net and the uer bound of kx (z)k. In ths ste, we stll fx z. Weusethe"-net method to derve the uer bound of kx (z)k, whch s gven by kx (z)k = su a T X (z)b ar,br From Lemma.5 n [39], we can fnd an "-net A n the unt shere of R,.e. for all a n the unt shere of R,thereexstsa A such that ka ak ale. Besdes, A ale( + /"). Smlarly, there exsts -net B of the unt ball of R such that B ale( + /"). By (.4), we have (.6) P a T X (z)b T,9a A, b B ale (+/") + ex( C( + )) Now we consder under the event that a T X (z)b ale T,8a A, b B. Suose µ = kx (z)k = max kak =kbk = a T X (z)b, (a,b ) = arg max a,b a T X (z)b, then we can fnd a A, b B such that ka a kale ", kb b kale". Then, µ = a T X (z)b = a T X (z)b + (a a ) T X (z)b + a T X (z)(b b ) alet +(ka a k + kb b k ) kx (z)k alet +"µ Ths means µ ale T/( "). Therefore, when a T X (z)b ale T,8a A, b B, wehavekx (z)k alet/( "). Fnally, we set " =/3, under the event that a T X (z)b ale T = kzk (C( + )) + kzk C( + ), 8a A, b B we have kx (z)k alet/( ") ale 3 kzk C( + )+kzk C( + ) By (.6), the robablty that all the event haen s at least ex( (C log 7)( + )) Ths fnshed the roof of the lemma.

11 Proof of Prooston 7... In the roof, we wll use to reresent P wthout any confuson. The deas of the roof of Prooston 7. follows from [36], [3]. Denote S kr = {X R : rank(x) ale kr, kxk F =}. By Lemma 3. n [3], for any ">, there exsts "-net Skr such that S kr ale(9/")( + +)kr. For gven C,C such that C < /(3 4 ), C >, we set C = C +/(3 4 ), C = C +. We can choose small enough such that /(3 + 4 ) C ale mn, C then by Lemma 7., for any gven A R,wehaveC ale kx (A)k /n ale C wth robablty at least ex( n). Hence, P C ale kx (A)k /n ale C, for all A S r (9/") kr( + +) ex( n) Next, we ll estmate the bound of kx (A)k /n on the whole set S kr rovded that C ale kx (A)k /n ale C for all A S kr.defne ale = nf kx (A)k /n and ale = su kx (A)k /n. AS kr AS kr For any A S kr,thereexstsa S kr such that ka A k F ale ". So kx (A)k /n ale kx (A )k /n+kx (A A )k /n ale C kak+ale ka A k F ale C +ale " kx (A)k /n kx (A )k /n kx (A A )k /n C kak ale ka A k F C ale " whch mean ale = su AS kr kx (A)k F ale C + "ale, ale = nf AS kr kx (A)k F C "ale namely, ale ale C /( "), ale C "ale. We choose C " ale mn, /(3 4 ) C, C C by some calculatons we can see ale C, ale ale C. To sum u, we can choose, " only deendng on C,C,,toensure C ale ale = nf AS kr kx (A)k /n ale su AS kr kx (A)k /n = ale ale C

12 wth robablty at least (9/") kr( + +) ex( n). The last ste s to estmate the robablty above. We choose D 8k log(9/")/, then for n Dr( + ), we have n/ 4 log(9/")kr( + ) log + log(9/")kr( + + ), (9/") kr( + +) e ex( n/). n = ex( n + log + kr( + + ) log(9/")) Fnally, we fnsh the roof of the Theorem by choosng ale /. Proof of Lemma The dea of the roof s orgnated from [4, 3]. We rovde the roof here for the comleteness of the aer. Note = mn(, ), suose for any matrx B, (B) sthe-th largest sngular value of B. By Lemma n [3], we have ka max(r) k + ka max(r) k X =kak ka k = ka ( R)k (A) ( R) rx ( (A) (R)) + X ( (R) (A)) =r+ =ka max(r) k ka max(r) k + kr max(r) k kr max(r) k whch mles (7.3). Proof of Lemma Suose R = A A,thenwehave (.7) kx (R)k /n ale Snce ka k ale kak. By Lemma 7.7, wemusthave(7.3). Suose = mn(, ) and R has the sngular value decomoston, R = P u v T = Udag(~ )V T,then~ max(kr) satsfes k~ max(kr) k ale kr, k~ max(kr) k =k~ max(r) k ( r+ + + kr ) Set alek~ max(r) k +ka max(r) k (k )r kr = max kr, (k~ max(r) k +ka max(r) k (k )r kr )/(kr),

13 then k~ max(kr) k ale, k~ max(kr) k ale kr. Smlarly to the roof of Theorem., aly Lemma 7.6, we can get b () R n,, =,,N such that ~ max(kr) = P N ub () and (7.). Hence, kb () k ale q kb () k kb () k ale q (k~ max(r) k +ka max(r) k (k )r kr ) If = kr, we can otmze over kr n the nequalty, q kb () k ale kr(k~ max(r) k +ka max(r) k (k )r kr ) ale k~ max(r)k +ka max(r) k ; r(k ) f =(k~ max(r) k +ka max(r) k (k )r kr )/(kr), we have 3 (.8) kb () k ale k~ max(r)k +ka max(r) k (k )r kr kr ale k~ max(r)k +ka max(r) k kr. Snce k, we always have (.8). Next, we defne B = Udag(b () )V T,then the rank of B are at most kr, P N B = R max(kr) and kb k F = kb () k. Then (.9) kx (R)k /n kx (R max(kr) )k /n kx (R max(kr) )k /n C kr max(kr) k F C kr max(kr) k F C kr max(kr) k F N X X N kx ( B )k /n C kb k F C kr max(r) k +ka max(r) k kr C kr max(kr) k F C k kr max(kr) k F C ka max(r) k kr, where the last nequalty s due to kr max(kr) k F kr max(r) k F rkrmax(r) k. Therefore, (.3) kr max(kr) k F ale C C / k + ka max(r)k r( kc /C )

14 4 Fnally, kr max(kr) k F =k~ q ale max(kr) k ale q k~ max(kr) k k~ max(kr) k kr (k~ max(r) k r(k ) kr ) ale k~ max(r)k r(k ) ale k~ max(r)k +ka max(r) k r(k ) ale kr max(r)k F k + ka max(r)k r(k ) Therefore, q krk F = kr max(kr) k F + kr max(kr)k F v u ale t krmax(kr) k F + kr max(r)k F k + ka! max(r)k r(k ) s (.3) ale + 4(k ) kr max(kr)k F + ka max(r)k r(k ) ale + kr 8(k ) max(kr) k F + ka max(r)k r(k ) ale C C / 3 k + kc /C + k ka max(r) k r Proof of Lemma The roof of ths theorem s smlar to the roof of Lemma 7.8. Suose R = A A. In ths case we have (.3) kx X (R)k ale nstead of (.7). Besdes, snce ka k alekak and Lemma 7.7, westllhave (7.3). Wth the smlar argument as (.9), we have (.33) krk hr, X X (R) = kx (R)k kx (R)k /n C C ka max(r) k n C kr max(kr) k F kr max(kr) k F k kr + Here (x) + means max(x, ). Besdes, krk ale kr max(r) k +(kr max(r) k +ka max(r) k ) ale ale rkrmax(r) k F + ka max(r) k rkrmax(kr) k F + ka max(r) k

15 Suose x = kr max(kr) k F, y = ka max(r) k / r. Based on the revous two nequaltes, we have When x (.34) n (C C / k)x C y ale r(x + y) k + C y k(c C /, the nequalty above leads to k) n(c C / k) x n(c C / k) C k y + r x r y ale. 5 Note that for second order nequalty ax bx c ale, a>, b, c, we have x ale b+ b +4ac a ale b/a + c/a. Hence we can get an uer bound of x from (.34). x ale ale r n(c C / k) + 4C / ky (C C / k) + r n(c C / k) + 4C / ky (C C / k) + r y n(c C / k) r n(c C / k) + y. Hence whenever x C y k(c C / k) or not, (.35) kr max(kr) k F = x ( ale max C y k(c C / k), 3 r 3 r ale n(c C / k) + 4C / (C C / k) + n(c C / k) + 4C / k k (C C / k) +! ka max(r) k. r Fnally, smlarly to (.3) n Lemma 7.8, we can get the uer bound of krk F, krk F ale + kr 8k 8 max(kr) k F + ka max(r)k r(k ) 4 r(" + ) ale (C C / k) 5 + n kc /C + ka max(r) k + k r! y ) whch fnshed the roof of lemma 7.9.

16 6.7. Proof of Proostons.,.3 and.4.. We frst show Prooston.. Denote X, X : R! R b n c such that (.36) [X ] (B) =( ( ) + () ) T B( ( ) () ), =,, b n c (.37) [X ](B) =( ( ) () ) T B( ( ) + () ), =,, b n c Note that ( ) + () and ( ) () are ndeendent..d. standard normal samles, so both X and X are ROP desgn (see (.6)). By Corollary., we know there exsts unform constant C such that whenever b n c Cr, X satsfes the followng roerty wth robablty at least ex( n ), (.38) 8A {A R : rank(a) ale r}, A = arg mn kbk subject to X (B) =X (A). BR Now we consder the event that (.38) holds. We note that for any symmetrc matrx B, [X ] (B) = ( )T B ( ) ()T B () = ([X ] (B) [X ] (B)) So X (B) =X (A) mlesx (B) =X (A) for symmetrc A and B. Also, snce A s feasble n rogrammng (.3), we have kak mn BS subject to X (B) =X (A) mn BR subject to X (B) =X (A) =kak, So we can conclude that A can be exactly recovered by (.3) gven (.38) holds. In summary, for n 6Cr, wth robablty at least ex( n ), X satsfes (.38), then rogrammng (.3) can recover all A S of rank at most r. Next, we consder Prooston.3. The dea of the roof s smlar to Prooston.. Defne z R b n c such that (.39) z = z z, =,, b n c.

17 7 Then z d N(, ). We shall also ont out two facts, z =ỹ X (A) and X = X + X (defned as (.36), (.37)). By Lemma 7.3, weknow P (kzk /n > ) ale 9 n Hence, P (k zk > n) ale bn/c P (k zk > log n) ale bn/c P (A s NOT n the feasble set of rogrammng (.7)) =P kzk /n >, or k X ( z)k > =P kzk /n >, or k X ( z)k > 4 n + 48 log n alep (kzk /n > )+P X (k zk > n)+p (k zk > log n) + P X k X ( z)k > 4k zk + k zk ale 9 n + bn/c + P X ale 5 n kx ( z)k > k zk +6 k zk + P X kx ( z)k > k zk +6 k zk +4ex( ( log 7)) When A s n the feasble set of rogrammng (.7), we have kâk alekak and (.4) k X (Â A)k = b n c X [X ] (Â A) [X ] (Â A) ale kx (Â A)k aleky X(Â)k + kx (A) yk aleky X(Â)k + kzk ale n (.4) k X X ( Â A)k alek X (ỹ X ( Â))k + k X ( X (A) ỹ)k =k X (ỹ X ( Â))k + k X ( z)k ale Smlarly as the roof to Prooston., by Theorem., there exsts constant D, such that f n Dr, X satsfes RUB of order r wth constants C,C such that C /C < wth robablty at least ex( n ). Now we suose the followng two events haen,

18 8. X satsfes RUB of order r and constants C,C satsfyng C /C <,. A s feasble n the rogrammng (.7). Snce X (B) =X (B) for any symmetrc matrx B, byx satsfes RUB condton, we have X satsfes RUB for symmetrc matrces of order r and constants C, C satsfyng (C )/(C ) <. We note that the roof of Lemmas 7.8 and 7.9 stll aly for X n the symmetrc matrces class, so we can get (.8) based on (.4) and (.4) under those two events haen. Fnally the robablty that these events haen s at least 5/n 4ex( ) ex( n ) for <mn(( log 7), ), whch fnshed the roof of Prooston.3. Fnally we consder Prooston.4. Denote = b/c, r = br/c. By r,, we have r r/3, /3. Defne a sub-class of the class rank-r symmetrc matrces, G = 8 >< A S : A = >: ale B B T,B R ( ), rank(b) ale r 9 >= >; we can see 8A G, [X (A)] = ()T A () = ( (),, () )B( () +,, so n G the SROP model becomes y =( (),, () )B( () +,, () ) T + z, z d N(, () ) T, /4) whch s an ROP model whch we already dscussed n secton. We omt the rest of the roof as t can be followed by the roof of Theorem.4. Proof of Prooston 3... The roof can follow from the roof of Prooston.3 and Theorem 3. once we can rove that n hgh robablty, X (defned n (.36)) satsfes RUB condton wth C,C such that C /C s bounded. Ths can be roved smlarly as Prooston 7., where we only need to edt the roof that we use the followng Lemma. nstead of Lemma 7.. Lemma.. Suose A R s a fxed matrx (not necessarly symmetrc) and X s gven by (.36). () s a set of -dmensonal vectors such

19 that d P, where P s some symmetrc varance sub-gaussan dstrbuton excet Rademacher ± dstrbuton. Then for >, we have! mn 3/ (Var(P )/, ) 3( P ) 4 8 P 4 P kak F ale kx (A)k (.4) bn/c ale 3/ P +8 P +4 P kak F. wth robablty at least ex( bn/c). The roof of Lemma. s n the Aendx rght after ths aragrah. Note that rovded P s symmetrc and wth varance, Var(P ) = f and only P s Rademacher ± and A s dagonal, n whch the lower bound of (.4) becomes meanngless. So we only exclude Rademacher ± dstrbuton from the result. Proof of Lemma... The roof of Lemma. s bascally the same to Lemma 7.. We only need to redo two arts of the roof, where there are major d erences.. Part. Ste. Even moments of ( () + () ) T A( () () ).. Frst, based on P s symmetrc and wth varance, we have EP k+ =, EP k ale k Ex k = k (k )!!.Then we can calculate that 9 E( () + () ) k+ = E( () () ) k+ = = E( () + () ) k kx l= k l = k kx l= Smlarly, E( () E( () ) l E( () ) (k l) ale k (k )!! k k! l l! k l (k l)! = k (k )!! kx l= kx l= k (l )!!((k l) )!! l k = k k (k )!! l () ) k ale k k (k )!!. Next, we can smlarly consder the exanson of E ( () + () )A( () () ) k,wherethe non-zero terms can be wrtten as k Y Y ky A l,j l E( () + () ) s l= j= E( () j () j ) t j

20 Here s + + s = t + + t = k. Ths term can be bounded as k ale k ale k Y ky A l,j l l= ky l= ky l= = 4k ky l= ale 4k ky l= A l,j l A l,j l A l,j l A l,j l E( () + () ) s ( () Y s s + t E () ) t () + () s +t + Y s +t (s +t ) ((s + t ) )!! Y ky = 4k 4 A l,j l l= ((s + t ))! s ale 4k +t (s + t )! ky l= A l,j l Y Y (s +t ) ((s + t )!) ky s ale 4k 4 A +t s l,j!t! l Y Ex s Y Ey t l= t () E s + t ((s +t )!!) s +t s!t! Y (s )!(t )! (s )!!(t )!! Here we assume that x,y d N(, ). The rght hand sde of the nequalty above s exactly the term n the exanson of ( ) 4k E(x T A abs y) k, where A abs s the element-wse absolute value of A. Therefore, we have E () + () T k A () () ale ( ) 4k E[x T A abs y] k Now we follow the same argument of the rest art of Ste n the roof of Lemma 7., we can rove that (.43) E () + () T k A () () ale ( ) 4k ((k )!!) kak k F () s +t. Part. The uer and lower bound of µ = E () + () T A () (). To follow the argument of the roof of Ste 3 n Lemma 7., weneedto derve a new bound for µ = E () + () T A () ().Frst,

21 we denote M = () + () T A () (),then µ =EM ale EM s = E ( () + () ) T A( () () ) s X = E( () + () ) 4 A j ( () j,j s X = ale A j + Var(P ) X 6=j r 3 PkAk F On the other hand, EM =E ( () + () ) T A( () () ) () j ) A ale s X 6=j A j P X A = X 6=j A j + Var(P ) X A mn( Var(P ), )kak F, By (.43), we also have EM 4 ale 9( ) 8 kak k F. By Hölder s nequalty, EM ale (EM) /3 (EM 4 ) /3. Hence, r (EM µ ) 3 mn 3/ (Var(P )/, )kak F EM 4 3( ) 4 To sum u, nstead of Lemma 7., we have the bound of µ as follows, (.44) mn 3/ (Var(P )/, )kak F 3( P ) 4 ale µ ale 3/ PkAk F The rest of the roof follows the roof of Lemma 7. wth modfcatons of constants, whch we do not go nto detals. Proof of Lemma 7... Note that () N(, ()T () ), we can assume that (.45) = ()T () Z, d where Z (N(, )) and Z, () are ndeendent. Based on the defnton of z n (4.3), we have (.46) z = y [X ( )] = () () = ()T () (Z ), =,,n.

22 We also denote Q = C n, Q = C n 4, Q 3 = C 3 log n max alealen. We ll frst consder the former art of (7.8), (7.9) and (7.). Suose Z (N(, )). It s well known that the non-central m-th moment of Z s (m )!!, so we have (.47) E (C Z Z ) C qe Z = C (.48) E (C Z Z ) ale E (C Z) + E (Z ) =3C + (.49) E CZ (Z ) =3C (.5) E CZ (Z ) ale C 4 E (Z) 4 +E (Z ) 4 = 5C 4 +6 Next we consder the random quadratc form of. Suose = (,, ) d N(, ), X, d N(, ), ( ),, ( ) are the egenvalues of. Snce s ostve defnte, we have (.5) E T =tr( )=k k E( T ) =E( X + X j j ) <j = X E 4 + X <j jj E j + X 4 je <j j =( X,j j)+( X ) =k k F + k k Hence, (.5) k k ale E( T ) ale 3k k (.53) X E( T ) 4 = E( ( )X ) 4 = X ( ) j ( ) s ( ) t ( )EX Xj Xs Xt ale ale,j,s,tale X ale,j,s,tale ( ) j ( ) s ( ) t ( )7!! = 5k k 4

23 Then we consder C z and C 4 z.by(.45) and (.46), we have C z = ()T () (C Z Z ), C 4 z = ()T () CZ (Z ), whle () and Z are ndeendent n the equaton above. By (.47)- (.53), we obtan an estmaton of the frst and second moment of these two quanttes as (.54) E C z (C )k k (.55) Var C z ale E C z =E (C Z Z ) E ()T () ale 9C +6 k k 3 (.56) E C 4 z 3C k k Var C 4 z ale E C 4 z (.57) ale 5(5C 4 + 6)k k 4 We note that Q kzk /n and Q kzk /n are the average of n..d. coy of C z and C 4 z. We can mmedately get an estmaton of the mean and varance of Q kzk /n and Q kzk /n based on (.54)-(.57). Fnally, by Chebyshev s nequalty, P (Q alekzk /n) ale Var(Q kzk /n) (E (Q kzk /n)) ale 9C +6 n(c ) P Q ale kzk ale Var Q kzk /n n E Q kzk /n ale 5 5C4 + 6 n(3c ) Snce C 3 > and n 3, we know C 3 log nz Z wth robablty. Suose = arg max ()T (),then P (Q 3 alekzk ) alep max C 3 log n( ()T () )Z ale max ( ()T () )(Z ) + P max C 3 log n( ()T () )Z ale max ( ()T () )( Z ) ale+p max C 3 log n( ()T () )Z ale max ()T () alep C 3 log n ()T () Z ale ()T ( ) alep Z ale = P N(, ) ale C 3 log n ale C3 log n C3 log n

24 4 Then we consder the latter art of (7.8)-(7.). We can do smlar calculatons as the frst art of the roof and get EQ = C E = C E( ()T () )=C k k Var(Q )= C n Var ale C n E 4 = C n E ()T () EZ ale 9C n k k EQ = C EZ E ()T () ale 9C k k Var(Q )= n Var C 4 ale C4 n E 8 = C4 n EZ4 E ()T () 4 5 C 4 ale k k 4 n So by Chebyshev s nequalty, P (Q M C k k ) ale P (Q EQ (M )C k k ) ale P Q M Ck k ale P Q EQ (M 9)Ck k Var(Q ) 5 ale (M 9) C 4 ale k k4 n(m 9) 9 (M ) n whch rovde the latter art of (7.8) and (7.9). Fnally we note that =( ()T () )Z. By Lemma n [8] and the fact that k k = P q P ( ), k k = max ( ), k k F = we have (.58) P =P T alen M 3 and ( )X k k + M 3 log nk k ( ), k k F ale k k k k, k k + M 3 log nk k k k +M 3 log nk k v u ( )(X ) t M3 log n ( ) + M 3 log n max P (Z M 3 log n) ale ex( M 3 log n/) = n M 3.! ( ) A

25 Hence, P C 3 log n C 3 M 3 log n k k + M 3 log nk k ale n M 3, 5 and consequently P C 3 log n max alealen C 3 M 3 log n k k + M 3 log nk k ale n M3+, whch gves the rght sde of (7.). Deartment of Statstcs The Wharton School Unversty of Pennsylvana Phladelha, PA, 94. E-mal: tca@wharton.uenn.edu URL: htt://www-stat.wharton.uenn.edu/ tca/ E-mal: anrzhang@wharton.uenn.edu URL: htt://www-stat.wharton.uenn.edu/ anrzhang/

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Exerments-I MODULE III LECTURE - 2 EXPERIMENTAL DESIGN MODELS Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 2 We consder the models

More information

Not-for-Publication Appendix to Optimal Asymptotic Least Aquares Estimation in a Singular Set-up

Not-for-Publication Appendix to Optimal Asymptotic Least Aquares Estimation in a Singular Set-up Not-for-Publcaton Aendx to Otmal Asymtotc Least Aquares Estmaton n a Sngular Set-u Antono Dez de los Ros Bank of Canada dezbankofcanada.ca December 214 A Proof of Proostons A.1 Proof of Prooston 1 Ts roof

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Supplementary material: Margin based PU Learning. Matrix Concentration Inequalities

Supplementary material: Margin based PU Learning. Matrix Concentration Inequalities Supplementary materal: Margn based PU Learnng We gve the complete proofs of Theorem and n Secton We frst ntroduce the well-known concentraton nequalty, so the covarance estmator can be bounded Then we

More information

CIS 700: algorithms for Big Data

CIS 700: algorithms for Big Data CIS 700: algorthms for Bg Data Lecture 5: Dmenson Reducton Sldes at htt://grgory.us/bg-data-class.html Grgory Yaroslavtsev htt://grgory.us Today Dmensonalty reducton AMS as dmensonalty reducton Johnson-Lndenstrauss

More information

Supplementary Material for Spectral Clustering based on the graph p-laplacian

Supplementary Material for Spectral Clustering based on the graph p-laplacian Sulementary Materal for Sectral Clusterng based on the grah -Lalacan Thomas Bühler and Matthas Hen Saarland Unversty, Saarbrücken, Germany {tb,hen}@csun-sbde May 009 Corrected verson, June 00 Abstract

More information

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD Matrx Approxmaton va Samplng, Subspace Embeddng Lecturer: Anup Rao Scrbe: Rashth Sharma, Peng Zhang 0/01/016 1 Solvng Lnear Systems Usng SVD Two applcatons of SVD have been covered so far. Today we loo

More information

Math 217 Fall 2013 Homework 2 Solutions

Math 217 Fall 2013 Homework 2 Solutions Math 17 Fall 013 Homework Solutons Due Thursday Sept. 6, 013 5pm Ths homework conssts of 6 problems of 5 ponts each. The total s 30. You need to fully justfy your answer prove that your functon ndeed has

More information

Solutions Homework 4 March 5, 2018

Solutions Homework 4 March 5, 2018 1 Solutons Homework 4 March 5, 018 Soluton to Exercse 5.1.8: Let a IR be a translaton and c > 0 be a re-scalng. ˆb1 (cx + a) cx n + a (cx 1 + a) c x n x 1 cˆb 1 (x), whch shows ˆb 1 s locaton nvarant and

More information

Lecture 9: Converse of Shannon s Capacity Theorem

Lecture 9: Converse of Shannon s Capacity Theorem Error Correctng Codes: Combnatorcs, Algorthms and Alcatons (Fall 2007) Lecture 9: Converse of Shannon s Caacty Theorem Setember 17, 2007 Lecturer: Atr Rudra Scrbe: Thanh-Nhan Nguyen & Atr Rudra In the

More information

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

Eigenvalues of Random Graphs

Eigenvalues of Random Graphs Spectral Graph Theory Lecture 2 Egenvalues of Random Graphs Danel A. Spelman November 4, 202 2. Introducton In ths lecture, we consder a random graph on n vertces n whch each edge s chosen to be n the

More information

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Solutions to exam in SF1811 Optimization, Jan 14, 2015 Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable

More information

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

More information

PARTIAL QUOTIENTS AND DISTRIBUTION OF SEQUENCES. Department of Mathematics University of California Riverside, CA

PARTIAL QUOTIENTS AND DISTRIBUTION OF SEQUENCES. Department of Mathematics University of California Riverside, CA PARTIAL QUOTIETS AD DISTRIBUTIO OF SEQUECES 1 Me-Chu Chang Deartment of Mathematcs Unversty of Calforna Rversde, CA 92521 mcc@math.ucr.edu Abstract. In ths aer we establsh average bounds on the artal quotents

More information

Managing Capacity Through Reward Programs. on-line companion page. Byung-Do Kim Seoul National University College of Business Administration

Managing Capacity Through Reward Programs. on-line companion page. Byung-Do Kim Seoul National University College of Business Administration Managng Caacty Through eward Programs on-lne comanon age Byung-Do Km Seoul Natonal Unversty College of Busness Admnstraton Mengze Sh Unversty of Toronto otman School of Management Toronto ON M5S E6 Canada

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Exerments-I MODULE II LECTURE - GENERAL LINEAR HYPOTHESIS AND ANALYSIS OF VARIANCE Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 3.

More information

Singular Value Decomposition: Theory and Applications

Singular Value Decomposition: Theory and Applications Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real

More information

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1 Problem Set 4 Suggested Solutons Problem (A) The market demand functon s the soluton to the followng utlty-maxmzaton roblem (UMP): The Lagrangean: ( x, x, x ) = + max U x, x, x x x x st.. x + x + x y x,

More information

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could

More information

Differentiating Gaussian Processes

Differentiating Gaussian Processes Dfferentatng Gaussan Processes Andrew McHutchon Aprl 17, 013 1 Frst Order Dervatve of the Posteror Mean The posteror mean of a GP s gven by, f = x, X KX, X 1 y x, X α 1 Only the x, X term depends on the

More information

Estimation: Part 2. Chapter GREG estimation

Estimation: Part 2. Chapter GREG estimation Chapter 9 Estmaton: Part 2 9. GREG estmaton In Chapter 8, we have seen that the regresson estmator s an effcent estmator when there s a lnear relatonshp between y and x. In ths chapter, we generalzed the

More information

Strong Markov property: Same assertion holds for stopping times τ.

Strong Markov property: Same assertion holds for stopping times τ. Brownan moton Let X ={X t : t R + } be a real-valued stochastc process: a famlty of real random varables all defned on the same probablty space. Defne F t = nformaton avalable by observng the process up

More information

First Year Examination Department of Statistics, University of Florida

First Year Examination Department of Statistics, University of Florida Frst Year Examnaton Department of Statstcs, Unversty of Florda May 7, 010, 8:00 am - 1:00 noon Instructons: 1. You have four hours to answer questons n ths examnaton.. You must show your work to receve

More information

MAT 578 Functional Analysis

MAT 578 Functional Analysis MAT 578 Functonal Analyss John Qugg Fall 2008 Locally convex spaces revsed September 6, 2008 Ths secton establshes the fundamental propertes of locally convex spaces. Acknowledgment: although I wrote these

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Modelli Clamfim Equazione del Calore Lezione ottobre 2014 CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g

More information

CSCE 790S Background Results

CSCE 790S Background Results CSCE 790S Background Results Stephen A. Fenner September 8, 011 Abstract These results are background to the course CSCE 790S/CSCE 790B, Quantum Computaton and Informaton (Sprng 007 and Fall 011). Each

More information

arxiv: v1 [math.co] 1 Mar 2014

arxiv: v1 [math.co] 1 Mar 2014 Unon-ntersectng set systems Gyula O.H. Katona and Dánel T. Nagy March 4, 014 arxv:1403.0088v1 [math.co] 1 Mar 014 Abstract Three ntersecton theorems are proved. Frst, we determne the sze of the largest

More information

1 Matrix representations of canonical matrices

1 Matrix representations of canonical matrices 1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:

More information

The Prncpal Component Transform The Prncpal Component Transform s also called Karhunen-Loeve Transform (KLT, Hotellng Transform, oregenvector Transfor

The Prncpal Component Transform The Prncpal Component Transform s also called Karhunen-Loeve Transform (KLT, Hotellng Transform, oregenvector Transfor Prncpal Component Transform Multvarate Random Sgnals A real tme sgnal x(t can be consdered as a random process and ts samples x m (m =0; ;N, 1 a random vector: The mean vector of X s X =[x0; ;x N,1] T

More information

An application of generalized Tsalli s-havrda-charvat entropy in coding theory through a generalization of Kraft inequality

An application of generalized Tsalli s-havrda-charvat entropy in coding theory through a generalization of Kraft inequality Internatonal Journal of Statstcs and Aled Mathematcs 206; (4): 0-05 ISS: 2456-452 Maths 206; (4): 0-05 206 Stats & Maths wwwmathsjournalcom Receved: 0-09-206 Acceted: 02-0-206 Maharsh Markendeshwar Unversty,

More information

ALGEBRA HW 7 CLAY SHONKWILER

ALGEBRA HW 7 CLAY SHONKWILER ALGEBRA HW 7 CLAY SHONKWILER 1 Whch of the followng rngs R are dscrete valuaton rngs? For those that are, fnd the fracton feld K = frac R, the resdue feld k = R/m (where m) s the maxmal deal), and a unformzer

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem.

Lecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem. prnceton u. sp 02 cos 598B: algorthms and complexty Lecture 20: Lft and Project, SDP Dualty Lecturer: Sanjeev Arora Scrbe:Yury Makarychev Today we wll study the Lft and Project method. Then we wll prove

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0 MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector

More information

First day August 1, Problems and Solutions

First day August 1, Problems and Solutions FOURTH INTERNATIONAL COMPETITION FOR UNIVERSITY STUDENTS IN MATHEMATICS July 30 August 4, 997, Plovdv, BULGARIA Frst day August, 997 Problems and Solutons Problem. Let {ε n } n= be a sequence of postve

More information

Lecture 4: Universal Hash Functions/Streaming Cont d

Lecture 4: Universal Hash Functions/Streaming Cont d CSE 5: Desgn and Analyss of Algorthms I Sprng 06 Lecture 4: Unversal Hash Functons/Streamng Cont d Lecturer: Shayan Oves Gharan Aprl 6th Scrbe: Jacob Schreber Dsclamer: These notes have not been subjected

More information

MATH Homework #2

MATH Homework #2 MATH609-601 Homework #2 September 27, 2012 1. Problems Ths contans a set of possble solutons to all problems of HW-2. Be vglant snce typos are possble (and nevtable). (1) Problem 1 (20 pts) For a matrx

More information

2.3 Nilpotent endomorphisms

2.3 Nilpotent endomorphisms s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms

More information

However, since P is a symmetric idempotent matrix, of P are either 0 or 1 [Eigen-values

However, since P is a symmetric idempotent matrix, of P are either 0 or 1 [Eigen-values Fall 007 Soluton to Mdterm Examnaton STAT 7 Dr. Goel. [0 ponts] For the general lnear model = X + ε, wth uncorrelated errors havng mean zero and varance σ, suppose that the desgn matrx X s not necessarly

More information

P exp(tx) = 1 + t 2k M 2k. k N

P exp(tx) = 1 + t 2k M 2k. k N 1. Subgaussan tals Defnton. Say that a random varable X has a subgaussan dstrbuton wth scale factor σ< f P exp(tx) exp(σ 2 t 2 /2) for all real t. For example, f X s dstrbuted N(,σ 2 ) then t s subgaussan.

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

arxiv: v2 [math.st] 17 Aug 2014

arxiv: v2 [math.st] 17 Aug 2014 ROP: Matrx Recovery va Rank-One Projectons arxv:1310.5791v [math.st] 17 Aug 014 T. Tony Ca and Anru Zhang Department of Statstcs The Wharton School Unversty of Pennsylvana Abstract Estmaton of low-rank

More information

Perron Vectors of an Irreducible Nonnegative Interval Matrix

Perron Vectors of an Irreducible Nonnegative Interval Matrix Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

Another converse of Jensen s inequality

Another converse of Jensen s inequality Another converse of Jensen s nequalty Slavko Smc Abstract. We gve the best possble global bounds for a form of dscrete Jensen s nequalty. By some examples ts frutfulness s shown. 1. Introducton Throughout

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

The lower and upper bounds on Perron root of nonnegative irreducible matrices

The lower and upper bounds on Perron root of nonnegative irreducible matrices Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College

More information

HANSON-WRIGHT INEQUALITY AND SUB-GAUSSIAN CONCENTRATION

HANSON-WRIGHT INEQUALITY AND SUB-GAUSSIAN CONCENTRATION HANSON-WRIGHT INEQUALITY AND SUB-GAUSSIAN CONCENTRATION MARK RUDELSON AND ROMAN VERSHYNIN Abstract. In ths expostory note, we gve a modern proof of Hanson-Wrght nequalty for quadratc forms n sub-gaussan

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

5 The Rational Canonical Form

5 The Rational Canonical Form 5 The Ratonal Canoncal Form Here p s a monc rreducble factor of the mnmum polynomal m T and s not necessarly of degree one Let F p denote the feld constructed earler n the course, consstng of all matrces

More information

Appendix B. Criterion of Riemann-Stieltjes Integrability

Appendix B. Criterion of Riemann-Stieltjes Integrability Appendx B. Crteron of Remann-Steltes Integrablty Ths note s complementary to [R, Ch. 6] and [T, Sec. 3.5]. The man result of ths note s Theorem B.3, whch provdes the necessary and suffcent condtons for

More information

Supplement to Clustering with Statistical Error Control

Supplement to Clustering with Statistical Error Control Supplement to Clusterng wth Statstcal Error Control Mchael Vogt Unversty of Bonn Matthas Schmd Unversty of Bonn In ths supplement, we provde the proofs that are omtted n the paper. In partcular, we derve

More information

Lecture 4: September 12

Lecture 4: September 12 36-755: Advanced Statstcal Theory Fall 016 Lecture 4: September 1 Lecturer: Alessandro Rnaldo Scrbe: Xao Hu Ta Note: LaTeX template courtesy of UC Berkeley EECS dept. Dsclamer: These notes have not been

More information

Edge Isoperimetric Inequalities

Edge Isoperimetric Inequalities November 7, 2005 Ross M. Rchardson Edge Isopermetrc Inequaltes 1 Four Questons Recall that n the last lecture we looked at the problem of sopermetrc nequaltes n the hypercube, Q n. Our noton of boundary

More information

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017 U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS

BOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all

More information

Complete subgraphs in multipartite graphs

Complete subgraphs in multipartite graphs Complete subgraphs n multpartte graphs FLORIAN PFENDER Unverstät Rostock, Insttut für Mathematk D-18057 Rostock, Germany Floran.Pfender@un-rostock.de Abstract Turán s Theorem states that every graph G

More information

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES

TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES TAIL BOUNDS FOR SUMS OF GEOMETRIC AND EXPONENTIAL VARIABLES SVANTE JANSON Abstract. We gve explct bounds for the tal probabltes for sums of ndependent geometrc or exponental varables, possbly wth dfferent

More information

Large Sample Properties of Matching Estimators for Average Treatment Effects by Alberto Abadie & Guido Imbens

Large Sample Properties of Matching Estimators for Average Treatment Effects by Alberto Abadie & Guido Imbens Addtonal Proofs for: Large Sample Propertes of atchng stmators for Average Treatment ffects by Alberto Abade & Gudo Imbens Remnder of Proof of Lemma : To get the result for U m U m, notce that U m U m

More information

REDUCTION MODULO p. We will prove the reduction modulo p theorem in the general form as given by exercise 4.12, p. 143, of [1].

REDUCTION MODULO p. We will prove the reduction modulo p theorem in the general form as given by exercise 4.12, p. 143, of [1]. REDUCTION MODULO p. IAN KIMING We wll prove the reducton modulo p theorem n the general form as gven by exercse 4.12, p. 143, of [1]. We consder an ellptc curve E defned over Q and gven by a Weerstraß

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

More information

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014)

10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014) 0-80: Advanced Optmzaton and Randomzed Methods Lecture : Convex functons (Jan 5, 04) Lecturer: Suvrt Sra Addr: Carnege Mellon Unversty, Sprng 04 Scrbes: Avnava Dubey, Ahmed Hefny Dsclamer: These notes

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

e - c o m p a n i o n

e - c o m p a n i o n OPERATIONS RESEARCH http://dxdoorg/0287/opre007ec e - c o m p a n o n ONLY AVAILABLE IN ELECTRONIC FORM 202 INFORMS Electronc Companon Generalzed Quantty Competton for Multple Products and Loss of Effcency

More information

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2 Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

Norms, Condition Numbers, Eigenvalues and Eigenvectors

Norms, Condition Numbers, Eigenvalues and Eigenvectors Norms, Condton Numbers, Egenvalues and Egenvectors 1 Norms A norm s a measure of the sze of a matrx or a vector For vectors the common norms are: N a 2 = ( x 2 1/2 the Eucldean Norm (1a b 1 = =1 N x (1b

More information

The Expectation-Maximization Algorithm

The Expectation-Maximization Algorithm The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.

More information

On the Connectedness of the Solution Set for the Weak Vector Variational Inequality 1

On the Connectedness of the Solution Set for the Weak Vector Variational Inequality 1 Journal of Mathematcal Analyss and Alcatons 260, 15 2001 do:10.1006jmaa.2000.7389, avalable onlne at htt:.dealbrary.com on On the Connectedness of the Soluton Set for the Weak Vector Varatonal Inequalty

More information

MATH 5707 HOMEWORK 4 SOLUTIONS 2. 2 i 2p i E(X i ) + E(Xi 2 ) ä i=1. i=1

MATH 5707 HOMEWORK 4 SOLUTIONS 2. 2 i 2p i E(X i ) + E(Xi 2 ) ä i=1. i=1 MATH 5707 HOMEWORK 4 SOLUTIONS CİHAN BAHRAN 1. Let v 1,..., v n R m, all lengths v are not larger than 1. Let p 1,..., p n [0, 1] be arbtrary and set w = p 1 v 1 + + p n v n. Then there exst ε 1,..., ε

More information

Algorithms for factoring

Algorithms for factoring CSA E0 235: Crytograhy Arl 9,2015 Instructor: Arta Patra Algorthms for factorng Submtted by: Jay Oza, Nranjan Sngh Introducton Factorsaton of large ntegers has been a wdely studed toc manly because of

More information

Some congruences related to harmonic numbers and the terms of the second order sequences

Some congruences related to harmonic numbers and the terms of the second order sequences Mathematca Moravca Vol. 0: 06, 3 37 Some congruences related to harmonc numbers the terms of the second order sequences Neşe Ömür Sbel Koaral Abstract. In ths aer, wth hels of some combnatoral denttes,

More information

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples

More information

ALGEBRA MID-TERM. 1 Suppose I is a principal ideal of the integral domain R. Prove that the R-module I R I has no non-zero torsion elements.

ALGEBRA MID-TERM. 1 Suppose I is a principal ideal of the integral domain R. Prove that the R-module I R I has no non-zero torsion elements. ALGEBRA MID-TERM CLAY SHONKWILER 1 Suppose I s a prncpal deal of the ntegral doman R. Prove that the R-module I R I has no non-zero torson elements. Proof. Note, frst, that f I R I has no non-zero torson

More information

A NOTE ON THE DISCRETE FOURIER RESTRICTION PROBLEM

A NOTE ON THE DISCRETE FOURIER RESTRICTION PROBLEM A NOTE ON THE DISCRETE FOURIER RESTRICTION PROBLEM XUDONG LAI AND YONG DING arxv:171001481v1 [mathap] 4 Oct 017 Abstract In ths aer we establsh a general dscrete Fourer restrcton theorem As an alcaton

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

THE WEIGHTED WEAK TYPE INEQUALITY FOR THE STRONG MAXIMAL FUNCTION

THE WEIGHTED WEAK TYPE INEQUALITY FOR THE STRONG MAXIMAL FUNCTION THE WEIGHTED WEAK TYPE INEQUALITY FO THE STONG MAXIMAL FUNCTION THEMIS MITSIS Abstract. We prove the natural Fefferman-Sten weak type nequalty for the strong maxmal functon n the plane, under the assumpton

More information

MATH 281A: Homework #6

MATH 281A: Homework #6 MATH 28A: Homework #6 Jongha Ryu Due date: November 8, 206 Problem. (Problem 2..2. Soluton. If X,..., X n Bern(p, then T = X s a complete suffcent statstc. Our target s g(p = p, and the nave guess suggested

More information

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede Fall 0 Analyss of Expermental easurements B. Esensten/rev. S. Errede We now reformulate the lnear Least Squares ethod n more general terms, sutable for (eventually extendng to the non-lnear case, and also

More information

DISCRIMINANTS AND RAMIFIED PRIMES. 1. Introduction A prime number p is said to be ramified in a number field K if the prime ideal factorization

DISCRIMINANTS AND RAMIFIED PRIMES. 1. Introduction A prime number p is said to be ramified in a number field K if the prime ideal factorization DISCRIMINANTS AND RAMIFIED PRIMES KEITH CONRAD 1. Introducton A prme number p s sad to be ramfed n a number feld K f the prme deal factorzaton (1.1) (p) = po K = p e 1 1 peg g has some e greater than 1.

More information

Slobodan Lakić. Communicated by R. Van Keer

Slobodan Lakić. Communicated by R. Van Keer Serdca Math. J. 21 (1995), 335-344 AN ITERATIVE METHOD FOR THE MATRIX PRINCIPAL n-th ROOT Slobodan Lakć Councated by R. Van Keer In ths paper we gve an teratve ethod to copute the prncpal n-th root and

More information

Affine transformations and convexity

Affine transformations and convexity Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/

More information

SL n (F ) Equals its Own Derived Group

SL n (F ) Equals its Own Derived Group Internatonal Journal of Algebra, Vol. 2, 2008, no. 12, 585-594 SL n (F ) Equals ts Own Derved Group Jorge Macel BMCC-The Cty Unversty of New York, CUNY 199 Chambers street, New York, NY 10007, USA macel@cms.nyu.edu

More information

Homework Notes Week 7

Homework Notes Week 7 Homework Notes Week 7 Math 4 Sprng 4 #4 (a Complete the proof n example 5 that s an nner product (the Frobenus nner product on M n n (F In the example propertes (a and (d have already been verfed so we

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1. 7636S ADVANCED QUANTUM MECHANICS Soluton Set 1 Sprng 013 1 Warm-up Show that the egenvalues of a Hermtan operator  are real and that the egenkets correspondng to dfferent egenvalues are orthogonal (b)

More information

Notes on Frequency Estimation in Data Streams

Notes on Frequency Estimation in Data Streams Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to

More information

Some basic inequalities. Definition. Let V be a vector space over the complex numbers. An inner product is given by a function, V V C

Some basic inequalities. Definition. Let V be a vector space over the complex numbers. An inner product is given by a function, V V C Some basc nequaltes Defnton. Let V be a vector space over the complex numbers. An nner product s gven by a functon, V V C (x, y) x, y satsfyng the followng propertes (for all x V, y V and c C) (1) x +

More information

A note on almost sure behavior of randomly weighted sums of φ-mixing random variables with φ-mixing weights

A note on almost sure behavior of randomly weighted sums of φ-mixing random variables with φ-mixing weights ACTA ET COMMENTATIONES UNIVERSITATIS TARTUENSIS DE MATHEMATICA Volume 7, Number 2, December 203 Avalable onlne at http://acutm.math.ut.ee A note on almost sure behavor of randomly weghted sums of φ-mxng

More information

STEINHAUS PROPERTY IN BANACH LATTICES

STEINHAUS PROPERTY IN BANACH LATTICES DEPARTMENT OF MATHEMATICS TECHNICAL REPORT STEINHAUS PROPERTY IN BANACH LATTICES DAMIAN KUBIAK AND DAVID TIDWELL SPRING 2015 No. 2015-1 TENNESSEE TECHNOLOGICAL UNIVERSITY Cookevlle, TN 38505 STEINHAUS

More information

Mixture of Gaussians Expectation Maximization (EM) Part 2

Mixture of Gaussians Expectation Maximization (EM) Part 2 Mture of Gaussans Eectaton Mamaton EM Part 2 Most of the sldes are due to Chrstoher Bsho BCS Summer School Eeter 2003. The rest of the sldes are based on lecture notes by A. Ng Lmtatons of K-means Hard

More information

Deriving the X-Z Identity from Auxiliary Space Method

Deriving the X-Z Identity from Auxiliary Space Method Dervng the X-Z Identty from Auxlary Space Method Long Chen Department of Mathematcs, Unversty of Calforna at Irvne, Irvne, CA 92697 chenlong@math.uc.edu 1 Iteratve Methods In ths paper we dscuss teratve

More information

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP C O L L O Q U I U M M A T H E M A T I C U M VOL. 80 1999 NO. 1 FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP BY FLORIAN K A I N R A T H (GRAZ) Abstract. Let H be a Krull monod wth nfnte class

More information

Chapter 7 Generalized and Weighted Least Squares Estimation. In this method, the deviation between the observed and expected values of

Chapter 7 Generalized and Weighted Least Squares Estimation. In this method, the deviation between the observed and expected values of Chapter 7 Generalzed and Weghted Least Squares Estmaton The usual lnear regresson model assumes that all the random error components are dentcally and ndependently dstrbuted wth constant varance. When

More information