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)

Hatchepsout
la source

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)
Dilip Sarwate
la source
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.

Glen_b -Reinstate Monica
la source
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)

Deep North
la source
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

Jonathan
la source