order statistics - university of torontofisher.utstat.toronto.edu/~hadas/sta257/lecture...

26
week 10 1 Order Statistics The order statistics of a set of random variables X 1 , X 2 ,…, X n are the same random variables arranged in increasing order. Denote by X (1) = smallest of X 1 , X 2 ,…, X n X (2) = 2 nd smallest of X 1 , X 2 ,…, X n X (n) = largest of X 1 , X 2 ,…, X n Note, even if X i ’s are independent, X (i) ’s can not be independent since X (1) X (2) X (n) Distribution of X i ’s and X (i) ’s are NOT the same. M

Upload: others

Post on 01-May-2020

2 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 1

Order Statistics• The order statistics of a set of random variables X1, X2,…, Xn are the same

random variables arranged in increasing order.

• Denote by X(1) = smallest of X1, X2,…, Xn

X(2) = 2nd smallest of X1, X2,…, Xn

X(n) = largest of X1, X2,…, Xn

• Note, even if Xi’s are independent, X(i)’s can not be independent sinceX(1) ≤ X(2) ≤ … ≤ X(n)

• Distribution of Xi’s and X(i)’s are NOT the same.

M

Page 2: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 2

Distribution of the Largest order statistic X(n)

• Suppose X1, X2,…, Xn are i.i.d random variables with common distribution function FX(x) and common density function fX(x).

• The CDF of the largest order statistic, X(n), is given by

• The density function of X(n) is then

( )( ) ( )( ) =≤= xXPxF nX n

( )( )

( )( ) == xF

dxdxf

nn XX

Page 3: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 3

Example

• Suppose X1, X2,…, Xn are i.i.d Uniform(0,1) random variables. Find the density function of X(n).

Page 4: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 4

Distribution of the Smallest order statistic X(1)

• Suppose X1, X2,…, Xn are i.i.d random variables with common distribution function FX(x) and common density function fX(x).

• The CDF of the smallest order statistic X(1) is given by

• The density function of X(1) is then

( )( ) ( )( ) ( )( ) =>−=≤= xXPxXPxFX 11 1

1

( )( )

( )( ) == xF

dxdxf XX 11

Page 5: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 5

Example

• Suppose X1, X2,…, Xn are i.i.d Uniform(0,1) random variables. Find the density function of X(1).

Page 6: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 6

Distribution of the kth order statistic X(k)

• Suppose X1, X2,…, Xn are i.i.d random variables with common distribution function FX(x) and common density function fX(x).

• The density function of X(k) is

( )( ) ( ) ( ) ( )( ) ( )( ) ( )xfxFxF

knknxf X

knX

kXX n

−− −−−

= 1!!1

! 1

Page 7: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 7

Example

• Suppose X1, X2,…, Xn are i.i.d Uniform(0,1) random variables. Find the density function of X(k).

Page 8: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 8

Some facts about Power Series∞

• Consider the power series with non-negative coefficients ak.

• If converges for any positive value of t, say for t = r, then it converges for all t in the interval [-r, r] and thus defines a function of t on that interval.

• For any t in (-r, r), this function is differentiable at t and the series converges to the derivatives.

• Example: For k = 0, 1, 2,… and -1< x < 1 we have that

(differentiating geometric series).

∑∞

=0k

kk ta

∑=0k

kkta

∑∞

=

0

1

k

kktka

( ) ∑∞

=

−−⎟⎟⎠

⎞⎜⎜⎝

⎛ +=−

0

11m

mk xk

mkx

Page 9: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 9

Generating Functions

• For a sequence of real numbers {aj} = a0, a1, a2 ,…, the generating function of {aj} is

if this converges for |t| < t0 for some t0 > 0.

( ) ∑∞

=

=0j

jjtatA

Page 10: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 10

Probability Generating Functions• Suppose X is a random variable taking the values 0, 1, 2, … (or a subset of

the non-negative integers).

• Let pj = P(X = j) , j = 0, 1, 2, …. This is in fact a sequence p0, p1, p2, …

• Definition: The probability generating function of X is

• Since if |t| < 1 and the pgf converges absolutely at leastfor |t| < 1.

• In general, πX(1) = p0 + p1 + p2 +… = 1.

• The pgf of X is expressible as an expectation:

( ) ∑∞

=

=+++=0

2210

j

jjX tptptppt Lπ

jj

j ptp ≤ ∑∞

=

=0

1j

jp

( ) ( )X

j

jjX tEtpt == ∑

=0π

Page 11: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 11

Examples• X ~ Binomial(n, p),

converges for all real t.

• X ~ Geometric(p),

converges for |qt| < 1 i.e.

Note: in this case pj = pqj for j = 1, 2, …

( ) ( )nn

j

jjnjX qpttqp

jn

t +=⎟⎟⎠

⎞⎜⎜⎝

⎛= ∑

=

( )qt

pttpqtj

jjX −

== ∑∞

=

11

pqt

−=<

111

Page 12: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 12

PGF for sums of independent random variables

• If X, Y are independent and Z = X+Y then,

• ExampleLet Y ~ Binomial(n, p). Then we can write Y = X1+X2+…+ Xn . Where Xi’sare i.i.d Bernoulli(p). The pgf of Xi is

The pgf of Y is then

( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )tttEtEttEtEtEt YXYXYXYXZ

Z πππ ===== +

( ) ( ) .1 10 qtpptpttiX +=+−=π

( ) ( ) ( ) ( ) ( ) ( ) .2121 nXXXXXXY qtptEtEtEtEt nn +=== +++ LLπ

Page 13: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 13

Use of PGF to find probabilities• Theorem

Let X be a discrete random variable, whose possible values are the nonnegative integers. Assume πX(t0) < ∞ for some t0 > 0. Then

πX(0) = P(X = 0),

etc. In general,

where is the kth derivative of πX with respect to t.• Proof:

( ) ( ),220'' == XPXπ

( )( ) ( ),!0 kXPkkX ==π

( )kXπ

( ) ( ),10' == XPXπ

Page 14: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 14

Example

• Suppose X ~ Poisson(λ). The pgf of X is given by

• Using this pgf we have that

( ) == ∑∞

=

0 !j

jj

X tj

et λπλ

Page 15: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 15

Finding Moments from PGFs• Theorem

Let X be a discrete random variable, whose possible values are the nonnegative integers. If πX(t) < ∞ for |t| < t0 for some t0 > 1. Then

etc. In general,

Where is the kth derivative of πX with respect to t.

• Note: E(X(X-1)···(X-k+1)) is called the kth factorial moment of X.• Proof:

( ) ( ),1' XEX =π

( ) ( )( ),11'' −= XXEXπ

( )( ) ( )( ) ( )( ),1211 +−−−= KXXXXEkX Lπ

( )kXπ

Page 16: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 16

Example• Suppose X ~ Binomial(n, p). The pgf of X is

πX(t) = (pt+q)n.

Find the mean and the variance of X using its pgf.

Page 17: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 17

Uniqueness Theorem for PGF

• Suppose X, Y have probability generating function πX and πY respectively. Then πX(t) = πY(t) if and only if P(X = k) = P(Y = k) for k = 0,1,2,…

• Proof:Follow immediately from calculus theorem:If a function is expressible as a power series at x=a, then there is only one such series.A pgf is a power series about the origin which we know exists with radius of convergence of at least 1.

Page 18: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 18

Moment Generating Functions

• The moment generating function of a random variable X is

mX(t) exists if mX(t) < ∞ for |t| < t0 >0

• If X is discrete

• If X is continuous

• Note: mX(t) = πX(et).

( ) ( )tXX eEtm =

( ) ( ).∑=x

Xtx

X xpetm

( ) ( ) .dxxfetm Xtx

X ∫∞

∞−=

Page 19: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 19

Examples

• X ~ Exponential(λ). The mgf of X is

• X ~ Uniform(0,1). The mgf of X is

( ) ( ) ∫∞ − ===0

dxeeeEtm xtxtXX

λλ

( ) ( ) ∫ ===1

0dxeeEtm txtX

X

Page 20: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 20

Generating Moments from MGFs• Theorem

Let X be any random variable. If mX(t) < ∞ for |t| < t0 for some t0 > 0. Then

mX(0) = 1

etc. In general,

Where is the kth derivative of mX with respect to t. • Proof:

( ) ( ),0' XEm X =

( ) ( ),0'' 2XEm X =

( ) ( ) ( ),0 kkX XEm =

( )kXm

Page 21: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 21

Example

• Suppose X ~ Exponential(λ). Find the mean and variance of X using its moment generating function.

Page 22: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 22

Example

• Suppose X ~ N(0,1). Find the mean and variance of X using its moment generating function.

Page 23: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 23

Example

• Suppose X ~ Binomial(n, p). Find the mean and variance of X using its moment generating function.

Page 24: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 24

Properties of Moment Generating Functions

• mX(0) = 1.

• If Y=a+bX, then the mgf of Y is given by

• If X,Y independent and Z = X+Y then,

( ) ( ) ( ) ( ) ( ).btmeeEeeEeEtm XatbtXatbtXattY

Y ==== +

( ) ( ) ( ) ( ) ( ) ( ) ( ) ( )tmtmeEeEeeEeEeEtm YXtYtXtYtXtYtXtZ

Z ===== +

Rba ∈,

Page 25: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 25

Uniqueness Theorem

• If a moment generating function mX(t) exists for t in an open interval containing 0, it uniquely determines the probability distribution.

Page 26: Order Statistics - University of Torontofisher.utstat.toronto.edu/~hadas/STA257/Lecture notes/week10.pdf · week 10 8 Some facts about Power Series ∞ • Consider the power series

week 10 26

Example• Find the mgf of X ~ N(μ,σ2) using the mgf of the standard normal random

variable.

• Suppose, , independent.

Find the distribution of X1+X2 using mgf approach.

( )2111 ,~ σμNX ( )2

222 ,~ σμNX