[R-br] funcao histogram

[Walmes Zeviani]

Uma sugestão seria

my.densityplot <- function(x, cols, ...){
  panel.densityplot(x, ...)
  dx <- density(x)
  breaks <- sort(c(range(dx$x), quantile(x, c(1,3)/4)))
  fx <- approxfun(dx$x, dx$y)
  do.polygon <- function(x){
    y <- fx(x)
    return(list(x=c(min(x), x, max(x)), y=c(0, y, 0)))
  }
  seqs <- lapply(1:(length(breaks)-1),
                 function(i){
                   x <- seq(breaks[i], breaks[i+1], l=30)
                   do.polygon(x)
                 })
  for(i in 1:length(seqs)){
    seqs[[i]]$col <- cols[i]
  }
  lapply(seqs, function(i) do.call(panel.polygon, i))
  panel.mathdensity(dmath = dnorm, col = "black",
                    args = list(mean=mean(x),sd=sd(x)), lty=3)
}

require(RColorBrewer)
cols <- brewer.pal(3, "Greens")

densityplot(~height|voice.part, data=singer, cols=cols,
panel=my.densityplot)

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================


2012/7/5 Samuel Carvalho <samukajm em yahoo.com.br>

> Falae Walmes,
> Cara é basicamente isso mesmo que queria quando voce gera o ultimo gráfico
> com as hachuras nos diferentes quantis.
> Mas voce teria alguma dica de como implementar isso dentro da função
> xyplot? Talvez a saída seja pelo panel.function() e depois pelo
> panel.histogram()
> Consegue dar um help?
> Valeu
>
> *====================================*
> *Samuel P. C. Carvalho
> *Mestre em Ciências Florestais [UFLA]
> Doutorando em Recursos Florestais [ESALQ/USP]
> =============================================
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120705/63e6c5aa/attachment.html>