Utilize a library fdth para estatísticas descritivas. É uma excelente biblioteca. Baixe o pdf do library para maiores detalhes. library(fdth) x <- rnorm(1e5, mean=5, sd=1) plot(fdt(x, k=100), type='d', col=heat.colors(100)) curve(dnorm(x, mean=5, sd=1), col='darkgreen', add=TRUE, lwd=2) legend(par('usr')[1],par('usr')[4]-0.1,'Valores observados',pch=1,bty='n',col=heat.colors(100)) legend(par('usr')[1],par('usr')[4]-0.15,'Valores estimados',lty=1,bty='n',col='darkgreen') Como já exposto por Lucas, consulte o help da função legend para maiores detalhes! (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5076 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}