I use to do this (It may help you) >boxplot(x) >bx<-boxplot(x,plot=F)
In bx object, you and object names out with yaxis cordinates. For the i-th boxplot in the figure, use >text(rep(i,length(bx$out)),bx$out,...) Justin BEM Elève Ingénieur Statisticien Economiste BP 294 Yaoundé. Tél (00237)9597295. ----- Message d'origine ---- De : Greg Snow <[EMAIL PROTECTED]> À : [EMAIL PROTECTED]; r-help@stat.math.ethz.ch Envoyé le : Mardi, 16 Janvier 2007, 21h35mn 43s Objet : Re: [R] labels outliers in boxplot Some information on the outliers is returned from the boxplot function. Try something like: set.seed(123) tmp <- data.frame( group=gl(3,10), y=rcauchy(30), sex=gl(2,5,30,c('M','F')) ) tmp2 <- boxplot( split(tmp$y,tmp$group) ) identify( tmp2$group, tmp2$out, tmp2$group ) Or if your grouping variable works out to the same values as the x axis (used integers) and you want to specify a 3rd variable to be the labels you can do: Identify( tmp$group, tmp$y, tmp$sex ) Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 > -----Original Message----- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of > [EMAIL PROTECTED] > Sent: Wednesday, January 10, 2007 2:08 PM > To: r-help@stat.math.ethz.ch > Subject: [R] labels outliers in boxplot > > Dear R-users, > > Following is part of my data, where slide has 36 levels and > block 48 levels. I have done boxplot for each slide on the > same graph. There are outliers for each slide and I tried to > use indentify functtion to identify outliers in such a way > that when I click on an outlier or point, the points will be > labelled by either their block or ID or by both but without > success. How can I make it work or are there other ways to do > it than using identify function? > > > Thanks in advance, > > Jenny, > > dat1[1:10,] > y Slide Block ID Control > 1 0.03147823 1 1 IgG-human 5 > 2 -0.23815974 1 1 LPPAANDVSVLTAAR 0 > 3 -0.71926359 1 1 HTKHYRVVSKPAALV 0 > 4 -0.14607826 1 1 FVALPAATADAYATT 0 > 5 0.89553073 1 1 NYPAMMAHAGDMAGY 0 > 6 -0.67587100 1 1 RRALRQIGVLERPVG 0 > 7 0.32636034 1 1 DCGTIRVGSFRGRWL 0 > 8 -1.44057259 1 1 MAKLSTDELLDAFKE 0 > 9 -0.37064338 1 1 LELSDFVKKFEETFE 0 > 10 -0.20387233 1 1 VSRRAKVDVLIVHTT 0 > > > tb_ncs<-subset(dat1,dat1$Control==1) ### this > data contains only negative controls > > par(las=2,mar=c(10.1,4.1,4.1,2.1)) > > boxplot(split(tb_ncs$y,tb_ncs$Slide),col="orange", > cex=.65, > outline=TRUE,main="Negative control > response of each patient", cex.main=1, font.main=1, > col.main="blue", > names=c(1:35,"B")) > grid(nx=NA, ny=NULL) > ### grid over boxplot > legend("bottomright", "B = Buffer + > sec",text.col="blue") > > > out.block<- > identify(tb_ncs$y,tb_ncs$Slide) > > > > > > > _________________________________________________________ > Flyger tiden iväg? Fånga dagen med Yahoo! Mails inbyggda > > ______________________________________________ > R-help@stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. ___________________________________________________________________________ Découvrez une nouvelle façon d'obtenir des réponses à toutes vos questions ! Profitez des connaissances, des opinions et des expériences des internaut [[alternative HTML version deleted]]
______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.