Variance des combinaisons linéaires de variables aléatoires corrélées


10

Je comprends la preuve que mais je ne comprends pas comment prouver la généralisation à des combinaisons linéaires arbitraires.

Var(aX+bY)=a2Var(X)+b2Var(Y)+2abCov(X,Y),

Soit scalaires pour donc nous avons un vecteur , et un vecteur de variables aléatoires corrélées. Alors Comment prouvons-nous cela? J'imagine qu'il y a des preuves dans la notation de sommation et dans la notation vectorielle?aii1,,na_X_=Xi,,Xn

Var(a1X1+anXn)=i=1nai2σi2+2i=1nj>inaiaj Cov(Xi,Xj)

Réponses:


18

Il s'agit simplement d'un exercice d'application des propriétés de base des sommes, de la linéarité des attentes et des définitions de la variance et de la covariance

var(i=1naiXi)=E[(i=1naiXi)2](E[i=1naiXi])2one definition of variance=E[i=1nj=1naiajXiXj](E[i=1naiXi])2basic properties of sums=i=1nj=1naiajE[XiXj](i=1naiE[Xi])2linearity of expectation=i=1nj=1naiajE[XiXj]i=1nj=1naiajE[Xi]E[Xj]basic properties of sums=i=1nj=1naiaj(E[XiXj]E[Xi]E[Xj])combine the sums=i=1nj=1naiajcov(Xi,Xj)apply a definition of covariance=i=1nai2var(Xi)+2i=1nj:j>inaiajcov(Xi,Xj)re-arrange sum
Notez que dans cette dernière étape, nous avons également identifié comme la variance .cov(Xi,Xi)var(Xi)

6

Vous pouvez réellement le faire par récursivité sans utiliser de matrices:

Prenez le résultat pour et laissez .Var(a1X1+Y1)Y1=a2X2+Y2

Var(a1X1+Y1)

=a12Var(X1)+2a1Cov(X1,Y1)+Var(Y1)

=a12Var(X1)+2a1Cov(X1,a2X2+Y2)+Var(a2X2+Y2)

=a12Var(X1)+2a1a2Cov(X1,X2)+2a1Cov(X1,Y2)+Var(a2X2+Y2)

Continuez ensuite à remplacer et en utilisant les mêmes résultats de base, puis à la dernière étape, utilisezYi1=aiXi+YiYn1=anXn

Avec des vecteurs (donc le résultat doit être scalaire):

Var(aX)=aVar(X)a

Ou avec une matrice (le résultat sera une matrice variance-covariance):

Var(AX)=AVar(X)A

Ceci a l'avantage de donner des covariances des différentes combinaisons linéaires dont les coefficients sont les rangées de sur les éléments hors diagonales du résultat.A

Même si vous ne connaissez que les résultats univariés, vous pouvez les confirmer en vérifiant élément par élément.


2

Fondamentalement, la preuve est la même que la première formule. Je vais prouver qu'il utilise une méthode très brutale.

Var(a1X1+...+anXn)=E[(a1X1+..anXn)2][E(a1X1+...+anXn)]2=E[(a1X1)2+...+(anXn)2+2a1a2X1X2+2a1a3X1X3+...+2a1anX1Xn+...+2an1anXn1Xn][a1E(X1)+...anE(Xn)]2

=a12E(X12)+...+an2E(Xn2)+2a1a2E(X1X2)+...+2an1anE(Xn1Xn)a12[E(X1)]2...an2[E(Xn)]22a1a2E(X1)E(X2)...2an1anE(Xn1)E(Xn)

=a12E(X12)a12[E(X1)]2+...+an2E(Xn2)an2[E(Xn)]2+2a1a2E(X1X2)2a1a2E(X1)E(X2)+...+2an1anE(Xn1Xn)2an1anE(Xn1)E(Xn)

Ensuite, notez simplement:

an2E(Xn2)an2[E(Xn)]2=anσn2

et

2unen-1unenE(Xn-1Xn)-2unen-1unenE(Xn-1)E(Xn)=2unen-1unenCov(Xn-1,Xn)


2

Juste pour le plaisir, preuve par induction!

Soit la déclaration queP(k)Vuner[je=1kunejeXje]=je=1kuneje2σje2+2je=1kj>jekunejeunejCov[Xje,Xj]

Alors est (trivialement) vrai (vous avez dit que vous en étiez satisfait dans la question).P(2)

Supposons que P (k) soit vrai. Donc,

Vuner[je=1k+1unejeXje]=Vuner[je=1kunejeXje+unek+1Xk+1]

=Vuner[je=1kunejeXje]+Vuner[unek+1Xk+1]+2Cov[je=1kunejeXje,unek+1Xk+1]

=je=1kuneje2σje2+2je=1kj>jekunejeunejCov[Xje,Xj]+unek+12σk+12+2Cov[je=1kunejeXje,unek+1Xk+1]

=je=1k+1uneje2σje2+2je=1kj>jekunejeunejCov[Xje,Xj]+2je=1kunejeunek+1Cov[Xje,Xk+1]

=je=1k+1uneje2σje2+2je=1k+1j>jek+1unejeunejCov[Xje,Xj]

Ainsi, est vrai.P(k+1)

Donc, par induction,

Vuner[je=1nunejeXje]=je=1nuneje2σje2+2je=1nj>jenunejeunejCov[Xje,Xj] pour tout entier .n2

En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.