Thank you, Riccardo! The script is working, I am only puzzled by one thing: in the second half the standard procedure is run once more and the HJC column is added. But between the first and the second run the result changes dramatically between the first three criterions. You can see the result and the aforementioned inconsistency in the attachment.
@Allin I noticed and I just assumed it was implicitly referring to known results with simulations on the other criterions(he quotes Lutkepol 1985 but my subscription doesn't allow me to read it),
http://onlinelibrary.wiley.com/doi/10.1111/j.1467-9892.1985.tb00396.x/abstract. It definetly warrants deeper investigation. It would be interesting to rerun his montecarlo simulations and see how different criteria compares.
Again thanks for your prompt responses.