Isa García Barón
2014-Nov-17 11:25 UTC
[R-es] Cómo aplicar weights a las observaciones en un GLM binomial
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]]