Hola a todos!
Tengo una duda acerca de una glm muy simple que estoy haciendo para
analizar bioensayos de mortalidad con
sedimento extra�do en cuatro estaciones de un puerto.
El modelo es el siguiente:
m1-glm (cbind (Mortalidad,Sobrev)~Muestra, family=binomial,data=dat2)
Hola, ¿qué tal?
Leyendo con más cuidado me he dado cuenta de lo de los cero casos en
control. Tu problema se puede reproducir con un ejemplo como
dat - data.frame(a = c(10, 10, 10, 10), b = c(0,2,3,4), grupo = letters[1:4])
res - glm(cbind(a,b) ~ grupo, data = dat, family = binomial)
Hola, ¿qué tal?
¿Cómo son tus datos? Para un modelo como el que planteas, deberían
tener cuatro filas (una para cada grupo) y en tal caso, sospecho, los
grados de libertad menores. ¿Tienes varias filas por grupo con
proporciones de mortalidad muy desiguales? ¿Existen otras variables