Medida apropriada para encontrar a menor matriz de covariância

10

No livro que estou lendo, eles usam definição positiva (definição semi-positiva) para comparar duas matrizes de covariância. A idéia é que, se é pd então é menor do que . Mas estou lutando para conseguir a intuição desse relacionamento?ABBA

Há um tópico semelhante aqui:

/math/239166/what-is-the-intuition-for-using-definiteness-to-compare-matrices

Qual é a intuição para usar a definição definitiva para comparar matrizes?

Embora as respostas sejam boas, elas realmente não abordam a intuição.

Aqui está um exemplo que eu acho confuso:

[1612129][1224]

Agora aqui o determinante da diferença é -25, de modo que a relação não é pd ou mesmo psd e, portanto, a primeira matriz não é maior que a primeira?

Eu simplesmente quero comparar duas matrizes de covariância 3 * 3 para ver qual é a menor? Parece-me mais intuitivo usar algo como a norma euclidiana para compará-los? No entanto, isso significa que a primeira matriz acima é maior que o segundo matix. Além disso, só vejo o critério pd / psd usado para comparar matrizes de covariância.

Alguém pode explicar por que o pd / psd é melhor do que usar outra medida, como a norma euclidiana?

Também publiquei essa pergunta no fórum de matemática (não sabia o que era melhor), espero que isso não viole nenhuma regra.

/math/628135/comparing-two-covariance-matrices

Baz
fonte
2
Você pode ler isso onde se considera a intuição por trás da (semi) definição positiva. Quando você compara 2 variações ae b, se a-bfor positivo, então diríamos que após a remoção variabilidade bfora do aresta alguma variabilidade "real" deixada no a. Da mesma forma, é um caso de variâncias multivariadas (= matrizes de covariância) Ae B. Se A-Bé definido positivamente, isso significa que a A-Bconfiguração de vetores é "real" no espaço euclidiano: em outras palavras, ao remover Bde A, o último ainda é uma variabilidade viável.
precisa saber é
2
O que você quer dizer com "menor" de duas matrizes de covariância?
whuber
Oi whuber, as matrizes de covariância referem-se a estimadores concorrentes, desejo selecionar o estimador que tem a menor variação. (Será que isso esclarecer as coisas?)
Baz
2
Baz: Então, por que não comparar diretamente as variações dos estimadores?
Glen_b -Reinstala Monica
Olá, o método está definido, a expressão para o que eles chamam de variação (que inclui covariâncias) é dada. No entanto, mesmo que eu comparasse apenas as variações, isso ainda envolveria a comparação de valores vetoriais que terão problemas semelhantes aos da matriz?
Baz

Respostas:

8

A ordem das matrizes a que você se refere é conhecida como ordem de Loewner e é uma ordem parcial muito usada no estudo de matrizes definidas positivas. Aqui está um tratamento de livro da geometria no coletor de matrizes positivas-definidas (posdef) .

Primeiro tentarei abordar sua pergunta sobre intuições . Uma matriz (simétrica) A é posdef se cTAc0 para todos os cRn . Se X é uma variável aleatória (rv) com matriz de covariância A , então cTX é (proporcional a) sua projeção em algum subespaço de uma dimensão e Var(cTX)=cTAc . Aplicando isso a ABem sua Q, em primeiro lugar: é uma matriz de covariância, segunda: uma variável aleatória com matriz COVAR B projectos em todas as direções com variação menor do que um rv com covariância matriz A . Isso deixa intuitivamente claro que essa ordem pode ser apenas parcial; existem muitos rvs que se projetam em direções diferentes com variações muito diferentes. Sua proposta de alguma norma euclidiana não tem uma interpretação estatística tão natural.

Seu "exemplo confuso" é confuso porque ambas as matrizes têm zero determinante. Portanto, para cada um, há uma direção (o vetor próprio com valor próprio zero) em que eles sempre se projetam para zero . Mas essa direção é diferente para as duas matrizes, portanto elas não podem ser comparadas.

A ordem Loewner é definida de forma que AB , B seja mais positivo definido que A , se BA for posdef. Esta é uma ordem parcial, para algumas matrizes posdef nem BA nem AB são posdef. Um exemplo é:

A=(10.50.51),B=(0.5001.5)
Uma maneira de mostrar isso graficamente é desenhar um gráfico com duas elipses, mas centralizado na origem, associado de maneira padrão às matrizes (então a distância radial em cada direção é proporcional à variação de projetando nessa direção):

Duas matrizes posdef mostradas como elipses

Nesse caso, as duas elipses são congruentes, mas giradas de maneira diferente (na verdade, o ângulo é de 45 graus). Isso corresponde ao fato de que as matrizes A e B têm os mesmos valores próprios, mas os vetores próprios são rotacionados.

Como essa resposta depende muito das propriedades das elipses, a seguir, qual é a intuição por trás das distribuições gaussianas condicionais? explicar elipses geometricamente, pode ser útil.

Agora vou explicar como as elipses associadas às matrizes são definidas. Uma matriz posdef A define uma forma quadrática QA(c)=cTAc . Isso pode ser plotado como uma função, o gráfico será quadrático. Se AB , em seguida, o gráfico de QB será sempre acima do gráfico de QA . Se recortarmos os gráficos com um plano horizontal na altura 1, os cortes descreverão elipses (que na verdade são uma maneira de definir elipses). As elipses deste corte são dadas pelas equações

QA(c)=1,QB(c)=1
e vemos queAB corresponde à elipse de B (agora com interior) está contida na elipse de A. Se não houver ordem, não haverá contenção. Observamos que a ordem de inclusão é oposta à ordem parcial de Loewner, se não gostarmos de que podemos desenhar elipses dos inversos. Isso porqueAB é equivalente aB1A1 . Mas vou ficar com as elipses, conforme definido aqui.

Uma elipse pode ser descrita com os semiaxos e seu comprimento. Discutiremos apenas 2×2 -matrizes aqui, pois são eles que podemos desenhar ... Portanto, precisamos dos dois eixos principais e seu comprimento. Isso pode ser encontrado, conforme explicado aqui, com uma composição independente da matriz posdef. Então os eixos principais são dados pelos vetores próprios, e seu comprimento uma,b pode ser calculado a partir dos valores próprios λ1 1,λ2 por

uma=1 1/λ1 1,b=1 1/λ2.
Também podemos ver que a área da elipse representandoUMAéπumab=π1 1/λ1 11 1/λ2=πdetUMA .

Vou dar um exemplo final em que as matrizes podem ser ordenadas:

Duas matrizes que podem ser ordenadas plotadas como elipses

As duas matrizes neste caso foram:

UMA=(2/31 1/51 1/53/4),B=(1 11 1/71 1/71 1)

kjetil b halvorsen
fonte
3

O @kjetil b halvorsen oferece uma boa discussão sobre a intuição geométrica por trás da semi-definição positiva como uma ordem parcial. Vou dar uma opinião mais suja da mesma intuição. Um que procede de que tipos de cálculos você gostaria de fazer com suas matrizes de variação.

xyV(x)V(y)V(x)=5V(y)=15xy

xy

V(x)=[1 10,50,51 1]V(y)=[8336]
x1 1y1 1V(x1 1)=1 1<8=V(y1 1)V(x2)=1 1<6=V(y2)xyxyV(x)V(y)V(x)V(y)

V(x)=[V(x1 1)0 00 0V(x2)]V(y)=[V(y1 1)0 00 0V(y2)]
V(y)-V(x)V(x)V(y)V(x1 1)V(y1 1)V(x2)V(y2)
V(x)=[1 10,10,11 1]V(y)=[1 10 00 01 1]
V(x)V(y)V(xk)V(yk)3x1 1+2x23y1 1+2y2V(3x1 1+2x2)>V(3y1 1+2y2)V(x)V(y)

xyV(x)V(y)umaV(umax)V(umay)

V(x)V(y)V(uma1 1x1 1+uma2x2)V(uma1 1y1 1+uma2y2)uma1 1uma2uma1 1=1 1,uma2=0 0V(x1 1)V(y1 1)uma1 1=0 0,uma2=1 1V(x2)V(y2)

V(x)V(y)V(umax)V(umay)uma

V(umay)-V(umax)=umaV(x)uma-umaV(y)uma=uma(V(x)-V(y))uma
V(x)V(y)V(umax)V(umay)uma(V(y)-V(x))

VWW-V

Conta
fonte