Kennzeichnung Ausreißer von Boxplots in R

Habe ich den code, der erzeugt einen boxplot mit ggplot in R, ich möchte mein label "Ausreißer" mit dem Jahr und der Schlacht.

Hier ist mein code, um meine 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")

Kann jemand helfen? Ich wusste, dass dies richtig ist, ich will nur den label der Ausreißer.

  • Was passiert mit den Daten seabattle kommen? Können Sie dput die Daten oder geben Sie Beispieldaten, um in diesem Beispiel reproduzierbar?
  • Alles, was du schon versucht hast?
InformationsquelleAutor Deborah_Watson | 2015-11-04
  1. 25

    Folgende ist eine reproduzierbare Lösung, verwendet dplyr und die built-in mtcars dataset.

    Fuß durch den code: erstellen Sie Zunächst eine Funktion, is_outlier wird, geben einen booleschen Wert zurück TRUE/FALSE wenn der übergebene Wert ist ein Ausreisser. Wir führen Sie dann den "Analyse/Kontrolle" und eine Darstellung der Daten -- zuerst haben wir group_by unsere variable (cyl in diesem Beispiel, in deinem Beispiel wäre dies PortugesOutcome) und fügen wir eine variable outlier im Aufruf mutate (wenn die drat variable ist ein Ausreißer [Hinweis: dies entspricht RatioPort2Dutch in deinem Beispiel], geben wir die drat Wert, sonst werden wir zurück NA so, dass Wert ist nicht gezeichnet). Endlich, wir zeichnen die Ergebnisse und zeichnen Sie die text-Werte über geom_text - und eine ästhetische label gleich unsere neue variable; darüber hinaus haben wir den offset der text (schieben Sie es etwas nach rechts) mit hjust so dass wir sehen können, die Werte weiter zu, anstatt auf der Oberseite, die Ausreißer-Punkte.

    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)

    Kennzeichnung Ausreißer von Boxplots in R

    InformationsquelleAutor JasonAizkalns
  2. 7

    Macht diese Arbeit für Sie?

    library(ggplot2)
library(data.table)

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

    ggplot definiert einen Ausreißer standardmäßig als etwas das ist > 1.5*IQR von den Grenzen der box.

    #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)

    Kennzeichnung Ausreißer von Boxplots in R

    InformationsquelleAutor Heroka
  3. 6

    Zur Beschriftung der Ausreißer mit rownames (basierend auf JasonAizkalns Antwort)

    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)

    Kennzeichnung Ausreißer von Boxplots in R

    InformationsquelleAutor qfazille
  4. 3

    Ähnliche Antwort zu oben, wird aber Ausreißer direkt aus ggplot2, damit Vermeidung von potenziellen Konflikten in der Methode:

    # calculate boxplot object
g <- ggplot(mtcars, aes(factor(cyl), drat)) + geom_boxplot()

# get list of outliers 
out <- ggplot_build(g)[["data"]][[1]][["outliers"]]

# label list elements with factor levels
names(out) <- levels(factor(mtcars$cyl))

# convert to tidy data
tidyout <- purrr::map_df(out, tibble::as_tibble, .id = "cyl")

# plot boxplots with labels
g + geom_text(data = tidyout, aes(cyl, value, label = value), 
              hjust = -.3)

    Kennzeichnung Ausreißer von Boxplots in R

    InformationsquelleAutor jtr13
  5. 0

    Können Sie dies einfach innerhalb von ggplot sich, mit einem geeigneten stat_summary nennen.

    ggplot(mtcars, aes(x = factor(cyl), y = drat, fill = factor(cyl))) + 
  geom_boxplot() +
  stat_summary(
    aes(label = round(stat(y), 1)),
    geom = "text", 
    fun.y = function(y) { o <- boxplot.stats(y)$out; if(length(o) == 0) NA else o },
    hjust = -1
  )

    Kennzeichnung Ausreißer von Boxplots in R

    InformationsquelleAutor Axeman
  6. -2

    Mit einem kleinen twist auf @JasonAizkalns Lösung können Sie label Ausreißer mit Ihrem Standort in Ihrem data frame.

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

    Lade ich die Daten-Frames in der R-Studio-Umgebung, kann ich dann also nehmen Sie einen genaueren Blick auf die Daten Ausreißer Zeilen.

    InformationsquelleAutor earcanal
Schreibe einen Kommentar