Gostaria de comparar as taxas de erro tipo I e o desempenho do poder de  dois testes de normalidade multivariada do R :  um é a função mvShapiroTest do pacote mvnormtest  e o outro teste é a função mvnorm.etest do pacote energy.  Essa comparação será feita por meio de simulação Monte Carlo.  Gostaria de saber se existe como armazenar  o valor-p sem ter que implementar essas funções uma vez que ela

s já estão disposta no R. Como fazer isso?

Tem como ver como esta função foi implementada no R?

Grato.

Fernando