[R-es] Comparación entre dos DL50's

[Argel Gastélum Arellánez]

     Hola compañeros de la lista.

     Tengo el análisis que muestro abajo, para la estimación de la Dosis 
Letal Media (DL50) por medio de regresión logística, de dos compuestos 
similares (denominados 1 y 2) aplicados a grupos de individuos 
similares, evaluando el número de individuos vivos y muertos al final 
del ensayo. Mi duda es qué forma me recomiendan para determinar si la 
diferencia es significativa entre las DL50's estimadas para cada compuesto.

     De antemano muchas gracias por su ayuda.

----------------------------------------------------------------------
Código:

datos.1 = data.frame(
     dosis.1 = c(1,10,10,10,50,50,50,100,100,100,1000,1000,1000),
     vivos.1 = c(20,20,20,20,19,14,15,5,3,2,0,0,0),
     muertos.1 = c(0,0,0,0,1,6,5,15,17,18,20,20,20)
)

datos.2 = data.frame(
     dosis.2 = c(1,10,10,10,50,50,50,100,100,100,1000,1000,1000),
     vivos.2 = c(20,18,18,18,7,6,6,1,1,2,0,0,0),
     muertos.2 = c(0,2,2,2,13,14,14,19,19,18,20,20,20)
)

attach(datos.1)

attach(datos.2)

y1 = cbind(vivos.1,muertos.1)

y2 = cbind(vivos.2,muertos.2)

modelo.logistico.1 = glm(y1~log(dosis.1),binomial)

modelo.logistico.2 = glm(y2~log(dosis.2),binomial)

library(MASS)

log.LD50.1 <- dose.p(modelo.logistico.1,p=0.5)

log.LD50.2 <- dose.p(modelo.logistico.2,p=0.5)

log.LD50.1

log.LD50.2

LD50.1 <- exp(log.LD50.1[1])

LD50.2 <- exp(log.LD50.2[1])

LD50.1

LD50.2

-- 

     Argel.