を示す構成例


12

E 1の確率分布の例を構築する方法X=1E(1X)=1E(X)PX0=1P(X0)=1と仮定して、 EX が成り立つ?

正の値RV用ジェンセンの不等式から以下の不平等XはX似ているE 1X1EX E(1X)1E(X)X<0のX<0場合、逆不等式)。マッピングがためですxは1xx1xは、x>0のx>0場合は凸で、x<0の場合はx<0凹です。ジェンセンの不等式の等号条件に従って、必要な等式が成立するためには、分布を縮退させる必要があると思います。等式が成り立つ些細なケースは、もちろんX=1X=1aeの場合です。問題の本で見つけた例は次のとおりです。PX=1=1のような離散確率変数XX9PX=12=PX=2=49P(X=1)=19,P(X=12)=P(X=2)=49。その後、E1X=1EX =1E(1X)=1E(X)=1

この例は、タイトルの平等が保持されるために、XXが正(または負)aeである必要がないことを示しています。ここでの分布も縮退していません。

本で見つけたようなサンプルを作成するにはどうすればよいですか?動機はありますか?


この例は、ゼロ以外の定数であるランダム変数を保持します。また、2番目の例は縮退ではありません。
マイケルR.チャーニック

3
不等式は、Xがほぼ確実に正であるとさらに仮定しない限り、ジェンセンの不等式から続きませX
whuber

@MichaelChernick私は、この例が退化した分布を持っていると言うつもりはありませんでした。
StubbornAtom

「ジェンセンの不等式の平等条件に従って、必要な平等が成り立つためには、分布を縮退させる必要があると思います」と述べました。しかし、あなたは非縮退の例を示しました。
マイケルR.チャーニック

1
@whuberタイトルの平等が真である例を見つける方法を知りたいだけです。
StubbornAtom

回答:


18

レッツ・構文すべての可能な例のランダム変数のXXのためのE [ X ] E [ 1 / X ] = 1E[X]E[1/X]=1。次に、それらのうち、いくつかのヒューリスティックに従って、最も単純な可能な例を取得します。 これらのヒューリスティックは、予備分析から脱落するすべての式に可能な限り単純な値を与えることで構成されています。 これが教科書の例であることがわかりました。

予備分析

これには、定義に基づいたほんの少しの分析が必要です。解決策は二次的な関心のみです。主な目的は、直感的に結果を理解するのに役立つ洞察を開発することです。

まず、Jensenの不等式(またはコーシー・シュワルツの不等式)はその正の確率変数のための暗示ことを観察XXE [ X ] E [ 1 / X ] 1E[X]E[1/X]1と場合にだけ、平等保持とXがX、それは次のとおりです。「縮退」であります、XXはほぼ確実に一定です。場合XはX負の確率変数であり、- XはX正であり、前述の結果は、不等号を逆にして成り立ちます。したがって、任意の実施例E [ 1 / X ] = 1 / E[ X ]E[1/X]=1/E[X]は、負になる正の確率と正になる正の確率を持つ必要があります。

ここでの洞察は、任意の点であるXXE [ X ] Eは、[ 1 / X ] = 1はE[X]E[1/X]=1何とかその負の部分から他の方向に不平等に対するその正の部分から不等式「バランス」しなければなりません。これは、進むにつれて明らかになります。

ゼロ以外のランダム変数Xを考えます。期待値の定義を策定する最初のステップ(少なくともこれがメジャー理論を使用して完全に一般化されている場合)は、Xを正と負の部分に分解することです。どちらも正のランダム変数です。XX

Y = 正の部分X = max 0 X ; Z = 負の部分X = min 0 X

YZ=Positive part(X)=max(0,X);=Negative part(X)=min(0,X).

レッツ考えるXとして混合物Y重量とP- Z量の1 - P P = のPr X > 0 1 - P = のPr X < 0 明らかに0 < p < 1。これにより、X1 / Xの期待値を記述できます。XYpZ1p 

p=Pr(X>0), 1p=Pr(X<0).
0<p<1.
X1/X正の変数YおよびZの期待値に関して。YZ

均一に再スケーリングすることを今後の代数少し、ノート簡単にするためにXを数によってσは変化しないEを[ X ] E [ 1 / X ] -しかしそれは多重しE [ Y ]E [ Z ]によってそれぞれσ。正のσの場合、これは単にXの測定単位を選択することになります。負のσYZの役割を切り替えます。σの符号の選択XσE[X]E[1/X]E[Y]E[Z]σσXσYZσ適切したがって、我々は仮定することができるE [ Z ] = 1  と  E [ Y ] E [ Zを]

E[Z]=1 and E[Y]E[Z].(1)

表記法

予備的な簡略化は以上です。いい表記を作成するために、次のように書きましょう。

μ = E [ Y ] ; ν = E [ 1 / Y ] ; λ = E [ 1 / Z ]  

μ=E[Y]; ν=E[1/Y]; λ=E[1/Z]

私たちがコントロールできない3つの期待のために。3つの量はすべて正です。ジェンセンの不平等は断言します

μ ν 1  と  λ 1。

μν1 and λ1.(2)

総確率の法則は、Xおよび1 / Xの期待値を、我々が命名した量の観点から表現しています。X1/X

E [ X ] = E [ X | X > 0 ] のPr X > 0 + E [ X | X < 0 ] のPr X < 0 = μ P - 1 - P = μ + 1 P 1

E[X]=E[XX>0]Pr(X>0)+E[XX<0]Pr(X<0)=μp(1p)=(μ+1)p1

そして、以降1 / Xは同じ符号有するXを1/XX

E [ 1X ]=E[1XX>0]PrX>0+E[1X |X<0]のPrX<0=νP-λ1-P=ν+λP-λ

E[1X]=E[1XX>0]Pr(X>0)+E[1XX<0]Pr(X<0)=νpλ(1p)=(ν+λ)pλ.

これら2つの式の積を1と等しくすると、変数間の本質的な関係が得られます。1

1=E[X]E[1X]=((μ+1)p1)((ν+λ)pλ).

1=E[X]E[1X]=((μ+1)p1)((ν+λ)pλ).(*)

Reformulation of the Problem

Suppose the parts of XX--YY and ZZ--are any positive random variables (degenerate or not). That determines μ,ν,μ,ν, and λλ. When can we find pp, with 0<p<10<p<1, for which ()() holds?

This clearly articulates the "balancing" insight previously stated only vaguely: we are going to hold YY and ZZ fixed and hope to find a value of pp that appropriately balances their relative contributions to XX. Although it's not immediately evident that such a pp need exist, what is clear is that it depends only on the moments E[Y]E[Y], E[1/Y]E[1/Y], E[Z]E[Z], and E[1/Z]E[1/Z]. The problem thereby is reduced to relatively simple algebra--all the analysis of random variables has been completed.

Solution

This algebraic problem isn't too hard to solve, because ()() is at worst a quadratic equation for pp and the governing inequalities (1)(1) and (2)(2) are relatively simple. Indeed, ()() tells us the product of its roots p1p1 and p2p2 is

p1p2=(λ1)1(μ+1)(ν+λ)0

p1p2=(λ1)1(μ+1)(ν+λ)0

and the sum is

p1+p2=(2λ+λμ+ν)1(μ+1)(ν+λ)>0.

p1+p2=(2λ+λμ+ν)1(μ+1)(ν+λ)>0.

Therefore both roots must be positive. Furthermore, their average is less than 11, because

1(p1+p2)2=λμ+ν+2μν2(μ+1)(ν+λ)>0.

1(p1+p2)2=λμ+ν+2μν2(μ+1)(ν+λ)>0.

(By doing a bit of algebra, it's not hard to show the larger of the two roots does not exceed 11, either.)

A Theorem

Here is what we have found:

Given any two positive random variables YY and ZZ (at least one of which is nondegenerate) for which E[Y]E[Y], E[1/Y]E[1/Y], E[Z]E[Z], and E[1/Z]E[1/Z] exist and are finite. Then there exist either one or two values pp, with 0<p<10<p<1, that determine a mixture variable XX with weight pp for YY and weight 1p1p for ZZ and for which E[X]E[1/X]=1E[X]E[1/X]=1. Every such instance of a random variable XX with E[X]E[1/X]=1E[X]E[1/X]=1 is of this form.

That gives us a rich set of examples indeed!


Constructing the Simplest Possible Example

Having characterized all examples, let's proceed to construct one that is as simple as possible.

  • For the negative part ZZ, let's choose a degenerate variable--the very simplest kind of random variable. It will be scaled to make its value 11, whence λ=1λ=1. The solution of ()() includes p1=0p1=0, reducing it to an easily solved linear equation: the only positive root is

    p=11+μ+11+ν.

    p=11+μ+11+ν.(3)
  • For the positive part YY, we obtain nothing useful if YY is degenerate, so let's give it some probability at just two distinct positive values a<ba<b, say Pr(X=b)=qPr(X=b)=q. In this case the definition of expectation gives

    μ=E[Y]=(1q)a+qb; ν=E[1/Y]=(1q)/a+q/b.

    μ=E[Y]=(1q)a+qb; ν=E[1/Y]=(1q)/a+q/b.
  • To make this even simpler, let's make YY and 1/Y1/Y identical: this forces q=1q=1/2q=1q=1/2 and a=1/ba=1/b. Now

    μ=ν=b+1/b2.

    μ=ν=b+1/b2.

    The solution (3)(3) simplifies to

    p=21+μ=42+b+1/b.

    p=21+μ=42+b+1/b.
  • How can we make this involve simple numbers? Since a<ba<b and ab=1ab=1, necessarily b>1b>1. Let's choose the simplest number greater than 11 for bb; namely, b=2b=2. The foregoing formula yields p=4/(2+2+1/2)=8/9p=4/(2+2+1/2)=8/9 and our candidate for the simplest possible example therefore is

    Pr(X=2)=Pr(X=b)=Pr(Y=b)p=qp=1289=49;Pr(X=1/2)=Pr(X=a)=Pr(Y=a)p=qp==49;Pr(X=1)=Pr(Z=1)(1p)=1p=19.

    Pr(X=2)=Pr(X=b)=Pr(Y=b)p=qp=1289=49;Pr(X=1/2)=Pr(X=a)=Pr(Y=a)p=qp==49;Pr(X=1)=Pr(Z=1)(1p)=1p=19.

This is the very example offered in the textbook.


2
Nice answer. Despite my initial scepticism it's easy to find an example with distinct solutions p(0,1)p(0,1).
P.Windridge

8

As you've mentioned, if XX is positive then E(1/X)=1/E(X)E(1/X)=1/E(X) occurs only when XX is almost surely constant. Otherwise you need XX to take both negative and positive values.

To construct such an example, first go as simple as possible. Assume X takes two values, a and b, with probabilities p and 1p respectively. Then E(X)=ap+b(1p)

and E(1/X)=1ap+1b(1p).
To have 1/E(X)=E(1/X) we require ap+b(1p)=11ap+1b(1p)
which rearranges to the requirement (ab)2p(1p)=0.
This means the only possible solution must have either a=b, or p=0, or p=1. In all cases we return to the degenerate case: X is constant.

Next try: a distribution with three possible values. Here there are many more choices. The example you cited tries an X such that 1/X has the same distribution. If we know X takes three values, it must be that one of the values is either 1 or 1, and the other two must be a and 1/a for some choice of a. For definiteness let's try P(X=a)=P(X=1/a)=p, and P(X=1)=12p. Then E(1/X)=E(X)=(a+1a)p(12p)=(2+a+1a)p1.

To meet the requirement 1/E(X)=E(1/X) we demand E(X)=1 or E(X)=1. Expression (1) is never 1 unless p=0, which returns us to the degenerate case again. So aim for E(X)=1, which gives (2+a+1a)p=2p=22+a+1a=2a(a+1)2.
Expression (2) gives an entire family of solutions that meet the requirement. The only constraint is that a must be positive. The example you cited takes a=2. Only the case a=1 is degenerate.

1
I assume your first line is intended to be "if X is positive then E[1/X]=1/E[X] occurs only when X is almost surely constant", like in the question. This drops out from (the proof of) Jensen's inequality, where we also use the fact that h(x)=1/x is not linear.
P.Windridge

@P.Windridge You are right! Fixed.
grand_chat
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.