Supplement to the papers on the polyserial correlation coefficients and discrepancy functions for general covariance structures

Size: px
Start display at page:

Download "Supplement to the papers on the polyserial correlation coefficients and discrepancy functions for general covariance structures"

Transcription

1 conomic Review (Otaru University of Commerce Vo.6 No March 0. Suppement to the papers on the poyseria correation coefficients and discrepancy functions for genera covariance structures Haruhio Ogasawara Parts A and B of this note are to suppement Ogasawara (0 00 respectivey. Part A 0. Derivation of the inverse expansion Let L / θ θ '. Then from (. the inverse expansion of ˆθ is 0 0 ˆ ˆ θθ L L θθ L θˆ θ 4 0 ( 0 ( 0 θ0 θ0( θ0 ' 6 θ0 ( θ0 ' O ( p N L L L θ0 θ0( θ0 ' θ0 L L L L θ0( θ0 ' θ0 θ0( θ0 ' θ 0 4 O p N L θ ( θ ' θ (. Let ( L Λ I where I is the information matrix per observation and et L Λ M where M O N ( / p. Then we obtain The ove resuts give / L Λ Λ M Λ Λ M Λ M Λ O N. p ( 5

2 ˆ θ θ0 Λ θ 0 / Op ( N Λ MΛ Λ Λ θ0 θ0( θ0 ' θ0 Λ MΛ MΛ Λ MΛ Λ θ0 θ0( θ0 ' θ0 Λ Λ MΛ Λ θ0( θ0 ' θ0 θ0 Λ Λ θ0( θ0 ' θ0( θ0 ' θ0 i 4 O p N Λ θ0 θ0 θ0 / O ( N ( i ( i 0 p p O ( N Λ Λ Λ Λ θ0( θ0 ' θ0 θ0( θ0 ' θ 0 Λ ( 6 ( ' Λ O ( N p where [ ] O p ( indicates that the sum of the terms in bracets is of order Op( for carity.. xpectations of the og ieihood derivatives Let Λ Bdiag( Λx Λ and I Bdiag( Ix I where Λ I Bdiag( denotes the boc diagona matrix with the diagona bocs being the matrices in parentheses and I is the information matrix per observation.. Information matrix ( μ and Σ 6

3 * x μ * x N Σ xi μ i ( * x N Σ μμ' * x μ N [ { Σ ( xi μ} a{ Σ ( xi μ} b] i N [( Σ a{ Σ ( xi μ} b ( Σ b{ Σ ( xi μ} a] i * N x ( ( cd cd 4 i where 4 7 ac db ad Σ x μ Σ x μ ac { ( i } d{ ( i } b Σ denotes the sum of terms with simiar patterns. Then we obtain N N N μμ' μμ' N( I 0 * * x x ( Ix μμ' ( Λx μμ' Σ x μ' * x N N( Ix ( ( ( N Λ cd x cd cd cd 4 that is ( Λx μμ' Σ ( Λx μ 0 ac db ad cb ( Λx ( ( ( cd cd 4 ( r a b ; r c d where e.g. ( uv ' cb ac db ad cb ( denotes the submatrix of the matrix in parentheses corresponding to the product pairs wi be used hereafter in simiar situations. u v ' and the ranges of (a b (c d and simiar

4 ( β( (... ' and ξ ( ( (... r ' * N i i ( i i i where i ( ( ( i ( d i i i i i ( exp i i ξ ' x ( 0 i i ξ * i N i x i i ( i i * N ( i i i ( i i i i i where i i * ξ ( ( i i ( i i i i ( i i i N xi ( i i i ( i i i i i ( ( i i ( i i i i i * ξξ' i i N xixi ' ( ( i i i i i i i 8

5 * N( I ( N Λ X Z X N X ( ( N where ( ( X {} {} p( d p X x x ( Σ ( x exp ( ' ( r/ / x μ Σ x μ * N( I ξ N( Λ ξ X Z X ξ N X x * N( I ξξ' N( Λ ξξ' XZ X N X xx' ( a a ξξ' a a that is ( ( ( Λ = X ( Λ ξξ' = X xx' ( a a. a a 9 ( Λ ξ ( ( X x

6 . xpectations of the products of the og ieihood derivatives * Reca L Λ M and N. Let M Bdiag( Mx M. ( μ and Σ ' ( ( N (- 0 N v( M' θ0 ' θ 0 ' (for The nonero resuts of n p( n p( N = x x X ( θx A( θx B ( θx C ( θx A( θx B ( θ x where C ( A denotes the A-th eement of the vector in parentheses are n p( x n p( x X Σ a Σ b Σ b Σ a μ μ' {( ( ( ( } 8 n p( x n p( x ac de fb X ( ( cd ( ef { }. cd ef 8 (- ( N (for ( 0 The nonero resuts of X p( x θ x n p( x n p( x n p( x X a b {( Σ ( Σ μ μ' are ( Σ b( Σ a} 0

7 n p( x n p( x n p( x X ( ( cd ( ef cd ef 8 cd ef cd ef cd ef ef cd ( 5 cd ef 8 ( ( ( ac de fb cd ef. 8 (- N ( ( ( N 0 0 (for ( θ ( θ ( θ ( θ x A x B x C x D n p( x n p( x n p( x n p( x = X X O( N ( θx A ( θx B ( θx C ( θx D N ( θ ( θ ( θ ( θ ( θ ( θ ( θ n p( x n p( x n p( x n p( x = X X O( N. ( θx A( θx B ( θx C ( θx D ( θx (-4 x A x B x A x B x C x D x ( ' ( ( N 0 0 (for Define { X ( } as { X ( X } N ( ( ( θx A( θx B ( θx C ( θx D ( θx ( θx F we derive the nonero resuts of. In

8 n p( x n p( x X X ( X ( ( θx A( θx B ( θx C ( θx D [{( L ( Λ }{( L ( Λ }] X x AB x AB x CD x CD [{( L ( Λ }( L ] where X x AB x AB x CD n p( x Lx ( L Bdiag( Lx L θ θ ' x x n p( x n p( x X X ( X ( c de f 4 n p( x n p( x ac bd ef = X ( ( de c de f 4 n p( x n p( x X X ( X ( cd ef gh ( ( cd ( ef ( gh 6 4 ac db ad cb ac X { ( } d{ ( } b Σ x μ Σ x μ 4 eg hf eh gf eg { Σ ( x μ} h{ Σ ( x μ} f 6 ( ( ( ( ac eg ( dh bf df bh cd ef gh. 6 (-5 For ( ' ( ( N 0 0 (for ( N ( 0 see (-. The nonero resuts of

9 n p( x n p( x X {( Lx AB ( Λx AB} ( θx C ( θx D n p( x n p( x ( L ( Λ ( Λ ( θx C ( θx D are X x AB x AB x CD 4 n p( x n p( x n p( x X = ( ( de c de f 4 n p( x n p( x n p( x X X ( cd e f 4 ac de bf df be = ( ( cd ( 4 n p( x n p( x n p( x X X ( cd ef gh 4 8 ac de fg hb = ( ( cd ( ef ( gh. 6 ca bd ef (-6 For ( ' ( ( N 0 0 (for ( N ( 4 0 ( θ x A( θ x B ( θ x C ( θ x D ( θ x and N ( N ( ( ( θx A( θx B ( θx C ( θx D ( θx ( θ x F see (- and (-4. The remaining resuts are

10 n p( x n p( x X X ( ( θx A( θx B( θx C ( θx D n p( x n p( x X ( θx A( θx B( θx C ( θx D in ( θ ( θ ( θ ( θ ( θ ( θ N ( x A x B x C x D x x F First we obtain the nonero n p( x μμ' n p( x ( θ ( θ ( θ x A x B x C {( Σ a( Σ b ( Σ b( Σ a } 4 n p( x ( ( ac [{ ( } ( cd Σ x μ b Σ d cd μ 4 and n p( x ( ( cd ( ef 8 cd ef as { Σ ( x μ} d ( Σ b] n p( x n p( x X ( θx A( θx B( θx C ( θx as D n p( x n p( x X cdμ μ' 4 ( ( ac {( ( ( ( cd Σ d Σ b Σ b Σ d } 4 which yied the nonero 8 4 ac de fb ac de { Σ ( x μ} f { Σ ( x μ} b 4

11 n p( x n p( x X ( ( cd ( ef ( gh cd ef gh 6 4 ac de fg hb fh gb (. (-7 ( 4 ( 0 0 ( N N {( N } (for 4 n ( x n ( x n ( x n ( x N N p p p p X X N ( θx A ( θx B ( θx C ( θx D n p( x n p( x n p( x n p( x X x AB x CD N( Λ ( Λ ( θx A ( θx B ( θx C ( θx D n p( x n p( x n p( x n p( x = X x AB x CD ( Λ ( Λ ( θx A ( θx B ( θx C ( θx D (the fourth mutivariate cumuant of n p( x / ( θ s. X The nonero expectations required for the fourth cumuant are n p( x n p( x n p( x n p( x cd e f 8 ( ( ef ( ac bd ad bc ae cf bd cd 4 n p( x n p( x n p( x n p( x X cd ef gh ( ( cd ( ef ( gh 6 x A 5

12 6 cd ef gh cd ( ef gh eg fh eh fg cd ef gh cd ef gh 64 ( ( cd ( ef ( gh 6 (note that X ac de fg hb ac bd eg fh ( cd ef gh ( c d ( g h n p( x n p( x n p( x n p( x a b c d nonero but the corresponding cumuant is ero. cd is (-8 ( ( ( N 0 0 (for 4 The expectations are expressed as 0 ( where each ( of two expectations which was shown in (- (- and (-5. is the product (-9 ( ( ( N 0 0 (for 4 The resuts are expressed as 5 ( where each ( expectations which was shown in (- and (-4. (-0 ( ( ( N 0 0 (for 4 The resuts are expressed as expectations which was shown in (-6. 5 ( where each ( is the product of three is the product of three (- ( θx A ( θx B ( θ x (for C and 6

13 The resuts are given by the Bartett identity as foows: n p( x n p( x n p( x ( θx A( θx B( θx C ( θx A( θx B ( θx C (- n p( x n p( x n p( x. ( θx A ( θx B ( θx C 4 ( θx A ( θx B ( θx C ( θ x (for D As in (- we obtain and n p( x n p( x n p( x ( θx A( θx B( θx C ( θx D ( θx A( θx B( θx C ( θx D n p( x n p( x n p( x ( θ ( θ ( θ ( θ ( θ ( θ 6 x A x B x C x D x A x B n p( x n p( x n p( x n p( x n p( x n p( x. ( θx C ( θx D ( θx A ( θx B ( θx C ( θx D ( β and ξ ' ( ( N (- 0 N v( M' θ0 ' θ 0 ' (for The resuts are 7

14 n Pr( Z x n Pr( Z x m ( ( ( ( m ( m m m m ( X ( m ( m ( ( m( m n Pr( Z x n Pr( Z x ξ ( ( ( ( a ( x a a a a ( X x ( ( ( 8

15 n Pr( Z x n Pr( Z x ξ ( ( x ( ( ( X x ( ( a ( a a a a ( ( ( ( n Pr( Z x n Pr( Z x ξ ξ' ( ( xx' ( ( X xx' ( ( ( ( 9

16 n Pr( Z x n Pr( Z x ξ ξ ' xx' ( mm m m ( m m m m m ( m X xx' ( ( n Pr( Z x n Pr( Z x m ( ( x x x ( ( m aa a a a a a a a a b ( b b b b X x x xm ( aa a a ( a a ( a a a a a (- ( N (for ( 0 n Pr( Z x n Pr( Z x n Pr( Z x m ( ( m ( m m m m. 0

17 X m ( m ( m ( ( m ( m n Pr( Z x n Pr( Z x n Pr( Z x ξ ( ( a x ( a a a a X x ( ( ( ( ( n Pr( Z x n Pr( Z x n Pr( Z x ' ξ ξ ( a xx' ( a a a a X xx' ( ( n Pr( Z x n Pr( Z x n Pr( Z x m a x xxm ( a a X x xxm ( a a. a a a a (- ( ( ( N 0 0 (for

18 N ( θ ( θ ( θ ( θ n Pr( Z x n Pr( Z x n Pr( Z x n Pr( Z x = ( θ A ( θ B ( θ C ( θ D O N ( (see the expression of N A B C D Λ in. ( ( ( θ ( θ ( θ ( θ ( θ n Pr( Z x n Pr( Z x n Pr( Z x n Pr( Z x = ( θ A( θ B ( θ C ( θ D ( θ O N ( A B C D (see (- and the expression of Λ in. (. (-4 For ( ' ( ( N 0 0 (for ( ( 0 0 ( N see (-. Reca L n Pr( Z x θ θ ' remaining resuts are [{( L ( Λ }{( L ( Λ }] in AB AB CD CD then the n Pr( Z x n Pr( Z x [{( L AB ( Λ AB}{( L CD ( Λ CD}]. ( θ ( θ F We obtain n Pr( Z x n Pr( Z x ( Λ ( Λ m n m n ( (

19 ( ( m ( m m ( m mn m m mn m m m m m m( n ( m n ( ( mn Λ Λ m m m n X mn 4 m( n ( mn ( m( n ( mn ( ( ( m n ( m( n ( m( n ( ( ( mn ( mn Λ Λ n Pr( Z x n Pr( Z x ( Λ ( Λ ξ m ξ x ( ( m m n ( ( m ( m m mm m( m m m m m

20 ( m ( ( ( Λ Λ m m m m ξ m m m X x m ( m m m m m mm ( m m m m m m m ( m ( m m m m m mm ( ( m ( m m m m m m m ( m ( m m m ( m m mm ( ( m ( m m m ( Λ ( Λ m ( m m ξ n Pr( Z x n Pr( Z x ( Λ ( Λ ' ξ ξ xx' ξξ' ( ( ( ( ( mm m m ( m m m ( Λ ( Λ ξ ξ' m m m 4

21 X xx' ( ( ( ( ( ( ( ( ( Λ ( Λ ξ ξ' n Pr( Z x n Pr( Z x ( Λ ( ' ξ Λ ξ ξ ( ( xx' ( ( ( ( ( ( ( Λ ( Λ ξ ξ' X xx' ( ( ( ( ( ( ( Λ ( Λ ξ ξ' ξ' 5

22 n Pr( Z x n Pr( Z x ( Λ ( Λ ' ξ ξ ( ( xxx' ( ( ( mm m m ( m m m ( Λ ( Λ ξ ξ' m m m ξ ξ' X xxx' ( ( ( ( ( ( ( Λ ( ' Λ ξξ n Pr( Z x n Pr( Z x ( Λ ( Λ m n x x x x ( ( ( Λ ( Λ m n aa a a a a a a a a m n X x x xmxn aa a a ( a a a a a ( Λ ( Λ m n m n (-5 For ( ' ( ( N 0 0 (for ( N ( 0 see (-. The remaining resuts are 6

23 n Pr( Z x n Pr( Z x {( L AB ( Λ AB} ( θ C ( θ D n Pr( Z x n Pr( Z x X ( L AB ( Λ AB( Λ CD. ( θ C ( θ D We obtain the first term on the right-hand side of the ove equation as foows: n Pr( Z x n Pr( Z x n Pr( Z x m n ( ( ( ( m ( m n ( n m n m m n n X mn m( n ( mn ( ( m( n ( mn ( m( n ( ( mn ( m( n ( ( mn ( ( m( n 7

24 n Pr( Z x n Pr( Z x n Pr( Z x m ξ ( ( x ( ( m ( m a m m m a a X x m ( ( a a ( ( m ( ( ( m ( ( m ( ( ( m ( ( n Pr( Z x n Pr( Z x n Pr( Z x ξ ξ' ( ( xx' ( ( a ( a a a a 8

25 X xx' ( ( ( ( ( n Pr( Z x n Pr( Z x n Pr( Z x ξ m ( ( x ( ( ( m ( m m m m ( X x m ( m ( ( m ( ( ( m ( ( ( m ( 9

26 n Pr( Z x n Pr( Z x n Pr( Z x ξ ξ' ( ( xx' ( ( ( ( a a a a a X xx' ( ( ( ( ( ( ( ( ( ( n Pr( Z x n Pr( Z x n Pr( Z x ξ m ( ( x x xm ( ( a a a ( a a 40

27 X x x xm ( ( n Pr( Z x n Pr( Z x n Pr( Z x ' ξ ξ ( ( xx' ( ( aa a a a a a a a a ( ( ' ( ( X xx ( ( ( ( ( ( n Pr( Z x n Pr( Z x n Pr( Z x ξ ξ ' xx' x ( ( aa a a a a a a a a ( ( b b b b b 4

28 X xx' x ( ( ( ( ( ( n Pr( Z x n Pr( Z x n Pr( Z x m n x x x x ( ( m n aa a a a a a a a a b b b ( b b X x x xmxn ( aa a a ( a a ( a a a a a. (-6 For ( ' ( ( N 0 0 (for ( N ( 4 0 ( θ A( θ B ( θ C ( θ D ( θ N ( and N ( ( ( θ ( θ ( θ ( θ ( θ ( θ see (- an (.4. The remaining resuts are A B C D F Pr( Z x Pr( Z x N ( ( θ A( θ B( θ C ( θ D Pr( Z x Pr( Z x ( θ A ( θ B ( θ C ( θ D 4

29 in N ( ( θ A( θ B( θ C ( θ D ( θ ( θ F. Pr( Z x First we derive ( θ ( θ ( θ as foows: A B C n Pr( Z x ( m ( m ( ( ( m ( m ( ( m n Pr( Z x ( x ( ξ ( ( ( ( ( ( ( ( ( ( ( ( ( ( 4

30 n Pr( Z x ( xx' ( ξ ξ' ( ( ( { ( } ( { ( } ( ( ( n Pr( Z x {( ( x x xm a a a a } m a a Then the expectations are ( aa a a( a a ( a a a. a a ( n Pr( Z x n Pr( Z x ( m m n ( ( ( m ( m n ( n ( ( m n n n 44

31 = ( X mn 4 ( mn ( m ( m( n ( ( mn ( ( m( n m n ( ( ( ( ( m n ( ( n Pr( Z x n Pr( Z x ( x m m ξ ( ( 45

32 ( m ( m a ( ( m ( a a a a = X x m ( ( ( ( ( m ( m ( ( m ( ( n Pr( Z x n Pr( Z x ( x ξ m ( ( ( ( ( ( ( m ( m ( ( ( m ( m m 46

33 = x ( X m ( ( ( ( ( m ( ( ( ( ( m ( ( ' ( n Pr( Z x n Pr( Z x ( xx ξ ξ' ( ( ( ( ( ( ( a ( ( ( ( a a ( a a = X xx' ( ( ( ( ( ( 47

34 ( ( ( ( ( ( ( n Pr( Z x n Pr( Z x ( xx' ( ' ξξ ( ( ( { ( } ( ( ( { ( } ( ( = xx' ( X ( ( { ( } ( { ( } ( 48

35 ( ( ( { ( } ( ( ( ( { ( } ( ' ( n Pr( Z x n Pr( Z x ( xx x ξ ξ ' ( ( ( { ( } ( ( ( { ( } ( a ( a a a a = X xx' x ( ( { ( } ( ( 49

36 ( ( { ( } ( ( n Pr( Z x n Pr( Z x m n x x xm {( a a ( a a } a a n ( n ( aa a a( a a ( a a a n a a n n X x x xm {( n n ( n n } n n ( ( ( n n n n n n n n n n n n {( n n ( n n } n n ( n n nn ( n n n ( n n n n n n Pr( Z x n Pr( Z x m n 50

37 x x xm xn {( a a ( a a } a a b ( aa a a( a a ( a a a ( b b a a b b X x x xm xn {( a a ( a a } a a ( aa a a( a a ( a a ( a a. a a (-7 ( 4 ( 0 0 ( N N {( N } (for 4 n Pr( n Pr( x ( θ A ( θ B N N Z x Z N n Pr( Z x n Pr( Z x n Pr( Z x n Pr( Z x ( θ C ( θ D ( θ A ( θ B n Pr( Z x n Pr( Z x N( Λ ( Λ ( θ C ( θ D n Pr( Z x n Pr( Z x n Pr( Z x n Pr( Z x = ( θ A ( θ B ( θ C ( θ D ( Λ ( Λ AB CD (the fourth mutivariate cumuant of n Pr( Z x / ( θ s. The expectations required in the ove cumuants are n Pr( Z x n Pr( Z x n Pr( Z x n Pr( Z x m n A AB CD 5

38 ( ( m ( m n ( n m n m m n n 4 4 X mn m( n ( mn 6 ( ( ( m n ( m( n ( mn n Pr( Z x n Pr( Z x n Pr( Z x n Pr( Z x m ξ ( ( x m ( m a m ( a a m m a a X x m ( ( ( m ( ( ( m ( ( m n Pr( Z x n Pr( Z x n Pr( Z x n Pr( Z x ξ ξ' ( ( a xx' ( a a a a 5

39 xx' ( ( X ( ( ( n Pr( Z x n Pr( Z x n Pr( Z x n Pr( Z x m n ( a x xm xn ( a a a a X x xm xn ( ( n Pr( Z x n Pr( Z x n Pr( Z x n Pr( Z x m n 4 a x x xm xn ( a a a a 4 X x x xm xn ( a a a a (-8 ( ( ( N 0 0 (for 4 The expectations are expressed as of two expectations (see (- (- and (-5. 0 ( where each ( is the product (-9 ( ( ( N 0 0 (for 4 5

40 The resuts are expressed as 5 ( where each ( is the product of three expectations (see (- and (-4. ( N (for 4 ( ( ( The resuts are simiar to those in (-9 (see (-6. (- (- ( θ A ( θ B ( θ (for C and The resuts are given by the Bartett identity (see (-. 4 ( θ A ( θ B ( θ C ( θ (for D and 4 The resuts are given by the Bartett identity (see (-. ( μ Σ β and ξ Some of the expectations of the og ieihood derivatives with respect to θx and θ e.g. N ( θx A ( θx B ( θ C ( θ are nonero D though the corresponding fourth cumuants are ero. The resuts are given as in (- (-4 (- and (-4.. The nonero partia derivatives of η with respect to θ. / We define S S( ξ Σ ( ξ' Σξ ( R then τ ( β ξ' μ S / and ρ (Diag Σ Σξ S.. First derivatives τ τ S S τ S τ τ S e S S τ S 54

41 where e is the vector of an appropriate dimension whose -th eement is and the remaining ones are 0. ρ ρ where Σξ S (Diag Σ ( ξ S ρs / / aa aa ba S S / (Diag Σ ( Σ S ρ S S S a b S ( Σξ S and is the matrix of an appropriate sie whose (a bth eement is and the remaining ones are 0.. Second derivatives τ S S τ τ S S S τ S cd ( cd cd cd τ S τ S es S S τ τ S τ S S S S τ S τ τ S τ τ S S S S τ S ( ρ 5/ / cd aa aaσξ S aa cd cd 4 ( cd 4 ( cd cd cd ( ρ S S ( acad adac ξ S S ρs 55

42 ρ ( Σ S (Diag Σ ( e e S / / aa aa a b b a ρ S ρ S S S S ρs ρ ρ S S S ρs ( where S cd ( ( cd cd S 4 S S S S S S S. ( ba S S. Third derivatives τ S S τ S τ S cd ( cd cd cd τ τ S τ S S S cd ef ( cd ef cd ef cd ef τ S S cd ef τ S S τ S τ S cd ( cd cd cd τ τ S τ S τ S S S S 56

43 τ τ S τ S S S S cd ( cd cd cd τ S S cd τ S S S S S τ S τ S cd cd cd τ τ S τ S τ S S S S τ S S τ S τ S m ( m m τ τ S τ S S S S ( τ S S τ S S S S S τ S τ S τ S S τ S τ S m ( m m τ τ S τ S S S m ( m m m τ S S m ρ 5 7/ 5/ cdef aa aaσξ S cd aa ef cd ef 8 ( cd ef 8 ρ S ρ S ( aeaf af ae ξ S S S cd ef ef cd ( 57

44 ρ S S cd ef ρ S 5/ / cd aa aa ( Σ S Σξ S aa cd 4 ( cd 4 S ( cd ( ac d a d c eas ( acad adac ξ S ρ S ρ S S ρ S S S S ( cd cd cd cd ρ S S S S S ρs ρs cd cd cd ρ ρ S ρ S S S S ( ρ S ρ S S S S S ρs ρs S ρ ρ S ρ S S S S ρs m ( m m m m where S cd ef ( ( ( 5 cd ef cde f S 8 S ( ( cd ( ba cds cd ( cd 4 cd S 4 4 S 58

45 S S ( ab b a S S S S S S S S S ( S S S S. m ( m m. Computation by Gaussian quadrature. Univariate case Stroud and Sechrest (966 Te Five pp.7-5 gave the foowing vaues of A and x : i i n exp( x f ( x dx Ai f ( xi i where n=(64(496(86. Let y x. Then ( y f ( y dy exp( y / f ( y dy n exp( x f ( x dx Ai f ( xi. The ove resut corresponds to Boc and Lieberman (970 quation (5. Boc and Lieberman (970 p.8 used n=64 where ony 40 vaues of 4 were empoyed since A i s for the remaining x i s are. Boc and Aitin (98 Te reported the resuts using n=0 and.. Bivariate case be i A i 0 Suppose that U ( X Y' N( μ Σ. Let the density of U at u ( x y' Σ ( x y exp ( ' ( / u μ Σ u μ x i 59

46 where Σ x yx xy y ( x y μ Σ. Define ( x0 ( x x that x x x and μ ( x y' ( U u ( x y and ( x ( x x exp x x. Then using the transformations ( x y g( x y dxdy ( x x x y y {( xy / x }( x x y y x y x / [ y {( xy / x}] y x [ y ( xy / x ( x y {( xy / x}] ( ( ( g( d d where i j x y x x x x y x y x y x x y x (note it foows y x g x y dxdy n ( A g( x d n x j x x x y x y x j x j n A A g( x x * i j x x i y x y x j. * y x ( / y xy x x y x i x x i 60

47 Part B. The partia derivatives of ULS with respect to ˆθ F F ˆ ˆ ˆ ˆ ULS ˆ Σ F ULS tr ( tr ( ˆ Σ Σ Σ ˆ Σ S ˆ ˆ ˆ Σ S ˆ ˆ ˆ ˆ i i i j i j i j F ˆ ULS Σ ( ˆ s ˆ i i F ˆ ˆ ˆ ULS tr ( ˆ Σ Σ Σ ˆ ˆ ˆ Σ S ˆ ˆ ˆ ˆ i j i j i ˆ ˆ j F ˆ ULS Σ ( ˆ ˆ ˆ ˆ i js i j 4 4 F ULS ˆ tr ( ˆ Σ ˆ ˆ ˆ ˆ Σ S ˆ ˆ ˆ ˆ i j i j 4 ˆ ˆ ˆ Σ Σ Σ Σˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ i j i j 4 F ˆ ULS Σ ( ˆ ˆ ˆ ˆ ˆ ˆ i js i j 5 5 F ULS ˆ tr ( ˆ Σ ˆ ˆ ˆ ˆ ˆ Σ S ˆ ˆ ˆ ˆ ˆ i j m i jm 5 ˆ 4 ˆ 0 ˆ Σ Σ Σ Σˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ i jm i j m 6

48 5 4 F ˆ ULS Σ ( ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ i js i j ( i j m... q; p a b. ˆF ULS Using the ove resuts with Lemma we have the partia derivatives of with respect to s as foows: Fˆ ˆ ˆ ULS tr ( ˆ Σ i Σ S ( ( ˆ ˆ Σ S s i s ˆ ˆ ˆ F ˆ ˆ ˆ ULS Σ Σ tr +( ˆ Σ i j Σ S s ˆ ˆ ˆ ˆ scd s i j i j scd ˆ ˆ ˆ ˆ tr ( ˆ Σ Σ S ( ˆ s i ( ac bd i i ˆ scd s i cd Fˆ ˆ ˆ ˆ ULS Σˆ Σˆ ˆ tr +( ˆ Σ i j Σ S ˆ ˆ ˆ ˆ ˆ ˆ sscd sef cd efi j s s s i j ˆ Σˆ Σˆ ˆ ˆ tr +( ˆ Σ i j Σ S ˆ ˆ ˆ ˆ s i j i j scd sef ˆ ˆ tr ( ˆ Σ i Σ S ˆ i s scd sef ˆ ˆ ˆ ˆ ˆ i j i ˆ ˆ s ˆ i j cd sef s i cd sef ( 4 4 Fˆ ˆ ˆ ˆ ˆ ULS tr Σ Σ ˆ ˆ ˆ ˆ Σ Σ s ˆ ˆ ˆ ˆ scd sef sgh i j i j 6

49 4 ˆ ˆ ˆ ˆ ˆ ˆ Σ i j +( Σ S ˆ ˆ ˆ ˆ i j s s s s ˆ ˆ ˆ 6 ˆ ˆ Σ Σ ˆ tr +( ˆ Σ i j Σ S ˆ ˆ ˆ ˆ ˆ ˆ s s s s i j i j cd ef gh cd ef gh Σˆ Σˆ ˆ ˆ ˆ tr +( ˆ Σ i j Σ S ˆ ˆ ˆ ˆ s i j i j scd sef sgh 4 ˆ ˆ ˆ 4 ˆ ii j tr ( ˆ Σ Σ S s s s s ˆ s s s s i cd ef gh cd ef gh 4 ˆ ˆ ˆ ˆ ˆ ˆ ˆ i j i j ( ˆ ˆ ˆ s s s ˆ ˆ s s s i j i j ˆ ˆ i ˆ s i cd sef s gh ( p a b p c d p e f p g h cd ef gh cd ef gh where instein s summation convention is to be used ony for subscripts i j and.. The chain rues gˆ gˆ Fˆ s Fˆ s ˆ ˆ gˆ gˆ F F gˆ Fˆ s s Fˆ s s Fˆ s s cd cd cd 6

50 gˆ gˆ Fˆ Fˆ Fˆ s s s Fˆ s s s cd ef cd ef gˆ Fˆ Fˆ gˆ Fˆ ˆ F s s s Fˆ s s s cd ef cd ef 4 4 gˆ gˆ Fˆ Fˆ Fˆ Fˆ s s s s Fˆ s s s s 4 cd ef gh cd ef gh 6 4 gˆ Fˆ Fˆ Fˆ gˆ Fˆ Fˆ ˆ F s s s s Fˆ s s s s cd ef gh cd ef gh 4 gˆ Fˆ Fˆ gˆ Fˆ ˆ F s s s s Fˆ s s s s ( p a b p cd ef gh cd ef gh c d p e f p g h. References Boc R. D. & Lieberman M. (970. Fitting a response mode for n dichotomousy scored items. Psychometria Boc R. D. & Aitin M. (98. Margina maximum ieihood estimation of item parameters: Appication of an M agorithm. Psychometria Ogasawara H. (0. Asymptotic expansions of the distributions of the poyseria correation coefficients. Behaviormetria 8 ( Ogasawara H. (00. Asymptotic expansions of the nu distributions of discrepancy functions for genera covariance structures under nonnormaity. American Journa of Mathematica and Management Sciences 0 ( & Stroud A. H. & Secrest D. (966. Gaussian quadrature formuas. ngewood Ciffs NJ: Prentice-Ha. 64

Parts Manual. EPIC II Critical Care Bed REF 2031

Parts Manual. EPIC II Critical Care Bed REF 2031 EPIC II Critical Care Bed REF 2031 Parts Manual For parts or technical assistance call: USA: 1-800-327-0770 2013/05 B.0 2031-109-006 REV B www.stryker.com Table of Contents English Product Labels... 4

More information

Graphing Square Roots - Class Work Graph the following equations by hand. State the domain and range of each using interval notation.

Graphing Square Roots - Class Work Graph the following equations by hand. State the domain and range of each using interval notation. Graphing quare Roots - lass Work Graph the following equations by hand. tate the domain and range of each using interval notation. 1. y = x + 2 2. f x = x 1. y = x + 4. g x = 2 x 1 5. y = x + 2 + 4 6.

More information

Multiple Random Variables

Multiple Random Variables Multiple Random Variables This Version: July 30, 2015 Multiple Random Variables 2 Now we consider models with more than one r.v. These are called multivariate models For instance: height and weight An

More information

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain

Alberto Maydeu Olivares Instituto de Empresa Marketing Dept. C/Maria de Molina Madrid Spain CORRECTIONS TO CLASSICAL PROCEDURES FOR ESTIMATING THURSTONE S CASE V MODEL FOR RANKING DATA Aberto Maydeu Oivares Instituto de Empresa Marketing Dept. C/Maria de Moina -5 28006 Madrid Spain Aberto.Maydeu@ie.edu

More information

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form

Fourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier

More information

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons

Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Brandon Maone Department of Computer Science University of Hesini February 18, 2014 Abstract This document derives, in excrutiating

More information

Math 124B January 17, 2012

Math 124B January 17, 2012 Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia

More information

CS229 Lecture notes. Andrew Ng

CS229 Lecture notes. Andrew Ng CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view

More information

Solutions for Chapter 3

Solutions for Chapter 3 Solutions for Chapter Solutions for exercises in section 0 a X b x, y 6, and z 0 a Neither b Sew symmetric c Symmetric d Neither The zero matrix trivially satisfies all conditions, and it is the only possible

More information

An Extension of Almost Sure Central Limit Theorem for Order Statistics

An Extension of Almost Sure Central Limit Theorem for Order Statistics An Extension of Amost Sure Centra Limit Theorem for Order Statistics T. Bin, P. Zuoxiang & S. Nadarajah First version: 6 December 2007 Research Report No. 9, 2007, Probabiity Statistics Group Schoo of

More information

Week 6: Differential geometry I

Week 6: Differential geometry I Week 6: Differential geometry I Tensor algebra Covariant and contravariant tensors Consider two n dimensional coordinate systems x and x and assume that we can express the x i as functions of the x i,

More information

A Probability Review

A Probability Review A Probability Review Outline: A probability review Shorthand notation: RV stands for random variable EE 527, Detection and Estimation Theory, # 0b 1 A Probability Review Reading: Go over handouts 2 5 in

More information

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS

SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS ISEE 1 SUPPLEMENTARY MATERIAL TO INNOVATED SCALABLE EFFICIENT ESTIMATION IN ULTRA-LARGE GAUSSIAN GRAPHICAL MODELS By Yingying Fan and Jinchi Lv University of Southern Caifornia This Suppementary Materia

More information

(This is a sample cover image for this issue. The actual cover is not yet available at this time.)

(This is a sample cover image for this issue. The actual cover is not yet available at this time.) (This is a sampe cover image for this issue The actua cover is not yet avaiabe at this time) This artice appeared in a journa pubished by Esevier The attached copy is furnished to the author for interna

More information

EEL 5544 Noise in Linear Systems Lecture 30. X (s) = E [ e sx] f X (x)e sx dx. Moments can be found from the Laplace transform as

EEL 5544 Noise in Linear Systems Lecture 30. X (s) = E [ e sx] f X (x)e sx dx. Moments can be found from the Laplace transform as L30-1 EEL 5544 Noise in Linear Systems Lecture 30 OTHER TRANSFORMS For a continuous, nonnegative RV X, the Laplace transform of X is X (s) = E [ e sx] = 0 f X (x)e sx dx. For a nonnegative RV, the Laplace

More information

Spanning tree invariants, loop systems and doubly stochastic matrices

Spanning tree invariants, loop systems and doubly stochastic matrices Spanning tree invariants, loop systems and doubly stochastic matrices Ricardo Gómez, José Miguel Salazar Instituto de Matemáticas de la Universidad Nacional Autónoma de México. Circuito Exterior, Ciudad

More information

Chapter 4 Multiple Random Variables

Chapter 4 Multiple Random Variables Review for the previous lecture Definition: n-dimensional random vector, joint pmf (pdf), marginal pmf (pdf) Theorem: How to calculate marginal pmf (pdf) given joint pmf (pdf) Example: How to calculate

More information

Homework 1/Solutions. Graded Exercises

Homework 1/Solutions. Graded Exercises MTH 310-3 Abstract Algebra I and Number Theory S18 Homework 1/Solutions Graded Exercises Exercise 1. Below are parts of the addition table and parts of the multiplication table of a ring. Complete both

More information

STRAIGHT LINES EXERCISE - 3

STRAIGHT LINES EXERCISE - 3 STRAIGHT LINES EXERCISE - 3 Q. D C (3,4) E A(, ) Mid point of A, C is B 3 E, Point D rotation of point C(3, 4) by angle 90 o about E. 3 o 3 3 i4 cis90 i 5i 3 i i 5 i 5 D, point E mid point of B & D. So

More information

CHAPTER 2 BOOLEAN ALGEBRA

CHAPTER 2 BOOLEAN ALGEBRA CHAPTER 2 BOOLEAN ALGEBRA This chapter in the book includes: Objectives Study Guide 2.1 Introduction 2.2 Basic Operations 2.3 Boolean Expressions and Truth Tables 2.4 Basic Theorems 2.5 Commutative, Associative,

More information

Research Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation

Research Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation Appied Mathematics and Stochastic Anaysis Voume 007, Artice ID 74191, 8 pages doi:10.1155/007/74191 Research Artice On the Lower Bound for the Number of Rea Roots of a Random Agebraic Equation Takashi

More information

Geometry 3 SIMILARITY & CONGRUENCY Congruency: When two figures have same shape and size, then they are said to be congruent figure. The phenomena between these two figures is said to be congruency. CONDITIONS

More information

OR MSc Maths Revision Course

OR MSc Maths Revision Course OR MSc Maths Revision Course Tom Byrne School of Mathematics University of Edinburgh t.m.byrne@sms.ed.ac.uk 15 September 2017 General Information Today JCMB Lecture Theatre A, 09:30-12:30 Mathematics revision

More information

Statistical Astronomy

Statistical Astronomy Lectures for the 7 th IAU ISYA Ifrane, nd 3 rd Juy 4 p ( x y, I) p( y x, I) p( x, I) p( y, I) Statistica Astronomy Martin Hendry, Dept of Physics and Astronomy University of Gasgow, UK http://www.astro.ga.ac.uk/users/martin/isya/

More information

Math 124B January 31, 2012

Math 124B January 31, 2012 Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat

More information

Pricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications

Pricing Multiple Products with the Multinomial Logit and Nested Logit Models: Concavity and Implications Pricing Mutipe Products with the Mutinomia Logit and Nested Logit Modes: Concavity and Impications Hongmin Li Woonghee Tim Huh WP Carey Schoo of Business Arizona State University Tempe Arizona 85287 USA

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

Brief introduction to groups and group theory

Brief introduction to groups and group theory Brief introduction to groups and group theory In physics, we often can learn a lot about a system based on its symmetries, even when we do not know how to make a quantitative calculation Relevant in particle

More information

Factorization of weighted EP elements in C -algebras

Factorization of weighted EP elements in C -algebras Factorization of weighted EP elements in C -algebras Dijana Mosić, Dragan S. Djordjević Abstract We present characterizations of weighted EP elements in C -algebras using different kinds of factorizations.

More information

Homework 5 Solutions

Homework 5 Solutions Stat 310B/Math 230B Theory of Probabiity Homework 5 Soutions Andrea Montanari Due on 2/19/2014 Exercise [5.3.20] 1. We caim that n 2 [ E[h F n ] = 2 n i=1 A i,n h(u)du ] I Ai,n (t). (1) Indeed, integrabiity

More information

Lecture Notes 4: Fourier Series and PDE s

Lecture Notes 4: Fourier Series and PDE s Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for

More information

( ) Chapter 6 ( ) ( ) ( ) ( ) Exercise Set The greatest common factor is x + 3.

( ) Chapter 6 ( ) ( ) ( ) ( ) Exercise Set The greatest common factor is x + 3. Chapter 6 Exercise Set 6.1 1. A prime number is an integer greater than 1 that has exactly two factors, itself and 1. 3. To factor an expression means to write the expression as the product of factors.

More information

Statistics for Applications. Chapter 7: Regression 1/43

Statistics for Applications. Chapter 7: Regression 1/43 Statistics for Appications Chapter 7: Regression 1/43 Heuristics of the inear regression (1) Consider a coud of i.i.d. random points (X i,y i ),i =1,...,n : 2/43 Heuristics of the inear regression (2)

More information

,,,,..,,., {. (, ),, {,.,.,..,,.,.,,....... {.. : N {, Z {, Q {, Q p { p{ {. 3, R {, C {. : ord p {. 8, (k) {.42,!() { {. 24, () { {. 24, () { {. 25,., () { {. 26,. 9, () { {. 27,. 23, '() { ( ) {. 28,

More information

Cross-sections. Kevin Cahill for 524 January 24, There are two standard ways of normalizing states. We have been using continuum normalization

Cross-sections. Kevin Cahill for 524 January 24, There are two standard ways of normalizing states. We have been using continuum normalization Cross-sections Kevin Cahill for 54 January 4, 013 1 Fermi s Golden Rule There are two standard ways of normalizing states. We have been using continuum normalization p p = δ (3) (p p). (1) The other way

More information

BASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set

BASIC NOTIONS AND RESULTS IN TOPOLOGY. 1. Metric spaces. Sets with finite diameter are called bounded sets. For x X and r > 0 the set BASIC NOTIONS AND RESULTS IN TOPOLOGY 1. Metric spaces A metric on a set X is a map d : X X R + with the properties: d(x, y) 0 and d(x, y) = 0 x = y, d(x, y) = d(y, x), d(x, y) d(x, z) + d(z, y), for a

More information

Final Exam Advanced Mathematics for Economics and Finance

Final Exam Advanced Mathematics for Economics and Finance Final Exam Advanced Mathematics for Economics and Finance Dr. Stefanie Flotho Winter Term /5 March 5 General Remarks: ˆ There are four questions in total. ˆ All problems are equally weighed. ˆ This is

More information

Generalized Bell polynomials and the combinatorics of Poisson central moments

Generalized Bell polynomials and the combinatorics of Poisson central moments Generaized Be poynomias and the combinatorics of Poisson centra moments Nicoas Privaut Division of Mathematica Sciences Schoo of Physica and Mathematica Sciences Nanyang Technoogica University SPMS-MAS-05-43,

More information

INSTITIÚID TEICNEOLAÍOCHTA CHEATHARLACH INSTITUTE OF TECHNOLOGY CARLOW MATRICES

INSTITIÚID TEICNEOLAÍOCHTA CHEATHARLACH INSTITUTE OF TECHNOLOGY CARLOW MATRICES 1 CHAPTER 4 MATRICES 1 INSTITIÚID TEICNEOLAÍOCHTA CHEATHARLACH INSTITUTE OF TECHNOLOGY CARLOW MATRICES 1 Matrices Matrices are of fundamental importance in 2-dimensional and 3-dimensional graphics programming

More information

Radar/ESM Tracking of Constant Velocity Target : Comparison of Batch (MLE) and EKF Performance

Radar/ESM Tracking of Constant Velocity Target : Comparison of Batch (MLE) and EKF Performance adar/ racing of Constant Veocity arget : Comparison of Batch (LE) and EKF Performance I. Leibowicz homson-csf Deteis/IISA La cef de Saint-Pierre 1 Bd Jean ouin 7885 Eancourt Cede France Isabee.Leibowicz

More information

(x x 0 )(x x 1 )... (x x n ) (x x 0 ) + y 0.

(x x 0 )(x x 1 )... (x x n ) (x x 0 ) + y 0. > 5. Numerical Integration Review of Interpolation Find p n (x) with p n (x j ) = y j, j = 0, 1,,..., n. Solution: p n (x) = y 0 l 0 (x) + y 1 l 1 (x) +... + y n l n (x), l k (x) = n j=1,j k Theorem Let

More information

3 Applications of partial differentiation

3 Applications of partial differentiation Advanced Calculus Chapter 3 Applications of partial differentiation 37 3 Applications of partial differentiation 3.1 Stationary points Higher derivatives Let U R 2 and f : U R. The partial derivatives

More information

Engg. Math. I. Unit-I. Differential Calculus

Engg. Math. I. Unit-I. Differential Calculus Dr. Satish Shukla 1 of 50 Engg. Math. I Unit-I Differential Calculus Syllabus: Limits of functions, continuous functions, uniform continuity, monotone and inverse functions. Differentiable functions, Rolle

More information

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES

MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is

More information

( ) is just a function of x, with

( ) is just a function of x, with II. MULTIVARIATE CALCULUS The first ecture covered functions where a singe input goes in, and a singe output comes out. Most economic appications aren t so simpe. In most cases, a number of variabes infuence

More information

ECE 238L Boolean Algebra - Part I

ECE 238L Boolean Algebra - Part I ECE 238L Boolean Algebra - Part I August 29, 2008 Typeset by FoilTEX Understand basic Boolean Algebra Boolean Algebra Objectives Relate Boolean Algebra to Logic Networks Prove Laws using Truth Tables Understand

More information

Algebraic Expressions

Algebraic Expressions Algebraic Expressions 1. Expressions are formed from variables and constants. 2. Terms are added to form expressions. Terms themselves are formed as product of factors. 3. Expressions that contain exactly

More information

MANY BILLS OF CONCERN TO PUBLIC

MANY BILLS OF CONCERN TO PUBLIC - 6 8 9-6 8 9 6 9 XXX 4 > -? - 8 9 x 4 z ) - -! x - x - - X - - - - - x 00 - - - - - x z - - - x x - x - - - - - ) x - - - - - - 0 > - 000-90 - - 4 0 x 00 - -? z 8 & x - - 8? > 9 - - - - 64 49 9 x - -

More information

PAIR OF LINES-SECOND DEGREE GENERAL EQUATION THEOREM If the equation then i) S ax + hxy + by + gx + fy + c represents a pair of straight lines abc + fgh af bg ch and (ii) h ab, g ac, f bc Proof: Let the

More information

4. CONTINUOUS RANDOM VARIABLES

4. CONTINUOUS RANDOM VARIABLES IA Probability Lent Term 4 CONTINUOUS RANDOM VARIABLES 4 Introduction Up to now we have restricted consideration to sample spaces Ω which are finite, or countable; we will now relax that assumption We

More information

Analysis of Emerson s Multiple Model Interpolation Estimation Algorithms: The MIMO Case

Analysis of Emerson s Multiple Model Interpolation Estimation Algorithms: The MIMO Case Technica Report PC-04-00 Anaysis of Emerson s Mutipe Mode Interpoation Estimation Agorithms: The MIMO Case João P. Hespanha Dae E. Seborg University of Caifornia, Santa Barbara February 0, 004 Anaysis

More information

SCHOOL OF MATHEMATICS AND STATISTICS. Mathematics II (Materials) Section A. Find the general solution of the equation

SCHOOL OF MATHEMATICS AND STATISTICS. Mathematics II (Materials) Section A. Find the general solution of the equation Data provided: Formua Sheet MAS250 SCHOOL OF MATHEMATICS AND STATISTICS Mathematics II (Materias Autumn Semester 204 5 2 hours Marks wi be awarded for answers to a questions in Section A, and for your

More information

Topic 2 [312 marks] The rectangle ABCD is inscribed in a circle. Sides [AD] and [AB] have lengths

Topic 2 [312 marks] The rectangle ABCD is inscribed in a circle. Sides [AD] and [AB] have lengths Topic 2 [312 marks] 1 The rectangle ABCD is inscribed in a circle Sides [AD] and [AB] have lengths [12 marks] 3 cm and (\9\) cm respectively E is a point on side [AB] such that AE is 3 cm Side [DE] is

More information

Topper Sample Paper- 3 CLASS XI MATHEMATICS Time Allowed: 3 Hrs Maximum Marks: 100

Topper Sample Paper- 3 CLASS XI MATHEMATICS Time Allowed: 3 Hrs Maximum Marks: 100 Topper Sample Paper- CLASS XI MATHEMATICS Time Allowed: Hrs Maximum Marks: 00. All questions are compulsory.. The question paper consist of 9 questions divided into three sections A, B and C. Section A

More information

We could express the left side as a sum of vectors and obtain the Vector Form of a Linear System: a 12 a x n. a m2

We could express the left side as a sum of vectors and obtain the Vector Form of a Linear System: a 12 a x n. a m2 Week 22 Equations, Matrices and Transformations Coefficient Matrix and Vector Forms of a Linear System Suppose we have a system of m linear equations in n unknowns a 11 x 1 + a 12 x 2 + + a 1n x n b 1

More information

Stat 206: Sampling theory, sample moments, mahalanobis

Stat 206: Sampling theory, sample moments, mahalanobis Stat 206: Sampling theory, sample moments, mahalanobis topology James Johndrow (adapted from Iain Johnstone s notes) 2016-11-02 Notation My notation is different from the book s. This is partly because

More information

Riemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,

Riemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E, Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of

More information

IE 230 Probability & Statistics in Engineering I. Closed book and notes. 120 minutes.

IE 230 Probability & Statistics in Engineering I. Closed book and notes. 120 minutes. Closed book and notes. 10 minutes. Two summary tables from the concise notes are attached: Discrete distributions and continuous distributions. Eight Pages. Score _ Final Exam, Fall 1999 Cover Sheet, Page

More information

Course 2BA1, Section 11: Periodic Functions and Fourier Series

Course 2BA1, Section 11: Periodic Functions and Fourier Series Course BA, 8 9 Section : Periodic Functions and Fourier Series David R. Wikins Copyright c David R. Wikins 9 Contents Periodic Functions and Fourier Series 74. Fourier Series of Even and Odd Functions...........

More information

Postulates and Theorems of Quantum Mechanics

Postulates and Theorems of Quantum Mechanics Postulates and Theorems of Quantum Mechanics Literally, a postulate is something taen as self-evident or assumed without proof as a basis for reasoning. It is simply is Postulate 1: State of a physical

More information

Note : This document might take a little longer time to print. more exam papers at : more exam papers at : more exam papers at : more exam papers at : more exam papers at : more exam papers at : more

More information

Chapter 2. Recurrence Relations. Divide and Conquer. Divide and Conquer Strategy. Another Example: Merge Sort. Merge Sort Example. Merge Sort Example

Chapter 2. Recurrence Relations. Divide and Conquer. Divide and Conquer Strategy. Another Example: Merge Sort. Merge Sort Example. Merge Sort Example Recurrence Relations Chapter 2 Divide and Conquer Equation or an inequality that describes a function by its values on smaller inputs. Recurrence relations arise when we analyze the running time of iterative

More information

Conjugacy in epigroups

Conjugacy in epigroups October 2018 Conjugacy in epigroups António Malheiro (CMA/FCT/NOVA University of Lisbon) @ the York semigroup (joint work with João Araújo, Michael Kinyon, Janusz Konieczny) Acknowledgement: This work

More information

Haruhiko Ogasawara. This article gives the first half of an expository supplement to Ogasawara (2015).

Haruhiko Ogasawara. This article gives the first half of an expository supplement to Ogasawara (2015). Economic Review (Otaru University of Commerce, Vol.66, No. & 3, 9-58. December, 5. Expository supplement I to the paper Asymptotic expansions for the estimators of Lagrange multipliers and associated parameters

More information

(X i,y, i i ), i =1, 2,,n, Y i = min (Y i, C i ), i = I (Y i C i ) C i. Y i = min (C, Xi t θ 0 + i ) E [ρ α (Y t )] t ρ α (y) =(α I (y < 0))y I (A )

(X i,y, i i ), i =1, 2,,n, Y i = min (Y i, C i ), i = I (Y i C i ) C i. Y i = min (C, Xi t θ 0 + i ) E [ρ α (Y t )] t ρ α (y) =(α I (y < 0))y I (A ) Y F Y ( ) α q α q α = if y : F Y (y ) α. α E [ρ α (Y t )] t ρ α (y) =(α I (y < 0))y I (A ) Y 1,Y 2,,Y q α 1 Σ i =1 ρ α (Y i t ) t X = x α q α (x ) Y q α (x )= if y : P (Y y X = x ) α. q α (x ) θ α q α

More information

Bourgain s Theorem. Computational and Metric Geometry. Instructor: Yury Makarychev. d(s 1, s 2 ).

Bourgain s Theorem. Computational and Metric Geometry. Instructor: Yury Makarychev. d(s 1, s 2 ). Bourgain s Theorem Computationa and Metric Geometry Instructor: Yury Makarychev 1 Notation Given a metric space (X, d) and S X, the distance from x X to S equas d(x, S) = inf d(x, s). s S The distance

More information

Conditional Statement: Statements in if-then form are called.

Conditional Statement: Statements in if-then form are called. Monday 9/21 2.2 and 2.4 Wednesday 9/23 2.5 and 2.6 Conditional and Algebraic Proofs Algebraic Properties and Geometric Proofs Unit 2 Angles and Proofs Packet pages 1-3 Textbook Pg 85 (14, 17, 20, 25, 27,

More information

Supplement to Ogasawara s papers on the ADF pivotal statistic, mean and covariance structure analysis, and maximal reliability. Haruhiko Ogasawara

Supplement to Ogasawara s papers on the ADF pivotal statistic, mean and covariance structure analysis, and maximal reliability. Haruhiko Ogasawara Economc Revew (Otaru Unversty of Commerce, Vol.60, o., -44, July, 009. Sulement to Ogasawara s aers on the ADF votal statstc, mean and covarance structure analyss, and maxmal rellty Haruhko Ogasawara In

More information

LINEAR SYSTEMS AND MATRICES

LINEAR SYSTEMS AND MATRICES CHAPTER 3 LINEAR SYSTEMS AND MATRICES SECTION 3. INTRODUCTION TO LINEAR SYSTEMS This initial section takes account of the fact that some students remember only hazily the method of elimination for and

More information

Characteristic Numbers of Matrix Lie Algebras

Characteristic Numbers of Matrix Lie Algebras Commun. Theor. Phys. (Beijing China) 49 (8) pp. 845 85 c Chinese Physical Society Vol. 49 No. 4 April 15 8 Characteristic Numbers of Matrix Lie Algebras ZHANG Yu-Feng 1 and FAN En-Gui 1 Mathematical School

More information

Part 1. For any A-module, let M[x] denote the set of all polynomials in x with coefficients in M, that is to say expressions of the form

Part 1. For any A-module, let M[x] denote the set of all polynomials in x with coefficients in M, that is to say expressions of the form Commutative Algebra Homework 3 David Nichols Part 1 Exercise 2.6 For any A-module, let M[x] denote the set of all polynomials in x with coefficients in M, that is to say expressions of the form m 0 + m

More information

Mathematical Preliminaries

Mathematical Preliminaries Mathematical Preliminaries Economics 3307 - Intermediate Macroeconomics Aaron Hedlund Baylor University Fall 2013 Econ 3307 (Baylor University) Mathematical Preliminaries Fall 2013 1 / 25 Outline I: Sequences

More information

AP Calculus BC Chapter 4 AP Exam Problems. Answers

AP Calculus BC Chapter 4 AP Exam Problems. Answers AP Calculus BC Chapter 4 AP Exam Problems Answers. A 988 AB # 48%. D 998 AB #4 5%. E 998 BC # % 5. C 99 AB # % 6. B 998 AB #80 48% 7. C 99 AB #7 65% 8. C 998 AB # 69% 9. B 99 BC # 75% 0. C 998 BC # 80%.

More information

T-2 In the equation A B C + D E F = G H I, each letter denotes a distinct non-zero digit. Compute the greatest possible value of G H I.

T-2 In the equation A B C + D E F = G H I, each letter denotes a distinct non-zero digit. Compute the greatest possible value of G H I. 2016 ARML Local Problems and Solutions Team Round Solutions T-1 All the shelves in a library are the same length. When filled, two shelves side-by-side can hold exactly 12 algebra books and 10 geometry

More information

Online Appendix. j=1. φ T (ω j ) vec (EI T (ω j ) f θ0 (ω j )). vec (EI T (ω) f θ0 (ω)) = O T β+1/2) = o(1), M 1. M T (s) exp ( isω)

Online Appendix. j=1. φ T (ω j ) vec (EI T (ω j ) f θ0 (ω j )). vec (EI T (ω) f θ0 (ω)) = O T β+1/2) = o(1), M 1. M T (s) exp ( isω) Online Appendix Proof of Lemma A.. he proof uses similar arguments as in Dunsmuir 979), but allowing for weak identification and selecting a subset of frequencies using W ω). It consists of two steps.

More information

Stochastic Variational Inference with Gradient Linearization

Stochastic Variational Inference with Gradient Linearization Stochastic Variationa Inference with Gradient Linearization Suppementa Materia Tobias Pötz * Anne S Wannenwetsch Stefan Roth Department of Computer Science, TU Darmstadt Preface In this suppementa materia,

More information

Answers Investigation 3

Answers Investigation 3 Answers Investigation Applications. a., b. s = n c. The numbers seem to be increasing b a greater amount each time. The square number increases b consecutive odd integers:,, 7,, c X X=. a.,,, b., X 7 X=

More information

STRAND J: TRANSFORMATIONS, VECTORS and MATRICES

STRAND J: TRANSFORMATIONS, VECTORS and MATRICES Mathematics SKE, Strand J STRAND J: TRANSFORMATIONS, VECTORS and MATRICES J4 Matrices Text Contents * * * * Section J4. Matrices: Addition and Subtraction J4.2 Matrices: Multiplication J4.3 Inverse Matrices:

More information

1-2 Measuring and Constructing Segments

1-2 Measuring and Constructing Segments 1-2 Measuring and Constructing Segments Warm Up Lesson Presentation Lesson Quiz Objectives Use length and midpoint of a segment. Construct midpoints and congruent segments. Vocabulary coordinate midpoint

More information

Partial permutation decoding for MacDonald codes

Partial permutation decoding for MacDonald codes Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics

More information

ALGORITHMIC SUMMATION OF RECIPROCALS OF PRODUCTS OF FIBONACCI NUMBERS. F. = I j. ^ = 1 ^ -, and K w = ^. 0) n=l r n «=1 -*/!

ALGORITHMIC SUMMATION OF RECIPROCALS OF PRODUCTS OF FIBONACCI NUMBERS. F. = I j. ^ = 1 ^ -, and K w = ^. 0) n=l r n «=1 -*/! ALGORITHMIC SUMMATIO OF RECIPROCALS OF PRODUCTS OF FIBOACCI UMBERS Staney Rabinowitz MathPro Press, 2 Vine Brook Road, Westford, MA 0886 staney@tiac.net (Submitted May 997). ITRODUCTIO There is no known

More information

Chapter 10: Moments of Inertia

Chapter 10: Moments of Inertia Chapter 10: Moments of Inertia Chapter Objectives To develop a method for determining the moment of inertia and product of inertia for an area with respect to given x- and y-axes. To develop a method for

More information

Research Article The Diagonally Dominant Degree and Disc Separation for the Schur Complement of Ostrowski Matrix

Research Article The Diagonally Dominant Degree and Disc Separation for the Schur Complement of Ostrowski Matrix Hindawi Pubishing Corporation Journa of Appied Mathematics Voume 2013, Artice ID 973152, 10 pages http://dxdoiorg/101155/2013/973152 Research Artice The Diagonay Dominant Degree and Disc Separation for

More information

Random Vectors, Random Matrices, and Matrix Expected Value

Random Vectors, Random Matrices, and Matrix Expected Value Random Vectors, Random Matrices, and Matrix Expected Value James H. Steiger Department of Psychology and Human Development Vanderbilt University James H. Steiger (Vanderbilt University) 1 / 16 Random Vectors,

More information

complicated calculations, similar to HKDSE difficult types ** Please answer any 30 questions ** Level 1 In the figure, sin θ =

complicated calculations, similar to HKDSE difficult types ** Please answer any 30 questions ** Level 1 In the figure, sin θ = Level : single concept, simple calculation, easier than HKDSE types Level : one or two concept, some calculations, similar to HKDSE easy types Level 3: involving high level, logical and abstract thinking

More information

Chapter 7 Introduction to vectors

Chapter 7 Introduction to vectors Introduction to ectors MC Qld-7 Chapter 7 Introduction to ectors Eercise 7A Vectors and scalars a i r + s ii r s iii s r b i r + s Same as a i ecept scaled by a factor of. ii r s Same as a ii ecept scaled

More information

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations)

Homework #04 Answers and Hints (MATH4052 Partial Differential Equations) Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =

More information

Matrix Mathematics. Errata and Addenda for the Second Edition. Dennis S. Bernstein February 13, 2014

Matrix Mathematics. Errata and Addenda for the Second Edition. Dennis S. Bernstein February 13, 2014 Matrix Mathematics Errata and Addenda for the Second Edition Dennis S. Bernstein dsbaero@umich.edu February 13, 2014 This document contains an updated list of errata for the second edition of Matrix Mathematics.

More information

A B CDE F B FD D A C AF DC A F

A B CDE F B FD D A C AF DC A F International Journal of Arts & Sciences, CD-ROM. ISSN: 1944-6934 :: 4(20):121 131 (2011) Copyright c 2011 by InternationalJournal.org A B CDE F B FD D A C A BC D EF C CE C A D ABC DEF B B C A E E C A

More information

Matrices and Determinants

Matrices and Determinants Chapter1 Matrices and Determinants 11 INTRODUCTION Matrix means an arrangement or array Matrices (plural of matrix) were introduced by Cayley in 1860 A matrix A is rectangular array of m n numbers (or

More information

Investigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l

Investigation on spectrum of the adjacency matrix and Laplacian matrix of graph G l Investigation on spectrum of the adjacency matrix and Lapacian matrix of graph G SHUHUA YIN Computer Science and Information Technoogy Coege Zhejiang Wani University Ningbo 3500 PEOPLE S REPUBLIC OF CHINA

More information

MA 201: Partial Differential Equations Lecture - 10

MA 201: Partial Differential Equations Lecture - 10 MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary

More information

chapter 1 vector geometry solutions V Consider the parallelogram shown alongside. Which of the following statements are true?

chapter 1 vector geometry solutions V Consider the parallelogram shown alongside. Which of the following statements are true? chapter vector geometry solutions V. Exercise A. For the shape shown, find a single vector which is equal to a)!!! " AB + BC AC b)! AD!!! " + DB AB c)! AC + CD AD d)! BC + CD!!! " + DA BA e) CD!!! " "

More information

Intermediate Tier - Algebra revision

Intermediate Tier - Algebra revision Intermediate Tier - Algebra revision Contents : Collecting like terms Multiplying terms together Indices Expanding single brackets Expanding double brackets Substitution Solving equations Finding nth term

More information

Homework 3/ Solutions

Homework 3/ Solutions MTH 310-3 Abstract Algebra I and Number Theory S17 Homework 3/ Solutions Exercise 1. Prove the following Theorem: Theorem Let R and S be rings. Define an addition and multiplication on R S by for all r,

More information

Multiple Integrals. Chapter 4. Section 7. Department of Mathematics, Kookmin Univerisity. Numerical Methods.

Multiple Integrals. Chapter 4. Section 7. Department of Mathematics, Kookmin Univerisity. Numerical Methods. 4.7.1 Multiple Integrals Chapter 4 Section 7 4.7.2 Double Integral R f ( x, y) da 4.7.3 Double Integral Apply Simpson s rule twice R [ a, b] [ c, d] a x, x,..., x b, c y, y,..., y d 0 1 n 0 1 h ( b a)

More information

Class-IX CBSE Latest Pattern Sample Paper {Mathematics}

Class-IX CBSE Latest Pattern Sample Paper {Mathematics} Class-IX CBSE Latest Pattern Sample Paper {Mathematics} Term-I Examination (SA I) Time: 3hours Max. Marks: 90 General Instructions (i) All questions are compulsory. (ii) The question paper consists of

More information

Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1

Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rules 1 Steepest Descent Adaptation of Min-Max Fuzzy If-Then Rues 1 R.J. Marks II, S. Oh, P. Arabshahi Λ, T.P. Caude, J.J. Choi, B.G. Song Λ Λ Dept. of Eectrica Engineering Boeing Computer Services University

More information

MATH 221: SOLUTIONS TO SELECTED HOMEWORK PROBLEMS

MATH 221: SOLUTIONS TO SELECTED HOMEWORK PROBLEMS MATH 221: SOLUTIONS TO SELECTED HOMEWORK PROBLEMS 1. HW 1: Due September 4 1.1.21. Suppose v, w R n and c is a scalar. Prove that Span(v + cw, w) = Span(v, w). We must prove two things: that every element

More information

8 APPENDIX. E[m M] = (n S )(1 exp( exp(s min + c M))) (19) E[m M] n exp(s min + c M) (20) 8.1 EMPIRICAL EVALUATION OF SAMPLING

8 APPENDIX. E[m M] = (n S )(1 exp( exp(s min + c M))) (19) E[m M] n exp(s min + c M) (20) 8.1 EMPIRICAL EVALUATION OF SAMPLING 8 APPENDIX 8.1 EMPIRICAL EVALUATION OF SAMPLING We wish to evauate the empirica accuracy of our samping technique on concrete exampes. We do this in two ways. First, we can sort the eements by probabiity

More information