Antes de mais nada vou recriar os dados com uma `seed` para que os resultados possamos ser reproduzidos. library(tidyverse) set.seed(123) dataset <- as_tibble(matrix(rnorm(6*1000,1500,200),ncol=6)) cluster <- kmeans(dataset,centers=3) dataset$kmeans <- as.factor(cluster[['cluster']]) e nomear o gráfico para que possamos usá-lo adiante. p <- dataset%>% gather(.,key='group',value='var',V1:V6) >% ggplot(aes(group,var,fill=kmeans))+ facet_grid(kmeans~.)+ geom_boxplot(outlier.color='yellow',outlier.shape=21, outlier.fill='black',outlier.size=1)+ stat_summary(fun.y=mean, geom="line", aes(group=1))+ stat_summary(fun.y=mean, geom="point")+ theme_dark() ## Obtendo _ouliers_ no `ggplot2` No [tag:ggplot2] os outliers são calculados no momento em que o gráfico vai ser desenhado (quando der `print(plot)`). A função responsável por isso será a função `StatBoxplot` ([código-fonte](https://github.com/tidyverse/ggplot2/blob/master/R/stat-boxplot.r)) que é executada dentro do "ambiente" do gráfico. É possível forçar a execução dos cálculos por meio da função `ggplot_build()`. Assim temos calculado <- ggplot_build(p) O objeto criado tem uma elemento chamado `data` que possui três `data.frames`. O primeiro deles se refere aos dados usados para a construção do __boxplot__. Neste `data.frame` há uma lista com os _outliers_ de cada boxplot chamada "outliers". retiramos dasd outliers <- calculado$data[[1]]$outliers str(outliers) List of 18 $ : num 2148 $ : num [1:4] 2059 2066 1191 2158 $ : num [1:5] 992 1914 930 973 1819 $ : num [1:3] 910 932 874 $ : num [1:7] 2094 1163 2154 2085 2189 ... $ : num [1:3] 966 2072 2117 $ : num(0) $ : num [1:2] 961 890 $ : num [1:5] 1036 998 1940 2047 988 $ : num [1:2] 916 2079 $ : num [1:6] 949 2000 922 1946 1909 ... $ : num [1:6] 1065 2038 2050 1031 1085 ... $ : num [1:2] 968 2015 $ : num [1:3] 2137 2178 2063 $ : num [1:4] 2104 1262 1259 2184 $ : num [1:3] 1003 1975 1931 $ : num 891 $ : num [1:3] 862 848 1962 # Obtendo os enésimos _outliers_ primeiros <- function(outliers, lower, upper) { primeiro_baixo <- max(outliers[outliers < lower]) primeiro_cima <- min(outliers[outliers > upper]) c(primeiro_baixo, primeiro_cima) %>% map_dbl(~ ifelse(is.infinite(.x), NA, .x)) } dados <- calculado$data[[1]] pmap(list(dados$outliers, dados$lower, dados$upper), primeiros) # Usando o pacote `graphics` O pacote `graphics` já vem no `R` (é o pacote que possui a função `plot`, por exemplo) e tem uma função que pode facilitar essa tarefa. res <- boxplot(var ~ group, gather(dataset, key = 'group', value = 'var', V1:V6)) res$out # [1] 2148.2080 971.3702 967.8154 938.0451 979.6601 2038.3428 2036.9718 ...