Eu tenho um dataframe onde todas as minhas 90 variáveis ​​têm dados inteiros, do tipo:

código | variável1 | variável2 | variável3 | ...

AB | 2 | 3 | 10 | ...

AH | 4 | 6 | 8 | ...

BC | 1 | 5 | 9 | ...

... | ... | ... | ...

Eu quero aplicar um teste de shapiro (shapiro.test {stats}) ao meu dataframe por variável e escrever os resultados em uma tabela como:

variable_name | W | valor p

Alguém tem uma pista?