Interpretieren von Schätzungen der logistischen Regression von Cloglogs


21

Könnte mir jemand raten, wie ich die Schätzungen aus einer logistischen Regression mithilfe eines Cloglog-Links interpretieren kann?

Ich habe folgendes Modell eingebaut lme4:

glm(cbind(dead, live) ~ time + factor(temp) * biomass,
    data=mussel, family=binomial(link=cloglog))

Beispielsweise beträgt die geschätzte Zeit 0,015. Ist es richtig zu sagen, dass die Wahrscheinlichkeit der Sterblichkeit pro Zeiteinheit mit exp (0,015) = 1,015113 multipliziert wird (~ 1,5% Zunahme pro Zeiteinheit)?
Mit anderen Worten, werden die in einem Cloglog erhaltenen Schätzungen in Log Odds ausgedrückt, wie dies bei einer logit logistic-Regression der Fall ist?


Bitte bearbeiten Sie den Code, um den RSyntaxregeln zu folgen . Sie können nicht haben (nach '
Frank Harrell

Bearbeiten Sie den ursprünglichen Beitrag und löschen Sie den Kommentar.
Frank Harrell

Antworten:


30

Bei einer Komplementär-Protokoll-Protokoll-Verknüpfungsfunktion handelt es sich nicht um eine logistische Regression, sondern um eine logistische Verknüpfung. Es ist natürlich immer noch eine binomische Regression.

Die geschätzte Zeit beträgt 0,015. Ist es richtig zu sagen, dass die Wahrscheinlichkeit der Sterblichkeit pro Zeiteinheit mit exp (0,015) = 1,015113 multipliziert wird (~ 1,5% Zunahme pro Zeiteinheit)?

Nein, weil es keine logarithmischen Gewinnchancen bietet. Das ist, was Sie mit einem Logit-Link haben würden; Wenn Sie ein Modell möchten, das in Bezug auf Log-Odds funktioniert, verwenden Sie einen Logit-Link.

Die Komplementär-Log-Log-Link-Funktion sagt das aus

η(x)=log(log(1πx))=xβ

πx=P(Y=1|X=x)

exp(η)exp(η)=log(1πx)

exp(exp(η))=(1πx)1exp(exp(η))=πxx

x


Wie Ben in seinen Kommentaren in seiner Frage sanft angedeutet hat:

Stimmt es, dass sich die Wahrscheinlichkeit der Sterblichkeit pro Zeiteinheit (dh die Gefahr) um 1,5% erhöht?

Die Parameter im komplementären Protokoll-Protokoll-Modell sind in Bezug auf das Gefährdungsverhältnis gut interpretiert. Wir haben das:

eη(x)=log(1πx)=log(Sx)S

(Das logarithmische Überleben sinkt also im Beispiel um ca. 1,5% pro Zeiteinheit.)

h(x)=ddxlog(Sx)=ddxeη(x)

P(Y=1)


7
Stimmt es, dass sich die Wahrscheinlichkeit der Sterblichkeit pro Zeiteinheit (dh die Gefahr) um 1,5% erhöht?
Ben Bolker
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.