Ok Ivan...

Obrigado pela dica....
Eu encontrei o erro aqui......Estava usando a função matrix() e na verdade no meu caso a função é as.matrix().

C=as.matrix(carvão[,(9:14)]) 
  
  mvShapiro.Test(C) #Teste de normalidade multivariado



Thiago de Paula Protásio
Acadêmico de Engenharia Florestal
Universidade Federal de Lavras
Laboratório de Energia da Biomassa Florestal
Laboratório de Biomateriais
(035) 9183-2246