Hello Sarah, Thanks for the suggestion of table(obs.group).
I checked it and found no problems of the obs.group. As for using dput(), please kindly share some examples to display the colored graph in r-help. In fact, some colors were assigned to all types of diets, like A=>red1 B=>red2 C=>green1 D=>green4 E=>darkgreen F=>orange3 G=>darkorange2 H=>navy I=>blue1 However, the boxplot displays as A=>red1 B=>red2 C=>green1 D=>green4 E=>red1 F=>red2 G=>green1 H=>green4 I=>red1. I am wondering if there is any bug in showing too many colors for boxplot. Elaine On Tue, Sep 25, 2012 at 6:52 AM, Sarah Goslee <sarah.gos...@gmail.com>wrote: > Hi Elaine, > > Without a reproducible example it's impossible to say, but I'd take a > hard look at: > table(obs.group) > > If that doesn't give you some insight, a small reproducible example > included in your email using dput() would allow us to answer you more > effectively. > > Sarah > > > On Mon, Sep 24, 2012 at 6:23 PM, Elaine Kuo <elaine.kuo...@gmail.com> > wrote: > > Hello, > > > > I am making a boxplot of 13 boxes. > > I tried to color the box using 13 colors but failed. > > Only red and brown were displayed. > > Green, blue, and grey disappeared. > > > > Please kindly advise modification after checking the code below. > > Thank you in advance. > > > > Elaine > > > > R code > > # data input > > dataN > > > <-read.csv("H:/a_mig_distance_B_NB/R_data/Mig_bird_586_20120925.csv",header=T, > > row.names=1) > > dim(dataN) > > dataN[1,] > > str(dataN) > > > > # graph > > par(mai=c(1,1.03,0.4,0.4)) > > > > obs.group<-dataN$Diet_B > > > > par(new=T) > > > > boxplot(GE_distance~Diet_B, data=dataN,xlab="Diet of Breeding > > Ground",ylab="", > > yaxt="n",type="p", > > pch=1,lwd=0.95, > > cex.lab=1.4, cex.axis=1.2, > > font.axis=2, > > cex=1.5, > > las=1, > > bty="l", > > col=c("forestgreen", > "darkgreen","chocolate1","darkorange2","sienna2", > > "red2","firebrick3", > > > > > "saddlebrown","coral4","chocolate4","darkblue","navy","grey38")[obs.group])) > > > > op = par(mar = c(5,5,4,2) + 0.1) > > title(ylab = "Distance between Centers of B and NB Range (1000 km)", > > cex.lab = 1.4,line = 3) > > > > axis(side=2,yaxp=c(0,20000,4),cex.lab=1.4, cex.axis=1.2,font.axis=2, > > las=1) > > -- > Sarah Goslee > http://www.functionaldiversity.org > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org 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.