Predicción inversa a partir de glm binomial: ¿cómo obtener intervalos de confianza?

Question

Estoy intentando utilizar un glm con una distribución binomial (enlace logit) para analizar los datos de una curva dosis-respuesta para los efectos letales de diferentes cepas bacterianas. Ahora quiero obtener estimaciones de la dosis letal 50 (LD50).

Mi modelo es simplemente de la forma

survival ~ a + b*dose

Así, como tengo un enlace logit, obtener las estimaciones de la DL50 (es decir, p=0,5) significa simplemente dividir -a/b. Sin embargo, también quiero conocer los intervalos de confianza para estas estimaciones y ahí es donde no estoy seguro de qué hacer. Si alguien tiene una sugerencia para esto, sería genial.

Answer 1

Hay un paquete "drc" para la estimación de la DL50 y sus intervalos de confianza.

De otra manera, sus resultados glm pueden ser tratados por una simple función "dose.p" en el paquete MASS que calcula la LD50 y su SE.