[R-es] Duda sobre GLM

[Matias Ledesma]

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]]