<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Diogo,<br><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Você vai ter que multiplicar a função densidade da normal pela área do histograma. Segue o código para você juntar as partes.<br><span style="font-family:monospace,monospace"><br># dados<br>res <- c(rnorm(1000, 20, 2),<br>         rnorm(1000, 15, 1.5),<br>         rnorm(1000, 25, 2.5),<br>         rnorm(1000, 25, 2.5))<br><br># fatores<br>fac <- factor(gl(4, 1000, labels = c("F1", "F2", "F3", "F4")))<br><br>library(lattice)<br><br>histogram(~res | fac,<br>          xlab = "Resultado",<br>          ylab = "Frequência",<br>          type = "density",<br>          layout = c(1, 4),<br>          endpoints = c(0, 40),<br>          nint = 30,<br>          scales = list(x = list(at = c(0, 10, 20, 30, 40))),<br>          panel = function(x, ...) {<br>              panel.histogram(x,<br>                              col = rgb(62, 88, 172, max = 255),<br>                              ...)<br>              panel.mathdensity(n = 301,<br>                                dmath = dnorm,<br>                                args = list(mean = mean(x), sd = sd(x)),<br>                                ...,<br>                                col = rgb(13, 141, 70, max = 255),<br>                                lwd = 2)<br>          })<br><br># Fazendo usando a hist().<br>y <- rnorm(100)<br>ht <- hist(y, breaks = seq(-5, 5, by = 0.5))<br><br># Área do histograma.<br>a <- sum(diff(ht$breaks) * ht$counts)<br>a<br><br>plot(ht)<br>curve(a * dnorm(x, mean(y), sd(y)),<br>      add = TRUE, col = 2)</span><br><br>À diposição.<br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Walmes.<br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div>​</div>