Re: [R] controlling the x axis of boxplots

From: Rich FitzJohn <rich.fitzjohn_at_gmail.com>
Date: Wed 20 Apr 2005 - 07:38:53 EST

Hi Chris,

You can get the desired effect by using the "at" argument to boxplot, and by setting up the plot dimensions manually.

men <- data.frame(grp=rep(letters[c(1, 2, 4, 6)], each=10),

                  response=rnorm(40))
women <- data.frame(grp=rep(letters[c(2:5, 7)], each=10),
                    response=rnorm(50))

## Determine all levels used, the number of levels, and appropriate ## xlim for the plot
grp.levs <- sort(unique(c(levels(men$grp), levels(women$grp)))) nlevs <- length(grp.levs)
xlim <- c(.5, nlevs+.5)

## Determine which of the levels are present in men and women; these
## are the x coordinates to draw the boxes (this will need tweaking if
## you are using ordered factors).

at.men <- match(levels(men$grp), grp.levs) at.women <- match(levels(women$grp), grp.levs)

par(mfrow=c(2,1))

## boxplot() does not take an xlim argument, so you'll have to set up
## the plot yourself (including the axis, since the default is not
## appropriate for a boxplot).
plot(NA, type="n", xlim=xlim, ylim=range(men$response), xlab="",
     ylab="Response", xaxt="n")

axis(1, 1:nlevs, grp.levs)
boxplot(response ~ grp, men, at=at.men, add=TRUE)

plot(NA, type="n", xlim=xlim, ylim=range(women$response),

     xlab="Group", ylab="Response", xaxt="n") axis(1, 1:nlevs, grp.levs)
boxplot(response ~ grp, women, at=at.women, add=TRUE)

Cheers,
Rich

On 4/20/05, Chris Evans <chris@psyctc.org> wrote:
> v 2.0.1 (sooooh old!) on Win2k
>
> I think I know the answer to this but I can hope ...
>
> I have data for continuous variables (measures of residents) by a
> categorical variable in range (1,22), the units in which they live.
>
> I want to plot these data with a pair of boxplots one above another
> with same x-axis (1,22) using par(mfrow=c(2,1)) and then plotting
> first for the women then for the men. My problem is that some units
> have only men, some have only women, and some have both. I'd like
> both plots to have the same x axis and the notched, varwidth boxplots
> to locate themselves with some gaps so that the units are in the same
> place on the x axis on each plot.
>
> I think that I can't do this with boxplot or bxp as both work out the
> x axis from the work that boxplot has done on the data. Although
> there also seem to be useful extensions or alternative boxplots in
> other packages, I can't see one that does what I want and I think
> that rolling my own from bxp is really beyond my skills.
>
> Am I right that it doesn't exist in the CRAN packages? If not,
> apologies and point me where I should look? If I am right (boo hoo!)
> I don't suppose anyone has written this or is feeling like
> demonstrating their personal genius with R coding?!!! If they were,
> I don't think I'd be the only one to end up owing them a great deal
> of gratitude!
>
> Thanks as ever to all who have made and continue to make R what it
> is: brilliant!
>
> Chris
>
> --
> Chris Evans <chris@psyctc.org>
> Consultant Psychiatrist in Psychotherapy, Rampton Hospital;
> Research Programmes Director, Nottinghamshire NHS Trust,
> Hon. SL Institute of Psychiatry
> *** My views are my own and not representative of those institutions
> ***
>
> ______________________________________________
> 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
>

-- 
Rich FitzJohn
rich.fitzjohn <at> gmail.com   |    http://homepages.paradise.net.nz/richa183
                      You are in a maze of twisty little functions, all alike

______________________________________________
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
Received on Wed Apr 20 07:48:05 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:31:17 EST