[R-es] Cómo aplicar weights a las observaciones en un GLM binomial

[Isa García Barón]

Muchas gracias por todas las respuestas ayudando a resolver el problema,
finalmente he encontrado la solución de la siguiente manera:

#El código asume que mi dataset es un data.frame llamado mydata
#y la variable respuesta (0 o 1) se encuentra en la columna presencia
#número de 1
n1s=sum(mydata$presence)
#número de datos
n=nrow(mydata)
#Obtención de los pesos inversamente proporcionales a los datos
#por ejemplo, mayores pesos a los 1 y menores pesos a los 0
mydata$myweights=ifelse(mydata$presence==1,(n-n1s)/n,n1s/n)
#realización del modelo
myglm=glm(presence~covariate1+...+covariatek,family="
binomial",weights=myweights,data=mydata)
#y listo, visualización de los resultados
summary(myglm)
plot(myglm)


Ahí lo dejo por si a alguien le sirviera en un futuro,

Un saludo

Isa
-- 
*-------------*
*Isabel García Barón*

	[[alternative HTML version deleted]]