climas = read_csv(here::here("data/clima_cg_jp-semanal.csv"),
col_types = "cTdddddd")
glimpse(climas)
## Observations: 2,748
## Variables: 8
## $ cidade <chr> "Campina Grande", "Campina Grande", "Campina Grande", "Campi...
## $ semana <dttm> 1992-12-27, 1993-01-03, 1993-01-10, 1993-01-31, 1993-02-07,...
## $ tmedia <dbl> 26.13333, 26.11905, 25.76667, 25.74000, 26.31429, 26.28571, ...
## $ tmax <dbl> 30.4, 32.4, 32.2, 32.0, 32.7, 32.7, 32.3, 32.3, 32.1, 31.2, ...
## $ tmin <dbl> 20.7, 19.3, 19.7, 19.9, 19.6, 20.0, 20.4, 21.2, 19.0, 19.0, ...
## $ chuva <dbl> 0.0, 0.0, 0.0, 0.4, 0.3, 0.0, 4.9, 0.0, 0.0, 6.1, 0.4, 1.2, ...
## $ mes <dbl> 12, 1, 1, 1, 2, 2, 2, 2, 10, 11, 11, 11, 11, 12, 12, 12, 12,...
## $ ano <dbl> 1992, 1993, 1993, 1993, 1993, 1993, 1993, 1993, 1993, 1993, ...
Temperatura
Onde é mais quente?
media_cidade = aggregate(climas$tmedia, list(climas$cidade), mean)
colnames(media_cidade) <- c("Cidade", "Temperatura média")
ggplot(media_cidade) +
geom_bar(
aes(x = `Cidade`, y = `Temperatura média`),
stat = "identity",
color = "red",
fill = "pink"
)

A média de temperatura é maior em João Pessoa, mas não difere tanto de Campina Grande.
As temperaturas máximas semanais (o pico do calor) são tão diferentes quanto as médias? Mais?
ggplot(climas) +
geom_point(aes(x = semana, y = tmax, colour = "Máxima")) +
geom_point(aes(x = semana, y = tmedia, colour = "Média")) +
xlab("Ano") + ylab("Temperatura")

A maioria da média de temperatura semanal está entre 25 e 30 °C intervalo esse que a média e a máxima mais se encontram. A temperatura máxima está majoritariamente entre 30 e 35 °C, com uma diferença de no máximo 10 °C para a média, mas que no ambiente fazem muita diferença (rsrs).
Quais os meses mais quentes e mais frios?
media_meses = aggregate(climas$tmedia, list(climas$mes), mean)
meses=c("JAN","FEV","MAR","ABR","MAI","JUN","JUL","AGO","SET","OUT","NOV","DEZ")
media_meses$mes=meses
colnames(media_meses) <- c("Número do mês", "Temperatura média", "Mês")
ggplot(media_meses)+
geom_bar(
aes(x = reorder(`Mês`, `Temperatura média`), y = `Temperatura média`),
stat = "identity",
color = "red",
fill = "pink"
) + xlab("Mês") + ylab("Temperatura média")

Calculou-se a média das temperaturas médias semanais de todos os anos agrupadas por mês. Apesar da diferença não ser tão gritante, março se torna o mês mais quente pois é fim de verão, já julho o mais frio por se encontrar no período de inverno.
Qual foi o São João mais frio de CG que está nos nossos dados?
filter_media_sj = filter(climas, cidade=="Campina Grande", mes=="6")
media_sj = aggregate(filter_media_sj$tmedia, list(filter_media_sj$ano), mean)
colnames(media_sj) <- c("Ano", "Temperatura média")
ggplot(media_sj) +
geom_point(aes(x = `Ano`, y = `Temperatura média`)) +
geom_line(
aes(x = `Ano`, y = `Temperatura média`)) +
theme(axis.text.x = element_text(angle = 45, vjust = 0.5)) +
scale_x_continuous("Ano", labels = as.character(media_sj$Ano), breaks = media_sj$Ano)

Calculando a média das médias de temperatura semanais, 2004 teve o São João mais frio nesses 23 anos.
Chuva
Quanto chove por semana em JP e CG?
media_chuva = aggregate(climas$chuva, list(climas$cidade), mean)
colnames(media_chuva) <- c("Cidade", "Média de chuva por semana")
ggplot(media_chuva) +
geom_bar(
aes(x = `Cidade`, y = `Média de chuva por semana`),
stat = "identity",
color = "red",
fill = "pink"
)

Talvez por ser litoral, ocorram mais chuvas em João Pessoa, já em Campina Grande choveu menos da metade do que na outra cidade.
No geral, onde chove mais?
soma_chuva = aggregate(climas$chuva, list(climas$cidade), sum)
colnames(soma_chuva) <- c("Cidade", "Chuva total")
ggplot(soma_chuva) +
geom_bar(
aes(x = `Cidade`, y = `Chuva total`),
stat = "identity",
color = "red",
fill = "pink"
)

Em João Pessoa chove bem mais que em Campina Grande, por ser litoral, as temperaturas são mais elevadas e as chuvas mais recorrentes.
A distribuição é diferente em meses específicos?
Campina Grande
chuva_cg = filter(climas, cidade=="Campina Grande")
ggplot(chuva_cg, aes(x = semana, y = chuva, color = mes)) +
geom_bar(stat = "identity") +
facet_wrap(~ mes, nrow=4, ncol=3) +
theme(legend.position = "none") +
xlab("Ano") + ylab("Chuva")

João Pessoa
chuva_jp = filter(climas, cidade=="João Pessoa")
ggplot(chuva_jp, aes(x = semana, y = chuva, color = mes)) +
geom_bar(stat = "identity") +
facet_wrap(~ mes, nrow=4, ncol=3) +
theme(legend.position = "none") +
xlab("Ano") + ylab("Chuva")

Nas duas cidades, os últimos 3 meses do ano (oututbro, novembro e dezembro) são os mais escassos de chuva, período esse que ocorre a primavera e o clima está mais ameno. Esse decaimento do volume de chuvas começa a ocorrer nos meses de inverno (julho, agosto, setembro).