[R] lattice histogram
Aimin Yan
aiminy at iastate.edu
Mon Mar 5 05:46:26 CET 2007
thank you very much. Your code almost solve my problem, but I have a
further question.
In my data, there is no observation in some group, I want to label that
panel by
n=0
hat(m)=NA
hat(s)=NA.
I try to modify your panel function, but it doesn't work out. Do you know
how to add something to your panel
function so that it can deal with some group that has 0 observation.
Aimin
At 02:54 AM 3/4/2007, 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")
> })
>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
More information about the R-help
mailing list