2015-04-05 124 views
8

极端值的指示,我在我的数据集的一些非常,非常少数离群作出困难的箱线图看:包括在ggplot

library(ggplot2) 
mtcars$mpg[1] <- 60 
p <- ggplot(mtcars, aes(factor(cyl), mpg)) 
p + geom_boxplot() 

enter image description here

因此,我想指出的极端值像这样:

enter image description here

任何想法如何做到这一点的ggplot2?改造轴不是我的选择...

回答

8

这是一个开始:

library("ggplot2") 
mtcars$mpg[1:2] <- c(50,60) 
p <- ggplot(mtcars, aes(factor(cyl), mpg)) 
p + geom_boxplot() 

定义最大值:

maxval <- 40 

使用dplyr(也可以在基础R完成或plyr)来提取异常值并将文本字符串放在一起:

library("dplyr") 
dd <- mtcars %>% filter(mpg>maxval) %>% 
    group_by(cyl) %>% 
     summarise(outlier_txt=paste(mpg,collapse=",")) 

Set ma X Y值和添加的箭头加标签:

library("grid") # needed for arrow() function 
p2 <- p + geom_boxplot() + 
    scale_y_continuous(limits=c(min(mtcars$mpg),maxval))+ 
     geom_text(data=dd,aes(y=maxval,label=outlier_txt), 
       size=3,vjust=1.5,hjust=-0.5)+ 
      geom_segment(data=dd,aes(y=maxval*0.95,yend=maxval, 
         xend=factor(cyl)), 
       arrow = arrow(length = unit(0.1,"cm"))) 
p2 

enter image description here