2015-04-05 40 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