Thank you for your comment. The difference is due to what I am calling the “sample”. In the app, “n” refers to the size of the block, i.e. max(rnorm(n)) which I replicate 20000 times. ie.

par(mfrow=c(2,1), las=1, cex.main=1)

hist(replicate(20000, max(rnorm(10))), xlim=c(0,5), col=’purple’)

hist(replicate(20000, max(rnorm(100))), xlim=c(0,5), col=’purple’)

Sorry for any confusion I have caused.

]]>However, there should be less values for smaller n in the histogram, right? I understand you compute max(rnorm(365)) n times, which would relate to the following R code:

par(mfrow=c(2,1), las=1, cex.main=1)

hist(replicate(10, max(rnorm(365))), xlim=c(0,5), col=”purple”)

hist(replicate(100, max(rnorm(365))), xlim=c(0,5), col=”purple”)

I also do not see such a strong sample size dependency as with the app…

Any thoughts on this? ]]>