hacer - Ponga puntas de bigotes en el cuadro de caja

graficos de cajas y bigotes en r

Question

Me gustaría poner líneas perpendiculares en los extremos de los bigotes como la función de boxplot da automáticamente.

Answer 1

Como se sugirió pero no implementado por @Roland, puede usar stat_boxplot para implementar esto. El truco es llamar a _boxplot dos veces y establecer el geom en la errorbar de errorbar para una de las llamadas.

Tenga en cuenta que a medida que R utiliza un enfoque de pluma y papel, es aconsejable implementar las barras de error primero, dibuje el cuadro de caja tradicional sobre la parte superior.

Usando los datos ficticios de @Roland df

ggplot(df, aes(x=cond, y = value)) +
stat_boxplot(geom =''errorbar'') +
geom_boxplot() # shorthand for stat_boxplot(geom=''boxplot'')

La ayuda para stat_boxplot ( ?stat_boxplot ) detalla los diversos valores calculados y guardados en un data.frame

Answer 2

Para cambiar el tamaño de las líneas de bigotes, podemos usar el argumento width = 0.5 dentro de la función: stat_boxplot

set.seed(42)
df <- data.frame(cond = factor(rep(c("A", "B"), each = 500)),
value = c(rnorm(500, mean = 1, sd = 0.2),
rnorm(500, mean = 1.5, sd = 0.1)))
library(ggplot2)
ggplot(df, aes(x = cond, y = value)) +
stat_boxplot(geom = "errorbar", width = 0.5) +
geom_boxplot()

Answer 3

Podría ser posible usar stat_boxplot para calcular los extremos del bigote, pero no soy lo suficientemente un asistente ggplot2 , entonces uso la función base para eso.

set.seed(42)
df <- data.frame(cond = factor( rep(c("A","B"), each=500) ),
value = c(rnorm(500,mean=1,sd=0.2),rnorm(500, mean=1.5,sd=0.1)))
whisk <- function(df,cond_col=1,val_col=2) {
require(reshape2)
condname <- names(df)[cond_col]
names(df)[cond_col] <- "cond"
names(df)[val_col] <- "value"
b <- boxplot(value~cond,data=df,plot=FALSE)
df2 <- cbind(as.data.frame(b$stats),c("min","lq","m","uq","max"))
names(df2) <- c(levels(df$cond),"pos")
df2 <- melt(df2,id="pos",variable.name="cond")
df2 <- dcast(df2,cond~pos)
names(df2)[1] <- condname
df2
}
library(ggplot2)
plot1 <- ggplot(df, aes(x=cond))
plot1 <- plot1 + geom_errorbar(aes(ymin=min,ymax=max),data=whisk(df),width = 0.5)
plot1 <- plot1 + geom_boxplot(aes(y=value))
plot1