Oznaczanie wartości odstających na boxplocie w R

Question

Chciałbym wykreślić każdą kolumnę macierzy jako planszę, a następnie oznaczyć wartości odstające w każdym polu w polu jako nazwę wiersza, do której należą w macierzy. By posłużyć się przykładem:

vv=matrix(c(1,2,3,4,8,15,30),nrow=7,ncol=4,byrow=F) 
rownames(vv)=c("one","two","three","four","five","six","seven") 
boxplot(vv) 

chciałbym oznaczyć odrębne stanowisko w każdej działki (w tym przypadku 30), co nazwa rzędzie należy, więc w tym przypadku 30 należący do rzędu 7. Czy istnieje prosty sposób to zrobić? Widziałem podobne pytania do tego pytania, ale żaden z nich nie działał tak, jak tego chcę.

Answer 1

W podanym przykładzie jest to trochę nudne, ponieważ wszystkie są w tym samym rzędzie. ale tutaj jest kod:

bxpdat <- boxplot(vv) 
text(bxpdat$group,            # the x locations 
    bxpdat$out,            # the y values 
    rownames(vv)[which(vv == bxpdat$out, arr.ind=TRUE)[, 1]], # the labels 
    pos = 4) 

ten wybiera się rownames które mają wartości równe z „zewnątrz” listy (tj odstających) w wyniku boxplot. Boxplot wywołuje i zwraca wartości z boxplot.stats. Spójrz na:

str(bxpdat) 

Answer 2

@ rozwiązanie Dwin za działa bardzo dobrze dla pojedynczego wykresu pudełkowego, ale nie na wszystko z dwóch wartości, jak zbiorze mam Utworzono:

#Create data 
set.seed(1) 
basenums <- c(1,2,3,4,8,15,30) 
vv=matrix(c(basenums, sample(basenums), 1-basenums, 
      c(0, 29, 30, 31, 32, 33, 60)),nrow=7,ncol=4,byrow=F) 
dimnames(vv)=list(c("one","two","three","four","five","six","seven"), 1:4) 

Na tym zbiorze, @ rozwiązanie Dwin daje:

Co jest fałszywe, bo w 4. Przykładowo, nie jest możliwe do minimum i maksimum, które należy w t w tym samym rzędzie.

To rozwiązanie jest monstrualne (i mam nadzieję, że można je uprościć), ale skuteczne.

#Reshape data 
vv_dat <- as.data.frame(vv) 
vv_dat$row <- row.names(vv_dat) 
library(reshape2) 
new_vv <- melt(vv_dat, id.vars="row") 

#Get boxplot data 
bxpdat <- as.data.frame(boxplot(value~variable, data=new_vv)[c("out", "group")]) 

#Get matches with boxplot data 
text_guide <- do.call(rbind, apply(bxpdat, 1, 
    function(x) new_vv[new_vv$value==x[1]&new_vv$variable==x[2], ])) 

#Add labels 
with(text_guide, text(x=as.numeric(variable)+0.2, y=value, labels=row)) 

Answer 3

Albo można po prostu uruchomić kod z this blog post:

source("https://raw.githubusercontent.com/talgalili/R-code-snippets/master/boxplot.with.outlier.label.r") # Load the function 
set.seed(6484) 
y <- rnorm(20) 
x1 <- sample(letters[1:2], 20,T) 
lab_y <- sample(letters, 20) 
# plot a boxplot with interactions: 
boxplot.with.outlier.label(y~x1, lab_y) 

(który obsługuje wiele odstających, które są blisko siebie)

Answer 4

@ sebastian-c Jest to niewielka modyfikacja rozwiązania Dwin które wydają się działać z większą ogólność

bx1<-boxplot(pb,las=2,cex.axis=.8) 
if(length(bx1$out)!=0){ 
    ## get the row of each outlier 
    out.rows<-sapply(1:length(bx1$out),function(i) which(vv[,bx1$group[i]]==bx1$out[i])) 
    text(bx1$group,bx1$out, 
    rownames(vv)[out.rows], 
    pos=4 
) 
} 

Answer 5

Lub alternatywnie, można użyć funkcji „boxplot” z samochodu} {opakowaniu, które etykiety odstających dla Ciebie.

Patrz poniższy link: http://www.inside-r.org/packages/cran/car/docs/Boxplot

wiele łatwiejsze!

Answer 6

Istnieje prosty sposób. Zauważ, że b w Boxplot w kolejnych wierszach jest wielką literą.

library(car) 

Boxplot(y ~ x, id.method="y")