Wie berechnet man die Erwartung von

12

Wenn Xi exponentiell verteilt ist (i=1,...,n) mit dem Parameter λ und Xi ‚s ist voneinander unabhängig, was die Erwartung

(i=1nXi)2

in Bezug auf n undλ und möglicherweise andere Konstanten?

Hinweis: Diese Frage wurde unter /math//q/12068/4051 mathematisch beantwortet . Die Leser würden es sich auch ansehen.

Isaac
quelle
5
Die beiden Kopien dieser Frage verweisen aufeinander und entsprechend hat die Statistikseite (hier) eine statistische Antwort und die Mathematikseite eine mathematische Antwort. Es scheint eine gute Einteilung zu sein: Lass es stehen!
Whuber

Antworten:

31

Wenn , dann (unter Unabhängigkeit), y = Σ x i ~ G a m m a ( n , 1 / λ ) , so dass y ist gamma verteilt (siehe wikipedia ). Also brauchen wir nur E [ y 2 ] . Da V a r [ y ] = E [ y 2 ] - ExiExp(λ)y=xiGamma(n,1/λ)yE[y2] , wissen wir, dass E [ y 2 ] = V a r [ Y ] + E [ y ] 2 . Daher ist E [ y 2 ] = n / λ 2 + n 2 / λ 2 = n ( 1 + n ) / λ 2 (sieheWikipediafür die Erwartung und Varianz der Gammaverteilung).Var[y]=E[y2]E[y]2E[y2]=Var[y]+E[y]2E[y2]=n/λ2+n2/λ2=n(1+n)/λ2

Wolfgang
quelle
Vielen Dank. Eine sehr gute Möglichkeit zur Beantwortung der Frage (die zur gleichen Antwort führte) wurde vor einigen Minuten auch auf math.stackexchange (Link oben in der Frage) bereitgestellt.
Wolfgang,
2
Die mathematische Antwort berechnet die Integrale unter Verwendung der Linearität der Erwartung. In gewisser Hinsicht ist es einfacher. Aber ich mag Ihre Lösung, weil sie statistisches Wissen ausnutzt : Weil Sie wissen, dass eine Summe unabhängiger Exponentialvariablen eine Gamma-Verteilung hat, sind Sie fertig.
Whuber
1
Ich habe es sehr genossen und bin keineswegs Statistiker oder Mathematiker.
Kortuk
sehr elegante Antwort.
Cyrus S
1
@Dilip Der Mathematiker neigt dazu, diese Frage als eine Frage nach einem Integral zu betrachten und sie direkt zu integrieren. Der Statistiker drückt es in bekannten statistischen Größen wie der Varianz und bekannten statistischen Beziehungen wie dem Exponential Gamma um und die Gamma-Familie wird unter Faltung geschlossen. Die Antworten sind die gleichen, aber die Ansätze sind völlig unterschiedlich. Dann stellt sich die Frage, was "eine Integration machen" wirklich bedeutet. Beispielsweise erfolgt dieses komplizierte Integral rein algebraisch.
Whuber
9

Die obige Antwort ist sehr schön und beantwortet die Frage vollständig, aber ich werde stattdessen eine allgemeine Formel für das erwartete Quadrat einer Summe bereitstellen und sie auf das hier erwähnte spezifische Beispiel anwenden.

Für jeden Satz von Konstanten ist es eine Tatsachedassa1,...,an

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

this is true by the Distributive property and becomes clear when you consider what you're doing when you calculate (a1+...+an)(a1+...+an) by hand.

Therefore, for a sample of random variables X1,...,Xn, regardless of the distributions,

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

provided that these expectations exist.

In the example from the problem, X1,...,Xn are iid exponential(λ) random variables, which tells us that E(Xi)=1/λ and 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.

Macro
quelle
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 :)

wolfies
quelle