Pacotes

library(tidyverse)  # manipulação de dados
library(metan)
library(lubridate)    # gráfico de radar
library(rio)        # importação/exportação de dados

clima <- import("https://bit.ly/inmet_fred_2020")

# gerar tabelas html
print_tbl <- function(table, digits = 3, n = NULL, ...){
  if(!missing(n)){
    knitr::kable(head(table, n = n), booktabs = TRUE, digits = digits, ...)
  } else{
  knitr::kable(table, booktabs = TRUE, digits = digits, ...)
  }
}

Radiação solar (média mensal)

df_rad <- 
clima %>% 
  sum_by(DIA, MES) %>% 
  means_by(MES) %>% 
  select(MES, RADIAC)

ggplot(df_rad, aes(factor(MES), RADIAC, group = 1)) + 
  stat_summary(geom = "point", 
               fun = mean) +
  stat_summary(geom = "line") + 
  theme(panel.grid.minor = element_blank(),
        legend.position = "bottom",
        legend.title = element_blank(),
        axis.title = element_text(size = 12),
        axis.text = element_text(size = 12)) + 
  labs(title = "Radiação solar",
       subtitle = "Estação INMET - UFSM-FW",
       caption = "Elaboração - Olivoto 2021",
       x = "Mês do ano",
       y = expression(paste("Radiação solar em média da soma diária (KJ m"^-2~")")))

Radiação solar em dois meses

O seguinte gráfico mostra a média horária da radiação solar em dois meses do ano (junho e dezembro).

#rad
df_rad <- 
  clima %>% 
  select(MES, HORA, RADIAC) %>% 
  subset(MES %in% c("6", "12")) %>% 
  mutate(Mês = ifelse(MES == "6", "Junho", "Dezembro"))


ggplot(df_rad, aes(HORA, RADIAC, color = factor(Mês), group = Mês)) + 
  stat_summary(geom = "point", 
               fun = mean) +
  stat_summary(geom = "line") + 
  stat_summary(geom = "errorbar", width = 0.5) +
  theme(panel.grid.minor = element_blank(),
        legend.position = "bottom",
        legend.title = element_blank(),
        axis.title = element_text(size = 12),
        axis.text = element_text(size = 12)) + 
  scale_x_continuous(breaks = seq(0,24, by = 2)) + 
  labs(title = "Radiação em dois meses do ano",
       subtitle = "Estação INMET - UFSM-FW",
       caption = "Elaboração - Olivoto 2021",
       x = "Hora do dia",
       y = expression(paste("Radiação média horária (KJ m"^-2~")")))