Intuicja stojąca za tym, dlaczego paradoks Stein'a dotyczy tylko wymiarów


46

Przykład Steina pokazuje, że oszacowanie maksymalnego prawdopodobieństwa n zmiennych o rozkładzie normalnym ze średnimi μ1,,μn i wariancjami 1 jest niedopuszczalne (pod funkcją straty kwadratowej) iff n3 . Aby uzyskać dobry dowód, zobacz pierwszy rozdział Wnioskowania na dużą skalę: empiryczne metody Bayesa do szacowania, testowania i przewidywania autorstwa Bradleya Effrona.

xN(μ,1)Ex2μ2+n

Moje pytanie brzmi raczej: jakiej właściwości wymiarowej przestrzeni (dla ) brakuje co ułatwia przykład Steina? Możliwe odpowiedzi mogą dotyczyć krzywizny sfery lub czegoś zupełnie innego.n 3 R 2 nnn3R2n

Innymi słowy, dlaczego MLE jest dopuszczalny w ?R2


Edycja 1: W odpowiedzi na obawy @mpiktas o 1.31 wynikające z 1.30:

Eμ(zμ^2)=Eμ(S(N2S)2)=Eμ((N2)2S).

μi^=(1N2S)zi
więc
Eμ(μi^zi)=Eμ(1N2S+2zi2S2).
Dlatego mamy:

2i=1NEμ(μi^zi)=2N2Eμ(N(N2)S)+4Eμ((N2)S)=2NEμ2(N2)2S.

Edycja 2 : W tym artykule Stein udowadnia, że ​​MLE jest dopuszczalny dla N=2 .


4
@mpiktas Nie jest tak niestosowne, jak się wydaje. Sytuacja jest podobna do ANOVA po zastosowaniu redukcji wystarczalności. Wskazuje to, że zwykłe oszacowania ANOVA średnich dla grupy są niedopuszczalne, pod warunkiem, że staramy się oszacować średnie dla więcej niż 3 grup (co okazuje się być prawdą). Poleciłbym przyjrzeć się dowodom, że MLE jest dopuszczalny dla i zobaczyć, gdzie zawodzą przy próbie rozszerzenia do zamiast po prostu spojrzeć na dowody, że estymator Stein robi to, co twierdzi, co jest łatwe raz. masz na myśli estymator. N = 3N=1,2N=3
facet

2
... i umiesz używać Lemmy Stein'a. Wydaje mi się, że jest to trochę mniej proste niż 6 minut temu.
facet

2
Zgadzam się. Czy masz jakieś dobre referencje (oprócz oryginału). Uważam, że oryginalny artykuł Stein'a był zbyt obliczeniowy i miałem nadzieję, że ktoś opracuje inną metodę w ciągu ostatnich pięćdziesięciu lat.
Har

2
Dowiedziałem się, że nauczono mnie Browna i Hwanga z 1983 r., W którym, jak sądzę, zastosowano metodę sugerowaną przez Blytha z wczesnych lat 50. Jest dość ogólny (bardziej ogólny niż wynik Stein, ponieważ działa dla wykładniczej rodziny) i, jak sądzę, całkiem różny od Steina. Ale to nie jest trywialne.
facet

2
@Har świetne pytanie! (+1)
suncoolsu

Odpowiedzi:


43

Dychotomia między przypadkami i dla dopuszczalności MLE średniej wymiarowej wielowymiarowej normalnej zmiennej losowej jest z pewnością szokująca.d 3 dd<3d3d

Istnieje inny bardzo znany przykład prawdopodobieństwa i statystyki, w którym istnieje dychotomia między przypadkami i . Jest to ponowny prosty losowy spacer po kratce . To znaczy, prosty wymiarowy wymiarowy jest powtarzalny w 1 lub 2 wymiarach, ale jest przejściowy w wymiarach . Obowiązuje również analog czasu ciągłego (w postaci ruchu Browna).d 3 Z d d d 3d<3d3Zddd3

Okazuje się, że oba są ze sobą ściśle powiązane.

Larry Brown udowodnił, że dwa pytania są zasadniczo równoważne. Oznacza to, że najlepszy niezmienny estymator o wymiarową wielowymiarowej normalnym średnim wektora jest dopuszczalna wtedy i tylko wtedy, gdy wymiarową Browna jest okresowe.ddμ^μ^(X)=Xdd

W rzeczywistości jego wyniki idą znacznie dalej. Dla każdego sensownego (tj. Uogólnionego estymatora Bayesa) z ograniczonym (uogólnionym) ryzykiem , istnieje wyraźne (!) Odpowiadające wymiarowe rozproszenie, tak że estymator jest dopuszczalny wtedy i tylko wtedy, gdy odpowiadająca mu dyfuzja jest powtarzalna.L2d ˜ μμ~μ~(X)L2dμ~

Lokalny Średni tej dyfuzji jest zasadniczo taka rozbieżność między tymi dwoma estymatorów, czyli i kowariancja dyfuzji jest . Z tego łatwo zauważyć, że w przypadku MLE odzyskujemy (przeskalowaliśmy) ruch Browna.μ~μ^2Iμ~=μ^=X

Tak więc, w pewnym sensie, możemy spojrzeć na kwestię dopuszczalności przez pryzmat procesów stochastycznych i użyć dobrze zbadanych właściwości dyfuzji, aby dojść do pożądanych wniosków.

Bibliografia

  1. L. Brown (1971). Dopuszczalne estymatory, powtarzające się dyfuzje i nierozpuszczalne problemy z wartością graniczną . Ann. Matematyka Stat. , vol. 42, nr 3, s. 855–903.
  2. RN Bhattacharya (1978). Kryteria ponownego wystąpienia i istnienia niezmiennych miar dla dyfuzji wielowymiarowych . Ann. Prob. , vol. 6, nr 4, 541–553.

2
Właściwie to coś takiego miałem nadzieję. Związek z inną dziedziną matematyki (geometria różniczkowa lub procesy stochastyczne), która pokazuje, że dopuszczalność dla nie była tylko przypadkiem. Świetna odpowiedź! n=2
Har

Zainspirowany twoją odpowiedzią, podałem
Henry.L

21

@cardinal dał świetną odpowiedź (+1), ale cały problem pozostaje tajemniczy, chyba że znamy dowody (a ja nie). Myślę więc, że pozostaje pytanie, co jest intuicyjnym powodem, dla którego paradoks Stein'a nie pojawia się w i .RR2

Bardzo pomocna jest dla mnie perspektywa regresji przedstawiona w Stephen Stigler, 1990, A Galtonian Perspective on Shrinkage Estimators . Rozważ niezależne pomiary , z których każdy mierzy pewne podstawowe (nieobserwowane) i pobiera próbki z . Gdybyśmy w jakiś sposób znali , moglibyśmy stworzyć wykres rozproszenia par :XiθiN(θi,1)θi(Xi,θi)

Paradoks Stein'a: perspektywa regresji

Linia ukośna odpowiada zerowemu szumowi i doskonałemu oszacowaniu; w rzeczywistości hałas jest niezerowy, a zatem punkty są przesunięte względem linii ukośnej w kierunku poziomym . Odpowiednio, może być postrzegane jako linia regresji na . Znamy jednak i chcemy oszacować , więc powinniśmy raczej rozważyć linię regresji na - która będzie miała inne nachylenie, odchylone w poziomie , jak pokazano na rysunku (linia przerywana).θ=Xθ=XXθXθθX

Cytując z pracy Stiglera:

Ta galtonowska perspektywa paradoksu Stein czyni z niej niemal przejrzystą. „Zwykłe” estymatory pochodzą od teoretycznej linii regresji on . Ta linia byłaby przydatna, gdyby naszym celem było przewidywanie z , ale nasz problem jest odwrotny, mianowicie przewidywanie z za pomocą sumy błędów kwadratowych jako kryterium. Dla tego kryterium optymalne estymatory liniowe są podane przez linię regresji najmniejszych kwadratów naXθXθθXĎ(θı - θ ı)2θXθ^i0=XiXθXθθX(θiθ^i)2θX, a estymatory Jamesa-Steina i Efrona-Morrisa same są estymatorami tego optymalnego estymatora liniowego. „Zwykłe” estymatory pochodzą z niewłaściwej linii regresji, estymatory Jamesa-Steina i Efrona-Morrisa pochodzą z aproksymacji do właściwej linii regresji.

A teraz pojawia się kluczowy bit (wyróżnienie dodane):

Możemy nawet zobaczyć dlaczego jest konieczne: jeśli lub , najmniej linia kwadratów na musi przechodzić przez punkty , a tym samym dla lub , w dwie linie regresji ( na i na ) muszą się zgadzać na każdym .k32 θ X ( X i , θ i ) k = 1 2 X θ θ X X ik=12θX(Xi,θi)k=12XθθXXi

Myślę, że to wyjaśnia, co jest specjalnego w i .k = 2k=1k=2

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.