Aplicar teste Kruskal Wallis a varias colunas

Olá pessoal, necessito de uma ajuda para a seguinte questão:

tenho uma planilha com algumas variáveis resposta, e queria executar o teste de Kruskal Wallis em todas ao mesmo tempo, isso é possível?

kruskal.test(coleta1$VARIAVEIS~coleta1$Tto)

Muito obrigado.

Não sei se entendi direito a pergunta, mas aqui vai a minha tentativa de resposta:

# Dados
respostas <- dplyr::select(airquality, Ozone:Temp)
grupo <- airquality$Month

# Iterar nas respostas
purrr::map(respostas, kruskal.test, grupo)
#> $Ozone
#> 
#>  Kruskal-Wallis rank sum test
#> 
#> data:  .x[[i]] and grupo
#> Kruskal-Wallis chi-squared = 29.267, df = 4, p-value = 6.901e-06
#> 
#> 
#> $Solar.R
#> 
#>  Kruskal-Wallis rank sum test
#> 
#> data:  .x[[i]] and grupo
#> Kruskal-Wallis chi-squared = 7.9246, df = 4, p-value = 0.09438
#> 
#> 
#> $Wind
#> 
#>  Kruskal-Wallis rank sum test
#> 
#> data:  .x[[i]] and grupo
#> Kruskal-Wallis chi-squared = 12.52, df = 4, p-value = 0.01388
#> 
#> 
#> $Temp
#> 
#>  Kruskal-Wallis rank sum test
#> 
#> data:  .x[[i]] and grupo
#> Kruskal-Wallis chi-squared = 73.328, df = 4, p-value = 4.496e-15

Created on 2021-03-03 by the reprex package (v1.0.0)

Excepcional!! muito obrigado!!