Dear Boris, Thanks very much, have a great week.
Best regards Antônio Olinto Fisheries Institute São Paulo, Brasil 2015-03-21 21:09 GMT-03:00 Boris Steipe <boris.ste...@utoronto.ca>: > ... just for completeness - the more concise way: (no need to go through > names()). > > boxplot(mydata[,order(apply(mydata,2,median))]) > > ... or descending > boxplot(mydata[,order(-apply(mydata,2,median))]) > > > > > B. > > > > On Mar 21, 2015, at 7:04 PM, Boris Steipe <boris.ste...@utoronto.ca> > wrote: > > > There may be more concise ways to do this - but you are 99% there with > your approach: > > try: > > > > boxplot(mydata[,names(sort(apply(mydata,2,median)))]) > > > > B. > > > > > > On Mar 21, 2015, at 6:49 PM, Antonio Silva <aolinto....@gmail.com> > wrote: > > > >> Thanks Bill and David > >> > >> Here goes an example > >> > >> SP1<-c(9,6,7,8,5,8,7,5,9,7) > >> SP2<-c(1,3,4,2,4,2,5,3,2,1) > >> SP3<-c(4,6,7,5,7,8,7,6,5,4) > >> SP4<-c(5,4,3,5,2,3,4,3,4,2) > >> mydata<-data.frame(SP1,SP2,SP3,SP4) > >> > rownames(mydata)<-c("ST1","ST2","ST3","ST4","ST5","ST6","ST7","ST8","ST9","ST10") > >> mydata > >> boxplot(mydata) > >> > >> Note that this data frame does not have the format Response ~ Group. In > my > >> real matrix I have up to 40 species. > >> > >> Is there any way to have the species ordered by their median abundance > (or > >> other parameter?) > >> > >> The desired order is given by names(sort(apply(mydata,2,median))) > >> > >> Thanks once more, > >> > >> Antonio Olinto > >> Fisheries Institute > >> Sao Paulo, Brazil > >> > >> 2015-03-21 15:13 GMT-03:00 William Dunlap <wdun...@tibco.com>: > >> > >>> You can use the reorder() function to reorder the grouping vector's > >>> factor levels according to a function of the data in each group. E.g., > >>> compare the following two plots: > >>> > >>> d <- data.frame(Response=cos(1:15), > Group=rep(c("A","B","C"),c(6,5,4))) > >>> par(mfrow=c(1,2)) > >>> boxplot(Response ~ Group, data=d) > >>> boxplot(Response ~ reorder(Group, X=Response, FUN=median), data=d) > >>> > >>> > >>> Bill Dunlap > >>> TIBCO Software > >>> wdunlap tibco.com > >>> > >>> On Fri, Mar 20, 2015 at 2:20 PM, Antonio Silva <aolinto....@gmail.com> > >>> wrote: > >>> > >>>> Hello > >>>> > >>>> I'm using a dataframe (mydata) where row names are sampling points and > >>>> column names are species in a multivariate analysis. > >>>> > >>>> If I write boxplot(mydata) I'll have boxplots for each species > abundance > >>>> in > >>>> alphabetical order. > >>>> > >>>> How to get the boxes orderer by the median? > >>>> > >>>> Usually for this I write > >>>> > >>>> > boxplot(value~category,at=rank(tapply(mydata$value,mydata$category,median))) > >>>> but for this to work I need a column for the categories and another > column > >>>> for the values. > >>>> > >>>> Thanks in advance, best regards. > >>>> > >>>> Antonio Olinto > >>>> > >>>> [[alternative HTML version deleted]] > >>>> > >>>> ______________________________________________ > >>>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > >>>> 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. > >>>> > >>> > >>> > >> > >> [[alternative HTML version deleted]] > >> > >> ______________________________________________ > >> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > >> 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@r-project.org mailing list -- To UNSUBSCRIBE and more, see > > 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. > > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.