para <- vec2par(c(10,100), type="nor")
# The two output values from these two lines should be
# very similar: (1) theoretical and (2) simulation
expect.max.ostat(10, para=para, cdf=cdfnor, pdf=pdfnor)
mean(sapply(1:1000, function(x) { max(rlmomco(10,para))}))
Run the code above in your browser using DataLab