R help - Mar 2007 - lattice histogram

How to add mean,sd, number of observation in each panel for lattice histogram?

Aimin

Here is an example using the grid package to annotate the graphs:

library(lattice)
library(grid)
resp  <- rnorm(200)
group <- sample(c("G1", "G2", "G3"), replace =
TRUE, size = 100)
histogram(~ resp | group,
  panel = function(x, ...){
    std <- round(sd(x), 2)
    n <- length(x)
    m <- round(mean(x), 2)
    panel.histogram(x, ...)
    x1 <- unit(1, "npc") - unit(2, "mm")
    y1 <- unit(1, "npc") - unit(2, "mm")
    grid.text(label = bquote(n == .(n)), x = x1, y = y1, just =
"right")
    grid.text(label = bquote(hat(m) == .(m)), x = x1, y = y1 - unit(1,
"lines"), just = "right")
    grid.text(label = bquote(hat(s) == .(std)), x = x1, y = y1 -
unit(2, "lines"), just = "right")
    })
Best,

Renaud

2007/3/4, Aimin Yan <aiminy at iastate.edu>:
> How to add mean,sd, number of observation in each panel for lattice
histogram?
>
> Aimin
>
> ______________________________________________
> R-help at 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.
>

-- 
Renaud LANCELOT
D?partement Syst?mes Biologiques du CIRAD
CIRAD, Biological Systems Department

Campus International de Baillarguet
TA 30 / B
F34398 Montpellier
Tel   +33 (0)4 67 59 37 17
Secr. +33 (0)4 67 59 37 37
Fax   +33 (0)4 67 59 37 95

Steve,

it says resp  <- rnorm(200) must be resp  <- rnorm(100) 

regards

Felipe

-----Original Message-----
From: r-help-bounces at stat.math.ethz.ch
[mailto:r-help-bounces at stat.math.ethz.ch]On Behalf Of steve
Sent: Sunday, March 04, 2007 8:58 AM
To: r-help at stat.math.ethz.ch
Subject: Re: [R] lattice histogram


Hi,

When I tried this the groups had sizes 36 30 34 as expected, but
the annotations "n = " were 68 72 60 - twice as large. I don't 
understand why.

Steve

Renaud Lancelot wrote:
> Here is an example using the grid package to annotate the graphs:
> 
> library(lattice)
> library(grid)
> resp  <- rnorm(200)
> group <- sample(c("G1", "G2", "G3"),
replace = TRUE, size = 100)
> histogram(~ resp | group,
>   panel = function(x, ...){
>     std <- round(sd(x), 2)
>     n <- length(x)
>     m <- round(mean(x), 2)
>     panel.histogram(x, ...)
>     x1 <- unit(1, "npc") - unit(2, "mm")
>     y1 <- unit(1, "npc") - unit(2, "mm")
>     grid.text(label = bquote(n == .(n)), x = x1, y = y1, just =
"right")
>     grid.text(label = bquote(hat(m) == .(m)), x = x1, y = y1 - unit(1,
> "lines"), just = "right")
>     grid.text(label = bquote(hat(s) == .(std)), x = x1, y = y1 -
> unit(2, "lines"), just = "right")
>     })
______________________________________________
R-help at 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.