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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2012/7/5 Samuel Carvalho <samukajm@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] =============================================