[R] Boxplot position on X-axis relative to it's value

[Felix Andrews]

Dennis, this works for me:

xyplot(nmgml ~ day, data = girafe, horizontal = FALSE,
        panel = function(...) {
          panel.bwplot(...)
          panel.average(..., fun = median, col.line = 'gray30')
         }
  )


On 13 May 2010 06:29, Dennis Murphy <djmuser at gmail.com> wrote:
> Hi:
>
> Here are a couple of ways to superimpose the boxplots with a line, using
> lattice and latticeExtra.
> Notice that the boxplots are located at the observed x values - the idea
> comes straight out of the
> example on p. 183 of the Lattice book.
>
> library(lattice)
> library(latticeExtra)
>
> # Connect the medians
> xyplot(nmgml ~ day, data = girafe,
>       panel = panel.bwplot, horizontal = FALSE) +
>  layer(panel.average(..., fun = median, col.line = 'gray30'))
>
> # Least squares regression line
> xyplot(nmgml ~ day, data = girafe,
>       panel = panel.bwplot, horizontal = FALSE) +
>       layer(panel.lmline(...))
>
> Caution: be careful if you have the ggplot2 package loaded, because the
> layer function
> in ggplot2 can mask the layer function in latticeExtra.
>
> ------------
> Aside:
> I couldn't figure out how to get the horizontal = FALSE argument to work
> within a
> single panel function that contains both panel.bwplot and panel.average. I
> tried
>
> xyplot(nmgml ~ day, data = girafe,
>         panel = function(...) {
>           panel.bwplot(..., horizontal = FALSE)
>           panel.average(..., fun = median, col.line = 'gray30')
>          }
>   )
>
> and a couple of variants thereof but couldn't get it to work. If there's an
> easy
> solution, I'd be interested in seeing it.
>
>
> HTH,
> Dennis
>
> On Wed, May 12, 2010 at 6:30 AM, Mächler Marc Jaques
> <marcm at student.ethz.ch>wrote:
>
>> Dear R-Experts.
>>
>> I collected different datas about Nitrogen content (mg/ml) in Dung. The
>> dung was eighter fresh (day=0) or had different ages (15,29,47) to observe
>> nutrient changes over time.
>>
>> Now I like to draw a boxplot.
>>
>> boxplot(nmgml~day)
>> abline((nmgml~day)
>>
>> The Problem is, that the boxplot considers the day values as groups and not
>> as time series (neighter when the days are numeric or factors).
>>
>> The result is, that the regression line looks different as when I draw a
>>  plot
>> plot(nmgml~day)
>> abline(nmgml~day)
>>
>> How can I reach that the boxplot's position on the x-axis is relative to
>> the time value? (and the days are not just treated as groups)
>>
>> Data sample (Girafe-Dung)
>>   day   nmgml
>> 1   29 2.72655
>> 2   29 2.48820
>> 3   15 2.85945
>> 4   15 2.58915
>> 5   15 2.88345
>> 6   15 2.66675
>> 7   47 3.29125
>> 8   15 2.44085
>> 9   29 2.43065
>> 10  29 2.43425
>> 11  29 2.42705
>> 12   0 1.12350
>> 13  15 5.10310
>> 14  15 0.99880
>> 15   0 2.22350
>> 16  47 2.18290
>> 17  47 2.21570
>> 18  47 2.46450
>> 19  29 2.53350
>> 20  47 2.78550
>> 21  47 3.06920
>> 22   0 2.65660
>> 23   0 2.16650
>> 24   0 2.20620
>> 25   0 1.91870
>> 26   0 2.45160
>> 27   0 2.35980
>> 28   0 1.99900
>> 29   0 2.16300
>> 30  15 2.00460
>> 31  15 2.56220
>> 32  47 2.38650
>> 33  29 2.65820
>> 34  29 2.47670
>> 35  47 2.23770
>> 36  29 2.40080
>> 37  47 2.12650
>> 38  29 2.35950
>>
>> I really apreciate your help
>>
>> Sincerly...
>>
>> MJMaechler (Zurich, Switzerland)
>> ______________________________________________
>> R-help at 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.
>>
>
>        [[alternative HTML version deleted]]
>
>
> ______________________________________________
> R-help at 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.
>
>



-- 
Felix Andrews / 安福立
Postdoctoral Fellow
Integrated Catchment Assessment and Management (iCAM) Centre
Fenner School of Environment and Society [Bldg 48a]
The Australian National University
Canberra ACT 0200 Australia
M: +61 410 400 963
T: + 61 2 6125 4670
E: felix.andrews at anu.edu.au
CRICOS Provider No. 00120C
-- 
http://www.neurofractal.org/felix/