
Bom dia pessoal, Estou precisando de uma ajuda para gerar um histograma de frequencias para dados agrupados e estou reproduzindo o exemplo da função ?histogram A seguir o #RCMR require(lattice) ?histogram histogram( ~ height | voice.part, prob=T, data = singer, xlab = "Height (inches)", type = "density", panel = function(x, ...) { panel.histogram(x, ...) panel.mathdensity(dmath = dnorm, col = "black", args = list(mean=mean(x),sd=sd(x))) panel.abline(v=qnorm(0.25, mean=mean(x), sd=sd(x)), lty=2, col='blue') panel.abline(v=qnorm(0.75, mean=mean(x), sd=sd(x)), lty=2, col='blue') } ) Fiz uma adaptação no exemplo e inseri as duas ultimas linhas para que pudesse ser retornado uma linha para os quantis da normal que me retornam 25% e 75% dos dados porem queria que ao invés de linhas eu possa fazer hachuras divindo a distribuição em tres pontos, sendo que fossem feitas tres hachuras, ou seja, uma primeira para a area abaixo de 0.25, uma segunda entre 0.25 e 0.75 e uma terceira para a area acima de 0.75 Espero que tenha conseguido expressar o que preciso. Obrigado a todos Samuel ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

Samuel, use o panel.rect(), veja histogram( ~ height | voice.part, prob=T, data = singer, xlab = "Height (inches)", type = "density", panel = function(x, ...) { * panel.rect(40, 0, qnorm(0.25, mean=mean(x), sd=sd(x)), 1, col="gray90") * panel.histogram(x, ...) panel.mathdensity(dmath = dnorm, col = "black", args = list(mean=mean(x),sd=sd(x))) panel.abline(v=qnorm(0.25, mean=mean(x), sd=sd(x)), lty=2, col='blue') panel.abline(v=qnorm(0.75, mean=mean(x), sd=sd(x)), lty=2, col='blue') } ) À 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 ==========================================================================
participantes (2)
-
Samuel Carvalho
-
Walmes Zeviani