R Binomialfamilie mit Identitätsverknüpfung

8

Ich möchte ein lineares Modell von R anpassen family=binomial(link="identity"), wobei jedoch die Binomialfamilie keine Identitätsverknüpfung hat. Was sollte ich tun?

David
quelle
Ich denke, hier liegt ein statistisches Problem zugrunde.
Glen_b -State Monica
Ja, die folgende Frage lautet: Passen Sie den Standardfehler für Überdispersion an.
David
Aber für den ersten muss ich den Identitätslink in der Binomialfamilie verwenden, aber R erlaubt es nicht.
David
1
Beachten Sie, dass Sie in Betracht ziehen sollten, direkt zu einem Quasi-Modell mit binomialer Varianzfunktion zu wechseln, wenn Sie sich sowohl mit Überdispersion als auch mit einem Identitätslink befassen müssen. Ein Nur-Intercept-Binomialmodell kann von Hand angepasst werden .
Glen_b -Rate State Monica
2
(Und wenn Sie wollten, dass der Dispersionsparameter auf eins festgelegt wird : family=binomial(link=make.link("identity")).
Scortchi - Reinstate Monica

Antworten:

9

Siehe Wikipedia zum linearen Wahrscheinlichkeitsmodell und CV-Beiträge hier und hier für den statistischen Hintergrund. Obwohl dies nicht "falsch" ist, möchten Sie einen guten Grund für die Verwendung eines Identitätslinks zur Modellierung einer Bernoulli-Wahrscheinlichkeit.

Nach dem familyHandbuch

die binomischen Familie [akzeptiert] die Verbindungen logit, probit, cauchit, (entsprechend logistischer, normal und Cauchy CDFs respectively) log und cloglog(complementary log-log)

Aber

Die Verknüpfungs- und Varianzargumente haben eine ziemlich umständliche Semantik für die Rückkompatibilität. Es wird empfohlen, sie als Zeichenfolgen in Anführungszeichen anzugeben. Sie können jedoch auch ohne Anführungszeichen (als Namen oder Ausdrücke) angegeben werden. Darüber hinaus können sie auch als Zeichenvektor der Länge 1 mit dem Namen einer der Optionen oder als Liste (für link, der Klasse "link-glm") angegeben werden. Die Einschränkungen gelten nur für Links, die als Namen angegeben werden: Wenn sie als Zeichenfolge angegeben werden, werden alle bekannten Links make.linkakzeptiert.

Funktioniert family=binomial(link="identity")aber family=binomial(link=identity)nicht. (Wenn Sie etwas anderes finden, hat dies möglicherweise mit der R-Version zu tun.) Um eine Überstreuung zu ermöglichen, verwenden Sie family=quasi(link="identity", variance = "mu(1-mu)").

Scortchi - Monica wieder einsetzen
quelle
Der Link = "Identität" vs. Link = Identitätsfix war eine große Hilfe. Dies ist ein Training in Agrestis CDA-Lehrbuch. Der Code, den er bereitstellt, ist der Quasi (Link ...), den Sie diskutieren. Die Einfachheit des Hinzufügens von "" ist jedoch eine elegante Lösung. Nach meinem Verständnis repräsentiert der Aufruf link = "identity" das Binom als lineares Modell.
Justin Peterson