Trazado de proporciones binomiales en un gráfico de cajas y bigotes con R

Question

Disculpas de antemano - soy un poco novato con R pero he estado buscando en Google durante semanas y no he encontrado una respuesta directa.

El problema

Básicamente, estoy buscando trazar proporciones con sus intervalos de confianza en un gráfico de caja y bigotes utilizando el entorno R.

es decir, trace los resultados de binom.test(100,1000) y binom.test(125,1000)

¿Cuál es el mejor paquete (no tengo ninguna preferencia) y la mejor manera de hacerlo?

Answer 1

Esto parece más apropiado para una pregunta etiquetada [r] en SO, pero parece haber un grado sorprendentemente amplio de tolerancia para este tipo de preguntas en CV, así que aquí va. La respuesta también asume que usted no quería realmente utilizar el existente boxplot() que es lo que se utilizaría para construir un verdadero gráfico de caja y bigotes (y que en realidad no utiliza los IC binomiales, sino que utiliza el Tukey definido por fivenum() función ). Los paquetes plotrix y gplots tienen lo que parece ser la misma función plotCI;

est1 <- binom.test(100,1000) ; est2 <- binom.test(125,1000)
estvals <- rbind( c(est1$estimate, est1$conf.int[1], est1$conf.int[2]),
                  c(est2$estimate, est2$conf.int[1], est2$conf.int[2]))
plotCI(x = c(1,2), estvals[,1], ui=estvals[,3], li=estvals[,2], 
         xlim=c(0.5,2.5), ylim=c(0, 0.15) , xaxt="n", xlab="Groups")

La función básica es la arrows() función de trazado con un argumento que divide los puntos de las flechas a 90 grados en ambas direcciones.