Supplement to the papers on the polyserial correlation coefficients and discrepancy functions for general covariance structures
|
|
- Anthony Fox
- 5 years ago
- Views:
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
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 informationGraphing 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 informationMultiple 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 informationAlberto 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 informationFourier 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 informationExpectation-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 informationMath 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 informationCS229 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 informationSolutions 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 informationAn 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 informationWeek 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 informationA 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 informationSUPPLEMENTARY 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 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 informationEEL 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 informationSpanning 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 informationChapter 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 informationHomework 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 informationSTRAIGHT 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 informationCHAPTER 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 informationResearch 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 informationGeometry 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 informationOR 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 informationStatistical 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 informationMath 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 informationPricing 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 informationA 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 informationBrief 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 informationFactorization 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 informationHomework 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 informationLecture 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 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 informationStatistics 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 informationCross-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 informationBASIC 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 informationFinal 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 informationGeneralized 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 informationINSTITIÚ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 informationRadar/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.
> 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 information3 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 informationEngg. 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 informationMATH 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
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 informationECE 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 informationAlgebraic 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 informationMANY 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 informationPAIR 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 information4. 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 informationAnalysis 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 informationSCHOOL 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 informationTopic 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 informationTopper 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 informationWe 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 informationStat 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 informationRiemann 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 informationIE 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 informationCourse 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 informationPostulates 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 informationNote : 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 informationChapter 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 informationConjugacy 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 informationHaruhiko 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 )
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 informationBourgain 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 informationConditional 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 informationSupplement 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 informationLINEAR 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 informationCharacteristic 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 informationPart 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 informationMathematical 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 informationAP 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 informationT-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 informationOnline 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 informationStochastic 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 informationAnswers 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 informationSTRAND 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 information1-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 informationPartial 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 informationALGORITHMIC 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 informationChapter 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 informationResearch 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 informationRandom 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 informationcomplicated 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 informationChapter 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 informationHomework #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 informationMatrix 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 informationA 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 informationMatrices 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 informationInvestigation 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 informationMA 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 informationchapter 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 informationIntermediate 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 informationHomework 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 informationMultiple 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 informationClass-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 informationSteepest 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 informationMATH 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 information8 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