R familia binomial con enlace de identidad

8

Quiero ajustar un modelo lineal por R con family=binomial(link="identity"), sin embargo, la familia binomial no tiene un enlace de identidad. ¿Qué tengo que hacer?

david
fuente
Creo que hay un problema estadístico subyacente aquí.
Glen_b -Reinstalar Monica
Sí, la siguiente pregunta le preguntará sobre Ajustar el error estándar para la sobredispersión.
David
Pero para el primero, necesito usar el enlace de identidad en la familia binomial, pero R no lo permite.
David
1
Tenga en cuenta que si necesita lidiar con la sobredispersión y utilizar un enlace de identidad, debería considerar ir directamente a un cuasimodelo con función de varianza binomial. Se puede instalar a mano un modelo binomial de solo intercepción .
Glen_b -Reinstalar Monica
2
(Y si se quería que el parámetro de dispersión fija en un: family=binomial(link=make.link("identity")).
Scortchi - Restablecer Mónica

Respuestas:

9

Consulte Wikipedia sobre el modelo de probabilidad lineal y publicaciones de CV aquí y aquí para obtener información estadística. Aunque no está "equivocado", querría una buena razón para usar un enlace de identidad para modelar una probabilidad de Bernoulli.

De acuerdo con el familymanual

la familia binomial [acepta] los enlaces logit, probit, cauchit, (correspondiente a la logística, normal y Cauchy CDF respectivamente) log y cloglog(complementario log-log)

Pero

Los argumentos de enlace y varianza tienen una semántica bastante incómoda para la compatibilidad con versiones anteriores. La forma recomendada es proporcionarlos como cadenas de caracteres entre comillas, pero también se pueden suministrar sin comillas (como nombres o expresiones). Además, también se pueden suministrar como un vector de caracteres de longitud uno con el nombre de una de las opciones, o como una lista (para link, de clase "link-glm"). Las restricciones se aplican solo a los enlaces dados como nombres: cuando se dan como una cadena de caracteres, make.linkse aceptan todos los enlaces conocidos .

Entonces family=binomial(link="identity")funciona pero family=binomial(link=identity)no lo hace. (Si encuentra algo diferente, podría tener que ver con la versión R). Para permitir una dispersión excesiva, utilice family=quasi(link="identity", variance = "mu(1-mu)").

Scortchi - Restablece a Monica
fuente
el enlace = "identidad" frente a enlace = corrección de identidad fue de gran ayuda. Este es un entrenamiento en el libro de texto CDA de Agresti. El código que proporciona es el cuasi (enlace ...) que discute, sin embargo, la simplicidad de agregar "" es una solución elegante. A mi entender, la llamada link = "identity" representa el binomio como un modelo lineal.
Justin Peterson el