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

shabbychef
źródło
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.

StasK
źródło
Nie jestem pewien, czy mogę z łatwością przetłumaczyć je na rozwiązanie mojego problemu, ale oba są fascynującymi lekturami. +1.
shabbychef