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)
summary(m1)
Call:
glm(formula = cbind (Mortalidad, Sobrev) ~ Muestra, family =
binomial,
data = dat2)
Deviance Residuals:
Min
1Q Median 3Q
Max
-2.80372 -0.44368 -0.00019 0.58488 1.78081
Coefficients:
Estimate Std. Error z value
Pr(>|z|)
(Intercept) -20.12 2587.87 -0.008
0.994
MuestraB1 19.39 2587.87 0.007
0.994
MuestraB2 19.27 2587.87 0.007
0.994
MuestraC
17.92 2587.87 0.007 0.994
(Dispersion parameter for
binomial family taken to be 1)
Null
deviance: 41.324 on 14 degrees of freedom
Residual deviance: 19.340
on 11 degrees of freedom
AIC: 51.564
Number of Fisher Scoring
iterations: 17
El control (intercept) tiene mortalidad cero, lo cual explica el resultado
err�neo. Me aconsejaron utilizar un aproximaci�n num�rica exacta. Alguien sabe
si hay alg�n paquete en R que lo haga?
SaludosMatias
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
R-help-es@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es
