[R-br] Como armazenar o valor-p de um teste de normalidade multivariada?

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 elas já estão disposta no R. Como fazer isso?Tem como ver como esta função foi implementada no R?Grato.Fernando