library(tidyverse)
library(here)
theme_set(theme_bw())

O que analizaremos são dados meteorológicos que foram coletadas de 3 cidades da paraiba: Campina Grande, João Pessoa e Patos. Os dados estão separados por semanas para melhorar a vizualização.

# SEMPRE read_csv NUNCA read.csv
clima_tudo = read_csv(
    here("data/tempo-jp-cg-pt.csv"),
    col_types = cols(
        .default = col_double(),
        cidade = col_character(),
        semana = col_date(format = ""),
        ano = col_integer(),
        mes = col_integer()
    )
)

clima = clima_tudo %>% 
    filter(ano == 2019)

clima %>% 
    slice(1:8)
## # A tibble: 8 x 11
##   cidade   semana     temp_max temp_media temp_min vento_medio vento_max umidade
##   <chr>    <date>        <dbl>      <dbl>    <dbl>       <dbl>     <dbl>   <dbl>
## 1 Campina… 2019-01-06     32.1       25.6     19.6        3.57       4.6    74.9
## 2 Campina… 2019-01-13     33.5       27.1     20          3.29       4.5    69  
## 3 Campina… 2019-01-20     33.3       26.7     19.6        4.01       6.5    72.5
## 4 Campina… 2019-01-27     32.4       26.4     19.6        3.53       5      75.0
## 5 Campina… 2019-02-03     32.2       25.3     19.5        2.95       4.8    83  
## 6 Campina… 2019-02-10     31.1       25.6     20.1        3.17       5      84.1
## 7 Campina… 2019-02-17     31.6       25.9     20.5        2.89       4.1    79.2
## 8 Campina… 2019-02-24     31.1       25.9     20.5        3.63       4.8    76.8
## # … with 3 more variables: chuva <dbl>, ano <int>, mes <int>

Primeiramente iremos gerar um gráfico simples da variação de temperatura media das 3 cidades que estamos analizando em 2019, que será necessario para compararmos com dados futuros.

clima %>% 
  ggplot(aes(x=temp_media, y = cidade))+
  geom_jitter(height = .1, color = "firebrick1")+
  labs(
         y = "", 
         x = "Temperatura Média Semanal (celcius)"
    )

A primeira analise a ser feita é sobre a distribuição de vento nas 3 cidades no ano de 2019, estamos analisando o parametro vento medio. Os resultados podem ser visto no grafico abaixo.

clima %>% 
  ggplot(aes(x = vento_medio, y = cidade))+
  geom_jitter(height = .1, color = "gold4")+
  labs(
         y = "", 
         x = "Velocidade media do vento em m/s"
    )

Como é possível ver nesse grafico de pontos:

A partir desse histograma abaixo, é possivel ver de outra maneira como se comporta a velocidade media do vento dessas cidades.

clima %>% 
  ggplot(aes(x=vento_medio))+
  facet_wrap(~ cidade)+
  geom_histogram(binwidth=.5, fill="deepskyblue", color="black")+
  labs(
         y = "contador de vezes", 
         x = "Velocidade do vento em m/s"
    )

Uma questão que pode ser analizada é se existe uma relação entre velocidade media do vento e temperatura media.

Podemos ver uma relação mais clara entre essas duas variaveis nas cidades de João Pessoa e Campina Grande.

Uma duvida que ficou sem explicação foi porque Joao Pessoa , uma cidade a beira mar possui uma velocidade media de vento bastante baixo. O que consegui entender é que pode dizer que seria por causa da quantidade de predios que existe a beira mar, que impede o fluxo de vento para dentro da cidade.

Agora, iremos analizar o calor dos meses de janeiro e fevereiro a partir de 2018

clima_ultimos_anos =clima_tudo %>% 
  filter(ano >= "2010" ,mes == "1" | mes =="2", !is.na(temp_media))
clima_ultimos_anos 
## # A tibble: 296 x 11
##    cidade  semana     temp_max temp_media temp_min vento_medio vento_max umidade
##    <chr>   <date>        <dbl>      <dbl>    <dbl>       <dbl>     <dbl>   <dbl>
##  1 Campin… 2010-01-03     31.4       25.5     21.2        4.08       5.3    76.2
##  2 Campin… 2010-01-10     32.1       25.6     21.3        3.95       6.6    76.0
##  3 Campin… 2010-01-17     31         24.2     20.5        2.86       5      85.3
##  4 Campin… 2010-01-24     31.2       26.1     21.2        3.69       5.5    77.8
##  5 Campin… 2010-01-31     32.1       26.0     21.4        3.97       6.5    76.0
##  6 Campin… 2010-02-07     32.2       26.0     20.4        3.40       6.3    77.8
##  7 Campin… 2010-02-14     31.3       25.8     21.5        3.40       5.5    80.4
##  8 Campin… 2010-02-21     32.9       27.0     22.1        3.54       5.5    76.0
##  9 Campin… 2010-02-28     34.1       28.0     22          3.06       4      71  
## 10 Campin… 2011-01-02     32.8       25.5     20.4        3.93       5      74.8
## # … with 286 more rows, and 3 more variables: chuva <dbl>, ano <int>, mes <int>
clima_ultimos_anos %>% 
  ggplot(aes(x = temp_media, y = cidade))+
    geom_point(alpha=.6)+
    facet_wrap(~ano, ncol=3)+
    labs(
         y = "", 
         x = "Temperatura media nos meses de janeiro e fevereiro"
    )

clima_ultimos_anos %>% 
  ggplot(aes(x = temp_media, y = cidade))+
    geom_jitter(height = .1, color = "firebrick1")+
    labs(
         y = "", 
         x = "Temperatura media nos meses de janeiro e fevereiro"
    )

Nos graficos gerado acima podemos tirar algumas conclusoes: