<< Chapter < Page Chapter >> Page >

A function with a compound definition

The pair { X , Y } has joint density f X Y ( t , u ) = 1 / 2 on the square region bounded by u = 1 + t , u = 1 - t , u = 3 - t , and u = t - 1 (see [link] ),

W = X for max { X , Y } 1 2 Y for max { X , Y } > 1 = I Q ( X , Y ) X + I Q c ( X , Y ) 2 Y

where Q = { ( t , u ) : max { t , u } 1 } = { ( t , u ) : t 1 , u 1 } . Determine E [ W ] .

Figure 2 is a density drawing, with a horizontal axis labeled as t, and a vertical axis labeled u. The drawing is a shaded square rotated 45 degrees to be sitting with one point on the horizontal axis.  The point sits on (1, 0) and a second point sits against the vertical axis, at (0, 1). In looking at the drawing it can be deduced that the third vertex is at (1, 2), and that the fourth vertex is at (2, 1). Each side of the square is labeled with an equation. Starting with the side between the vertices that are sitting on the axes, an reading them clockwise, the equations are listed as u= 1 - t, u = 1 + t, u= 3 - t, and u = t - 1. There is also an equation inside the shaded square, reading f_xy (t, u) = 1/2. Figure 2 is a density drawing, with a horizontal axis labeled as t, and a vertical axis labeled u. The drawing is a shaded square rotated 45 degrees to be sitting with one point on the horizontal axis.  The point sits on (1, 0) and a second point sits against the vertical axis, at (0, 1). In looking at the drawing it can be deduced that the third vertex is at (1, 2), and that the fourth vertex is at (2, 1). Each side of the square is labeled with an equation. Starting with the side between the vertices that are sitting on the axes, an reading them clockwise, the equations are listed as u= 1 - t, u = 1 + t, u= 3 - t, and u = t - 1. There is also an equation inside the shaded square, reading f_xy (t, u) = 1/2.
The density for [link] .

ANALYTIC SOLUTION

The intersection of the region Q and the square is the set for which 0 t 1 and 1 - t u 1 . Reference to the figure shows three regions of integration.

E [ W ] = 1 2 0 1 1 - t 1 t d u d t + 1 2 0 1 1 1 + t 2 u d u d t + 1 2 1 2 t - 1 3 - t 2 u d u d t = 11 / 6 1 . 8333

APPROXIMATION

tuappr Enter matrix [a b]of X-range endpoints [0 2] Enter matrix [c d]of Y-range endpoints [0 2] Enter number of X approximation points 200Enter number of Y approximation points 200 Enter expression for joint density ((u<=min(t+1,3-t))&... (u>=max(1-t,t-1)))/2 Use array operations on X, Y, PX, PY, t, u, and PM = max(t,u)<=1; G = t.*M + 2*u.*(1 - M); % Z = g(X,Y)EG = total(G.*P) % E[g(X,Y)] EG = 1.8340 % Theoretical 11/6 = 1.8333[Z,PZ] = csort(G,P); % Distribution for ZEZ = dot(Z,PZ) % E[Z] from distributionEZ = 1.8340
Got questions? Get instant answers now!

Special forms for expectation

The various special forms related to property (E20a) are often useful. The general result, which we do not need, is usually derived by an argument which employs a generalform of what is known as Fubini's theorem. The special form (E20b)

E [ X ] = - [ u ( t ) - F X ( t ) ] d t

may be derived from (E20a) by use of integration by parts for Stieltjes integrals. However, we use the relationship between the graph of the distribution function and the graph ofthe quantile function to show the equivalence of (E20b) and (E20f) . The latter property is readily established by elementary arguments.

The property (E20f)

If Q is the quantile function for the distribution function F X , then

E [ g ( X ) ] = 0 1 g [ Q ( u ) ] d u

VERIFICATION

If Y = Q ( U ) , where U uniform on ( 0 , 1 ) , then Y has the same distribution as X . Hence,

E [ g ( X ) ] = E [ g ( Q ( U ) ) ] = g ( Q ( u ) ) f U ( u ) d u = 0 1 g ( Q ( u ) ) d u
Got questions? Get instant answers now!

Reliability and expectation

In reliability, if X is the life duration (time to failure) for a device, the reliability function is the probability at any time t the device is still operative. Thus

R ( t ) = P ( X > t ) = 1 - F X ( t )

According to property (E20b)

E [ X ] = 0 R ( t ) d t
Got questions? Get instant answers now!

Use of the quantile function

Suppose F X ( t ) = t a , a > 0 , 0 t 1 . Then Q ( u ) = u 1 / a , 0 u a .

E [ X ] = 0 1 u 1 / a d u = 1 1 + 1 / a = a a + 1

The same result could be obtained by using f X ( t ) = F X ' ( t ) and evaluating t f X ( t ) d t .

Got questions? Get instant answers now!

Equivalence of (E20b) And (E20f)

For the special case g ( X ) = X , Figure 3(a) shows 0 1 Q ( u ) d u is the difference in the shaded areas

0 1 Q ( u ) d u = Area A - Area B

The corresponding graph of the distribution function F is shown in Figure 3(b) . Because of the construction, the areas of the regions marked A and B are the same in the two figures. As may be seen,

Area A = 0 [ 1 - F ( t ) ] d t and Area B = - 0 F ( t ) d t

Use of the unit step function u ( t ) = 1 for t > 0 and 0 for t < 0 (defined arbitrarily at t = 0 ) enables us to combine the two expressions to get

0 1 Q ( u ) d u = Area A - Area B = - [ u ( t ) - F ( t ) ] d t
Figure three contains two graphs. The first graph has a horizontal axis labeled t, and a vertical axis labeled u. The large label of the graph reads,  u = Q(t). A dashed vertical line along t = 1 bounds an increasing curved plot. The curve starts with a vertical asymptote along the vertical axis below the horizontal axis, and as it approaches the horizontal axis, the slope becomes more shallow. The curve's slope shallows until it is midway in horizontal distance between the vertical axis and the dashed vertical line. At this point, the slope begins to increase again, until it reaches a vertical asymptote along the dashed line at t = 1. The horizontal and vertical axes, along with the curve itself, create a bounded shape. A small right triangle loosely fits this bounded shape, and is labeled as B. The dashed line, horizontal axis, and the segment of the curve above the horizontal axis create a larger bounded shape, and a larger right triangle loosely fits this bounded shape, labeled A. The second graph is roughly similar. The axes are in the same place, but with this figure, s dashed line is now drawn horizontally along u = 1. A curve of the same shape now begins as a horizontal asymptote along the t - axis. It increases in slope at an increasing rate for half of the vertical distance and then decreases in slope back to a horizontal asymptote at u = 1. The same triangles fitting the same bounded regions as in the first figure are used in the second figures, only because of the rotated nature of the new curve, these triangles are rotated in the same fashion. Figure three contains two graphs. The first graph has a horizontal axis labeled t, and a vertical axis labeled u. The large label of the graph reads,  u = Q(t). A dashed vertical line along t = 1 bounds an increasing curved plot. The curve starts with a vertical asymptote along the vertical axis below the horizontal axis, and as it approaches the horizontal axis, the slope becomes more shallow. The curve's slope shallows until it is midway in horizontal distance between the vertical axis and the dashed vertical line. At this point, the slope begins to increase again, until it reaches a vertical asymptote along the dashed line at t = 1. The horizontal and vertical axes, along with the curve itself, create a bounded shape. A small right triangle loosely fits this bounded shape, and is labeled as B. The dashed line, horizontal axis, and the segment of the curve above the horizontal axis create a larger bounded shape, and a larger right triangle loosely fits this bounded shape, labeled A. The second graph is roughly similar. The axes are in the same place, but with this figure, s dashed line is now drawn horizontally along u = 1. A curve of the same shape now begins as a horizontal asymptote along the t - axis. It increases in slope at an increasing rate for half of the vertical distance and then decreases in slope back to a horizontal asymptote at u = 1. The same triangles fitting the same bounded regions as in the first figure are used in the second figures, only because of the rotated nature of the new curve, these triangles are rotated in the same fashion.
Equivalence of properties (E20b) and (E20f) .
Got questions? Get instant answers now!

Property (E20c) is a direct result of linearity and (E20b) , with the unit step functions cancelling out.

Property (E20d) Useful inequalities

Suppose X 0 . Then

n = 0 P ( X n + 1 ) E [ X ] n = 0 P ( X n ) N k = 0 P ( X k N ) , for all N 1

VERIFICATION

For X 0 , by (E20b)

E [ X ] = 0 [ 1 - F ( t ) ] d t = 0 P ( X > t ) d t

Since F can have only a countable number of jumps on any interval and P ( X > t ) and P ( X t ) differ only at jump points, we may assert

a b P ( X > t ) d t = a b P ( X t ) d t

For each nonnegative integer n , let E n = [ n , n + 1 ) . By the countable additivity of expectation

E [ X ] = n = 0 E [ I E n X ] = n = 0 E n P ( X t ) d t

Since P ( X t ) is decreasing with t and each E n has unit length, we have by the mean value theorem

P ( X n + 1 ) E [ I E n X ] P ( X n )

The third inequality follows from the fact that

k N ( k + 1 ) N P ( X t ) d t N E k N P ( X t ) d t N P ( X k N )
Got questions? Get instant answers now!

Remark . Property (E20d) is used primarily for theoretical purposes. The special case (E20e) is more frequently used.

Property (E20e)

If X is nonnegative, integer valued, then

E [ X ] = k = 1 P ( X k ) = k = 0 P ( X > k )

VERIFICATION

The result follows as a special case of (E20d) . For integer valued random variables,

P ( X t ) = P ( X n ) on E n and P ( X t ) = P ( X > n ) = P ( X n + 1 ) on E n + 1
Got questions? Get instant answers now!

An elementary derivation of (E20e) can be constructed as follows.

(E20e) For integer-valued random variables

By definition

E [ X ] = k = 1 k P ( X = k ) = lim n k = 1 n k P ( X = k )

Now for each finite n ,

k = 1 n k P ( X = k ) = k = 1 n j = 1 k P ( X = k ) = j = 1 n k = j n P ( X = k ) = j = 1 n P ( X j )

Taking limits as n yields the desired result.

Got questions? Get instant answers now!

The geometric distribution

Suppose X geometric ( p ) . Then P ( X k ) = q k . Use of (E20e) gives

E [ X ] = k = 1 q k = q k = 0 q k = q 1 - q = q / p
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