Olá a todos! O objetivo desse notebook é responder algumas perguntas relacionadas ao clima de algumas cidades da Paraiba. Dentre elas, selecionamos os dados das seguinte cidades:

Tal análise é importante para que possamos nos aprofundar nos questionamentos a respeito dos esteriótipos comumente falados no dia a dia sobre cada uma dessas cidades. Por exemplo, por ser uma cidade situada no sertão, Patos leva a fama de uma cidade muito quente, já Campina Grande, cidade situada no planalto da Borborema, é uma cidade com o clima mais fresco e assim por diante.

Para realizar essa análise utilizamos os dados do Instituto Nacional de Meteorologia (INMET) desde o ano de 2010.

Primeiramente vamos carregar nossos dados

clima = 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()
    )
)

1. Podemos começar olhando como foram os ventos nesses 3 cidades no ano de 2019

Pelo gráfico podemos entender que João Pessoa é uma cidade com poucos indícios de ventos comparado a Campina Grande e Patos, oque é bem interessante, visto que a mesma é uma cidade litorânea. Para Campina Grande, tal análise era esperada devido a sua localização ser no planalto da Borborema, ou seja, uma posição mais alta se comparada as outras duas cidades.

Podemos explorar mais tentando observar agora como é a concentração de ventos ao longo dos meses

Aqui podemos ver mais claramente que João Pessoa tem ventos médios similares a Patos apenas de Janeiro a Junho, depois disso os ventos em Patos aumentam até valores proximos aos de Campina Grande, tal fator pode ser explicado devido a altitude de cada cidade, após uma rápida googlada temos os seguintes valores de altitude para cada cidade:

2. Se olharmos apenas para as semanas no período de janeiro e fevereiro (pico de verão), como foi o calor das 3 cidades nos últimos anos ?

clima %>% 
    filter(mes == 1 | mes == 2) %>% 
    ggplot(aes(x = temp_media)) + 
    facet_wrap(~ cidade, ncol = 1) + 
    geom_histogram(binwidth = .5, fill = "coral", color = "black") + 
    geom_rug()
Warning: Removed 10 rows containing non-finite values (stat_bin).

Para responder a essa pergunta, o gráfico acima é certeiro em nos mostrar a distribuição das temperaturas durantes janeiro e fevereiro nos últimos anos.

Podemos chegar a conclusão de que Campina Grande possui histórico de variação de temperatura de 22° até 28° no verão, com maior frequência de 26°.

Já em Patos a temperatura varia muito durante o verão, mas é sempre bem quente, estando entre 27° até 32°.

Em João Pessoa temos um concentração de valores de temperatura próximos a 28°~29°.

3. Que tal agora olharmos como a umidade das 3 cidades varia ao longo do ano ?

clima %>% 
    filter(mes == 1 | mes == 2) %>% 
    ggplot(aes(x = umidade)) + 
    facet_wrap(~ cidade, ncol = 1) + 
    geom_histogram(binwidth = .5, fill = "steelblue", color = "black") + 
    geom_rug()
Warning: Removed 10 rows containing non-finite values (stat_bin).

Essa visualização nos dá um poder comparativo muito grande, no qual podemos enxergar nítidamente que Patos dentre as 3 cidades é a que possui menores índices de umidade. Já Campina Grande e João Pessoa possuem valores mais altos de umidade. Dentre elas João Pessoa possui valores mais concentrados na faixa de 70~80 umidade. Tal fator pode ser devido a localização das duas cidades, seja a altitude ou a litoralidade.

4. Pra fechar com chave de ouro, vamos investigar como é a temperatura durante as festas juninas nas 3 cidades nos últimos anos ?

clima %>% 
    filter(mes == 6) %>% 
    ggplot(aes(x = temp_media)) + 
    facet_wrap(~ cidade, ncol = 1) + 
    geom_density(binwidth = .5, fill = "coral", color = "black") + 
    geom_rug()
Warning: Ignoring unknown parameters: binwidth
Warning: Removed 8 rows containing non-finite values (stat_density).

Para mudar um poucos os gráficos, que tal analizarmos agora usando um gráfico de densidade ? Durante as tão amadas festas juninas, podemos compreender que ao longo dos anos Campina Grande é cidade com o clima mais frio/agradável para se comemorar esse mês incrível, mais um fator para se levar em consideração onde você deve passar o seu São João :)

LS0tDQp0aXRsZTogIkVEQSBubyBjbGltYSBkZSBjaWRhZGVzIGRhIFBhcmFpYmEiDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQpPbMOhIGEgdG9kb3MhIE8gb2JqZXRpdm8gZGVzc2Ugbm90ZWJvb2sgw6kgcmVzcG9uZGVyIGFsZ3VtYXMgcGVyZ3VudGFzIHJlbGFjaW9uYWRhcyBhbyBjbGltYSBkZSBhbGd1bWFzIGNpZGFkZXMgZGEgUGFyYWliYS4gRGVudHJlIGVsYXMsIHNlbGVjaW9uYW1vcyBvcyBkYWRvcyBkYXMgc2VndWludGUgY2lkYWRlczoNCg0KKiBDYW1waW5hIEdyYW5kZQ0KKiBQYXRvcw0KKiBKb8OjbyBQZXNzb2ENCg0KVGFsIGFuw6FsaXNlIMOpIGltcG9ydGFudGUgcGFyYSBxdWUgcG9zc2Ftb3Mgbm9zIGFwcm9mdW5kYXIgbm9zIHF1ZXN0aW9uYW1lbnRvcyBhIHJlc3BlaXRvIGRvcyBlc3RlcmnDs3RpcG9zIGNvbXVtZW50ZSBmYWxhZG9zIG5vIGRpYSBhIGRpYSBzb2JyZSBjYWRhIHVtYSBkZXNzYXMgY2lkYWRlcy4gUG9yIGV4ZW1wbG8sIHBvciBzZXIgdW1hIGNpZGFkZSBzaXR1YWRhIG5vIHNlcnTDo28sIFBhdG9zIGxldmEgYSBmYW1hIGRlIHVtYSBjaWRhZGUgbXVpdG8gcXVlbnRlLCBqw6EgQ2FtcGluYSBHcmFuZGUsIGNpZGFkZSBzaXR1YWRhIG5vIHBsYW5hbHRvIGRhIEJvcmJvcmVtYSwgw6kgdW1hIGNpZGFkZSBjb20gbyBjbGltYSBtYWlzIGZyZXNjbyBlIGFzc2ltIHBvciBkaWFudGUuIA0KDQpQYXJhIHJlYWxpemFyIGVzc2EgYW7DoWxpc2UgdXRpbGl6YW1vcyBvcyBkYWRvcyBkbyBJbnN0aXR1dG8gTmFjaW9uYWwgZGUgTWV0ZW9yb2xvZ2lhIFsoSU5NRVQpXShodHRwczovL3RlbXBvLmlubWV0Lmdvdi5ici8pIGRlc2RlIG8gYW5vIGRlIDIwMTAuDQoNCiMjIyMgUHJpbWVpcmFtZW50ZSB2YW1vcyBjYXJyZWdhciBub3Nzb3MgZGFkb3MNCg0KYGBge3J9DQpjbGltYSA9IHJlYWRfY3N2KA0KICAgIGhlcmUoImRhdGEvdGVtcG8tanAtY2ctcHQuY3N2IiksDQogICAgY29sX3R5cGVzID0gY29scygNCiAgICAgICAgLmRlZmF1bHQgPSBjb2xfZG91YmxlKCksDQogICAgICAgIGNpZGFkZSA9IGNvbF9jaGFyYWN0ZXIoKSwNCiAgICAgICAgc2VtYW5hID0gY29sX2RhdGUoZm9ybWF0ID0gIiIpLA0KICAgICAgICBhbm8gPSBjb2xfaW50ZWdlcigpLA0KICAgICAgICBtZXMgPSBjb2xfaW50ZWdlcigpDQogICAgKQ0KKQ0KYGBgDQoNCg0KIyMjIyAxLiBQb2RlbW9zIGNvbWXDp2FyIG9saGFuZG8gY29tbyBmb3JhbSBvcyB2ZW50b3MgbmVzc2VzIDMgY2lkYWRlcyBubyBhbm8gZGUgMjAxOQ0KDQpgYGB7cn0NCmNsaW1hICU+JSANCiAgICBmaWx0ZXIoYW5vID09IDIwMTkpICU+JSANCiAgICBnZ3Bsb3QobWFwcGluZyA9IGFlcyh5PWNpZGFkZSwgeD12ZW50b19tZWRpbykpICsgDQogICAgZ2VvbV9qaXR0ZXIoaGVpZ2h0PS4xLCBjb2xvcj0ic3RlZWxibHVlIiwgYWxwaGE9LjUpICsNCiAgICBsYWJzKA0KICAgICAgICAgeSA9ICIiLCANCiAgICAgICAgIHggPSAiVmVudG8gTcOpZGlvIChtL3MpIg0KICAgICkNCmBgYA0KDQpQZWxvIGdyw6FmaWNvIHBvZGVtb3MgZW50ZW5kZXIgcXVlIEpvw6NvIFBlc3NvYSDDqSB1bWEgY2lkYWRlIGNvbSBwb3Vjb3MgaW5kw61jaW9zIGRlIHZlbnRvcyBjb21wYXJhZG8gYSBDYW1waW5hIEdyYW5kZSBlIFBhdG9zLCBvcXVlIMOpIGJlbSBpbnRlcmVzc2FudGUsIHZpc3RvIHF1ZSBhIG1lc21hIMOpIHVtYSBjaWRhZGUgbGl0b3LDom5lYS4gUGFyYSBDYW1waW5hIEdyYW5kZSwgdGFsIGFuw6FsaXNlIGVyYSBlc3BlcmFkYSBkZXZpZG8gYSBzdWEgbG9jYWxpemHDp8OjbyBzZXIgbm8gcGxhbmFsdG8gZGEgQm9yYm9yZW1hLCBvdSBzZWphLCB1bWEgcG9zacOnw6NvIG1haXMgYWx0YSBzZSBjb21wYXJhZGEgYXMgb3V0cmFzIGR1YXMgY2lkYWRlcy4NCg0KUG9kZW1vcyBleHBsb3JhciBtYWlzIHRlbnRhbmRvIG9ic2VydmFyIGFnb3JhIGNvbW8gw6kgYSBjb25jZW50cmHDp8OjbyBkZSB2ZW50b3MgYW8gbG9uZ28gZG9zIG1lc2VzDQoNCmBgYHtyfQ0KY2xpbWEgJT4lIA0KICAgIGZpbHRlcihhbm8gPT0gMjAxOSkgJT4lIA0KICAgIGdncGxvdChhZXMoeCA9IHNlbWFuYSwgeSA9IHZlbnRvX21lZGlvLCBjb2xvciA9IGNpZGFkZSkpICsgDQogICAgZ2VvbV9wb2ludCgpDQpgYGANCkFxdWkgcG9kZW1vcyB2ZXIgbWFpcyBjbGFyYW1lbnRlIHF1ZSBKb8OjbyBQZXNzb2EgdGVtIHZlbnRvcyBtw6lkaW9zIHNpbWlsYXJlcyBhIFBhdG9zIGFwZW5hcyBkZSBKYW5laXJvIGEgSnVuaG8sIGRlcG9pcyBkaXNzbyBvcyB2ZW50b3MgZW0gUGF0b3MgYXVtZW50YW0gYXTDqSB2YWxvcmVzIHByb3hpbW9zIGFvcyBkZSBDYW1waW5hIEdyYW5kZSwgdGFsIGZhdG9yIHBvZGUgc2VyIGV4cGxpY2FkbyBkZXZpZG8gYSBhbHRpdHVkZSBkZSBjYWRhIGNpZGFkZSwgYXDDs3MgdW1hIHLDoXBpZGEgZ29vZ2xhZGEgdGVtb3Mgb3Mgc2VndWludGVzIHZhbG9yZXMgZGUgYWx0aXR1ZGUgcGFyYSBjYWRhIGNpZGFkZToNCg0KKiBDYW1waW5hIEdyYW5kZTogNTUxbQ0KKiBQYXRvczogMjQwbQ0KKiBKb8OjbyBQZXNzb2E6IDM3bQ0KDQojIyMjIDIuIFNlIG9saGFybW9zIGFwZW5hcyBwYXJhIGFzIHNlbWFuYXMgbm8gcGVyw61vZG8gZGUgamFuZWlybyBlIGZldmVyZWlybyAocGljbyBkZSB2ZXLDo28pLCBjb21vIGZvaSBvIGNhbG9yIGRhcyAzIGNpZGFkZXMgbm9zIMO6bHRpbW9zIGFub3MgPw0KDQpgYGB7cn0NCmNsaW1hICU+JSANCiAgICBmaWx0ZXIobWVzID09IDEgfCBtZXMgPT0gMikgJT4lIA0KICAgIGdncGxvdChhZXMoeCA9IHRlbXBfbWVkaWEpKSArIA0KICAgIGZhY2V0X3dyYXAofiBjaWRhZGUsIG5jb2wgPSAxKSArIA0KICAgIGdlb21faGlzdG9ncmFtKGJpbndpZHRoID0gLjUsIGZpbGwgPSAiY29yYWwiLCBjb2xvciA9ICJibGFjayIpICsgDQogICAgZ2VvbV9ydWcoKQ0KYGBgDQoNClBhcmEgcmVzcG9uZGVyIGEgZXNzYSBwZXJndW50YSwgbyBncsOhZmljbyBhY2ltYSDDqSBjZXJ0ZWlybyBlbSBub3MgbW9zdHJhciBhIGRpc3RyaWJ1acOnw6NvIGRhcyB0ZW1wZXJhdHVyYXMgZHVyYW50ZXMgamFuZWlybyBlIGZldmVyZWlybyBub3Mgw7psdGltb3MgYW5vcy4gDQoNClBvZGVtb3MgY2hlZ2FyIGEgY29uY2x1c8OjbyBkZSBxdWUgQ2FtcGluYSBHcmFuZGUgcG9zc3VpIGhpc3TDs3JpY28gZGUgdmFyaWHDp8OjbyBkZSB0ZW1wZXJhdHVyYSBkZSAyMsKwIGF0w6kgMjjCsCBubyB2ZXLDo28sIGNvbSBtYWlvciBmcmVxdcOqbmNpYSBkZSAyNsKwLiANCg0KSsOhIGVtIFBhdG9zIGEgdGVtcGVyYXR1cmEgdmFyaWEgbXVpdG8gZHVyYW50ZSBvIHZlcsOjbywgbWFzIMOpIHNlbXByZSBiZW0gcXVlbnRlLCBlc3RhbmRvIGVudHJlIDI3wrAgYXTDqSAzMsKwLg0KDQpFbSBKb8OjbyBQZXNzb2EgdGVtb3MgdW0gY29uY2VudHJhw6fDo28gZGUgdmFsb3JlcyBkZSB0ZW1wZXJhdHVyYSBwcsOzeGltb3MgYSAyOMKwfjI5wrAuDQoNCiMjIyMgMy4gUXVlIHRhbCBhZ29yYSBvbGhhcm1vcyBjb21vIGEgdW1pZGFkZSBkYXMgMyBjaWRhZGVzIHZhcmlhIGFvIGxvbmdvIGRvIGFubyA/DQoNCmBgYHtyfQ0KY2xpbWEgJT4lIA0KICAgIGZpbHRlcihtZXMgPT0gMSB8IG1lcyA9PSAyKSAlPiUgDQogICAgZ2dwbG90KGFlcyh4ID0gdW1pZGFkZSkpICsgDQogICAgZmFjZXRfd3JhcCh+IGNpZGFkZSwgbmNvbCA9IDEpICsgDQogICAgZ2VvbV9oaXN0b2dyYW0oYmlud2lkdGggPSAuNSwgZmlsbCA9ICJzdGVlbGJsdWUiLCBjb2xvciA9ICJibGFjayIpICsgDQogICAgZ2VvbV9ydWcoKQ0KYGBgDQoNCkVzc2EgdmlzdWFsaXphw6fDo28gbm9zIGTDoSB1bSBwb2RlciBjb21wYXJhdGl2byBtdWl0byBncmFuZGUsIG5vIHF1YWwgcG9kZW1vcyBlbnhlcmdhciBuw610aWRhbWVudGUgcXVlIFBhdG9zIGRlbnRyZSBhcyAzIGNpZGFkZXMgw6kgYSBxdWUgcG9zc3VpIG1lbm9yZXMgw61uZGljZXMgZGUgdW1pZGFkZS4NCkrDoSBDYW1waW5hIEdyYW5kZSBlIEpvw6NvIFBlc3NvYSBwb3NzdWVtIHZhbG9yZXMgbWFpcyBhbHRvcyBkZSB1bWlkYWRlLiBEZW50cmUgZWxhcyBKb8OjbyBQZXNzb2EgcG9zc3VpIHZhbG9yZXMgbWFpcyBjb25jZW50cmFkb3MgbmEgZmFpeGEgZGUgNzB+ODAgdW1pZGFkZS4gVGFsIGZhdG9yIHBvZGUgc2VyIGRldmlkbyBhIGxvY2FsaXphw6fDo28gZGFzIGR1YXMgY2lkYWRlcywgc2VqYSBhIGFsdGl0dWRlIG91IGEgbGl0b3JhbGlkYWRlLg0KDQojIyMjIDQuIFByYSBmZWNoYXIgY29tIGNoYXZlIGRlIG91cm8sIHZhbW9zIGludmVzdGlnYXIgY29tbyDDqSBhIHRlbXBlcmF0dXJhIGR1cmFudGUgYXMgZmVzdGFzIGp1bmluYXMgbmFzIDMgY2lkYWRlcyBub3Mgw7psdGltb3MgYW5vcyA/DQoNCmBgYHtyfQ0KY2xpbWEgJT4lIA0KICAgIGZpbHRlcihtZXMgPT0gNikgJT4lIA0KICAgIGdncGxvdChhZXMoeCA9IHRlbXBfbWVkaWEpKSArIA0KICAgIGZhY2V0X3dyYXAofiBjaWRhZGUsIG5jb2wgPSAxKSArIA0KICAgIGdlb21fZGVuc2l0eShiaW53aWR0aCA9IC41LCBmaWxsID0gImNvcmFsIiwgY29sb3IgPSAiYmxhY2siKSArIA0KICAgIGdlb21fcnVnKCkNCmBgYA0KDQpQYXJhIG11ZGFyIHVtIHBvdWNvcyBvcyBncsOhZmljb3MsIHF1ZSB0YWwgYW5hbGl6YXJtb3MgYWdvcmEgdXNhbmRvIHVtIGdyw6FmaWNvIGRlIGRlbnNpZGFkZSA/IA0KRHVyYW50ZSBhcyB0w6NvIGFtYWRhcyBmZXN0YXMganVuaW5hcywgcG9kZW1vcyBjb21wcmVlbmRlciBxdWUgYW8gbG9uZ28gZG9zIGFub3MgQ2FtcGluYSBHcmFuZGUgw6kgY2lkYWRlIGNvbSBvIGNsaW1hIG1haXMgZnJpby9hZ3JhZMOhdmVsIHBhcmEgc2UgY29tZW1vcmFyIGVzc2UgbcOqcyBpbmNyw612ZWwsIG1haXMgdW0gZmF0b3IgcGFyYSBzZSBsZXZhciBlbSBjb25zaWRlcmHDp8OjbyBvbmRlIHZvY8OqIGRldmUgcGFzc2FyIG8gc2V1IFPDo28gSm/Do28gOik=