On 7/22/06, Nair, Murlidharan T <[EMAIL PROTECTED]> wrote: > -----Original Message----- > From: Gabor Grothendieck [mailto:[EMAIL PROTECTED] > Sent: Saturday, July 22, 2006 1:37 PM > To: Nair, Murlidharan T > Cc: R-help@stat.math.ethz.ch > Subject: Re: [R] Multcomp > > On 7/22/06, Nair, Murlidharan T <[EMAIL PROTECTED]> wrote: > > Here it is again, hope this is more clear > > > > I am using the following data (only a small subset is given): > > > > Habitat Fungus.yield > > Birch 20.83829053 > > Birch 22.9718181 > > Birch 22.28216829 > > Birch 24.23136797 > > Birch 22.32147961 > > Birch 20.30783598 > > Oak 27.24047258 > > Oak 29.7730014 > > Oak 30.12608508 > > Oak 25.76088669 > > Oak 30.14750974 > > Hornbeam 17.05307949 > > Hornbeam 15.32805111 > > Hornbeam 18.26920177 > > Hornbeam 21.30987049 > > Hornbeam 21.7173223 > > > > I am using the multcomp package to do multiple comparisons as follows > > > > library(multcomp) # loads the package > > > > fungus<-read.table("fungi.txt", Header=T) # Reads the data from > file saved as fungi.txt > > > > > > fungus.cirec<-simint(Fungus.yield~Habitat, > data=fungus,conf.level=0.95,type =c("Tukey")) # Computes cimultaneous > intervals using Tukey's method > > > > > > plot(fungus.cirec) # plots the data > > > > The plot function plots all the comparisons, I want to plot only part > of the data since it clutters the graph. > > > > How do I plot only part of the data ? > > Don't understand what "part of the data" means. Use data = > fungus[1:10,] > in the simint call to just process the first 10 data rows. To > eliminate a portion of > the plot note in ?plot.hmtest that there is a ... argument and its > description > is that its passed to plot which in turn passes them to plot.default so > you > could use ylim = 1:2, say, to show only part of the plot vertically. > > > >>I have to use all the data for my computation. Since the number of > >>comparisons are many I want to plot it in different graphs so that the > >>graph does not look cluttered. So "part of the data" means a subset of > the >>comparisons in one graph and another subset in another and so > on.... > >> Can you give me an example of the "ylim" parameter. I tried > >>plot(fungus.cirec, ylim=1:2) and it was not happy with it. I am not > >>completely comfortable with the ... argument. I am trying to read up > on >>it. > > > > > > How do I tell it to mark the significant comparisons? > > # after your plot statement: > pp <- locator() > # now click on a spot on the plot > # and then right click and choose stop > text(pp$x, pp$y, "some text") > > >> This is very useful > > > > > How do I get rid of the field names in the plot? For eg. The plot > labels are HabitatBirch-HabitatOak, I want it to be labeled as > Birch-Oak. > > # change rownames of the estimates which has the effect > # of changing the y axis labels > rownames(fungus.cirec$estimate) <- LETTERS[1:3] > plot(fungus.cirec) > > >>This only labels it a A,B or C. What I wanted was to remove the field > name >>"Habitat" in this case and out put the label as Birch-Oak. >
See ?gsub ______________________________________________ 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.