Czy ktoś mógłby mi doradzić, jak interpretować szacunki z regresji logistycznej za pomocą linku cloglog?
Zamontowałem następujący model lme4
:
glm(cbind(dead, live) ~ time + factor(temp) * biomass,
data=mussel, family=binomial(link=cloglog))
Na przykład szacunkowy czas wynosi 0,015. Czy słuszne jest twierdzenie, że iloraz śmiertelności na jednostkę czasu jest mnożony przez exp (0,015) = 1,015113 (~ 1,5% wzrost na jednostkę czasu).
Innymi słowy, czy szacunki uzyskane w cloglog wyrażone są w ilorazach logarytmicznych, jak w przypadku regresji logistycznej logit?
logistic
regression-coefficients
jatalah
źródło
źródło
R
reguł składni. Nie możesz mieć (po 'Odpowiedzi:
Dzięki funkcji uzupełniającej log-log nie jest to regresja logistyczna - termin „logistic” oznacza link logit. Oczywiście nadal jest to regresja dwumianowa.
Nie, ponieważ nie modeluje się w kategoriach logarytmicznych szans. To właśnie miałbyś z linkiem logit; jeśli chcesz model działający pod względem logarytmicznych szans, użyj linku logit.
Mówi to funkcja łącza uzupełniającego dziennika-dziennika
gdzie .πx= P( Y= 1 | X= x )
Zatem nie jest ilorazem szans; faktycznie .exp( η) exp( η) = - log( 1 - πx)
Stąd i . W rezultacie, jeśli potrzebujesz ilorazu szans dla określonego , możesz go obliczyć, ale parametry nie mają bezpośredniej prostej interpretacji pod względem udziału w log-odds.exp( - exp( η) ) = ( 1 - πx) 1 - exp( - exp( η) ) = πx x
Zamiast tego (co nie dziwi) parametr pokazuje (dla zmiany jednostki ) wkład do logu komplementarnego.x
Jak Ben delikatnie wskazał w swoim pytaniu w komentarzach:
Parametry w uzupełniającym modelu log-log mają zgrabną interpretację pod względem współczynnika ryzyka. Mamy to:
(W tym przykładzie przeżycie dziennika spadnie o około 1,5% na jednostkę czasu).
Teraz zagrożenie, , więc rzeczywiście wydaje się, że w przykładzie podane w pytaniu prawdopodobieństwo śmierci * na jednostkę czasu wzrasta o około 1,5%h ( x ) = - drexlog( Sx) = drexmiη( x )
* (lub bardziej ogólnie dla modeli dwumianowych z linkiem cloglog, )P.(Y= 1 )
źródło