の期待値をどのように計算します


12

場合Xi指数関数的に分布している(i=1,...,n)パラメータとλXi「sは、互いに独立している、の期待何

(i=1nXi)2

面でnλとおそらく他の定数は?

注:この質問は、https : //math.stackexchange.com/q/12068/4051で数学的な答えを得ています。読者もそれを見るでしょう。


5
この質問の2つのコピーは相互に参照しており、適切には、統計サイト(ここ)には統計的な答えがあり、数学サイトには数学的な答えがあります。それは良い分裂のようです:それを立ててください!
whuber

回答:


31

もしは、(独立下)、Y = Σ X IG A M M A nは1 / λ ので、yはガンマ分布する(参照あるウィキペディア)。したがって、必要なのはE [ y 2 ]だけです。以降V R [ Y ] = E [ Y 2 ] - ExiExp(λ)y=xiGamma(n,1/λ)yE[y2]、我々は知っている E [ Y 2 ] = V R [ Y ] + E [ Y ] 2。したがって、 E [ Y 2 ] = N / λ 2 + N 2 / λ 2 = N 1 + N / λ 2(参照ウィキペディアガンマ分布の期待値と分散のために)。Var[y]=E[y2]E[y]2E[y2]=Var[y]+E[y]2E[y2]=n/λ2+n2/λ2=n(1+n)/λ2


ありがとう。数分前にmath.stackexchange(質問の上のリンク)で質問に答える(同じ答えにつながる)非常にきちんとした方法も提供されました。
ヴォルフガング

2
数学の答えは、期待の線形性を使用して積分を計算します。いくつかの点で簡単です。しかし、統計的知識を活用するため、ソリューションが気に入っています。独立した指数変数の合計がガンマ分布を持っていることがわかっているので、完了です。
whuber

1
私はそれをかなり楽しんだし、統計学者でも数学者でもない。
コルトゥク

非常にエレガントな答え。
サイラスS

1
@Dilip数学者は、この質問を積分を求めていると見なす傾向があり、直接積分を進めます。統計学者は、分散などのなじみのある統計量、および指数がガンマであり、たたみ込みの下でガンマ族が閉じているなどのなじみのある統計的関係の観点からそれを再表現します。答えは同じですが、アプローチは完全に異なります。次に、「統合を行う」とはどういう意味かという疑問があります。たとえば、この複雑な積分は純粋に代数的に行われます。
whuber

9

上記の答えは非常に素晴らしく、質問に完全に答えますが、代わりに、予想される合計の二乗の一般式を提供し、ここで言及した特定の例に適用します。

定数のいずれかのセットについて1nはそれが事実でありますa1,...,an

(i=1nai)2=i=1nj=1naiaj

これはtrueです分配的性質、あなたが計算するとき、あなたがやっているかを考えたときに明確になる手で。(a1+...+an)(a1+...+an)

したがって、ランダム変数 …のサンプルに対してX n、分布に関係なく、X1,...,Xn

E([i=1nXi]2)=E(i=1nj=1nXiXj)=i=1nj=1nE(XiXj)

これらの期待が存在する場合。

X1,...,Xnexponential(λ)E(Xi)=1/λvar(Xi)=1/λ2 for each i. By independence, for ij, we have

E(XiXj)=E(Xi)E(Xj)=1λ2

There are n2n of these terms in the sum. When i=j, we have

E(XiXj)=E(Xi2)=var(Xi)+E(Xi)2=2λ2

and there are n of these term in the sum. Therefore, using the formula above,

E(i=1nXi)2=i=1nj=1nE(XiXj)=(n2n)1λ2+n2λ2=n2+nλ2

is your answer.


3

This problem is just a special case of the much more general problem of 'moments of moments' which are usually defined in terms of power sum notation. In particular, in power sum notation:

s1=i=1nXi

Then, irrespective of the distribution, the original poster seeks E[s12] (provided the moments exist). Since the expectations operator is just the 1st Raw Moment, the solution is given in the mathStatica software by:

enter image description here

[ The '___ToRaw' means that we want the solution presented in terms of raw moments of the population (rather than say central moments or cumulants). ]

Finally, if X ~ Exponential(λ) with pdf f(x):

f = Exp[-x/λ]/λ;      domain[f] = {x, 0, ∞} &&  > 0};

then we can replace the moments μi in the general solution sol with the actual values for an Exponential random variable, like so:

enter image description here

All done.


P.S. The reason the other solutions posted here yield an answer with λ2 in the denominator rather than the numerator is, of course, because they are using a different parameterisation of the Exponential distribution. Since the OP didn't state which version he was using, I decided to use the standard distribution theory textbook definition Johnson Kotz et al … just to balance things out :)

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