Histogram, estimated normal pdf (red) and density estimation (green). Source code:
datu1=rnorm(40,22,4)
datu1
[1] 12.07930 24.23205 24.22420 18.61956 25.47271 24.76057 17.13244 27.69723
[9] 20.90116 21.33277 20.60566 17.56045 21.60253 19.78497 18.54969 25.88454
[17] 23.29435 19.53775 30.82505 22.27685 13.24867 27.61367 20.69942 20.68274
[25] 20.66159 19.21159 20.08311 25.95496 22.80454 18.89508 28.19135 27.58037
[33] 24.54323 18.76870 23.44803 23.30925 17.82521 22.78719 19.97626 15.28610
datu2=rnorm(40,32,5)
datu2
[1] 29.07175 24.00086 39.36396 29.60436 36.79717 35.55459 22.48502 35.60534
[9] 34.06736 43.97922 32.64646 33.80424 36.51349 35.00986 25.48278 26.98489
[17] 31.81420 38.03397 41.53816 33.46317 25.45158 25.45902 33.25066 24.30571
[25] 34.59699 31.87765 35.56252 39.12881 30.74369 26.54277 34.49255 36.97809
[33] 32.23456 35.01120 35.63466 28.65729 30.54426 33.27269 27.59977 26.73819
denak=stack(list(bat=datu1,bi=datu2))
denak$values
mean(denak$values)
sd(denak$values)
png("dentsitate_histograma_0001.png",width=1400,height=500)
par(mfrow=c(1,2))
hist(denak$values,col=c("gray"),prob=T,ylim=c(0,0.06),yaxt="n",ylab="",xaxt="n",xlab="",main="")
axis(1,at=seq(10,45,by=5),pos=0)
titulu1=expression(atop(Histograma~"eta",datuetatik~zenbatetsitako~banakuntza~normala))
curve(dnorm(x, mean =27.02, sd =7.04), add = T,col="red",lwd=5)
title(titulu1)
hist(denak$values,col=c("gray"),prob=T,yaxt="n",ylim=c(0,0.06),ylab="",xaxt="n",xlab="",main="")
axis(1,at=seq(10,45,by=5),pos=0)
titulu2=expression(atop(Histograma~"eta",dentsitate~zenbatespen~lerro~bat))
title(titulu2)
lines(density(denak$values, na.rm = T, from = 10, to = 45),col=c("green"),lwd=5)
dev.off()