Noncentral Hypergeometric Distribution: J S M J M P
Noncentral Hypergeometric Distribution: J S M J M P
Noncentral Hypergeometric Distribution: J S M J M P
(i) Univariate noncentral hypergeometric distribution There are several ways to derive the univariate noncentral hypergeometric distribution. Source 1 The noncentral hypergeometric distribution with odd ratio is an
exponentially weighted version of the central hypergeometric distribution. Let Y has a noncentral hypergeometric distribution. Then,
m1 m 2 y s y 1 P (Y = y | ) = P0 ( )
where
m1 m 2 y y s y 1 = b m1 m 2 j s j j=a 1
P0 ( ) =
m1 m 2 j j=a s1
b
and
a = max
(0 , s 1
m 2 ) y min
(m 1 , s 1 ) =
b .
. We denote
Y ~ H (m , s ;
or Y ~ H (s , m ;
1
),
where
m = (m 1 , m 2
and
s = (s 1 , s 2 ) . Note that
H (m , s ;1 ) = H (m , s )
is the central hypergeometric distribution. Source 2 Let Y ~ B (m 1 ,
),
X ~ B (m
and
1 1
2
.
2
H (m , s ; ) .
Y ~ H (m , s ;
can be
= log (
) . The distribution of
is
y f ( y , ) exp log ( ) P 0
= exp = exp
k (t ) = log P 0 e t + E (Y
where
[ (
)] log [P (e )] P (e )
0 t +
)=
'
(0 ) =
P0 e
( )
j
t=0
P1 ( P 0 (
) )
P0 e t + t t=0
b t + e j=a =
m1 m 2 j s1 t
j t=0
= P1 (
Similarly,
j t=0
2 P0 e t + P2 ( ) = = 2 t t=0
j=a
m1 m 2 j 2 j j s1
j .
1 = ( 11 , 12 , K , 1 k ) and
= (
21
22
,K ,
2k
).
is
m 1 m 2 y1 y 2 yk L k y s y 1 2 P (Y = y | X + Y = s ) = m1 m 2 j1 j 2 jk j s j 1 2 L k j
where
1 j 1 j
2 j 2 j
The conditional expectation and variance of Y can be derived as follows. The conditional distribution of Y is
(y
| s ,
)=
P (Y = y | s , log
k exp y j j =1
where
[P 0 ( )]
= log (
) and
m1 m 2 j = a j s1
b
P0 ( ) = P0 ( 1 , 2 , K , k ) =
j1 1 L kj k j .
Since it is the standard form of exponential family, the cumulant generating function is
k (t ) = k (t1 , t 2 , K , t k ) = log P0 e t 1 + 1 , e t 2 + 2 , K , e t k + k
Then,
[ (
)] log [P (e
0
, e 2 , K , e k
)]
E (Y l
)=
P0 e t 1 + 1 , e t 2 + 2 , K , e t k + k tl tl = 0 = P0 e 1 , e 2 , K , e k =
where
k (t ) tl t=0
P1 , l (
) P 0 ( )
P0 e t 1 + 1 , e t 2 + 2 , K , e t k + k t l tl = 0 j = m 1 m 2 (t 1 + 1 ) j 1 L e (t k + k ) j k j s j e tl t=0
= j = =
m 1 m 2 (t k + k ) j k (t 1 + 1 ) j 1 L j e e j s j l t=0 m 1 m 2 k jk 1 j1 L j e e l j s j
= P1 , l (
m 1 m 2 jk j1 j2 j L k j s j l 1 2
Similarly,
P2 , lm ( P0 (
P ( ) P1 , m ( ) 1,l ( ) ( ) P P 0 0
P2 , lm ( ) =
m1 m 2 jk j1 j2 j j L 1 2 l m k j s j