<< Chapter < Page Chapter >> Page >

Remark . It is easy to show that the function λ ( · ) is nondecreasing. This fact is used in establishing Jensen's inequality for conditional expectation.

The product rule for expectations of independent random variables

Product rule for simple random variables

Consider an independent pair { X , Y } of simple random variables

X = i = 1 n t i I A i Y = j = 1 m u j I B j (both in canonical form)

We know that each pair { A i , B j } is independent, so that P ( A i B j ) = P ( A i ) P ( B j ) . Consider the product X Y . According to the pattern described after Example 9 from "Mathematical Expectation: Simple Random Variables."

X Y = i = 1 n t i I A i j = 1 m u j I B j = i = 1 n j = 1 m t i u j I A i B j

The latter double sum is a primitive form, so that

E [ X Y ] = i = 1 n j = 1 m t i u j P ( A i B j ) = i = 1 n j = 1 m t i u j P ( A i ) P ( B j ) = i = 1 n t i P ( A i ) j = 1 m u j P ( B j ) = E [ X ] E [ Y ]

Thus the product rule holds for independent simple random variables.

Got questions? Get instant answers now!

Approximating simple functions for an independent pair

Suppose { X , Y } is an independent pair, with an approximating simple pair { X s , Y s } . As functions of X and Y , respectively, the pair { X s , Y s } is independent. According to [link] , above, the product rule E [ X s Y s ] = E [ X s ] E [ Y s ] must hold.

Got questions? Get instant answers now!

Product rule for an independent pair

For X 0 , Y 0 , there exist nondecreasing sequences { X n : 1 n } and { Y n : 1 n } of simple random variables increasing to X and Y , respectively. The sequence { X n Y n : 1 n } is also a nondecreasing sequence of simple random variables, increasing to X Y . By the monotone convergence theorem (MC)

E [ X n ] E [ X ] , E [ Y n ] E [ Y ] , and E [ X n Y n ] E [ X Y ]

Since E [ X n Y n ] = E [ X n ] E [ Y n ] for each n , we conclude E [ X Y ] = E [ X ] E [ Y ]

In the general case,

X Y = ( X + - X - ) ( Y + - Y - ) = X + Y + - X + Y - - X - Y + + X - Y -

Application of the product rule to each nonnegative pair and the use of linearity gives the product rule for the pair { X , Y }

Got questions? Get instant answers now!

Remark . It should be apparent that the product rule can be extended to any finite independent class.

The joint distribution of three random variables

The class { X , Y , Z } is independent, with the marginal distributions shown below. Let W = g ( X , Y , Z ) = 3 X 2 + 2 X Y - 3 X Y Z . Determine E [ W ] .

X = 0:4; Y = 1:2:7;Z = 0:3:12; PX = 0.1*[1 3 2 3 1]; PY = 0.1*[2 2 3 3]; PZ = 0.1*[2 2 1 3 2]; icalc3 % Setup for joint dbn for{X,Y,Z} Enter row matrix of X-values XEnter row matrix of Y-values Y Enter row matrix of Z-values ZEnter X probabilities PX Enter Y probabilities PYEnter Z probabilities PZ Use array operations on matrices X, Y, Z,PX, PY, PZ, t, u, v, and P EX = X*PX'% E[X] EX = 2EX2 = (X.^2)*PX' % E[X^2] EX2 = 5.4000EY = Y*PY' % E[Y] EY = 4.4000EZ = Z*PZ' % E[Z] EZ = 6.3000G = 3*t.^2 + 2*t.*u - 3*t.*u.*v; % W = g(X,Y,Z) = 3X^2 + 2XY - 3XYZ EG = total(G.*P) % E[g(X,Y,Z)]EG = -132.5200 [W,PW]= csort(G,P); % Distribution for W = g(X,Y,Z) EW = W*PW' % E[W]EW = -132.5200 ew = 3*EX2 + 2*EX*EY - 3*EX*EY*EZ % Use of linearity and product ruleew = -132.5200
Got questions? Get instant answers now!

A function with a compound definition: truncated exponential

Suppose X exponential (0.3). Let

Z = X 2 for X 4 16 for X > 4 = I [ 0 , 4 ] ( X ) X 2 + I ( 4 , ] ( X ) 16

Determine E [ Z ] .

ANALYTIC SOLUTION

E [ g ( X ) ] = g ( t ) f X ( t ) d t = 0 I [ 0 , 4 ] ( t ) t 2 0 . 3 e - 0 . 3 t d t + 16 E [ I ( 4 , ] ( X ) ]
= 0 4 t 2 0 . 3 e - 0 . 3 t d t + 16 P ( X > 4 ) 7 . 4972 (by Maple)

APPROXIMATION

To obtain a simple aproximation, we must approximate the exponential by a bounded random variable. Since P ( X > 50 ) = e - 15 3 · 10 - 7 we may safely truncate X at 50.

tappr Enter matrix [a b]of x-range endpoints [0 50] Enter number of x approximation points 1000Enter density as a function of t 0.3*exp(-0.3*t) Use row matrices X and PX as in the simple caseM = X<= 4; G = M.*X.^2 + 16*(1 - M); % g(X)EG = G*PX' % E[g(X)] EG = 7.4972[Z,PZ] = csort(G,PX); % Distribution for Z = g(X)EZ = Z*PZ' % E[Z] from distributionEZ = 7.4972

Because of the large number of approximation points, the results agree quite closely with the theoretical value.

Got questions? Get instant answers now!

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Applied probability. OpenStax CNX. Aug 31, 2009 Download for free at http://cnx.org/content/col10708/1.6
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Applied probability' conversation and receive update notifications?

Ask