Bayesowski model logitów - intuicyjne wyjaśnienie?

11

Muszę wyznać, że wcześniej nie słyszałem o tym terminie na żadnej z moich zajęć, na studiach licencjackich ani na studiach.

Co to znaczy, że regresja logistyczna jest Bayesowska? Szukam wyjaśnienia z przejściem od logistyki zwykłej do logistyki bayesowskiej podobnej do następującej:

Jest to równanie w modelu regresji liniowej: .E(y)=β0+β1x1+...+βnxn

Jest to równanie w modelu regresji logistycznej: . Odbywa się to, gdy y jest kategoryczny.ln(E(y)1E(y))=β0+β1x1+...+βnxn

To, co zrobiliśmy, to zmiana na .ln ( E ( y )E(y)ln(E(y)1E(y))

Co zatem zrobiono z modelem regresji logistycznej w bayesowskiej regresji logistycznej? Zgaduję, że to nie ma nic wspólnego z równaniem.

Ten podgląd książki wydaje się definiować, ale tak naprawdę nie rozumiem. Czym są te wszystkie wcześniejsze rzeczy dotyczące prawdopodobieństwa? Co to jest ? Czy ktoś mógłby wyjaśnić tę część książki lub model logitów bayesowskich w inny sposób?α

Uwaga: Myślę, że zostało to już wcześniej zadane, ale nie zostało udzielone zbyt dobrej odpowiedzi.

BCLC
źródło
1
Nie chcę odpowiedzieć na to pytanie, ponieważ myślę, że @Tim ma na to większość odpowiedzi. Jedyną rzeczą, której brakuje w tej świetnej odpowiedzi, jest to, że w bayesowskiej regresji logistycznej i bayesowskich uogólnionych modelach liniowych (GLM) bardziej ogólnie, wcześniejsze rozkłady nie są umieszczane tylko nad współczynnikami, ale nad wariancjami i kowariancją tych współczynników. Jest to niezwykle ważne, aby o tym wspomnieć, ponieważ jedną z kluczowych zalet bayesowskiego podejścia do GLM jest większa łatwość określania, aw wielu przypadkach także dopasowania złożonych modeli kowariancji współczynników.
Brash Equilibrium
2
@BrashEquilibrium: wspominasz o możliwym hierarchicznym rozszerzeniu standardowego modelowania bayesowskiego dla modelu logit. W naszej książce używamy na przykład przed G na „s, przed ustalający macierzy kowariancji pochodzi od współzmiennych X . βX
Xi'an
1
Wystarczająco uczciwy na g Prior.
Brash Equilibrium
1
To powiedziawszy, nadal istnieje przejęcie kowariancji !!!!!! Jeśli nie dyskutujesz, nie opisujesz, jak regresja logistyczna działa całkowicie.
Brash Equilibrium

Odpowiedzi:

19

Regresję logistyczną można opisać jako kombinację liniową

η=β0+β1X1+...+βkXk

który jest przekazywany przez funkcję link :g

g(E(Y))=η

gdzie funkcja link jest funkcją logowania

E(Y|X,β)=p=logit1(η)

gdzie przyjmuje tylko wartości z { 0 , 1 }, a odwrotne funkcje logitowe przekształcają kombinację liniową η na ten zakres. Tu kończy się klasyczna regresja logistyczna.Y{0,1}η

Jeśli jednak przypomnisz sobie, że dla zmiennych, które przyjmują tylko wartości w { 0 , 1 } , to E ( Y | X , β ) można uznać za P ( Y = 1 | X , β ) . W takim przypadku dane wyjściowe funkcji logit można uznać za warunkowe prawdopodobieństwo „sukcesu”, tj. P ( Y = 1 | X ,E(Y)=P(Y=1){0,1}E(Y|X,β)P(Y=1|X,β) . Rozkład Bernoulliegojest rozkładem opisującym prawdopodobieństwo zaobserwowania wyniku binarnego, z pewnymparametrem p , więc możemy opisać Y jakoP(Y=1|X,β)pY

yiBernoulli(p)

Zatem przy regresji logistycznej szukamy niektórych parametrów które razem z niezależnymi zmiennymi X tworzą kombinację liniową η . W regresji klasycznej E ( Y | X , β ) = η (zakładamy, że funkcja link jest funkcją tożsamości), jednak w modelu Y, który przyjmuje wartości w { 0 , 1 } , musimy przekształcić η , aby dopasować [ 0 , 1 ] zasięg.βXηE(Y|X,β)=ηY{0,1}η[0,1]

Teraz, aby oszacować regresję logistyczną w sposób bayesowski, wyłapujesz niektóre priorytety dla parametrów jak w przypadku regresji liniowej (patrz Kruschke i in., 2012 ), a następnie użyj funkcji logit do przekształcenia kombinacji liniowej η , aby użyć jej wyniku jako p parametr rozkładu Bernoulliego opisujący zmienną Y. Tak, tak, faktycznie używasz równania i funkcji logit w taki sam sposób, jak w przypadku częstotliwości, a reszta działa (np. Wybierając priory), podobnie jak szacowanie regresji liniowej metodą bayesowską.βiηpY

Prostym podejściem do wyboru priorów jest wybranie rozkładów normalnych (ale możesz również użyć innych rozkładów, np. Rozkład - lub Laplace'a dla bardziej wytrzymałego modelu) dla β i z parametrami μ i i σ 2 i, które są ustawione lub wzięte z hierarchiczne priory . Teraz, mając definicję modelu, możesz użyć oprogramowania takiego jak JAGS, aby przeprowadzić symulację Markov Chain Monte Carlo , aby oszacować model. Poniżej kod pocztowy Jags dla prostego modelu logistycznego (sprawdź tutaj po więcej przykładów).tβiμiσi2

model {
   # setting up priors
   a ~ dnorm(0, .0001)
   b ~ dnorm(0, .0001)

   for (i in 1:N) {
      # passing the linear combination through logit function
      logit(p[i]) <- a + b * x[i]

      # likelihood function
      y[i] ~ dbern(p[i])
   }
}

Jak widać, kod bezpośrednio przekłada się na definicję modelu. Jakie oprogramowanie robi to rysuje pewne wartości z normalnego priors na ai b, a następnie wykorzystuje te wartości do oszacowania pi wreszcie, wykorzystuje funkcję prawdopodobieństwa ocenić na ile prawdopodobne jest dane podane te parametry (to jest, gdy używasz Twierdzenie Bayesa, patrz tutaj dla bardziej szczegółowy opis).

Podstawowy model regresji logistycznej można rozszerzyć w celu modelowania zależności między predyktorami przy użyciu modelu hierarchicznego (w tym hiperpriorów ). W takim przypadku możesz narysować z wielowymiarowego rozkładu normalnego, który pozwala nam zawrzeć informację o kowariancji Σ między zmiennymi niezależnymiβiΣ

(β0β1βk)MVN([μ0μ1μk],[σ02σ0,1σ0,kσ1,0σ12σ1,kσk,0σk,1σk2])

... ale chodzi tu o szczegóły, więc zatrzymajmy się tutaj.

Część „bayesowska” tutaj wybiera priorytety, używając twierdzenia Bayesa i definiując model w kategoriach probabilistycznych. Zobacz tutaj definicję „modelu bayesowskiego”, a tutaj ogólną intuicję dotyczącą podejścia bayesowskiego . Można także zauważyć, że dzięki temu podejściu definiowanie modeli jest dość proste i elastyczne.


Kruschke, JK, Aguinis, H., i Joo, H. (2012). Nadszedł czas: Bayesowskie metody analizy danych w naukach organizacyjnych. Metody badań organizacyjnych, 15 (4), 722-752.

Gelman, A., Jakulin, A., Pittau, GM, i Su, Y.-S. (2008). Słabo informacyjna domyślna wcześniejsza dystrybucja modeli logistycznych i innych modeli regresji. The Annals of Applied Statistics, 2 (4), 1360–1383.

Tim
źródło
1
Potrzebujesz dowodów na rozbieżności, a nie tylko współczynniki.
Brash Equilibrium,
3
gη = β 0 + β 1 X 1 g E ( Y ) = ηηη=β0+β1X1gE(Y)=η
1
@BCLC sprawdź linki w mojej odpowiedzi, zawierają one ogólne wprowadzenie do statystyki bayesowskiej. Jest to znacznie szerszy temat niż ten wymieniony w początkowym pytaniu, ale można znaleźć miłe wprowadzenie w odniesieniach, które podałem w mojej odpowiedzi.
Tim
1
@Tim Zrobiłem tam literówkę. Dowody mają czytać priory. Zasadniczo współczynniki to nie jedyne nieznane parametry. Rozkład wielomianowy ma również macierz kowariancji wariancji i zazwyczaj nie zakładamy, że jest znany.
Brash Equilibrium,
3
„Część„ bayesowska ”tutaj wybiera priorytety, używając twierdzenia Bayesa i definiując model w kategoriach probabilistycznych.” Dobrym odniesieniem jest tutaj Gelman i in. SŁABO INFORMACYJNA DOMYŚLNA DOMYŚLNA DYSTRYBUCJA DLA MODELI LOGISTYCZNYCH I INNYCH REJESTRACJI stat.columbia.edu/~gelman/research/published/priors11.pdf
Dalton Hance
6

Czym są te wszystkie wcześniejsze rzeczy dotyczące prawdopodobieństwa?

Właśnie dlatego Bayesian. Generatywny model danych jest taki sam; różnica polega na tym, że analiza bayesowska wybiera wcześniejszy rozkład parametrów będących przedmiotem zainteresowania i oblicza lub aproksymuje rozkład tylny , na którym opiera się wszystkie wnioski. Zasada Bayesa odnosi się do dwóch: Tylna część ciała jest proporcjonalna do czasów prawdopodobieństwa wcześniejszych.

β

Niektóre modele częstych mogą być powiązane z odpowiednikiem bayesowskim z konkretnym przeorem, choć nie jestem pewien, co odpowiada w tym przypadku.

Sean Easter
źródło
βββ1,β2,...,βnX1X2Xnβ
1
@BCLC Aby odpowiedzieć na te pytania, zacznę od samego procesu wnioskowania bayesowskiego i definiuję terminy w miarę upływu czasu: Bayesianie traktują wszystkie interesujące parametry jako zmienne losowe i aktualizują swoje przekonania na temat tych parametrów w świetle danych. Przed dystrybucja wyraża swoje przekonanie o parametrach Przed analizą danych; * rozkład tylny * - zgodnie z regułą Bayesa, znormalizowanym iloczynem wcześniejszego i prawdopodobieństwa - podsumowuje niepewne przekonanie o parametrach w świetle wcześniejszego i danych. Obliczanie tyłu jest miejscem, w którym ma miejsce dopasowanie.
Sean Easter
1
βpp
Okej, myślę, że rozumiem cię lepiej po przeczytaniu Eseju na temat rozwiązania problemu w Doktrynie szans . Dzięki SeanEster
BCLC
1
P(B)