Estou tentando estender a resposta de uma perguntaR: filtrando dados e calculando a correlação.

Para obter a correlação de temperatura e umidade para cada mês do ano (1 = janeiro), teríamos que fazer o mesmo para cada mês (12 vezes).

cor(airquality[airquality$Month == 1, c("Temp", "Humidity")])

Existe alguma maneira de fazer todos os meses automaticamente?

No meu caso, tenho mais de 30 grupos (não meses, mas espécies) para os quais gostaria de testar correlações, só queria saber se existe uma maneira mais rápida do que fazê-lo um por um.

Obrigado!