Testowanie hipotez na macierzy odwrotnej kowariancji


10

Załóżmy, że obserwuję iid i chcę przetestować H 0 : A vech ( Σ - 1 ) = a dla zgodnej macierzy A i wektora a . Czy znane są prace nad tym problemem?xiN(μ,Σ)H0:A (Σ1)=aAa

Oczywistą (do mnie) próba byłaby za pomocą testu współczynnika prawdopodobieństwa, ale wydaje się, zwiększając prawdopodobieństwo zastrzeżeniem ograniczeń wymagałoby SDP solver i może być dość owłosione.H0


1
Czy masz jakieś dodatkowe ograniczenia dotyczące ? Jeśli A jest odwracalny, wówczas H 0 = v e c h ( Σ - 1 ) = A - 1 a . Wówczas problem sprowadza się do dobrze znany problem: że bada się czy Ď - 1 = B . Tutaj v e c h ( B ) = A - 1 a (pamiętaj, że v e c h ( B )AAH0=vech(Σ1)=A1aΣ1=Bvech(B)=A1avech(B)określa jednoznacznie). B
MånsT

@ MånsT; niestety jestem zainteresowany ogólną sprawą. Zazwyczaj ma około 10 wierszy i około 400 kolumn. A
shabbychef

1
(A,a)

Σ1A

Odpowiedzi:


3

Beran i Srivastava (1985, Annals of Statistics) napisali artykuł, w którym zaproponowali ogólne podejście bootstrap w celu zastosowania rotacji do macierzy kowariancji, która dopasowuje rozkład do wartości zerowej. @ Kardynał spostrzeżenie o istnieniu takiej macierzy jest tutaj jednak bardzo istotne. Musisz być w stanie wymyślić przynajmniej jakieś przybliżenie macierzy, która spełnia ograniczenia narzucone pod wartością zerową.

Chen, Variyath i Bovas napisali artykuł o skorygowanym prawdopodobieństwie empirycznym, w którym zademonstrowali, jak można go wykorzystać do przetestowania dość dziwnej struktury na macierzy kowariancji. Myślę, że ten artykuł ostatecznie ukazał się w CJS.


Nie jestem pewien, czy mogę z łatwością przetłumaczyć je na rozwiązanie mojego problemu, ale oba są fascynującymi lekturami. +1.
shabbychef
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.