Comment calculez-vous l'espérance de ?

12

Si est distribué de façon exponentielle (i = 1, ..., n) avec le paramètre \ lambda et que les X_i sont mutuellement indépendants, quelle est l'attente de ( i = 1 , . . . , N ) λ X iXi(i=1,...,n)λXi

(i=1nXi)2

en termes de n et λ et éventuellement d'autres constantes?

Remarque: Cette question a obtenu une réponse mathématique sur /math//q/12068/4051 . Les lecteurs y jetteraient également un coup d'œil.

Isaac
la source
5
Les deux copies de cette question se référencent et, de manière appropriée, le site de statistiques (ici) a une réponse statistique et le site de mathématiques a une réponse mathématique. Cela semble être une bonne division: laissez-le reposer!
whuber

Réponses:

31

Si , alors (sous indépendance), , donc est distribué gamma (voir wikipedia ). Donc, nous avons juste besoin de . Puisque , nous savons que . Par conséquent, (voir wikipedia pour l'espérance et la variance de la distribution gamma).y = x iG a m m a ( n , 1 / λ ) y E [ y 2 ] V a r [ y ] = E [ y 2 ] - E [ y ] 2 E [ y 2 ] = V a r [xiExp(λ)y=xiGamma(n,1/λ)yE[y2]Var[y]=E[y2]E[y]2 E [ y 2 ] = n / λ 2 + n 2 / λ 2 = n ( 1 + n ) / λ 2E[y2]=Var[y]+E[y]2E[y2]=n/λ2+n2/λ2=n(1+n)/λ2

Wolfgang
la source
Merci. Une manière très nette de répondre à la question (conduisant à la même réponse) a également été fournie sur math.stackexchange (lien ci-dessus dans la question) il y a quelques minutes.
Wolfgang
2
La réponse mathématique calcule les intégrales en utilisant la linéarité de l'attente. À certains égards, c'est plus simple. Mais j'aime votre solution car elle exploite les connaissances statistiques : parce que vous savez qu'une somme de variables exponentielles indépendantes a une distribution Gamma, vous avez terminé.
whuber
1
J'ai beaucoup aimé et je ne suis en aucun cas un statisticien ou un mathématicien.
Kortuk
réponse très élégante.
Cyrus S
1
@Dilip Le mathématicien a tendance à voir cette question comme demandant une intégrale et procède directement à son intégration. Le statisticien le ré-exprime en termes de quantités statistiques familières, telles que la variance, et de relations statistiques familières, telles que l'exponentielle est Gamma et la famille Gamma est fermée par convolution. Les réponses sont les mêmes mais les approches sont complètement différentes. Il y a ensuite la question de savoir ce que signifie réellement «faire une intégration». Par exemple, cette intégrale compliquée se fait de manière purement algébrique.
whuber
9

La réponse ci-dessus est très agréable et répond complètement à la question, mais je vais plutôt fournir une formule générale pour le carré attendu d'une somme et l'appliquer à l'exemple spécifique mentionné ici.

Pour tout ensemble de constantes c'est un fait quea1,...,an

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

cela est vrai par la propriété Distributive et devient clair lorsque vous considérez ce que vous faites lorsque vous calculez à la main.(a1+...+an)(a1+...+an)

Par conséquent, pour un échantillon de variables aléatoires , quelles que soient les distributions,X1,...,Xn

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

à condition que ces attentes existent.

Dans l'exemple du problème, sont des variables aléatoires iid , ce qui nous dit que et pour chaque . Par indépendance, pour , nous avonsX1,...,Xnexponential(λ)E(Xi)=1/λvar(Xi)=1/λ2iij

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

Il y a de ces termes dans la somme. Lorsque , nous avonsn2ni=j

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

et il y a de ces termes dans la somme. Par conséquent, en utilisant la formule ci-dessus,n

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

est votre réponse.

Macro
la source
3

Ce problème n'est qu'un cas particulier du problème beaucoup plus général des «moments des moments» qui sont généralement définis en termes de notation de somme de puissance. En particulier, en notation de somme de puissance:

s1=i=1nXi

Ensuite, quelle que soit la distribution , l'affiche originale cherche (à condition que les moments existent). L'opérateur des attentes n'étant que le 1er moment brut, la solution est donnée dans le logiciel mathStatica par:E[s12]

entrez la description de l'image ici

[Le '___ToRaw' signifie que nous voulons que la solution soit présentée en termes de moments bruts de la population (plutôt que de moments centraux ou cumulants). ]

Enfin, si ~ Exponentiel ( ) avec pdf :λ f ( x )Xλf(x)

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

alors nous pouvons remplacer les moments dans la solution générale par les valeurs réelles pour une variable aléatoire exponentielle, comme ceci:μisol

entrez la description de l'image ici

Terminé.


PS La raison pour laquelle les autres solutions affichées ici donnent une réponse avec dans le dénominateur plutôt que dans le numérateur est, bien sûr, parce qu'elles utilisent un paramétrage différent de la distribution exponentielle. Étant donné que l'OP n'a pas indiqué la version qu'il utilisait, j'ai décidé d'utiliser la définition standard du manuel de théorie de la distribution Johnson Kotz et al… juste pour équilibrer les choses :)λ2

Wolfies
la source