Ярлыки Ярлыки ящиков в R

Question

У меня есть код, который создает коробку, используя ggplot в R, я хочу обозначить мои выбросы годом и битвой.

Вот мой код, чтобы создать мой boxplot

require(ggplot2) 
ggplot(seabattle, aes(x=PortugesOutcome,y=RatioPort2Dutch),xlim="OutCome", 
y="Ratio of Portuguese to Dutch/British ships") + 
geom_boxplot(outlier.size=2,outlier.colour="green") + 
stat_summary(fun.y="mean", geom = "point", shape=23, size =3, fill="pink") + 
ggtitle("Portugese Sea Battles") 

Может кто-нибудь помочь? Я знал, что это правильно, я просто хочу обозначить выбросы.

Answer 1

Это работает для вас?

library(ggplot2) 
library(data.table) 

#generate some data 
set.seed(123) 
n=500 
dat <- data.table(group=c("A","B"),value=rnorm(n)) 

ggplot определяет выброс по умолчанию, как нечто такое, что> 1.5 * МКР от границ окна.

#function that takes in vector of data and a coefficient, 
#returns boolean vector if a certain point is an outlier or not 
check_outlier <- function(v, coef=1.5){ 
    quantiles <- quantile(v,probs=c(0.25,0.75)) 
    IQR <- quantiles[2]-quantiles[1] 
    res <- v < (quantiles[1]-coef*IQR)|v > (quantiles[2]+coef*IQR) 
    return(res) 
} 

#apply this to our data 
dat[,outlier:=check_outlier(value),by=group] 
dat[,label:=ifelse(outlier,"label","")] 

#plot 
ggplot(dat,aes(x=group,y=value))+geom_boxplot()+geom_text(aes(label=label),hjust=-0.3) 

Answer 2

Ниже воспроизводимое решение, которое использует dplyr и встроенный mtcars данных.

Прогулка по коду: сначала создайте функцию is_outlier, которая вернет логическое значение TRUE/FALSE, если переданное ему значение является выбросом. Затем мы выполняем «анализ/проверку» и построим данные - сначала мы перечислим нашу переменную (cyl в этом примере, в вашем примере это будет PortugesOutcome) и добавим переменную outlier в звонок mutate (если drat переменная является outlier [обратите внимание, что это соответствует RatioPort2Dutch в вашем примере], мы передадим значение drat, в противном случае мы вернем NA, чтобы значение не отображалось). Наконец, мы построим результаты и построим текстовые значения через geom_text и эстетическую метку, равную нашей новой переменной; кроме того, мы компенсируем текст (сдвиньте его немного вправо) с hjust, чтобы мы могли видеть значения рядом с точками выброса, а не сверху.

library(dplyr) 
library(ggplot2) 

is_outlier <- function(x) { 
    return(x < quantile(x, 0.25) - 1.5 * IQR(x) | x > quantile(x, 0.75) + 1.5 * IQR(x)) 
} 

mtcars %>% 
    group_by(cyl) %>% 
    mutate(outlier = ifelse(is_outlier(drat), drat, as.numeric(NA))) %>% 
    ggplot(., aes(x = factor(cyl), y = drat)) + 
    geom_boxplot() + 
    geom_text(aes(label = outlier), na.rm = TRUE, hjust = -0.3) 

Answer 3

С небольшим поворотом на @JasonAizkalns решении, которое вы можете маркировать выбросы с их расположением в кадре данных.

mtcars[,'row'] <- row(mtcars)[,1] 
... 
mutate(outlier = ifelse(is_outlier(drat), row, as.numeric(NA))) 
... 

загружает кадр данных в R студии окружающей среду, поэтому я могу взять более близкий взгляд на данных в строках резко отклоняющихся значения.

Answer 4

Чтобы пометить выбросы с rownames (на основе JasonAizkalns ответа)

library(dplyr) 
library(ggplot2) 
library(tibble) 

is_outlier <- function(x) { 
    return(x < quantile(x, 0.25) - 1.5 * IQR(x) | x > quantile(x, 0.75) + 1.5 * IQR(x)) 
} 

dat <- mtcars %>% tibble::rownames_to_column(var="outlier") %>% group_by(cyl) %>% mutate(is_outlier=ifelse(is_outlier(drat), drat, as.numeric(NA))) 
dat$outlier[which(is.na(dat$is_outlier))] <- as.numeric(NA) 

ggplot(dat, aes(y=drat, x=factor(cyl))) + geom_boxplot() + geom_text(aes(label=outlier),na.rm=TRUE,nudge_y=0.05)