Combinação linear de duas variáveis ​​aleatórias normais multivariadas dependentes

13

Suponha que temos dois vetores de variáveis ​​aleatórias, ambos são normais, ou seja, e . Estamos interessados ​​na distribuição de sua combinação linear , onde A e B são matrizes, C é um vetor. Se X e Y forem independentes, Z \ sim N (A \ mu_X + B \ mu_Y + C, A \ Sigma_X A ^ T + B \ Sigma_Y B ^ T) . A questão está no caso dependente, assumindo que sabemos a correlação de qualquer par (X_i, Y_i) . Obrigado.XN(μX,ΣX)YN(μY,ΣY)Z=AX+BY+CABCXYZN(AμX+BμY+C,AΣXAT+BΣYBT)(Xi,Yi)

Muitas felicidades, Ivan

Ivan
fonte

Respostas:

7

Nesse caso, você deve escrever (com notações claras) ( editado: assumindo a normalidade conjunta de ) Então e ie (X,Y)AX+BY= ( A B ) ( X Y ) AX+BY+C N [ ( A B ) ( μ X μ Y ) + C ,

(XY)N[(μXμY),ΣX,Y]
(X,Y)
AX+BY=(AB)(XY)
AX+BY+CN [ A μ X +B μ Y +C,A Σ X X A T +B Σ T X Y A T +A Σ X Y B T +B Σ Y Y
AX+BY+CN[(AB)(μXμY)+C,(AB)ΣX,Y(ATBT)]
AX+BY+CN[AμX+BμY+C,AΣXXAT+BΣXYTAT+AΣXYBT+BΣYYBT]
Xi'an
fonte
3
Caso seja ignorado, observe que o tópico de comentário para outra resposta indica (a) esses cálculos de covariância são bons (entendendo que envolvem uma notação de matriz de bloco natural, mas não declarada), mas (b) não podemos concluir validamente que as combinações lineares são normalmente distribuído até que façamos uma suposição adicional; ou seja, que e têm uma distribuição normal multivariada conjunta . YXY
whuber
2
Você poderia explicar como chegou de a na última linha? Eu teria pensado que e o resultado não simplifica mais. Aqui não é uma matriz simétrica, pois seu -ésimo elemento é enquanto seu -ésimo elemento é , e não há razão para que essas covariâncias sejam iguais. BΣXYTAT+AΣXYBT2AΣXYBT
BΣXYTAT+AΣXYBT=(AΣXYBT)T+AΣXYBT
ΣXY(i,j)cov(Xi,Yj)(j,i)cov(Xj,Yi)
Dilip Sarwate
1
@DilipSarwate: (+1) você está certo, no caso geral, não há razão para esses dois termos serem iguais.
Xian
3

Sua pergunta não tem uma resposta única, como é atualmente colocado, a menos que você assuma que e são normalmente distribuídos em conjunto com o bloco superior direito . Eu acho que você quer dizer isso porque diz que tem cada covariância entre X e Y. Nesse caso, podemos escrever que também é normal multivariada. então é dado em termos de como:XYΣXYW=(XT,YT)TZW

Z=(A,B)W+C

Então você usa sua fórmula usual para combinação linear. Observe que a média é inalterada, mas a matriz de covariância possui dois termos extras adicionadosAΣXYBT+BΣXYTAT

probabilityislogic
fonte
Obrigado por apontar esse problema, na verdade, eu nem pensei nisso, mas parece que as variáveis ​​podem realmente ser vistas, no meu caso, como normalmente distribuídas em conjunto, mesmo que seus componentes estejam correlacionados.
1526 Ivan Ivan
Concordo que a questão não pode ser resolvida como foi colocada. Ele pode ser resolvido de uma maneira simples , se um assume, como @ A resposta de Xian faz, que e , em conjunto, normalmente distribuídos. Poderia ser solucionável, presumivelmente com mais dificuldade, se a distribuição da articulação fosse especificada como algo diferente da normal da articulação. Mas apenas conhecer para todos , não significa que é normal multivariada . Quaisquer duas variáveis ​​aleatórias com variações finitas têm uma covariância. A covariância não é definida apenas para variáveis ​​aleatórias normais ou conjuntamente normais.XYcov(Xi,Yj)i,jW=(XT,YT)T
Dilip Sarwate
No meu caso, X e Y são normais em conjunto; tentarei explicar o porquê; por favor, corrija-me se estiver errado. Suponha que exista um conjunto de RVs normais univariados independentes. Cada elemento de X e Y é uma combinação linear arbitrária dessas variáveis ​​univariadas do conjunto. Portanto, uma vez que as variáveis ​​iniciais são independentes e apenas transformações lineares são envolvidas, os vetores resultantes X, Y e Z são todos rv normais normais multivariados. Segue a definição de um rv normal multivariado, em que deve ser um rv normal univariado para qualquer vetor . Isso faz sentido? aTXa
1526 Ivan Ivan
1
@Ivan Sua explicação faz sentido, mas a reclamação é sobre a afirmação "Suponha que tenhamos dois vetores de variáveis ​​aleatórias, ambos são normais, ou seja, e ", o que não significa que e sejam conjuntamente normais . Tampouco dizer que "sabemos a correlação de qualquer par " significa que e são conjuntamente normais , embora, como você declara corretamente, implique que seja normal (e da mesma forma para .) Normalidade univariadaXN(μX,ΣX)YN(μY,ΣY)XY(Xi,Yi)XiYiXN(μX,ΣX)XiYinão implica normalidade das articulações. Veja a referência abaixo.
precisa
@Ivan Veja a discussão após esta pergunta
Dilip Sarwate