Tests of Hypo PDF
Tests of Hypo PDF
HYPOTHESES
NET/JRF/IAS/ISS/JAM/GATE/STATISTICAL INFERENCE
A. SANTHAKUMARAN
About the Author
A. Santhakumaran received his Ph.D. in Mathematics - Statistics from
the Ramanujan Institute for Advanced Study in Mathematics, Univer-
sity of Madras.He has a rich experience in teaching and research. He
had positions as Associate Professor and Head of the Department of
Statistics at Salem Sowdeswari College,Salem, and Professor of Math-
ematics at the Indian Institute of Food Processing Technology, Than-
javur, Tamil Nadu. He has published research papers in Queuing
Theory, Statistical Quality Control,Neural Networks, Fuzzy
Statistics and Food Processing. He is the author of the book
Mathematical Statistics of Probability Models
A FEW REVIEWS OF THE FIRST EDITION
1.1 Introduction 1
Problems 24
2.1 Introduction 37
Problems 64
Problems 109
Problems 143
Notation 161
Bibliography 163
1.1 Introduction
α = EH0 [φ(X)]
1
= 1 × PH0 {X = 3} + × PH0 {X = 2}
8
+ 0 × PH0 {X = 0 or 1}
1 1 3 11
= 1× + × +0=
8 8 8 64
Note The class of all randomized tests is a larger class that
contains all non - randomized tests
Eθ [φ(X)] ≤ α ∀ θ ∈ ΩH0
i.e., Pθ {X ∈ C} ≤ α ∀ θ ∈ ΩH0
βφ (θ) = Eθ [φ(X)] ∀ θ ∈ Ω
= Pθ {X ∈ C} ∀ θ ∈ Ω
H0 : P0 ∈ P0 = ∪(0, 2) vs H1 : P1 ∈ P1 = ∪(3, 4)
Governing Basic Principles 6
and q 2
2 − x2
π e 0<x<∞
p1 (x) =
0 otherwise
Let H0 : X ∼ p0 vs H1 : X ∼ p1 . Here one who chooses
C = (0, ∞) and C 0 = (−∞, 0). Thus the probability of type
1 error is
Z ∞
PH0 {X ∈ C} = p0 (x)dx = 0
0
test function
0.3 if x = 0
φ(x) = 0.2 if x = 1
0 otherwise
Find the power and the probability of type 1 error of the test
function φ. The probability of type 1 error is
α = EH0 [φ(X)]
β = EH1 [φ(X)]
3
The curve of x1 x2 = 4 is given in Figure 1.1 The shaded
portion is the critical region. The size of the test is
Governing Basic Principles 8
X2
1 Critical
region
x1 x2 = 3/4
0 1 X1
FIgure 1.1 Critical region
α = EH0 [φ(X1 , X2 )]
3
= PH0 X2 ≥
4X1
Z 1Z 1
= dx1 dx2
3 3
4 4x1
Z 1
3
= 1− dx1
3 4x1
4
Z 1
3 3
= 1− − dx1
4 3 4x1
4
1 3 3
= + log
4 4 4
β = EH1 [φ(X1 , X2 )]
Z 1Z 1
= 4x1 x2 dx1 dx2
3 3
4 4x1
1
32
Z
= 2 1− dx1
3
4
16x21
1
9 1 1
Z Z
= 2 x1 dx1 − dx1
3 8 3 x1
4 4
7 9 3
= + log
16 8 4
Governing Basic Principles 9
i.e., α = PH0 {X ∈ C}
= PH0 {X < C} =
6 0
Z c
θ0
= 2
dx
θ0 x
θ0
⇒c =
1−α
θ0
The desirable critical region C = {x | x < 1−α }
β = PH1 {X ∈ C}
θ0
= PH 1 X <
1−α
Z θ0
1−α θ1
= dx
θ1 x2
θ1 θ 1
= 1− + α 6= 1
θ0 θ 0
Also β ≥ α if θ1 ≤ θ0
Governing Basic Principles 10
( i) Eθ [φ(X)] ≤ α θ ∈ ΩH0
Note If ΩH1 contains only one point, then the test is the
MP test. If ΩH1 contains more than one point, then each
θ ∈ ΩH1 the test is the MP test
x 0 1 2
p0 0.1 0.4 0.5
p1 0.02 0.08 0.9
The MP test φ1 is
1 if x = 0, 2
φ1 (x) =
0 otherwise
x 0 1 2
p0 0.1 0.2 0.7
p1 0.3 0.4 0.3
Then α = EH0 [φ1 (X)] = 0.2 and βφ1 = EH1 [φ1 (X)] = 0.4
.. . φ1 is a test function with level α = 0.2. Again define a
test
1 if x = 0
φ2 = 1
2 if x = 1
0 otherwise
so that α = EH0 [φ2 (X)] = 0.2 and βφ2 = EH1 [φ2 (X)] =
1
0.3 + 2 × 0.4 = 0.5. Thus φ2 is also a test function with
level α = 0.2. One can easily see that φ1 is a sub class of φ2
with same level α = 0.2, i.e., φ1 ⊆ φ2 and power of φ1 is less
than the power of φ2 . The test φ2 has more power, since the
critical region of φ2 is larger as the power is larger
β Best test of α0
(1,1)
Cq
1 1
(2, 2)
α0 α
Figure 1.2 Convex set critical region
α = EH0 [φ(X)]
Z 1
3γ
= γdx = 0≤γ≤1
1 4
4
β
(1,1)
C
1
4
1
0 4 α
Figure 1.3 Convex set critical region
x 0 1 x 0 1
H0 : vs H1 :
p0 0.4 0.6 p1 0.2 0.8
and γ1 , γ2 ≥ 0
γ1
1/30
"
"
"
0 1/20 γ2
Figure 1..4 Solution region
(0, 0) 0
( 1 , 1/30) 0.027
( 1/20 , 0) 0.01
• If H0 : θ < θ0 vs H1 : θ ≥ θ0 , or H0 : θ > θ0 vs
H1 : θ ≤ θ0 , then for any θ < θ0 or θ > θ0 , the UMP test
Governing Basic Principles 20
= 0.4γ1 + 0.6γ2
where γ1 ≥ 0, γ2 ≥ 0
x p0 p1 p2 p3
-1 0.2 0.3 0.4 0.7
1 0.8 0.7 0.6 0.3
subject to
3
X 3!
γi pi (1 − p0 )n−i ≤ 0.05
i!(3 − i)! 0
i=0
where p0 = 0.1, 0.3, 0.5
3
X 3!
also maximize β = γi (0.7)i (0.3)n−i
i!(3 − i)!
i=0
subject to
3
X 3!
γi pi (1 − p0 )n−i ≤ 0.05
i!(3 − i)! 0
i=0
where p0 = 0.1, 0.3, 0.5
where 0 ≤ γ0 , γ1 , γ2 , γ3 ≤ 1
Using simplex method to solve the LPP, there are two opti-
mum tests
(i) The optimum test under the alternative H1 : p = 0.6 is
given by
0.4 if x = 3
φ1 (x) =
0 otherwise
The power of the test φ1 (x) is β = 0.0864 under the alterna-
tive H1 : p = 0.6
(ii) The optimum test under the alternative H1 : p = 0.7 is
given by
0.1333 if x = 2
φ2 (x) =
0 otherwise
Governing Basic Principles 23
3
x2 − x =
4
3 1
⇒x = or x = −
2 2
Problems
hypothesis?
1
Ans: α = 4 + 43 log( 43 )
Governing Basic Principles 26
R1
1.7 Find a test function φ which maximizes 0 φ(x)3x2 dx
R1
subject to the two conditions (i) 0 φ(x)dx = 0.1 and
R1
(ii) 0 φ(x)2xdx = 0.19
9
1 if 10 ≤x≤1
Ans: φ(x) =
0 otherwise
1.8 Define the terms (a) size and power of a test and (b)
level of significance
1.9 Define the two kind of errors and the power of a test.
Which error is minimized in a statistical test? why not
both the errors?
1.11 Define (a) critical function (b) power function and (c)
randomized test
1
Let H0 : θ = 2 vs H1 : θ = 34 . Let the critical region
be C = {1, 2}. Find the probability of type 1 error and
power associated with C Ans: α = 34 , β = 15
16
x 0 1 x 0 1
H0 : H1 :
p 0.5 0.5 p 0.4 0.6
0.02 if x = 0 0.02 if x = 1
φ1 (x) = φ2 (x) =
0 otherwise 0 otherwise
x 1 2 3 x 1 2 3
H0 : or H0 :
p1 0.2 0.6 0.2 p2 0.4 0.6 0
x 1 2 3
vs H1 :
p3 0.5 0.5 0
Governing Basic Principles 28
β = 0.0625
x 0 1
H0 :
p0 0.2 0.8
x 0 1 x 0 1
H1 : or H1 :
p1 0.7 0.3 p2 0.1 0.9
β = 0.08
x p0 p1 p2 p3
1 0.05 0.15 0.20 0.25
2 0.95 0.85 0.80 0.75
Let H0 : X ∼ p0 or p1 vs H1 : X ∼ p2 or p3 . Determine
the MP level α = 0.05 test by using the LPP
0.3333 if x = 1
Ans: φ(x) =
0 if x = 2
α = 0.15
Governing Basic Principles 29
5x1 + 3x2 ≤ 15
−x1 + x2 ≤ 1
2x1 + 5x2 ≤ 10 x1 , x2 ≥ 0
y1 + y2 ≤ 1
y2 + y3 ≤ 1
y3 + y4 ≤ 1
y4 + y1 ≤ 1 yi ≥ 0 ∀ i = 1, 2, 3, 4
x1 + x2 ≥ 300
x2 + x3 ≥ 500
x3 + x4 ≥ 400
x4 + x1 ≥ 200, ∀ xi ≥ 0, i = 1, 2, 3, 4
x ≥ 0, y ≥ 0, x ≤ 3
1
x+y ≤ 4
2
x+y ≤ 5
x + 5y ≤ 10
2x + 2y ≤ 5, x ≥ 0, y ≥ 0
composite alternative
(d) MP test is not unique, if it exists
Ans:(a) , (b) and (d)
testing H0 : θ = 0 vs H1 : θ > 0
(a) The test function φ which rejects H0 if and only if
X1 + X2 > C is UMP of its size
(b) The test function φ which rejects H0 if and only if
X1 + X2 > (0.05)2 is of size 0.05
(c) The test function φ which rejects H0 if and only if
X1 + X2 > 0.1 is of size 0.05
(d) None of these Ans: (d)
1 if x > c
φ(x) = γ if x = c where 0 < γ < 1
0 otherwise
then the test is
(a) randomized test
(b) non- randomized
(c) neither randomized nor non- randomized
Governing Basic Principles 35
γ if x > c where 0 < γ < 1
φ(x) =
0 otherwise
as a non-randomized procedure
(b) Every randomized test procedure can be thought of
as a randomized procedure
(c) Every randomized test procedure can be thought of
as equivalent to a non-randomized procedure
(d) All the above Ans: (b)
2.1 Introduction
Commodity 1 2 3 4 5 6 7 8 9 10
Cost ct Rs. 5 10 6 20 25 30 40 50 35 12
The hawker has more profit 1.3 with the given capacity
of Rs. 25.00. for purchasing the commodity set { 1, 2, 3 }
only. Similarly the method of finding the best test for which
the critical region should include the value x of X for the ratio
p1 (x)
p0 (x) is large, where p1 (x) is the pdf of X under H1 and p0 (x)
is the pdf of X under H0 . First consider the case in which
both H0 and H1 are simple, i.e., ΩH0 = {θ0 }, ΩH1 = {θ1 }. For
testing H0 : θ = θ0 vs H1 : θ = θ1 , now states the Neyman
Pearson Lemma(NPL)
(a) Sufficiency
EH0 [φ(X)] = α
then φ is the MP level α test for testing H0 vs H1
(b) Existence
(c) Necessary
If x ∈ C −
We have
EH1 [φ(X)] − EH1 [φ∗ (X)] ≥ kEH0 [φ(X)] − kEH0 [φ∗ (X)] ≥ 0
⇒ β φ ≥ β φ∗
(b) Existence part
n o
Define α(c) = PH0 pp01 (X)
(X)
> c , the probability of the ran-
p1 (X)
dom variable p0 (X) under H0
p1 (X)
.. . 1 − α(c) is the distribution function of the following p0 (X)
under H0
The distribution function 1−α(c) has the following properties
• It is non-negative
• It is non-decreasing
• It is right continuous
Most powerful test 41
= α(c0 )
α − α(c0 )
+ [α(c0 − 0) − α(c0 )]
[α(c0 − 0) − α(c0 )]
= α
(φ − φ∗ )[p1 − kp0 ] = 0
Time x x1 = 0 x2 = 5 x3 = 10 x4 = 15 x5 = 20
Ho p0 0.60 0.20 0.15 0.05 0.00
H1 p1 0.00 0.10 0.20 0.30 0.40
x x1 = 0 x2 =5 x3 = 10 x4 = 15 x5 = 20
p1 (x) 1 4
p0 (x) 0 2 3 6 ∞
p1 (x)
It can be seen that p0 (x) ↑ x ⇒ x > c. so the first choice is
to select the constant c = ∞ for which φ is a test function
For α = 0 or α = 1 is ruled out in practical situations.If α = 0
means 100% , he is in good mood and if α = 1 means 100%
he is not in good mood. Thus if α = 0, the test function is
1 if p0 (x) = 0
φ(x) =
0 if p (x) > 0
1
and α = 1
1 if p1 (x) > 0
φ(x) =
0 if p0 (x) = 0
Now one can choose c = 6 and define the test function
p1 (x)
1 if p0 (x) >6
p1 (x)
φ(x) = γ if p0 (x) =6
0 p1 (x)
if p0 (x) <6
x n
p1 (x) θ1 (1 − θ0 ) (1 − θ1 )
=
p0 (x) θ0 (1 − θ1 ) (1 − θ0 )
= ax b
h i h i
θ1 (1−θ0 ) (1−θ1 )
where a = θ0 (1−θ1 ) > 1 if θ1 > θ0 and b = (1−θ0 ) >0
ax b ↑ x ⇒ x > c
Most powerful test 47
n
X n! n!
θ0x (1 − θ0 )n−x + γ θ0c (1 − θ0 )n−c = α
x=c+1
x!(n − x)! c!(n − c)!
The power of the test is
n
X n! n!
βφ (θ1 ) = θ1x (1 − θ1 )n−x + γ θ1c (1 − θ1 )n−c
x=c+1
x!(n − x)! c!(n − c)!
p1 (x)
For determining the constant c, p0 (x) is given in Table 2.6
p1 (x)
x p0 (x) p1 (x) p0 (x)
p1 (x)
The MP level α = 0.05 is obtained by NPL. Since p0 (x) ↑x
choose c = 32 ⇔ x = 5. The MP test is
1 if x > 5
φ1 (x) = γ1 if x = 5
0 otherwise
= 0.28
consider
1 if t > c
φ(t) = γ if t = c
0 otherwise
β = EH1 [φ(T )]
pθ1 (t)
t pθ0 (t) pθ1 (t) pθ0 (t)
pθ1 (t)
Since pθ0 (t) ↑ t, choose the constant c = 69 ⇔ x = 9 and
define the test function
1 if t > 9
φ1 (t) = γ1 if t = 9
0 otherwise
0.019
Here γ2 = 0.050 = 0.38 < 1, so φ2 is a test function. Thus
the MP level α = 0.05 test is
1 if t > 6
φ2 (t) = 0.38 if t = 6
0 otherwise
nx̄ n 2 2
= e σ2 (θ1 −θ0 ) e− 2σ2 (θ1 −θ0 )
= aebx̄ ↑ x̄ ⇒ x̄ > c
c − θ0 σ
zα = ⇒ c = zα √ + θ0
√σ n
n
where c is given by
β = EH1 [φ(X̄)]
4303
Decision The sample mean x̄ = 25 = 172.12psi is greater
than 171.32psi. The hypothesis H0 : θ = 170 is rejected in
favour of H1 : θ = 172 . The experimenter concludes that
the mean pressure strength of the bottles exceed 170psi, the
pressure strength of the bottles are 172psi with probability
80%
Example 2.8 Let X1 , X2 , X3 , · · · , Xn be iid random sample
drawn from normal distribution with known mean θ and un-
Most powerful test 56
EH0 [φ(X)] = α
PH0 {T > c} = α
(Xi − θ)2
P
c
PH 0 > 2 = α
σ02 σ0
Z ∞
pn (χ2 )dχ2 = α
c
σ02
β = EH1 [φ(T )]
nX o
= PH 1 (Xi − θ)2 > σ02 χ2α
(Xi − θ)2 σ02 2
P
= PH 1 > 2 χα
σ12 σ1
Z ∞
= σ2
pn (χ2 )dχ2
0 χ2
2 α
σ1
Specimen 1 2 3 4 5 6
Strength psi 48.89 52.07 49.29 51.66 52.16 49.72
Speciman 7 8 9 10 11 12
Srength psi 48.00 49.96 49.20 48.10 47.90 46.94
Speciman 13 14 15 16 - -
Strength psi 51.76 50.75 49.86 51.57 - -
and mean θ = 50
The MP level α = 0.05 test is
1 if P(xi − 50)2 > 72.325
φ(t) =
0 otherwise
R∞ c
Since c p16 (χ2 )dχ2 = 0.05 ⇒ 2.75 = 26.3, i.e., c = 72.325
σ02
R∞
where 26.3 p16 (χ2 )dχ2 = 0.05
The power of the test is
Z ∞
β= p16 (χ2 )dχ2 = 0.33
18.08
p1 (x) 1
Consider = ↓ x ∈ (0, 1) ⇒ x < c
p0 (x) 2x
= PH0 {X < c}
Z c
= 2xdx = c2
0
√
⇒c = α
β = EH1 [φ(X)]
Z √α
√
= dx = α
o
Most powerful test 60
p1 (x)
Consider >k
p0 (x)
2
ex √
> c where c = πk
1 + x2
ey
f (y) = where y = x2
1+y
yey
f 0 (y) = > 0 if y > 0
(1 + y)2
c is determined by
β = EH1 [φ(X)]
Z −1.396 Z ∞
= p1 (x)dx + p1 (x)dx
−∞ 1.386
Z −1.386 Z ∞
1 1 1 1
= 2
dx + dx
π −∞ 1+x π 1.386 1 + x2
1 1.386
Z
1
= 1− dx
π −1.386 1 + x2
2
= 1 − tan−1 (1.386) since tan(−θ) = − tan(θ)
π
2 π
= 1− (540 .18950 )
π 1800
= 1 − 0.6021 = 0.3979
Most powerful test 62
Example 2.12
6x for 0 ≤ x ≤ 12
Let p0 (x) =
2 − 2x for 1 ≤ x ≤ 1
2
1
2 if 0 ≤ x ≤ 2
and p1 (x) =
0 otherwise
H0 specifies the function p0 (x) and H1 specifies the function
p1 (x). Find out the MP test of size α = 0.05 for testing H0
vs H1 based on a single observation. Also obtain the power
of the test
p( x) 1
Consider p0 (x) in the interval 0 < x < 2
p1 (x) 1
i.e., = ↓ x ⇒x<c
p0 (x) 12x
1 p1 (x) 1
When 2 < x < 1, = ↑x ⇒x>c
p0 (x) 4(1 − x)
p1 (x)
Thus p0 (x) is monotonic in x
.. . The MP test is
1 x < c or x > c
φ(x) =
0 otherwise
Most powerful test 63
β = EH1 [φ(X)]
Z 0.0913 Z 2
1 1
= dx + dx = 0.6249
0 2 0.8419 2
1
Example 2.13 Find the MP test of size 10 based on a
single observation x of X for testing H0 : X ∼ p0 (x) vs
H1 : X ∼ p1 (x) where
2 1
if 0 < x < 4
p0 (x) =
2 if 1
≤x<1
3 4
1 1
2 if 0 < x < 2
and p1 (x) =
3 1
2 if 2 ≤x<1
Most powerful test 64
p1 (x) 1
Consider p0 (x) in the interval 0 < x < 4
p1 (x)
i.e., = 0.25
p0 (x)
p1 (x)
when 14 ≤ x < 12 ⇒ = 0.75
p0 (x)
p1 (x)
when 12 ≤ x < 1 ⇒ = 2.25
p0 (x)
p1 (x)
Thus ↑ x ⇒x>c
p0 (x)
The MP level α = 0.1 is
1 if x > c
φ(x) =
0 otherwise
R1 2
where 0.1 = c 3 dx ⇒ c = 0.85. The MP level α = 0.1 test is
1 if x > 0.85
φ(x) =
0 otherwiase
2.8
1
6x for 0 ≤ x < 2
Let p(x) =
2 − 2x 1
for 2 ≤x≤1
1
2 for 0 ≤ x ≤ 2
and g(x) =
0 otherwise
H0 specifies the function p(x) and H1 specifies the fre-
quency function g(x). Find the best test of size α = 0.1
for testing H0 vs H1 , based on a single observation x
on X. Also obtain the power of the test
1 if x < 0.1291 or x > 0.7764
Ans: φ(x) =
0 otherwise β = 0.6764
X ∼ p1 (x) where
2 1
if 0 < x < 4
p0 (x) =
2 if 1
≤x<1
3 4
1 1
2 if 0 < x < 2
and p1 (x) =
3 1
4 if 2 ≤x<1
1 if x > 0.925
Ans: φ(x) =
0 otherwise β = 0.1125
1 if yn > θ0
and φ2 (x) =
α if yn ≤ θ0
Explain why the necessity part of the NPL is not con-
tradicted
(a) β ≥ β 0 (c) β = β 0
(b) β ≤ β 0 (d) β 6= β 0 Ans: (c)
Number
of goals 0 1 2 3 4 5
Frequency 92 121 91 50 19 7
x 0 1 2 3
pθ0 (x) 0.01 0.04 0.5 0.45
pθ1 (x) 0.02 0.08 0.4 0.5
3.1 Introduction
pθ1 (x)
function t(x) = t. If for θ0 < θ1 , pθ0 (x) is non-decreasing function
in t
pθ1 (x)
i.e., = g(t) ↑ t
pθ0 (x)
g(t1 ) ≤ g(t2 ) whenever t1 ≤ t2
or
1 if t(x) > c
φ(t) = γ if t(x) = c where −c0 = c
0 otherwise
f (x1 − θ2 ) f (x2 − θ2 )
≤ for all x1 < x2 , θ1 < θ2
f (x1 − θ1 ) f (x2 − θ1 )
i.e., log f (x2 −θ1 )+log f (x1 −θ2 ) ≤ log f (x1 −θ1 )+log f (x2 −θ2 ) (3.1)
(x2 − x1 )
Let t =
(x2 − x1 + θ2 − θ1 )
and x1 − θ1 = t(x1 − θ2 ) + (1 − t)(x2 − θ1 )
x2 − θ2 = (1 − t)(x1 − θ2 ) + t(x2 − θ1 )
log f (x1 −θ1 ) ≥ t log f (x1 −θ2 )+(1−t) log f (x2 −θ1 ) (3.2)
and log f (x2 −θ2 ) ≥ (1−t) log f (x1 −θ2 )+t log f (x2 −θ1 ) (3.3)
Adding (3.2) and (3.3) give (3.1). Thus pθ (x) has MLR prop-
erty
Assume pθ (x) has MLR property, to prove that f (x − θ) is
log - concave. Let a < b be any real number, and set
x1 − θ2 = a, x2 − θ1 = b
and x1 − θ1 = x2 − θ2
a+b
But x1 − θ1 = x2 − θ2 = and (3.1) becomes
2
a+b
log f (a) + log f (b) ≤ 2 log f ∀ a, b ∈ <
2
⇒ f (x − θ) is concave
Example 3.1 Show that there is no UMP test exist for test-
ing H0 : θ = θ0 vs H1 : θ 6= θ0 in one parameter exponential
family of distributions
The one parameter exponential family of densities is
θ n!
i.e., pθ (x) = (1 − θ)n ex log [ 1−θ ]
x!(n − x)!
= c(θ)eQ(θ)t(x) h(x)
θ
where c(θ) = (1 − θ)n , Q(θ) = log[ 1−θ ], t(x) = x and
n! θ −1
h(x) = x!(n−x)! . Further Q(θ) = − log( 1−θ ) = − log 1−θ
θ
= c(θ)eQ(θ)t(x) h(x)
Pn
where c(θ) = θn , Q(θ) = log(1 − θ), and t(x) = i=1 xi
−1
Q0 (θ) = 1−θ < 0, ⇒ Q(θ) ↓ θ (θ ≥ 0)
.. . The UMP test is
1 if t(x) < c
φ(t) = γ if t(x) = c
0 otherwise
(n+t−1) θn (1 − θ)t if t = 0, 1, 2, 3, · · ·
(n−1)!t!
pθ (x) =
0 otherwise
Eθ0 [φ(T )] = α
0.05
⇒γ= 0.05265 = 0.95003 and 0 ≤ γ ≤ 1
.. . φ is a test function. The UMP level α = 0.05 is
0.95003 if x = 2
φ(x) =
0 otherwise
pθ (x) = pθ (xi , x2 , · · · , xn )
P
= θn e−θ xi
= c(θ)eQ(θ)t(x) h(x)
Pn
Let T = i=1 Xi ∼ G(n, 1θ ) and its pdf is
θn −θt n−1
pθ (t) = e t , 0<t<∞
Γn
1
Let 2θt = y, dt = 2θ dy
θ n − y y n 1
.. . p(y) = e 2
Γn 2θ 2θ
1 y 2n
= 2n e− 2 y 2 −1
2 2 Γ 2n
2
Z ∞
i.e., p(χ2 )dχ2 = 0.05 at 30 degrees of freedom
18.493
where n = 15 and θ0 = 2
.. . 2cθ0 = 18.493 ⇒ c = 4.6232. The UMP test is
1 if t < 4.6232
φ(t) =
0 otherwise
= Pθ {T < 4.6232} ∀ θ ≥ 2
Z 4.6232 n
θ −θt n−1
= e t dt ∀ θ ≥ 2
0 Γn
Z 2×4.6232θ
1 y 2n
= 2n e− 2 y 2 −1
2n
0 22Γ2
Neyman-Pearson Lemma 88
.5
.1
0 2 4 6 8 θ
Figure 3.1 UMP test power curve
n
θ1
n θ1 Observations θ2 θ2
4 2 {2, 3, 4, 5} 5 0.0256
3 2 {2, 3, 4} 4 0.1250
2 2 {2, 3} 3 0.4444
α = EH0 [φ(T )]
= PH0 {T > 8}
= PH1 {T > 8}
1 2
where t(x) = x̄, c(θ) = e 8 θ and
25 1 P 2
h(x) = 10√1 2π e− 200 xi
1
Q(θ) = θ ↑ θ if θ > 75
4
where c is determined by
α = EH0 [φ(T )]
Figure 3.2 is the power curve of the UMP level α = 0.1 test
for the hypothesis H0 : θ = 75 vs H1 : θ > 75
1.0
β(θ)
.50
.10
75 77 79 81 θ
The UMP test power function values are given in Table 3.3
1 1 P 2
pθ (x) = √ e− 32 (xi −θ)
4 2π
= c(θ)eQ(θ)t(x) h(x)
nθ 2
where Q(θ) = nθ
16 , t(x) = x̄, c(θ) = e− 32 and
Neyman-Pearson Lemma 94
1
x2i
P
h(x) = √1 e− 32 , Q(θ) ↓ θ if θ < 25
4 2π
The UMP test is given by
1 if x̄ < c
φ(x̄) =
0 otherwise
β(θ) = Eθ [φ(X̄)]
= Pθ {X̄ < c}
X̄ − θ c−θ
= Pθ √ < √
σ/ n σ/ n
c−θ
= Pθ Z < √
σ/ n
X̄−θ
where Z = √
σ/ n
c−θ
= φ √
σ/ n
Given β(25) = 0.1
c−25
i.e., φ σ/ √
n
= 0.10
Z c−25
√
σ/ n
p(z)dz = 0.10 (3.6)
−∞
Z −1.28
But p(z)dz = 0.10 (3.7)
−∞
From equations (3.6) and (3.7)
√
(c − 25) n
⇒ = −1.28 (3.8)
σ
c − 23
and β(23) = 0.90, i.e., φ √ = 0.90
σ/ n
√
(c − 23) n
⇒ = 1.28 (3.9)
σ
Solve the equations (3.8) and (3.9) with σ = 4, give n =
26.265, i.e., n = 26 or 27 and c = 24.001. The UMP test is
1 if x̄ < 24.001
φ(x̄) =
0 otherwise
Neyman-Pearson Lemma 95
c
Comparing (3.10) and (3.11) give σ02
= 31.41
.. . c = 5 × 31.41 = 157.05. The UMP level α = 0.05 test is
1 if P20 (xi − θ)2 > 157.05
i=1
φ(t) =
0 otherwise
1
β(σ)
.5
.1
5 7 9 11 13
σ2
Figure 3.3 UMP test power curve
where c1 is determined by
α = EH0 [φ(X̄)]
X̄ − θ0 c1 − θ0
= Pθ0 √ > √
σ/ n σ/ n
Z ∞
= √ p(z)dz (3.14)
n
(c1 −θ0 ) σ
X̄−θ
where Z = √0
σ/ n
∼ N (0, 1)
c3 − θ0
From equation (3.16) ⇒ − z α2 = √
σ/ n
Z −z α
2 α
where −z α is obtained from p(z)dz =
2
−∞ 2
σ
⇒ c3 = −z α2 √ + θ0
n
c4 − θ0
From equation (3.17) z α2 = √
σ/ n
Z ∞
where z α2 is obtained from p(z)dz = α/2
zα
2
σ
⇒ c4 = z α2 √ + θ0
n
The power function of θ ∈ ΩH1 is
√
Z −zα/2 −(θ−θ0 ) n/σ Z ∞
βφ3 (θ) = p(z)dz + √
p(z)dz
−∞ zα/2 −(θ−θ0 ) n/σ
(i) H0 : θ = 10 vs H1 : θ > 10
(ii) H0 : θ = 10 vs H1 : θ < 10
(iii) H0 : θ = 10 vs H1 : θ 6= 10
As in Example 3.12 , the UMP level α = 0.05 test
for the hypothesis H) : θ = 10 vs H1 : θ > 10 is given by
1 if x̄ > 10.66
φ1 (x̄) =
0 otherwise
Figure 3.4 shows the UMP right tail test power curve
Neyman-Pearson Lemma 103
1.0
βφ1 (θ)
,50
,10
As in Example 3.12, the UMP level α = 0.05 test for the hypothesis
H0 : θ = 10 vs H1 : θ < 10 is given by
1 if x̄ < 9.34
φ2 (x̄) =
0 otherwise
1.0
βφ2 (θ)
.50
.10
6 8
10 θ 9
Figure 3.5 UMP left tail test power curve
βφ3 (θ) 0.05 0.08 0.17 0.32 0.52 0.80 0.85 0.94 0.98
θ 10 9.8 9.6 9.4 9.2 9.0 8.8 8.6 8.4
βφ3 (θ) 0.05 0.08 0.17 0.32 0.52 0.80 0.85 0.94 0.98
1.0
βφ3 (θ)
.5
.1
0
8.4 9.4 10 10.6 11.6 θ
Figure 3.6 Two tails test power curve
Example 3.14 Sixteen items put on test and the failure times(
in hours) are 13.4, 14.2 , 28.8, 29.0, 29.8, 33.0, 37.8, 39.6, 43.4,
49.8, 54.8, 58.2, 67.4, 70.2, 91.2, 44.04. Assuming the failure time
distribution to be given by N (θ, σ 2 = 16). Obtain the UMP level
α = 0.05 test for the failure times coming from the normal popu-
lation with mean failure time H0 : θ ≤ 40 vs H1 : θ > 40
The hypothesis is to test that H0 : θ ≤ 40 vs H1 : θ > 40
The joint pdf of the random sample of X1 , X2 , X3 , · · · , Xn from
N (θ, 16) is
n 1
(xi −θ)2
P
√1
2πσ
e− 2σ2 if −∞ < x < ∞
pθ (x) =
0 otherwise
Table 3.8 shows the UMP test power function values for testing the
hypothesis H0 : θ ≤ 40 vs H1 : θ > 40
Consider
or
∂
1 if ∂θ log pθ (x) > k
φ0 (x) = ∂
γ if ∂θ log pθ (x) = k
0 ∂
if ∂θ log pθ (x) < k
Example 3.15 Find the LMP level α test for testing the hypoth-
esis H0 : θ ≤ 0 vs H1 : θ > 0 of the pdf
1 1
if − ∞ < x < ∞
π 1+(x−θ)
pθ (x) =
0 otherwise
∂ 2x
where ∂θ log pθ (x) |θ=θ0 = 1+x2 and k is determined by
Eθ0 [φ0 (X)] = α
Example 3.16 Construct LMP level α test, if any for testing H0 :
σ = σ0 vs H: σ > σ0 on the basis of n observations drawn from
N (θ, σ 2 ) where θ is known.
The joint pdf of the random sample X1 , X2 , X3 , · · · , Xn is
n
1 1
P 2
pσ (x) = √ e− 2σ2 (xi −θ)
2πσ 2
n n 1 X
log pσ (x) = − log 2π − log σ 2 − 2 (xi − θ)2
2 2 2σ
∂ n 1 X
log pσ (x) = − + 3 (xi − θ)2
∂σ σ σ
n
Let c = σ03 (k + σ0 ), then the LMP level α test is given by
(xi − θ)2 > c
1 P
if
φ0 (t) =
0 otherwise
where c is determined by
EH0 [φ0 (T )] = α
nX o
PH 0 (Xi − θ)2 > c = α
Xi − θ)2
P
c
PH 0 > 2 = α
σ2 σ0
Z 0∞
pn (χ2 )dχ2 = α (3.18)
c
σ02
Neyman-Pearson Lemma 109
R∞
and χ2α
pn (χ2 )dχ2 = α( using χ2 table) (3.19)
c
From equations (3.18) and (3.19) give σ02
= χ2α ⇒ c = σ02 χ2α . The
LMP level α test is
1 if P(x − θ)2 > σ 2 χ2
i 0 α
φ0 (t) =
0 otherwise
Problems
3.1 Stating the results used, obtain the UMP level α test for
testing H0 : θ < θ0 vs H1 : θ > θ0 , using a random sample
of size n from pθ (x) = θe−θx , x > 0, θ > 0
3.4 Starting the necessary results to obtain the UMP test for
testing θ ≤ θ0 vs θ > θ0 , using a random sample of size n
from a Poisson distribution with parameter θ
3.5 Show that UMP tests for testing one sided hypothesis for the
family of distributions with MLR property
3.6 Define UMP test, State and establish a sufficient condition for
existence of the UMP test for testing one sided hypothesis
against an one sided alternative in the case of exponential
families
3.7 Show that for a family having MLR property, there exists
UMP test for testing one sided hypothesis against one sided
alternative
Neyman-Pearson Lemma 110
3.8 Define the MLR property. Verify whether the following family
has the MLR property
1
2(1 − θ) if 0 < x < 2
pθ (x) = 1
2θ if 2 ≤ x < 1, 0 < θ < 1
0 otherwise
3.9 When do you say that the real parameter family of densities
pθ (x) to have MLR ? Explain
3.12 Obtain the UMP test, if one such exists, for testing the hy-
1 1
pothesis H0 : θ ≤ 4 vs H1 : θ > 4 using a random sample
from the uniform distribution over (0, 1)
1
3.13 Examine whether a UMP test exists for testing H0 : θ ≤ 2
1
vs H1 : θ > 2 given n independent observations on a random
variable X with pdf
θ if 0 < x < 1
pθ (x) = 1−θ if 1 < x < 2, 0 < θ < 1
0 otherwise
3.14 Show that the family of Cauchy densities does not possess
MLR property
Neyman-Pearson Lemma 111
3.15 If pθ (x) is a family with MLR property, show that there exists
a UMP test of H0 : θ ≤ θ0 vs H1 : θ > θ0
3.20 What are LMP tests? When are they preferred? Let
X1 , X2 , X3 , · · · , Xn be a random sample from a distribution
with density
(x−1)2 x2
e− 2 e− 2
pθ (x) = θ √ + (1 − θ) √ , 0 ≤ θ < 1, x ∈ <
2π 2π
Derive the LMP test for testing H0 : θ = 0 vs H1 : θ > 0
3.21 Describe LMP test. How will you derive the same?
3.31 A data set gave a 95% confidence interval (2.5, 3.6) for the
mean µ of a normal population with known variance . Let
µ0 < 25 be a fixed number. If we use the same data to
test H0 : µ = µ0 va H1 : µ 6= µ0 . Which of the following
statements are true?
(a) H0 would be necessarily rejected at α = 0.1
(b) H0 would be necessarily rejected at α = 0.025
(c) For α = 0.1, the information is not enough to draw a
conclusion
(d) For α = 0.025, the information is not enough to draw a
conclusion Ans: (a) and (d)
dn ∀ n Ans: (b)
PEARSON LEMMA
4.1 Introduction
R
0 ≤ φ(x) ≤ 1 and φ(x)pi (x)dx = ci , i = 1, 2, 3, · · · , m a
sufficient condition of the form is
1 if pm+1 (x) > Pm ki pi (x)
i=1
φ(x) = Pm
0 if p
m+1 (x) < i=1 ki pi (x)
R
for some k1 , k2 , · · · , km such that φ(x)pi (x)dx = ci , i =
1, 2, · · · , m
Theorem 4.1 For testing the hypothesis H0 : θ ≤ θ1 or
θ ≥ θ2 (θ1 < θ2 ) vs H1 : θ1 < θ < θ2 in an one parameter
exponential family, there exists an UMP test given by
1 if ci < t < c2
φ(t) = γi if t = ci , i = 1, 2
0 otherwise
where c0 s and γ 0 s are determined by Eθ1 [φ(T )] = Eθ2 [φ(T )] = α
Consider
k1 c(θ1 ) [Q(θ1 )−Q(θ0 )]t k2 c(θ2 ) [Q(θ2 )−Q(θ0 )]t
1 > e + e
c(θ0 ) c(θ0 )
1 > a1 eb1 t + a2 eb2 t
Case(i) If a1 ≤ 0 and a2 ≤ 0
then g(t) = a1 eb1 t + a2 eb2 t < 1 is true ∀ t
In this case H0 : θ ≤ θ1 or θ ≤ θ2 is always rejected
since the entire sample space is the critical region which is
not our interest for testing the hypothesis
Case (ii) If a1 ≥ 0 and a2 ≤ 0
then g(t) = a1 eb1 t + a2 eb2 t ↓ t
since g(t)0 = a1 b1 eb1 t + a2 b2 eb2 t ≤ 0 ∀ t
⇒ g(t) is monotonic decreasing function in t
.. . The test function has critical region on one sided (left tail)
of the form as in Figure 4.1
g(t)
g(t) > 1
1
g(t) < 1
c1 t
Figure 4.1 Left tail critical region
The test function is
1 if t < c1
φ(t) =
0 if t > c
1
For such test the power function is increasing and hence the
test does not satisfy the two sided conditions
g(t)
g(t) > 1
1
g(t) < 1
c2
t
Figure 4.2 Right tail critical region
The power function is increasing and hence, the test does not
satisfy the two sided conditions
log c
say t0 = b2 −b1 . The test function φ for the hypothesis
H00 : θ = θ1 or θ = θ2 vs H10 : θ = θ0 ∈ (θ1 , θ2 ) has the form
as in Figure 4.3
g(t)
t0 c2
c1 t
Figure 4.3 Two sided critical region
Eθ2 [φ(T )] = α
1 if pθ0 (x) > k ⇒ t(x) < c1
Pθ1 (x)
i.e., φ(t) =
0 otherwise
≤ α for θ < θ1
Hence, the test φ is the UMP for testing the two sided hy-
pothesis H0 : θ ≤ θ1 or θ ≥ θ2 vs H1 : θ1 < θ < θ2
Example 4.1 Let X be a random sample drawn from bino-
mial distribution with parameters θ and n = 2. Derive the
1 1
UMP level α = 27 for testing the hypothesis H0 : θ ≤ 3 or
1 1
θ≥ 2 vs H1 : 3 < θ < 12 . Also draw the power curve
The pmf of the random variable X is
2!
pθ (x) = θ2 (1 − θ)2−x , x = 0, 1, 2
x!(2 − x)!
x
2! θ
= (1 − θ)2
x!(2 − x)! 1 − θ
θ
= c(θ)ex log 1−θ h(x)
Applications of GNPL 122
θ 2!
where c(θ) = (1 − θ)2 , Q(θ) = log 1−θ and h(x) = x!(2−x)!
test φ1 is
1 if c1 < x < c2
γ1 if x = c1
φ1 (x) =
γ2 if x = c2
0 otherwise
where c1 , c2 and 0 ≤ γ1 , γ2 ≤ 1 are found by using the
1
Eθ1 [φ1 (X)] = Eθ2 [φ1 (X)] = 27 . By trial and er-
ror method choose c1 = 1 and c2 = 2, define the test function
φ1 as
1 if 1 < x < 2
γ1 if x = 1
φ1 (x) =
γ2 if x = 2
0 otherwise
1 1
so that Eθ1 = 1 [φ1 (X)] = 27 and Eθ2 = 1 [φ1 (X)] = 27
3 2
1
γ1 Pθ1 = 1 {X = 1} + γ2 Pθ1 = 1 {X = 2} =
3 3 27
1
and γ1 Pθ2 = 1 {X = 1} + γ2 Pθ2 = 1 {X = 2} =
2 2 27
1
4γ1 + γ2 =
3
8
and 4γ1 + 2γ2 =
27
1
Solving the above two equations, one gets γ2 = − 27 < 0 and
5
γ1 = 54 . Since γ2 < 0, φ1 is not a test function
Applications of GNPL 123
1 1
so that Eθ1 = 1 [φ2 (X)] = 27 and Eθ2 = 1 [φ2 [(X)] = 27
3 2
1
γ1 Pθ1 = 1 {X = 0} + γ2 Pθ1 = 1 {X = 1} =
3 3 27
1
and γ1 Pθ2 = 1 {X = 0} + γ2 Pθ2 = 1 {X = 1} =
2 2 27
1
4γ1 + 4γ2 =
3
4
and γ1 + 2γ2 =
27
7
Solving the above two equations, one gets γ2 = 108 and γ1 =
1
54 . Since 0 ≤ γ1 , γ2 ≤ 1, so the φ2 is a test function.Thus the
1
UMP level α = 27 test is
1
54 if x = 0
φ2 (x) = 7
108 if x = 1
0 otherwise
β(θ)
.0380
.0375
.0370
.33 .40 .45 .50
θ
Figure 4,4 UMP two sided test power curve
• Eθ [φ(X)] ≤ α, θ ∈ ΩH0
• Eθ [φ(X)] ≥ α, θ ∈ ΩH1
α = EH0 [φ(X)]
1 9
= PH0 X < or X >
10 10
Z 1 Z 1
10
= p(x)dx + p(x)dx
9
0 10
1
Z Z 1
10 1 7
= dx + dx
0 4 9 4
10
= 0.20
β = EH1 [φ(X)]
Z 1 Z 1
10
= g(x)dx + g(x)d
9
0 10
1
Z
10
Z 1
θ−1
= θx dx + θxθ−1 dx
9
0 10
θ
1
= 2 > 0.20 ∀ θ > 0
10
where c is determined by
α = EH0 [φ(X)]
Z 1
1
= 1+ dx
c 2
3 14
0.10 = [1 − c] ⇒ c =
2 15
The test φ is
1 14
if x > 15
φ(x) =
0 otherwise
β = EH1 [φ(X)]
Z 1
1+θ 1
= (1 + θ)dx = > 0.1 ∀ θ ≥
14 15 2
15
1 x
pθ (x) = e− θ , θ > 0, x > 0
θ
α = EH0 [φ(X)]
= PH0 {X > c}
Z ∞
= e−x dx = e−x dx = e−c ⇒ c = − loge α
c
Applications of GNPL 128
β = EH1 [φ(X)]
has a minimum at θ = θ0
c0 (θ0 )
.. . Eθ0 [ψ(T )] + Eθ [T ψ(T )] = 0
c(θ)
c0 (θ0 )
⇒ α + Eθ0 [T α] = 0
c(θ0 )
c0 (θ0 )
⇒ Eθ0 [T ] = −
c(θ0 )
c0 (θ0 )
.. . Eθ0 [T ψ(T )] = − Eθ [ψ(T )]
c(θ0 ) 0
= Eθ0 [T ]Eθ0 [ψ(T )]
= αEθ0 [T ]
Hence ψ ∈ D ⇒ Uα ⊂ Dα
Now to find the UMP test in Uα which will be the required
UMPU test. For θ0 6= θ0 , Eθ0 [ψ(T )] is maximized in Dα
subject to (iii) and (iv) is given by using GNPL
1 if p 0 (t) > k1 p (t) + k2 tp (t)
θ θ0 θ0
i.e., ψ(t) =
0 otherwise
0
c(θ0 )eθ t > k1 c(θ0 )eθ0 t + k2 tc(θ0 )eθ0 t
c(θ0 ) (θ0 −θ0 )t
⇒ e > k1 + tk2
c(θ0 )
c(θ0 )
⇒ aebt > k1 + tk2 where b = θ0 − θ0 , a = c(θ0 )
One can show that the inequality aebt > k1 + tk2 implies t
must be either outside certain interval as in Figure 4.5 and
Applications of GNPL 131
y = aebt , b < 0
y
@ y = k1 + k2 t
@
@
@
@
c1 c2 t
t∈
/ (c1 , c2 )
Figure 4.5 Two sided critical regions
y = aebt , b > 0
y
y = k1 + k2 t
c1 c2
t∈/ (c1 , c2 ) t
Figure 4.6 Two sided critical regions
1 if t < c1 or t > c2
ψ(t) =
0 otherwise
Applications of GNPL 132
y = aebt , b < 0
y y = k1 + k2 t
c2 t > c2 t
Figure 4.7 One sided right tail critical region
y = aebt , b > 0
y = k1 + k2 t
t < c1 c1 t
Figure 4.8 One sided left tail critical region
or
1 if t > c2
ψ(t) =
0 otherwise
If there is only one point of intersection the critical region
will be either left tail or right tail of the distribution of t (
depending upon the sign of a). But a single tail test will
lead to a monotonic power curve by virtue of the NPL which
against the property of unbiasedness . Thus the one sided
tests are ruled out. But unbiasedness ⇒ β(θ) has a turning
Applications of GNPL 133
n!
pθ (x) = θx (1 − θ)n−x , x = 0, 1, 2 · · · , n
x!(n − x)!
+γ2 PH0 {X = c2 }
1 −1
cX n
n x n−x
X n x
= θ0 (1 − θ0 ) + θ (1 − θ0 )n−x
x x 0
x=0 x=c2 +1
n ci
+ γ1 θ (1 − θ0 )n−ci , i = 1, 2
ci 0
2 −1
cX
n x
1−α = θ (1 − θ0 )n−x
x 0
x=c1 +1
2
X n ci
+ (1 − γi ) θ (1 − θ0 )n−ci (4.3)
ci 0
i=1
1 −1
cX n
n x n−x
X n x
αnθ0 = x θ0 (1 − θ0 ) + x θ (1 − θ0 )n−x
x x 0
x=0 x=c2 +1
n c1 n−c1 n c2
+ γ1 c1 θ0 (1 − θ0 ) + γ2 c2 θ (1 − θ0 )n−c2
c1 c2 0
Again consider
n x n!
x θ (1 − θ)n−x = x θx (1 − θ)n−x
x x!(n − x)!
n(n − 1)!
= θx (1 − θ)n−x
x(x − 1)!(n − x)!
n − 1 x−1
= nθ θ (1 − θ)(n−1)−(x−1)
x−1
1 −1
cX
n − 1 x−1
nθ0 α = nθ0 θ (1 − θ0 )n−x
x−1 0
x=0
n
X n − 1 x−1
+ nθ0 θ (1 − θ0 )n−x
x−1 0
x=c2 +1
n − 1 c1 −1
+ γ1 nθ0 θ (1 − θ0 )(n−1)−(c1 −1)
c1 − 1 0
n − 1 c2 −1
+ γ2 nθ0 θ (1 − θ0 )(n−1)−(c2 −1)
c2 − 1 0
Applications of GNPL 135
2 −1
cX
n − 1 x−1
1−α = θ (1 − θ0 )n−x
x−1 0
x=c2 +1
2
X n − 1 ci −1
+ (1 − γi ) θ (1 − θ0 )(n−1)−(ci −1)
ci − 1 0
i=1
(4.4)
α
PH0 {X < c1 } + γ1 PH0 {X = c1 } =
2
α
PH0 {X > c2 } + γ2 PH0 {X = c2 } =
2
= PH1 {X = 0} + PH1 {X = 8}
+ 0.675PH1 {X = 1} + 0.675PH1 {X = 7}
θ ∈ [0, 1]
1
β(θ)
.5
.1
0 .3 .5 .8 1 θ
Figure 4.9 UMPU test power curve
Example 4.8 Obtain the UMPU level α test for testing
the hypothesis H0 : θ = θ0 vs H1 : θ 6= θ0 based on normal
distribution with mean θ and known variance σ 2 for a sample
of size n
The joint pdf of the random sample X1 , X2 , · · · , Xn is
n
1 1 P 2 n 2 θ
pθ (x) = √ e− 2σ2 xi − 2σ2 θ + σ2 nx̄
2πσ
= c(θ)eQ(θ)t h(x)
n 2
where c(θ) = e− 2σ2 θ , Q(θ) = nθ σ2
, t = x̄ and h(x) =
n 1 P 2
√1
2πσ
e− 2σ2 xi . Thus pθ (x) is an one parameter expo-
nential family. The UMPU level α test is
1 if t < c1 or t > c2
φ(t) =
0 otherwise
Applications of GNPL 138
From equation(4.5)
Z c1 Z ∞ Z ∞
(t−θ0 )p(t−θ0 )dt+ (t−θ0 )p(t−θ0 )dt = α (t−θ0 )p(t−θ0 )dt
−∞ c2 −∞
Z ∞ Z ∞
=α tp(t − θ0 )dt − αθ0 p(t − θ0 )dt = α[θ0 − θ0 ] = 0
−∞ −∞
Z −(c1 −θ0 ) Z ∞
tp(t)dt = − tp(t)dt
−∞ (c2 −θ0 )
Z (c2 −θ0 )
= tp(t)dt
−∞
Example 4.9 Draw the power curve and testing the hy-
pothesis H0 : θ = 2 vs H1 : θ 6= 2 with level α = 0.05 based
on N (0, 1) by taking a sample of size n = 25
As in Example 4.8, the UMPU level 0.05 test is given by
1 if x̄ < c1 or x̄ > 2θ0 − c1
φ(x̄) =
0 otherwise
Applications of GNPL 139
Table 4.3 shows the UMPU test power function values. Fig-
ure 4.10 visualizes the UMPU test power curve
β(θ)
1.8 2θ 2.2
Figure 4.10 UMPU test power curve
Xi2
P
Let Y = σ2
∼ χ2 with n degrees of freedom
1 1 n
pn (y) = n
n
e− 2 y y 2 −1
2 Γ2
2
1 1 n
ypn (y) = n
n
e− 2 y y 2 +1−1
22 Γ2
n 1 n+2
npn+2 (y) = n+2 e− 2 y y 2 −1
n+2
2 2 Γ 2
n − 12 y n
= n
n
e y 2 +1−1
2 Γ2
2
= ypn (y)
Z c1 Z ∞
From (4.7) pn (y)dy + pn (y)dy = α
0 c2
Z c2
⇒ pn (y)dy = 1 − α (4.9)
c1
Z c1 Z ∞ Z ∞
From (4.8) ypn (y)dy + ypn (y)dy = α ypn (y)dy
0 c 0
Z 1∞
Since ypn (y)dy = n
0
Z c1 Z ∞
⇒ npn+2 (y)dy + npn+2 (y)dy = αn
0 c2
Z c2
.. . pn+2 (y)dy = 1−α (4.10)
c1
Solve (4.9) and (4.10) for c1 and c2 the UMPU level α can
be obtained
two. Thus
Z
∂
βφ0 (θ) = φ(x)
pθ (x)dx
∂θ
∂2
Z
βφ00 (θ) = φ(x) 2 pθ (x)dx
∂θ
∂ 1 ∂
log pθ (x) = pθ (x)
∂θ pθ (x) ∂θ
2
∂2 ∂2
1 1 ∂
log pθ (x) = pθ (x) − pθ (x)
∂θ2 ∂θ2 pθ (x) p2θ (x) ∂θ
2
∂2 ∂2
1 ∂ 1
2
pθ (x) = 2
log pθ (x) + pθ (x)
∂θ pθ (x) ∂θ ∂θ [pθ (x)]2
2
∂2 ∂
1 if ∂θ 2 log p θ (x) + ∂θ log p θ (x)
∂
> k1 + k2 ∂θ log pθ (x)
γ if ∂ log p (x) + ∂ log p (x)2
2
∂θ2 θ ∂θ θ
i.e., φ0 (x) =
∂
= k1 + k2 ∂θ log pθ (x)
∂ 2 ∂
2
0 if ∂θ 2 log pθ (x) + ∂θ log pθ (x)
∂
< k1 + k2 ∂θ log pθ (x)
Problems
4.12 Obtain the UMP test, if one such exists for testing
1 1
H0 : θ ≤ 4 vs H1 : θ > 4 using a random sample from
the uniform distribution over (0, θ)
Applications of GNPL 145
4.14 Derive the UMPU level α test for testing the hypoth-
esis equality of variances of two independent normal
distributions
4.17 Obtain the LMPU test for testing the statistical hy-
pothesis H0 : θ = 0 vs H1 : θ 6= 0 using a random
sample from N (0, 1)
5.1 Introduction
1 3
To test H0 : 4 ≤θ≤ 4 the critical function φ suggested is
1 1 3
8 if 4 ≤x≤ 4
φ(x) =
0 otherwise
α = sup Eθ [φ(X)]
θ∈ΩH0
1
= sup Pθ {X ∈ C | θ}
8 θ∈ΩH0
1 1 3 3
= PH ≤X≤ |θ=
8 0 4 4 4
Z 3
1 4 1
= xdx =
3 1 12
4
Now ω = {0}, since the closure of Ω̄H0 = {0} and Ω̄H1 = {0}
1 − 12 (2y−2µ)+ θ1
pθ (y, z) = e θ
θ2
1 if z < − loge 100 or z > loge 100
i.e., φ(z) =
0 otherwise
The power function of the test φ is
β(θ) = Eθ [φ(Z)] ∀ θ ≥ 1
Z − loge 100 Z ∞
1 z 1 −z
= e θ dz + e θ dz
−∞ 2 loge 100 2
1 − 1 loge 100 1 − 1 loge 100
= e θ + e θ
2 2
1
1
1
log( 100 )θ 1 θ
= e =
100
1
Here ω = Ω̄H0 ∩ Ω̄H1 = { 1 } and β(ω) = 100 = 0.01 = α the
size of the test function φ. .. . The test φ is similar on ω
Theorem 5.1 Suppose φ is a level α unbiased for testing
H0 : θ ∈ ΩH0 vs H1 : θ ∈ ΩH1 . If the power function βφ (θ) be
continuous function of θ , then the test function φ is similar
on ω
Proof Suppose βφ (θ) is a continuous function in θ
i.e., limn→∞ βφ (θn ) = θ exists finitely
⇒ θ is a limit point
Also φ is a level α unbiased test
i.e., EH0 [φ(X)] ≤ α ∀ θ ∈ ΩH0 and βφ (θ) ≥ α ∀ θ ∈ ΩH1
When θn0 ∈ ΩH0 , consider a sequence {θn0 , n = 1, 2, 3 · · · }
⇒ limn→∞ θn0 = θ
i.e., limn→∞ βφ (θn0 ) = βφ (limn→∞ θn0 )
= βφ (θ) ≤ α ∀ θ ∈ ΩH0 (5.3)
Again, consider a sequence {θn1 , n = 1, 2, · · · } when θn1 ∈ ΩH1
⇒ limn→∞ θn1 = θ
i.e., limn→∞ βφ (θn1 ) = βφ (limn→∞ θn1 )
= βφ (θ) ≥ α ∀ θ ∈ ΩH1 (5.4)
Neyman - Structure 153
Eθ [φ(X)] = α ∀ θ ∈ ω (5.5)
Eθ [E[φ(X)] | T = t] = Eθ (α) ∀ θ ∈ ω
= E[φ(X) − α] = α − α = 0, ∀ P T ∈ P T
M
Again φ(t) ≥ −c +α
M
≥ −c min(α, 1 − α) + α
0 if α ≤ 0.5
Now φ(t) =
≥ 0 if α > 0.5
.. . φ is a test function
and E[φ(T ] = cE[g(T )] + α = α ∀ P T ∈ P T
Thus φ is a similar test
This test φ has not the Neyman - Structure
since P {g(T ) 6= 0} > 0 for some P T ∈ P T
⇒ P {φ(T ) 6= α} > 0 for some P T ∈ P T
So the contradiction arises. To avoid the contradiction of the
assumption,P T must be boundedly complete. Hence φ has
the Neyman - Structure
Problems
1 1
pθ (x) = e− θ (x−µ) , µ < x < ∞, θ > 0
θ
5.6 Define (a) UMP unbiased test (b) UMP similar test
When will they coincide?
is similar on ω = {0}
(d) The test function φ is depending on the sufficient
statistic T = X̄ Ans: (a), (b), (c) and (d)
X - Sample space
Ω - Parameter space
H - Hypotheses
H0 - Null hypothesis
H1 - Alternative hypothesis
X - Random variable
x - Random observation
θ - Parameter
ω - Boundary point
P- Family of distributions
Pθ - Distribution with θ
Pθ {X = x} - Probability of x of X with θ
PH {X = x} - Probability of x of X under H
A.Santhakumaran 162
T (X) - Function in X
t(x) - Function in x
φ - Test function
ψ - Test function
MP - Most Powerful
LR - Likelihood Ratio
11. Nelson, W., Accelerated Testing Statistical Models, Tests Plan and Data Analysis,
John Wiley and Sons, Inc., 2004.
12. Nelson,W., Accelerated life testing step stress models and data analysis, IEEE
Trans. Reliability , 29, 103 - 108, 1980.
13. Neyman,J. and E.S. Pearson, On the problem of the most efficient tests of statis-
tical hypotheses, Phil. Trans.Roy. Soc., A, 231, 289 - 337,1933.
14. Rao,C.R., Linear Statistical Inference and Its Applications, Wiley and Sons, 1984
A.Santhakumaran 164
16. Santhakumaran,A.,Decision making tools in real life problem, The Hindu, 27th
Oct. 1997.
18. Zacks,S., Theory of Statistical Inference, John Wiley and Sons, New York, 1971.
SUBJECT INDEX