Un nadal aproximadament normal

x <- seq(-1,1, length.out = 500)
plot(0,0, type="n", xlim=c(-1,1), ylim=c(0,.6), cex.axis=.5, 
     xlab="i que tinguem un nadal aproximadament normal",
     ylab="Bon 2022")
curve(.6*dnorm(x, sd=.4), lwd=2, add=T) # cova
curve(.02*dnorm(x, mean=-.6, sd=.1), col="brown", add=T,n=500) # bou
curve(.02*dnorm(x, mean=.6, sd=.12), col="gray20", add=T,n=500) # mula
polygon(x,.04*dnorm(x, mean=.17, sd=.055/1.1), col="yellow", border=NA) # corona Maria
polygon(x,.04*dnorm(x, mean=.17, sd=.055), col="blue", border=NA) # Maria
polygon(x,.04*dnorm(x, mean=-.17, sd=.045/1.1), col="yellow", border=NA) # corona Josep
polygon(x,.04*dnorm(x, mean=-.17, sd=.045), col="gray40", border=NA) # Josep
polygon(x,.012*dnorm(x, sd=.045/1.2), col="yellow", border=NA) # corona Jesús
polygon(x,.012*dnorm(x, sd=.045), col="lightpink2", border=NA) # Jesús
text(.82,-.01,labels="PLB Nadal 2021", cex=.5)

Vegeu també https://rpubs.com/perelopez/nadala2018 i https://rpubs.com/perelopez/nadala2019