[R-br] funcao histogram

[Walmes Zeviani]

Samuel, tente adaptar isso aqui

set.seed(123)
x <- rnorm(1000, 5, 2)
plot(density(x))
breaks <- sort(c(range(x), quantile(x, c(1,3)/4)))

dx <- density(x)
fx <- approxfun(dx$x, dx$y)
fx(0.5)

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)
               })
str(seqs)

plot(density(x))
polygon(seqs[[1]], density=2)

seqs[[1]]$density <- 3
seqs[[2]]$density <- 5
seqs[[3]]$density <- 9

plot(density(x))
lapply(seqs, function(i) do.call(polygon, i))

À disposição.
Walmes.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120704/8a37d0ef/attachment.html>